refactor(source_io): replace std::exception with WARNING_QUIT

Replace standard exception handling with project-specific error handling
in source_io directory to comply with cpp-code-style skill guidelines.

Changes:
- Replace all instances of throw with ModuleBase::WARNING_QUIT
  - throw std::runtime_error -> WARNING_QUIT
  - throw std::logic_error -> WARNING_QUIT
  - throw std::invalid_argument -> WARNING_QUIT
  - throw std::domain_error -> WARNING_QUIT
- Remove unnecessary #include <stdexcept>
- Add #include "source_base/tool_quit.h" where needed

Files modified:
- module_restart/restart.cpp (4 instances)
- module_chgpot/write_libxc_r.cpp (9 instances)
- module_wannier/to_wannier90_lcao.h (1 instance)
- module_wannier/to_wannier90_lcao_in_pw.h (1 instance)
- module_wannier/to_wannier90_pw.h (1 instance)
- module_wannier/to_wannier90.cpp (1 instance)
- module_unk/berryphase.h (1 instance)
- module_parameter/read_input_tool.h (1 instance)
- module_parameter/input_conv.cpp (1 instance)

This follows the ABACUS code style rule that requires using
ModuleBase::WARNING/WARNING_QUIT instead of standard exceptions.

Author: abacus_fixer

source/source_io/module_parameter/input_conv.cpp

@@ -2,6 +2,7 @@
 
 #include "source_base/global_function.h"
 #include "source_base/global_variable.h"
+#include "source_base/tool_quit.h"
 #include "source_cell/module_symmetry/symmetry.h"
 #include "source_cell/unitcell.h"
 #include "source_estate/occupy.h"
@@ -398,7 +399,7 @@ void Input_Conv::Convert()
             }
             else
             {
-                throw std::invalid_argument(std::string(__FILE__)+" line "+std::to_string(__LINE__));
+                ModuleBase::WARNING_QUIT("input_conv", "Invalid exx_hybrid_alpha type");
             }
         }
         if(!erfc_alpha.empty())

source/source_io/module_parameter/read_input_tool.h

@@ -1,9 +1,9 @@
 #include <string>
 #include <vector>
-#include <stdexcept>
 #ifdef __MPI
 #include "source_base/parallel_common.h"
 #endif
+#include "source_base/tool_quit.h"
 
 #define strvalue item.str_values[0]
 #define intvalue std::stoi(item.str_values[0])
@@ -200,7 +200,7 @@ void parse_expression(const std::vector<std::string>& expressions, std::vector<T
         // If more than one '*' found, output an error message
         else
         {
-            throw std::runtime_error("Invalid expression: " + expr + " - More than one '*' found.");
+            ModuleBase::WARNING_QUIT("split_expr", "Invalid expression: " + expr + " - More than one '*' found.");
         }
     }
 }

source/source_io/module_unk/berryphase.h

@@ -8,6 +8,7 @@
 #include "source_cell/klist.h"
 #include "source_psi/psi.h"
 #include "source_basis/module_ao/parallel_orbitals.h"
+#include "source_base/tool_quit.h"
 
 class berryphase
 {
@@ -69,7 +70,7 @@ class berryphase
                                   const ModulePW::PW_Basis_K* wfcpw,
                                   const K_Vectors& kv)
     {
-        throw std::logic_error("berry phase supports only multi-k");
+        ModuleBase::WARNING_QUIT("berryphase", "berry phase supports only multi-k");
     };
     void Macroscopic_polarization(const UnitCell& ucell,
                                   const int npwx,

source/source_io/module_wannier/to_wannier90.cpp

@@ -5,6 +5,7 @@
 #include "source_base/math_polyint.h"
 #include "source_base/math_sphbes.h"
 #include "source_base/math_ylmreal.h"
+#include "source_base/tool_quit.h"
 
 toWannier90::toWannier90()
 {
@@ -425,8 +426,7 @@ bool toWannier90::try_read_nnkp(const UnitCell& ucell, const K_Vectors& kv)
             }
             else
             {
-                throw std::runtime_error("numkpt_nnkp uninitialized in " + std::string(__FILE__)
-                                         + " line " + std::to_string(__LINE__));
+                ModuleBase::WARNING_QUIT("toWannier90::read_nnkp", "numkpt_nnkp uninitialized, nspin must be 1 or 2");
             }
 
             for (int ik = 0; ik < numkpt_nnkp; ik++)

source/source_io/module_wannier/to_wannier90_lcao.h

@@ -10,6 +10,7 @@
 #include "source_base/parallel_reduce.h"
 #include "source_base/sph_bessel_recursive.h"
 #include "source_base/timer.h"
+#include "source_base/tool_quit.h"
 #include "source_base/vector3.h"
 #include "source_base/ylm.h"
 #include "source_basis/module_ao/ORB_atomic_lm.h"
@@ -90,7 +91,7 @@ class toWannier90_LCAO : public toWannier90
                    const psi::Psi<double>& psi,
                    const Parallel_Orbitals* pv)
     {
-        throw std::logic_error("The wave function of toWannier90_LCAO_IN_PW is generally a std::complex<double> type.");
+        ModuleBase::WARNING_QUIT("toWannier90_LCAO", "The wave function of toWannier90_LCAO is generally a std::complex<double> type.");
     }
 
     void cal_Amn(const UnitCell& ucell, const K_Vectors& kv, const psi::Psi<std::complex<double>>& psi);

source/source_io/module_wannier/to_wannier90_lcao_in_pw.h

@@ -9,6 +9,7 @@
 #include "source_base/matrix3.h"
 #include "source_base/parallel_reduce.h"
 #include "source_base/timer.h"
+#include "source_base/tool_quit.h"
 #include "source_base/vector3.h"
 #include "source_base/ylm.h"
 #include "source_basis/module_ao/ORB_atomic_lm.h"
@@ -61,7 +62,7 @@ class toWannier90_LCAO_IN_PW : public toWannier90_PW
                    const psi::Psi<double>* psi,
                    const Parallel_Orbitals* pv)
     {
-        throw std::logic_error("The wave function of toWannier90_LCAO_IN_PW is generally a std::complex<double> type.");
+        ModuleBase::WARNING_QUIT("toWannier90_LCAO_IN_PW", "The wave function of toWannier90_LCAO_IN_PW is generally a std::complex<double> type.");
     }
 
   protected:

source/source_io/module_wannier/to_wannier90_pw.h

@@ -13,6 +13,7 @@
 #include "source_base/global_function.h"
 #include "source_base/global_variable.h"
 #include "source_base/matrix.h"
+#include "source_base/tool_quit.h"
 #include "source_base/matrix3.h"
 #include "source_cell/klist.h"
 #include "source_lcao/wavefunc_in_pw.h"
@@ -50,7 +51,7 @@ class toWannier90_PW : public toWannier90
       const psi::Psi<double>* psi
     )
     {
-      throw std::logic_error("The wave function of toWannier90_PW is generally a std::complex<double> type.");
+      ModuleBase::WARNING_QUIT("toWannier90_PW", "The wave function of toWannier90_PW is generally a std::complex<double> type.");
     }
 
     void cal_Amn(const psi::Psi<std::complex<double>>& psi_pw, const ModulePW::PW_Basis_K* wfcpw);