/
m_scheme_test.cpp
59 lines (43 loc) · 1.53 KB
/
m_scheme_test.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
/****************************************************************
m_scheme_test.cpp
Patrick J. Fasano
University of Notre Dame
****************************************************************/
#include <iostream>
#include "m_scheme.h"
#include "nlj_orbital.h"
////////////////////////////////////////////////////////////////
// test code
////////////////////////////////////////////////////////////////
void TestStatesNmax()
{
////////////////////////////////////////////////////////////////
// orbital tests (should come after suborbital tests are passed)
////////////////////////////////////////////////////////////////
std::cout << "Orbitals -- Nmax scheme" << std::endl;
// set up space
std::cout << "Space" << std::endl;
int Nmax = 4;
basis::OrbitalSpacePN orbital_space(Nmax);
basis::SingleParticleSpacePN space(orbital_space);
std::cout << space.DebugStr();
// check subspaces
std::cout << "Subspaces" << std::endl;
for (std::size_t subspace_index=0; subspace_index<space.size(); ++subspace_index)
{
const basis::SingleParticleSubspacePN& subspace = space.GetSubspace(subspace_index);
std::cout << " index " << subspace_index
<< " species " << int(subspace.orbital_species())
<< std::endl;
std::cout << subspace.DebugStr();
}
}
////////////////////////////////////////////////////////////////
// main
////////////////////////////////////////////////////////////////
int main(int argc, char **argv)
{
TestStatesNmax();
// termination
return 0;
}