-
Notifications
You must be signed in to change notification settings - Fork 0
/
example.log
34 lines (33 loc) · 1.59 KB
/
example.log
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
lambda=1.014
d =15.1
k_x =-0.68023412906043
r =(1.11,0.12)
eps =1e-15
-- normal sum --
qgf=sum_{l=-inf}^inf G(x+ld,y) exp(il k_x d)
|l|=100000 qgf=-1.40658671867873e-02 +6.45496053340515e-02 I
|l|=200000 qgf=-1.39812247031168e-02 +6.44780160455837e-02 I
|l|=300000 qgf=-1.39811184485477e-02 +6.45414749780242e-02 I
|l|=400000 qgf=-1.40040163809374e-02 +6.45239930518603e-02 I
|l|=500000 qgf=-1.40043768858218e-02 +6.45188362852444e-02 I
|l|=600000 qgf=-1.39844833174258e-02 +6.45059159183048e-02 I
|l|=700000 qgf=-1.39881483479676e-02 +6.45449502298631e-02 I
|l|=800000 qgf=-1.40202399173256e-02 +6.45102812409256e-02 I
|l|=900000 qgf=-1.39724855042906e-02 +6.45101290566453e-02 I
qgf =-1.39741869773027e-02 +6.45099399077704e-02 I, |l|=1000000
dqgf/dx=-3.97401252089164e-01 -1.35142565747076e-01 I
dqgf/dy=-5.00380051883253e-02 +3.25617528750411e-03 I
-- Ewald method --
qgf =-1.39953725950511e-02 +6.45210925705291e-02 I, erc=104
dqgf/dx=-4.51671901640698e-01 +3.53438797565453e-02 I
dqgf/dy=-5.00380051888735e-02 +3.25617528646067e-03 I
dqgf/dr=-4.54433564815398e-01 +3.54891140631323e-02 I
1/r*dqgf/dtheta=-1.20160899424966e-03 -5.61512916391293e-04 I
cylindrical coordinate ( d2hm_qpgf_d1_ew_cs() )
qgf =-1.39953725950512e-02 +6.45210925705290e-02 I, erc=108
dqgf/dr=-4.54433564815398e-01 +3.54891140631322e-02 I
1/r*dqgf/dtheta=-1.20160899424966e-03 -5.61512916391194e-04 I
-- integral representation --
qgf =-1.39953725950511e-02 +6.45210925705291e-02 I, erc=795 ,err=5.67643e-16
dqgf/dx=-4.51671901640697e-01 +3.53438797565451e-02 I
dqgf/dy=-5.00380051888736e-02 +3.25617528646068e-03 I