Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1,022 changes: 522 additions & 500 deletions tools/01_NAO_generation/SIAB/src_spillage/Calculate_C4.cpp

Large diffs are not rendered by default.

65 changes: 36 additions & 29 deletions tools/01_NAO_generation/SIAB/src_spillage/Calculate_C4.h
Original file line number Diff line number Diff line change
Expand Up @@ -5,35 +5,42 @@

class Calculate_C4
{
public:

Calculate_C4();
~Calculate_C4();

void init( ifstream &ifs, const double &ecut, const double &rcut, const double &tolerence, const int &lmax );
void Test( ifstream &ifs, const double &ecut, const double &tolerence);
void plot( ifstream &ifs );

// Peize Lin add accept_rate 2015-12-24
static void norm_ic( realArray &C4, realArray &accept_rate, const int &ntype, const int &lmax,
const int &nmax, const int &enumber, const double &tolerence,
const double &rcut, const double &dr);

static void Find_Eigenvalues(const int &Enumber, matrix &en,
const int &lmax, const double &rcut, const double &tolerence);

private:

void Test_Orthogonal(const double &rcut, const double &dr,
const int &L, const matrix &en, const int &Enumber, const int &test_ei);

void key_solver(const string &name, const int &L,
const matrix &en, const int &Enumber, double *C4);

static int test;



public:
Calculate_C4();
~Calculate_C4();

void init(ifstream& ifs, const double& ecut, const double& rcut, const double& tolerence, const int& lmax);
void Test(ifstream& ifs, const double& ecut, const double& tolerence);
void plot(ifstream& ifs);

// Peize Lin add accept_rate 2015-12-24
static void norm_ic(realArray& C4,
realArray& accept_rate,
const int& ntype,
const int& lmax,
const int& nmax,
const int& enumber,
const double& tolerence,
const double& rcut,
const double& dr);

static void Find_Eigenvalues(const int& Enumber,
matrix& en,
const int& lmax,
const double& rcut,
const double& tolerence);

private:
void Test_Orthogonal(const double& rcut,
const double& dr,
const int& L,
const matrix& en,
const int& Enumber,
const int& test_ei);

void key_solver(const string& name, const int& L, const matrix& en, const int& Enumber, double* C4);

static int test;
};

#endif
Loading
Loading