mathzone_add1.h

Go to the documentation of this file.

#ifndef MATHZONE_ADD1_H
#define MATHZONE_ADD1_H
 
#include <cassert>
#include <cmath>
 
namespace ModuleBase
{
 
/************************************************************************
LiaoChen add @ 2010/03/09 to add efficient functions in LCAO calculation
************************************************************************/
// Only used in source_basis/module_ao
class Mathzone_Add1
{
  public:
    Mathzone_Add1();
    ~Mathzone_Add1();
 
    static double dualfac(const int& l);
    static double factorial(const int& l);
    static void SplineD2(const double* rad,
                         const double* rad_f,
                         const int& mesh,
                         const double& yp1,
                         const double& ypn,
                         double* y2);
 
    static void Cubic_Spline_Interpolation(const double* const rad,
                                           const double* const rad_f,
                                           const double* const y2,
                                           const int& mesh,
                                           const double* const r,
                                           const int& rsize,
                                           double* const y,
                                           double* const dy);
 
    static double Uni_RadialF(const double* rad_f, const int& msh, const double& dr, const double& R);
 
    static void Uni_Deriv_Phi(const double* radf, const int& mesh, const double& dr, const int& nd, double* phind);
 
  private:
    const static int sph_lmax = 20;
    static double** c_ln_c;
    static double** c_ln_s;
};
 
} // namespace ModuleBase
 
#endif