forked from JeffersonLab/simc_gfortran
-
Notifications
You must be signed in to change notification settings - Fork 6
/
h2.npst
142 lines (109 loc) · 8.01 KB
/
h2.npst
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
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
********************************************************************************
Theory IBM RS/6000-397 AIX 4.3.2 XL Fortran 5.1 10:30:16 18-Nov-99
2H
Hamiltonian H = Ti + Vij
Vij : Argonne v18 + Vem
r(1)= .7812500E-02, r(25600)= 200.00000; dr= .244140625E-03
|psi>=|psiD>
|psi_move>=|psiD>
alphas
1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
ebind eta dstate qm
2.2246 .0250 .0576 .2696
edeut ti vij v6 vb vqbb vem vope
-2.2249 19.8100 -22.0349 -21.0321 -.0828 -.9376 .0176 -21.2830
nucleon momentum distribution calculated for 51 q values
from .00 to 5.00
using 250 gauss points in r out to 50.00 fm
random number seed 9724109.
initial position:
particle 1 : ( .00, .00, .00)
particle 2 : ( 3.20, .00, .00)
wave function set-up time = .480 seconds
random walk step size = 3.2000
weight = 1.0000 + .0000 * sum(ri**2)
of 1000 initial moves 547 were accepted
randomized position:
particle 1 : ( -.84, .57, .30)
particle 2 : ( .84, -.57, -.30)
antisymmetry test
n1 n2 1 + PxPsPt
1 2 .33307E-15
initial move time = .010 seconds
metropolis walk attempts 10 steps per sample
********************************************************************************
total average after 50000 points
dst = .0575 ( .0002)
mu = .0000 ( .0000) qm = .0000 ( .0000) hex = .0000 ( .0000)
rmsp = 1.9799 ( .0074) rmsn = 1.9799 ( .0074)
norm = 1.0000 (.00000E+00)
.00000E+00 (.00000E+00)
.00000E+00 (.00000E+00)
k(fm-1) rhond(k) rhonu(k) rhopd(k) rhopu(k)
.00 .8591E-01 (.8940E-01) 3243. (37.73) .8591E-01 (.8940E-01) 3243. (37.73)
.10 .7506E-01 (.5278E-01) 2260. (19.95) .7506E-01 (.5278E-01) 2260. (19.95)
.20 .2710 (.3488E-01) 1004. (7.319) .2710 (.3488E-01) 1004. (7.319)
.30 .5231 (.2582E-01) 404.5 (2.826) .5231 (.2582E-01) 404.5 (2.826)
.40 .6623 (.2013E-01) 169.2 (1.230) .6623 (.2013E-01) 169.2 (1.230)
.50 .7124 (.1578E-01) 75.84 (.5843) .7124 (.1578E-01) 75.84 (.5843)
.60 .6902 (.1244E-01) 36.60 (.3038) .6902 (.1244E-01) 36.60 (.3038)
.70 .6368 (.9898E-02) 18.60 (.1682) .6368 (.9898E-02) 18.60 (.1682)
.80 .5766 (.7917E-02) 9.735 (.9822E-01) .5766 (.7917E-02) 9.735 (.9822E-01)
.90 .5029 (.6381E-02) 5.278 (.6143E-01) .5029 (.6381E-02) 5.278 (.6143E-01)
1.00 .4321 (.5227E-02) 2.944 (.4261E-01) .4321 (.5227E-02) 2.944 (.4261E-01)
1.10 .3698 (.4345E-02) 1.662 (.3417E-01) .3698 (.4345E-02) 1.662 (.3417E-01)
1.20 .3139 (.3651E-02) .9969 (.2743E-01) .3139 (.3651E-02) .9969 (.2743E-01)
1.30 .2623 (.3101E-02) .5416 (.2425E-01) .2623 (.3101E-02) .5416 (.2425E-01)
1.40 .2201 (.2654E-02) .3335 (.2153E-01) .2201 (.2654E-02) .3335 (.2153E-01)
1.50 .1847 (.2264E-02) .2067 (.1787E-01) .1847 (.2264E-02) .2067 (.1787E-01)
1.60 .1551 (.1952E-02) .1543 (.1592E-01) .1551 (.1952E-02) .1543 (.1592E-01)
1.70 .1285 (.1673E-02) .5921E-01 (.1332E-01) .1285 (.1673E-02) .5921E-01 (.1332E-01)
1.80 .1065 (.1431E-02) .1895E-01 (.1132E-01) .1065 (.1431E-02) .1895E-01 (.1132E-01)
1.90 .8876E-01 (.1221E-02) .1011E-01 (.9264E-02) .8876E-01 (.1221E-02) .1011E-01 (.9264E-02)
2.00 .7437E-01 (.1039E-02) .2325E-01 (.7879E-02) .7437E-01 (.1039E-02) .2325E-01 (.7879E-02)
2.10 .6193E-01 (.8780E-03) .1630E-01 (.6251E-02) .6193E-01 (.8780E-03) .1630E-01 (.6251E-02)
2.20 .5159E-01 (.7417E-03) .1719E-01 (.5215E-02) .5159E-01 (.7417E-03) .1719E-01 (.5215E-02)
2.30 .4288E-01 (.6231E-03) .1725E-01 (.4094E-02) .4288E-01 (.6231E-03) .1725E-01 (.4094E-02)
2.40 .3557E-01 (.5230E-03) .1372E-01 (.3380E-02) .3557E-01 (.5230E-03) .1372E-01 (.3380E-02)
2.50 .2950E-01 (.4389E-03) .1155E-01 (.2609E-02) .2950E-01 (.4389E-03) .1155E-01 (.2609E-02)
2.60 .2449E-01 (.3684E-03) .1229E-01 (.2120E-02) .2449E-01 (.3684E-03) .1229E-01 (.2120E-02)
2.70 .2030E-01 (.3099E-03) .1143E-01 (.1610E-02) .2030E-01 (.3099E-03) .1143E-01 (.1610E-02)
2.80 .1683E-01 (.2618E-03) .1161E-01 (.1234E-02) .1683E-01 (.2618E-03) .1161E-01 (.1234E-02)
2.90 .1393E-01 (.2226E-03) .1126E-01 (.9635E-03) .1393E-01 (.2226E-03) .1126E-01 (.9635E-03)
3.00 .1146E-01 (.1902E-03) .7793E-02 (.7019E-03) .1146E-01 (.1902E-03) .7793E-02 (.7019E-03)
3.10 .9453E-02 (.1639E-03) .7517E-02 (.5317E-03) .9453E-02 (.1639E-03) .7517E-02 (.5317E-03)
3.20 .7802E-02 (.1422E-03) .7319E-02 (.4061E-03) .7802E-02 (.1422E-03) .7319E-02 (.4061E-03)
3.30 .6423E-02 (.1240E-03) .6229E-02 (.3358E-03) .6423E-02 (.1240E-03) .6229E-02 (.3358E-03)
3.40 .5282E-02 (.1086E-03) .5664E-02 (.2987E-03) .5282E-02 (.1086E-03) .5664E-02 (.2987E-03)
3.50 .4339E-02 (.9530E-04) .4879E-02 (.2811E-03) .4339E-02 (.9530E-04) .4879E-02 (.2811E-03)
3.60 .3566E-02 (.8369E-04) .4531E-02 (.2829E-03) .3566E-02 (.8369E-04) .4531E-02 (.2829E-03)
3.70 .2921E-02 (.7340E-04) .3941E-02 (.2783E-03) .2921E-02 (.7340E-04) .3941E-02 (.2783E-03)
3.80 .2392E-02 (.6422E-04) .3721E-02 (.2689E-03) .2392E-02 (.6422E-04) .3721E-02 (.2689E-03)
3.90 .1951E-02 (.5596E-04) .2916E-02 (.2591E-03) .1951E-02 (.5596E-04) .2916E-02 (.2591E-03)
4.00 .1591E-02 (.4852E-04) .2654E-02 (.2413E-03) .1591E-02 (.4852E-04) .2654E-02 (.2413E-03)
4.10 .1290E-02 (.4184E-04) .2076E-02 (.2217E-03) .1290E-02 (.4184E-04) .2076E-02 (.2217E-03)
4.20 .1052E-02 (.3584E-04) .2082E-02 (.2015E-03) .1052E-02 (.3584E-04) .2082E-02 (.2015E-03)
4.30 .8508E-03 (.3047E-04) .1642E-02 (.1774E-03) .8508E-03 (.3047E-04) .1642E-02 (.1774E-03)
4.40 .6875E-03 (.2571E-04) .1482E-02 (.1555E-03) .6875E-03 (.2571E-04) .1482E-02 (.1555E-03)
4.50 .5500E-03 (.2151E-04) .1074E-02 (.1339E-03) .5500E-03 (.2151E-04) .1074E-02 (.1339E-03)
4.60 .4429E-03 (.1782E-04) .1029E-02 (.1145E-03) .4429E-03 (.1782E-04) .1029E-02 (.1145E-03)
4.70 .3548E-03 (.1462E-04) .8868E-03 (.9470E-04) .3548E-03 (.1462E-04) .8868E-03 (.9470E-04)
4.80 .2841E-03 (.1187E-04) .8142E-03 (.8080E-04) .2841E-03 (.1187E-04) .8142E-03 (.8080E-04)
4.90 .2245E-03 (.9507E-05) .5828E-03 (.6442E-04) .2245E-03 (.9507E-05) .5828E-03 (.6442E-04)
5.00 .1765E-03 (.7519E-05) .4864E-03 (.5486E-04) .1765E-03 (.7519E-05) .4864E-03 (.5486E-04)
nd nu pd pu
.043 .957 .043 .957
mok = 273266 nwarn = 0 0 0
********************************************************************************
total average for 5 blocks
dst = .0575 ( .0002)
mu = .0000 ( .0000) qm = .0000 ( .0000) hex = .0000 ( .0000)
rmsp = 1.9799 ( .0094) rmsn = 1.9799 ( .0094)
norm = 1.0000 (.00000E+00)
.00000E+00 (.00000E+00)
.00000E+00 (.00000E+00)
********************************************************************************
energy calculation time = 839.700 seconds
total flops = 23687016795.
energy calculation speed = 28.209 Mflops
total job time = 840.190 seconds