ABACUS develop
Atomic-orbital Based Ab-initio Computation at UStc
Loading...
Searching...
No Matches
make_gaunt_table-unittest.h
Go to the documentation of this file.
1//==========================================================
2// AUTHOR : Peize Lin
3// DATE : 2016-09-09
4//==========================================================
5
6#ifndef MAKE_GAUNT_TABLE_UNITTEST_H
7#define MAKE_GAUNT_TABLE_UNITTEST_H
8
10
11static void cout_MGT ( const ORB_gaunt_table & MGT, const int Lmax )
12{
13 for( int LA=0; LA<=Lmax; ++LA )
14 {
15 for( int LB=0; LB<=Lmax; ++LB )
16 {
17 for( int LAB=std::abs(LA-LB); LAB<=LA+LB; ++LAB )
18 {
19 for( int mA=0; mA!=2*LA+1; ++mA )
20 {
21 for( int mB=0; mB!=2*LB+1; ++mB )
22 {
23 for( int mAB=0; mAB!=2*LAB+1; ++mAB )
24 {
25 std::cout<<LA<<"\t"
26 <<LB<<"\t"
27 <<LAB<<"\t"
28 <<mA<<"\t"
29 <<mB<<"\t"
30 <<mAB<<"\t"
31 <<ORB_gaunt_table::Index_M(mA)<<"\t"
32 <<ORB_gaunt_table::Index_M(mB)<<"\t"
33 <<ORB_gaunt_table::Index_M(mAB)<<"\t"
34 <<MGT.Gaunt_Coefficients (
35 MGT.get_lm_index(LA,mA),
36 MGT.get_lm_index(LB,mB),
37 MGT.get_lm_index(LAB,mAB))<<std::endl;
38 }
39 }
40 }
41 }
42
43 }
44 }
45}
46#endif
Definition ORB_gaunt_table.h:9
static int get_lm_index(const int l, const int m)
Definition ORB_gaunt_table.cpp:157
static int Index_M(const int &m)
Definition ORB_gaunt_table.cpp:454
ModuleBase::realArray Gaunt_Coefficients
Definition ORB_gaunt_table.h:56