|
@@ -22,9 +22,9 @@ int main()
|
22
|
22
|
{
|
23
|
23
|
clock_t StartTime = clock();
|
24
|
24
|
|
25
|
|
- DP c_int = 2.0;
|
26
|
|
- DP L = 3.0;
|
27
|
|
- int N = 3;
|
|
25
|
+ DP c_int = 25.0;
|
|
26
|
+ DP L = 50.0;
|
|
27
|
+ int N = 10;
|
28
|
28
|
DP nbar_required = 1.0;
|
29
|
29
|
DP kBT = 4.0;
|
30
|
30
|
int Npts = 4*N;
|
|
@@ -34,9 +34,9 @@ int main()
|
34
|
34
|
|
35
|
35
|
if (true) { // State-by-state checks
|
36
|
36
|
|
37
|
|
- DP c_int = 4.0;
|
38
|
|
- DP L = 16.0;
|
39
|
|
- int N = 16;
|
|
37
|
+ DP c_int = 25.0;
|
|
38
|
+ DP L = 50.0;
|
|
39
|
+ int N = 10;
|
40
|
40
|
|
41
|
41
|
LiebLin_Bethe_State gstate (c_int, L, N);
|
42
|
42
|
gstate.Compute_All(true);
|
|
@@ -71,9 +71,9 @@ int main()
|
71
|
71
|
// << exp(real(ln_Density_ME(gstate, estate))) << "\t" << exp(real(ln_Density_ME(estate, gstate))) << endl;
|
72
|
72
|
|
73
|
73
|
//Field operator matrix elements:
|
74
|
|
- cout << "omega\tiK\t< estate | Psi | gstate > matrix element:" << endl;
|
75
|
|
- cout << setprecision(16) << estate.E - gstate.E << "\t" << estate.iK - gstate.iK << "\t"
|
76
|
|
- << exp(ln_Psi_ME(estate, gstate)) << endl;
|
|
74
|
+ // cout << "omega\tiK\t< estate | Psi | gstate > matrix element:" << endl;
|
|
75
|
+ // cout << setprecision(16) << estate.E - gstate.E << "\t" << estate.iK - gstate.iK << "\t"
|
|
76
|
+ // << exp(ln_Psi_ME(estate, gstate)) << endl;
|
77
|
77
|
|
78
|
78
|
//LiebLin_Bethe_State flipstate = estate;
|
79
|
79
|
//flipstate.Parity_Flip();
|