git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@15149 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp 2016-06-07 18:17:05 +00:00
parent b5abf9342b
commit 0129ff9696
3 changed files with 873 additions and 0 deletions

View File

@ -0,0 +1,225 @@
LAMMPS (7 Jun 2016)
# sample LAMMPS input script for thermal conductivity of liquid LJ
# use fix ehex to add/subtract energy from 2 regions
# settings
variable x equal 10
variable y equal 10
variable z equal 20
variable rho equal 0.6
variable t equal 1.35
variable rc equal 2.5
#variable rho equal 0.85
#variable t equal 0.7
#variable rc equal 3.0
# setup problem
units lj
atom_style atomic
lattice fcc ${rho}
lattice fcc 0.6
Lattice spacing in x,y,z = 1.88207 1.88207 1.88207
region box block 0 $x 0 $y 0 $z
region box block 0 10 0 $y 0 $z
region box block 0 10 0 10 0 $z
region box block 0 10 0 10 0 20
create_box 1 box
Created orthogonal box = (0 0 0) to (18.8207 18.8207 37.6414)
2 by 1 by 4 MPI processor grid
create_atoms 1 box
Created 8000 atoms
mass 1 1.0
velocity all create $t 87287
velocity all create 1.35 87287
pair_style lj/cut ${rc}
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
neighbor 0.3 bin
neigh_modify delay 0 every 1
# heat layers
region hot block INF INF INF INF 0 1
region cold block INF INF INF INF 10 11
compute Thot all temp/region hot
compute Tcold all temp/region cold
# 1st equilibration run
fix 1 all nvt temp $t $t 0.5
fix 1 all nvt temp 1.35 $t 0.5
fix 1 all nvt temp 1.35 1.35 0.5
thermo 100
run 1000
Neighbor list info ...
1 neighbor list requests
update every 1 steps, delay 0 steps, check yes
max neighbors/atom: 2000, page size: 100000
master list distance cutoff = 2.8
ghost atom cutoff = 2.8
binsize = 1.4 -> bins = 14 14 27
Memory usage per processor = 2.55761 Mbytes
Step Temp E_pair E_mol TotEng Press
0 1.35 -4.1241917 0 -2.0994448 -3.1961612
100 1.1819832 -3.7640881 0 -1.991335 0.53985757
200 1.2578365 -3.7395333 0 -1.8530144 0.69591862
300 1.3282971 -3.7215427 0 -1.7293461 0.79036065
400 1.3714367 -3.7043826 0 -1.6474847 0.85873226
500 1.3590952 -3.6707735 0 -1.6323855 0.99602024
600 1.3575117 -3.7118244 0 -1.6758114 0.81454305
700 1.3284444 -3.7075488 0 -1.7151313 0.81136596
800 1.3419995 -3.7155648 0 -1.7028172 0.82925676
900 1.3562214 -3.6965609 0 -1.6624831 0.88908117
1000 1.3732017 -3.7100044 0 -1.6504594 0.83982701
Loop time of 0.878884 on 8 procs for 1000 steps with 8000 atoms
Performance: 491532.534 tau/day, 1137.807 timesteps/s
99.3% CPU use with 8 MPI tasks x no OpenMP threads
MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 0.48108 | 0.4925 | 0.50422 | 1.0 | 56.04
Neigh | 0.23355 | 0.23565 | 0.23888 | 0.4 | 26.81
Comm | 0.088565 | 0.10044 | 0.11391 | 2.7 | 11.43
Output | 0.00029421 | 0.00031444 | 0.00041795 | 0.2 | 0.04
Modify | 0.027427 | 0.032625 | 0.035225 | 1.6 | 3.71
Other | | 0.01735 | | | 1.97
Nlocal: 1000 ave 1020 max 982 min
Histogram: 1 0 2 1 0 1 1 1 0 1
Nghost: 2299.5 ave 2331 max 2268 min
Histogram: 1 1 1 1 0 0 0 3 0 1
Neighs: 27122 ave 28382 max 26337 min
Histogram: 2 0 2 1 1 0 0 1 0 1
Total # of neighbors = 216976
Ave neighs/atom = 27.122
Neighbor list builds = 162
Dangerous builds = 0
velocity all scale $t
velocity all scale 1.35
unfix 1
# 2nd equilibration run
fix 1 all nve
fix hot all ehex 1 100.0 region hot
fix cold all ehex 1 -100.0 region cold
thermo_style custom step temp c_Thot c_Tcold
thermo 1000
run 10000
Memory usage per processor = 2.80761 Mbytes
Step Temp Thot Tcold
1000 1.35 1.431295 1.2955644
2000 1.3537291 1.6418772 1.1875127
3000 1.3615152 1.6451299 1.1769094
4000 1.3612129 1.5281727 1.2022419
5000 1.3552182 1.6672955 1.2212864
6000 1.3643442 1.6072213 1.2390567
7000 1.3665773 1.6909819 1.1466611
8000 1.375741 1.6144274 1.1691231
9000 1.3701136 1.8238424 1.136342
10000 1.3563004 1.8059065 1.1547129
11000 1.3794051 1.692299 1.0515688
Loop time of 10.292 on 8 procs for 10000 steps with 8000 atoms
Performance: 419743.444 tau/day, 971.628 timesteps/s
99.5% CPU use with 8 MPI tasks x no OpenMP threads
MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 4.1797 | 4.902 | 5.5136 | 25.3 | 47.63
Neigh | 2.1913 | 2.4616 | 2.6938 | 14.2 | 23.92
Comm | 0.87831 | 1.64 | 2.611 | 56.5 | 15.94
Output | 0.00043011 | 0.00063166 | 0.00075626 | 0.5 | 0.01
Modify | 0.78531 | 0.92487 | 1.0647 | 10.6 | 8.99
Other | | 0.3628 | | | 3.53
Nlocal: 1000 ave 1105 max 883 min
Histogram: 1 1 2 0 0 0 0 0 2 2
Nghost: 2319.38 ave 2502 max 2114 min
Histogram: 1 3 0 0 0 0 0 0 0 4
Neighs: 27387.9 ave 32453 max 21803 min
Histogram: 2 2 0 0 0 0 0 0 0 4
Total # of neighbors = 219103
Ave neighs/atom = 27.3879
Neighbor list builds = 1696
Dangerous builds = 0
# thermal conductivity calculation
compute ke all ke/atom
variable temp atom c_ke/1.5
compute layers all chunk/atom bin/1d z lower 0.05 units reduced
fix 2 all ave/chunk 10 100 1000 layers v_temp file profile.heat
variable tdiff equal f_2[11][3]-f_2[1][3]
fix ave all ave/time 1 1 1000 v_tdiff ave running start 13000
thermo_style custom step temp c_Thot c_Tcold v_tdiff f_ave
run 20000
Memory usage per processor = 3.0578 Mbytes
Step Temp Thot Tcold tdiff ave
11000 1.3794051 1.6903393 1.0515688 0 0
12000 1.3799777 1.8004888 1.1032219 -0.63860014 0
13000 1.3733605 1.7823094 1.0553582 -0.65827891 -0.65827891
14000 1.3749743 1.7852256 1.1674016 -0.68463005 -0.67145448
15000 1.3863795 1.8538127 1.0056247 -0.73524813 -0.69271903
16000 1.3731955 1.7518546 1.0741458 -0.74810775 -0.70656621
17000 1.3771856 1.9016869 1.0090502 -0.73999567 -0.7132521
18000 1.3766032 1.7616195 1.1142155 -0.73769104 -0.71732526
19000 1.3815934 1.7791247 1.1406987 -0.73617832 -0.72001855
20000 1.3725543 1.8637436 1.0799364 -0.73435569 -0.7218107
21000 1.3817369 1.8808771 1.0642524 -0.76702329 -0.72683432
22000 1.3968704 1.840287 1.072304 -0.82496419 -0.7366473
23000 1.3895558 1.9427293 1.0766665 -0.75363908 -0.73819201
24000 1.3900493 1.9883976 1.1081017 -0.86394774 -0.74867166
25000 1.3838912 1.8853041 1.0795751 -0.83043902 -0.75496145
26000 1.3912105 1.9330259 1.1070335 -0.79880182 -0.75809291
27000 1.3891151 1.8548451 1.0676153 -0.81856523 -0.7621244
28000 1.3942624 1.9796706 1.1251407 -0.81762456 -0.76559316
29000 1.3819302 1.8619138 1.0495292 -0.78627491 -0.76680973
30000 1.3968366 1.883107 1.1004588 -0.83902548 -0.77082172
31000 1.3822489 1.8220413 1.0322271 -0.7550338 -0.76999077
Loop time of 22.3319 on 8 procs for 20000 steps with 8000 atoms
Performance: 386890.255 tau/day, 895.579 timesteps/s
99.3% CPU use with 8 MPI tasks x no OpenMP threads
MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 7.4818 | 9.8886 | 12.008 | 66.1 | 44.28
Neigh | 4.1157 | 5.0693 | 5.9657 | 38.0 | 22.70
Comm | 1.8343 | 4.6586 | 7.8747 | 129.4 | 20.86
Output | 0.00095606 | 0.0017744 | 0.0019104 | 0.7 | 0.01
Modify | 1.6812 | 1.9319 | 2.2204 | 13.5 | 8.65
Other | | 0.7817 | | | 3.50
Nlocal: 1000 ave 1121 max 857 min
Histogram: 2 0 1 1 0 0 0 0 1 3
Nghost: 2299.75 ave 2541 max 2067 min
Histogram: 3 1 0 0 0 0 0 0 2 2
Neighs: 27487.2 ave 33361 max 20651 min
Histogram: 2 1 1 0 0 0 0 0 0 4
Total # of neighbors = 219898
Ave neighs/atom = 27.4872
Neighbor list builds = 3474
Dangerous builds = 0
Total wall time: 0:00:33

View File

@ -0,0 +1,225 @@
LAMMPS (7 Jun 2016)
# sample LAMMPS input script for thermal conductivity of liquid LJ
# use fix heat to add/subtract energy from 2 regions
# settings
variable x equal 10
variable y equal 10
variable z equal 20
variable rho equal 0.6
variable t equal 1.35
variable rc equal 2.5
#variable rho equal 0.85
#variable t equal 0.7
#variable rc equal 3.0
# setup problem
units lj
atom_style atomic
lattice fcc ${rho}
lattice fcc 0.6
Lattice spacing in x,y,z = 1.88207 1.88207 1.88207
region box block 0 $x 0 $y 0 $z
region box block 0 10 0 $y 0 $z
region box block 0 10 0 10 0 $z
region box block 0 10 0 10 0 20
create_box 1 box
Created orthogonal box = (0 0 0) to (18.8207 18.8207 37.6414)
2 by 1 by 4 MPI processor grid
create_atoms 1 box
Created 8000 atoms
mass 1 1.0
velocity all create $t 87287
velocity all create 1.35 87287
pair_style lj/cut ${rc}
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
neighbor 0.3 bin
neigh_modify delay 0 every 1
# heat layers
region hot block INF INF INF INF 0 1
region cold block INF INF INF INF 10 11
compute Thot all temp/region hot
compute Tcold all temp/region cold
# 1st equilibration run
fix 1 all nvt temp $t $t 0.5
fix 1 all nvt temp 1.35 $t 0.5
fix 1 all nvt temp 1.35 1.35 0.5
thermo 100
run 1000
Neighbor list info ...
1 neighbor list requests
update every 1 steps, delay 0 steps, check yes
max neighbors/atom: 2000, page size: 100000
master list distance cutoff = 2.8
ghost atom cutoff = 2.8
binsize = 1.4 -> bins = 14 14 27
Memory usage per processor = 2.55761 Mbytes
Step Temp E_pair E_mol TotEng Press
0 1.35 -4.1241917 0 -2.0994448 -3.1961612
100 1.1819832 -3.7640881 0 -1.991335 0.53985757
200 1.2578365 -3.7395333 0 -1.8530144 0.69591862
300 1.3282971 -3.7215427 0 -1.7293461 0.79036065
400 1.3714367 -3.7043826 0 -1.6474847 0.85873226
500 1.3590952 -3.6707735 0 -1.6323855 0.99602024
600 1.3575117 -3.7118244 0 -1.6758114 0.81454305
700 1.3284444 -3.7075488 0 -1.7151313 0.81136596
800 1.3419995 -3.7155648 0 -1.7028172 0.82925676
900 1.3562214 -3.6965609 0 -1.6624831 0.88908117
1000 1.3732017 -3.7100044 0 -1.6504594 0.83982701
Loop time of 0.889036 on 8 procs for 1000 steps with 8000 atoms
Performance: 485919.613 tau/day, 1124.814 timesteps/s
99.4% CPU use with 8 MPI tasks x no OpenMP threads
MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 0.47867 | 0.4928 | 0.5006 | 1.0 | 55.43
Neigh | 0.23302 | 0.23557 | 0.23833 | 0.3 | 26.50
Comm | 0.1001 | 0.11264 | 0.12817 | 2.6 | 12.67
Output | 0.00032806 | 0.00034043 | 0.00040531 | 0.1 | 0.04
Modify | 0.028788 | 0.031511 | 0.033365 | 0.9 | 3.54
Other | | 0.01618 | | | 1.82
Nlocal: 1000 ave 1020 max 982 min
Histogram: 1 0 2 1 0 1 1 1 0 1
Nghost: 2299.5 ave 2331 max 2268 min
Histogram: 1 1 1 1 0 0 0 3 0 1
Neighs: 27122 ave 28382 max 26337 min
Histogram: 2 0 2 1 1 0 0 1 0 1
Total # of neighbors = 216976
Ave neighs/atom = 27.122
Neighbor list builds = 162
Dangerous builds = 0
velocity all scale $t
velocity all scale 1.35
unfix 1
# 2nd equilibration run
fix 1 all nve
fix hot all heat 1 100.0 region hot
fix cold all heat 1 -100.0 region cold
thermo_style custom step temp c_Thot c_Tcold
thermo 1000
run 10000
Memory usage per processor = 2.55761 Mbytes
Step Temp Thot Tcold
1000 1.35 1.431295 1.2955644
2000 1.3518468 1.5562602 1.154905
3000 1.3477229 1.5890075 1.2395414
4000 1.3487175 1.5491615 1.2019696
5000 1.3594394 1.5780597 1.1824492
6000 1.3583923 1.541735 1.1675586
7000 1.3700321 1.6735877 1.1279114
8000 1.3631993 1.6367675 1.0697225
9000 1.3739201 1.6846211 1.1138829
10000 1.3751455 1.8039471 1.1500399
11000 1.3716416 1.833336 1.1267278
Loop time of 10.488 on 8 procs for 10000 steps with 8000 atoms
Performance: 411900.216 tau/day, 953.473 timesteps/s
99.5% CPU use with 8 MPI tasks x no OpenMP threads
MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 4.138 | 4.9118 | 5.6168 | 26.2 | 46.83
Neigh | 2.1715 | 2.4541 | 2.7147 | 14.7 | 23.40
Comm | 0.94395 | 1.8449 | 2.9212 | 55.3 | 17.59
Output | 0.00043416 | 0.00046757 | 0.00054884 | 0.2 | 0.00
Modify | 1.0234 | 1.0946 | 1.2539 | 8.7 | 10.44
Other | | 0.1821 | | | 1.74
Nlocal: 1000 ave 1131 max 878 min
Histogram: 3 1 0 0 0 0 0 1 1 2
Nghost: 2312.88 ave 2525 max 2114 min
Histogram: 2 2 0 0 0 0 0 1 1 2
Neighs: 27457 ave 33797 max 21031 min
Histogram: 3 1 0 0 0 0 0 0 1 3
Total # of neighbors = 219656
Ave neighs/atom = 27.457
Neighbor list builds = 1691
Dangerous builds = 0
# thermal conductivity calculation
compute ke all ke/atom
variable temp atom c_ke/1.5
compute layers all chunk/atom bin/1d z lower 0.05 units reduced
fix 2 all ave/chunk 10 100 1000 layers v_temp file profile.heat
variable tdiff equal f_2[11][3]-f_2[1][3]
fix ave all ave/time 1 1 1000 v_tdiff ave running start 13000
thermo_style custom step temp c_Thot c_Tcold v_tdiff f_ave
run 20000
Memory usage per processor = 2.8078 Mbytes
Step Temp Thot Tcold tdiff ave
11000 1.3716416 1.833336 1.1267278 0 0
12000 1.3703433 1.7829467 1.1194444 -0.66044316 0
13000 1.3686734 1.8334366 1.1193477 -0.71431978 -0.71431978
14000 1.3856987 1.8048077 1.1052708 -0.73112558 -0.72272268
15000 1.3820117 1.7460559 1.110725 -0.72927647 -0.72490728
16000 1.3911309 1.923603 1.1161499 -0.77407515 -0.73719925
17000 1.3841301 1.7276486 1.0283807 -0.77278638 -0.74431667
18000 1.3888918 1.7944951 1.0532944 -0.75665895 -0.74637372
19000 1.3876032 1.838174 1.059715 -0.71342263 -0.74166642
20000 1.3724644 1.8297128 1.1439176 -0.77352223 -0.7456484
21000 1.3798921 1.7968403 1.0288381 -0.70077132 -0.74066206
22000 1.3763952 1.8202225 1.0658157 -0.75629111 -0.74222496
23000 1.3911378 1.8691478 1.018589 -0.76094865 -0.74392711
24000 1.3867754 1.7826523 1.09347 -0.80367344 -0.74890597
25000 1.385877 1.9029313 1.0815131 -0.73559505 -0.74788206
26000 1.3791773 1.8904022 1.0151678 -0.7729123 -0.74966993
27000 1.3800063 1.729283 1.127594 -0.71473941 -0.74734123
28000 1.3757197 1.7823772 1.084523 -0.73849831 -0.74678855
29000 1.3777555 1.8287284 1.0715132 -0.70375514 -0.74425717
30000 1.3821118 1.7382856 1.1078333 -0.79892499 -0.74729427
31000 1.3870476 1.8410063 1.1235958 -0.76218423 -0.74807795
Loop time of 22.3813 on 8 procs for 20000 steps with 8000 atoms
Performance: 386036.860 tau/day, 893.604 timesteps/s
99.4% CPU use with 8 MPI tasks x no OpenMP threads
MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 7.4578 | 9.8912 | 11.96 | 64.5 | 44.19
Neigh | 4.0721 | 5.0195 | 5.8898 | 36.9 | 22.43
Comm | 1.9554 | 4.7393 | 7.8917 | 125.4 | 21.18
Output | 0.00099301 | 0.0017841 | 0.0019252 | 0.7 | 0.01
Modify | 2.2325 | 2.321 | 2.4866 | 6.5 | 10.37
Other | | 0.4085 | | | 1.83
Nlocal: 1000 ave 1134 max 850 min
Histogram: 2 1 0 1 0 0 0 1 0 3
Nghost: 2307.75 ave 2561 max 2083 min
Histogram: 4 0 0 0 0 0 0 0 2 2
Neighs: 27561.1 ave 34071 max 19891 min
Histogram: 2 1 1 0 0 0 0 0 1 3
Total # of neighbors = 220489
Ave neighs/atom = 27.5611
Neighbor list builds = 3442
Dangerous builds = 0
Total wall time: 0:00:33

View File

@ -0,0 +1,423 @@
# Chunk-averaged data for fix 2 and group all
# Timestep Number-of-chunks Total-count
# Chunk Coord1 Ncount v_temp
12000 20 8000
1 0.025 326.75 1.76011
2 0.075 330.61 1.63635
3 0.125 350.4 1.54636
4 0.175 366.28 1.52687
5 0.225 387.83 1.42461
6 0.275 401.89 1.3514
7 0.325 419.9 1.32411
8 0.375 432.38 1.26508
9 0.425 442.47 1.21776
10 0.475 460.36 1.13062
11 0.525 473.38 1.09967
12 0.575 467.33 1.14472
13 0.625 443.61 1.19478
14 0.675 439.36 1.26366
15 0.725 425.36 1.32373
16 0.775 400.31 1.41125
17 0.825 366.52 1.53434
18 0.875 362.38 1.58774
19 0.925 356.68 1.58572
20 0.975 346.2 1.63678
13000 20 8000
1 0.025 318.55 1.78446
2 0.075 339.57 1.6873
3 0.125 352.13 1.60227
4 0.175 371.22 1.50409
5 0.225 391.7 1.43235
6 0.275 403.96 1.37436
7 0.325 411.42 1.37362
8 0.375 433.7 1.27575
9 0.425 454.72 1.17224
10 0.475 465.31 1.11717
11 0.525 469.77 1.07014
12 0.575 463.38 1.1317
13 0.625 450.42 1.19345
14 0.675 431.46 1.2564
15 0.725 416.59 1.30352
16 0.775 401.69 1.36855
17 0.825 375.83 1.50988
18 0.875 364.35 1.54229
19 0.925 348.01 1.57123
20 0.975 336.22 1.68615
14000 20 8000
1 0.025 319.25 1.81747
2 0.075 330.82 1.71274
3 0.125 347.34 1.60095
4 0.175 376.89 1.49722
5 0.225 388.07 1.42925
6 0.275 403.57 1.37239
7 0.325 413.7 1.34564
8 0.375 432.5 1.24231
9 0.425 454.45 1.18194
10 0.475 469.33 1.13752
11 0.525 469.94 1.08635
12 0.575 469.52 1.13378
13 0.625 451.41 1.19001
14 0.675 429.75 1.27268
15 0.725 420.06 1.32601
16 0.775 407.93 1.37822
17 0.825 379.47 1.50346
18 0.875 360.53 1.5774
19 0.925 345.93 1.6271
20 0.975 329.54 1.70301
15000 20 8000
1 0.025 322.43 1.79905
2 0.075 340.11 1.69857
3 0.125 350.03 1.62191
4 0.175 364.41 1.58429
5 0.225 390.55 1.45217
6 0.275 406.7 1.40234
7 0.325 424.92 1.30617
8 0.375 441.54 1.25242
9 0.425 448.57 1.21286
10 0.475 465.2 1.12839
11 0.525 476.34 1.06977
12 0.575 466.6 1.09034
13 0.625 445.15 1.18607
14 0.675 426.93 1.24344
15 0.725 420.07 1.29251
16 0.775 398.46 1.37042
17 0.825 373.88 1.50286
18 0.875 361.55 1.55703
19 0.925 346.41 1.62695
20 0.975 330.15 1.72757
16000 20 8000
1 0.025 306.33 1.8508
2 0.075 334.86 1.68695
3 0.125 356.87 1.59134
4 0.175 364.2 1.52883
5 0.225 383.43 1.45981
6 0.275 406.97 1.35464
7 0.325 420.1 1.28596
8 0.375 434.18 1.26101
9 0.425 446.62 1.20148
10 0.475 466.11 1.12534
11 0.525 483.91 1.07673
12 0.575 475.77 1.12051
13 0.625 449.05 1.20793
14 0.675 428.75 1.29484
15 0.725 418.25 1.34228
16 0.775 400.21 1.40828
17 0.825 392.3 1.44318
18 0.875 365.76 1.54801
19 0.925 343.74 1.64342
20 0.975 322.59 1.75857
17000 20 8000
1 0.025 310.93 1.82241
2 0.075 329.92 1.75204
3 0.125 350.13 1.60837
4 0.175 373.94 1.52827
5 0.225 388.63 1.47066
6 0.275 401.74 1.42567
7 0.325 428.12 1.30522
8 0.375 427.11 1.28586
9 0.425 453.97 1.20607
10 0.475 467.21 1.10887
11 0.525 484.58 1.04962
12 0.575 463.16 1.12469
13 0.625 448.74 1.18792
14 0.675 424.58 1.28373
15 0.725 419.94 1.31273
16 0.775 404.58 1.37026
17 0.825 388.97 1.44704
18 0.875 362.68 1.5453
19 0.925 346.29 1.6109
20 0.975 324.78 1.75026
18000 20 8000
1 0.025 317.81 1.81751
2 0.075 334.67 1.69028
3 0.125 347.24 1.65479
4 0.175 370.6 1.53592
5 0.225 384.31 1.49562
6 0.275 407.08 1.38014
7 0.325 421.73 1.29883
8 0.375 433.41 1.27181
9 0.425 447.79 1.18853
10 0.475 471.17 1.10893
11 0.525 478.19 1.06085
12 0.575 471.56 1.12156
13 0.625 452.58 1.19406
14 0.675 432.91 1.25675
15 0.725 408.06 1.3457
16 0.775 397.54 1.37446
17 0.825 381.4 1.46095
18 0.875 369.25 1.54246
19 0.925 345.61 1.63271
20 0.975 327.09 1.75129
19000 20 8000
1 0.025 315.61 1.79908
2 0.075 329.82 1.70308
3 0.125 343.96 1.63278
4 0.175 372.3 1.53323
5 0.225 384.82 1.45132
6 0.275 402.05 1.37734
7 0.325 426.21 1.2951
8 0.375 437.2 1.25041
9 0.425 448.66 1.19301
10 0.475 460.07 1.15046
11 0.525 477.48 1.08566
12 0.575 466.92 1.12916
13 0.625 450 1.21173
14 0.675 432.43 1.26854
15 0.725 422.77 1.30946
16 0.775 404.66 1.38438
17 0.825 392.04 1.44156
18 0.875 360.65 1.57203
19 0.925 346.34 1.63047
20 0.975 326.01 1.71263
20000 20 8000
1 0.025 317.33 1.87602
2 0.075 333.96 1.71854
3 0.125 345.85 1.68538
4 0.175 370.59 1.55097
5 0.225 380.61 1.45728
6 0.275 403.48 1.36314
7 0.325 423.92 1.29723
8 0.375 437.69 1.21691
9 0.425 445.13 1.20635
10 0.475 457.07 1.16643
11 0.525 470.23 1.10249
12 0.575 460.96 1.11255
13 0.625 452.73 1.15892
14 0.675 438.27 1.23227
15 0.725 416.93 1.29087
16 0.775 405.97 1.35751
17 0.825 393.18 1.44766
18 0.875 375.9 1.50591
19 0.925 346.29 1.61577
20 0.975 323.91 1.77921
21000 20 8000
1 0.025 307.58 1.79848
2 0.075 325.01 1.71065
3 0.125 350.24 1.61243
4 0.175 365.42 1.54134
5 0.225 386.17 1.47103
6 0.275 403.83 1.36486
7 0.325 426.21 1.30362
8 0.375 436.69 1.24808
9 0.425 444.23 1.25963
10 0.475 462.32 1.16034
11 0.525 472.1 1.09771
12 0.575 470.26 1.13665
13 0.625 454.16 1.19576
14 0.675 445.95 1.22857
15 0.725 427.33 1.32477
16 0.775 404.5 1.38979
17 0.825 382.25 1.47878
18 0.875 372.91 1.50944
19 0.925 340.86 1.61304
20 0.975 321.98 1.74402
22000 20 8000
1 0.025 318.14 1.82641
2 0.075 330.97 1.77621
3 0.125 339 1.67839
4 0.175 366.72 1.56607
5 0.225 390.88 1.45595
6 0.275 405.22 1.35061
7 0.325 424.26 1.29091
8 0.375 430.23 1.24841
9 0.425 449.76 1.18526
10 0.475 459.29 1.11561
11 0.525 472.2 1.07012
12 0.575 460.38 1.11959
13 0.625 451.95 1.16339
14 0.675 438.03 1.22565
15 0.725 424.58 1.30695
16 0.775 406.55 1.40099
17 0.825 390.12 1.46227
18 0.875 362.32 1.54729
19 0.925 346.31 1.64909
20 0.975 333.09 1.71265
23000 20 8000
1 0.025 305.23 1.85561
2 0.075 316.24 1.7503
3 0.125 342.29 1.63776
4 0.175 361.32 1.51963
5 0.225 387.91 1.46099
6 0.275 408.48 1.3798
7 0.325 424.77 1.35815
8 0.375 435.94 1.25237
9 0.425 461.41 1.16863
10 0.475 470.9 1.13807
11 0.525 474.68 1.09466
12 0.575 464.69 1.14285
13 0.625 450.09 1.20136
14 0.675 433.29 1.26397
15 0.725 421.43 1.33433
16 0.775 410.74 1.39007
17 0.825 379.41 1.49507
18 0.875 360.62 1.55625
19 0.925 356.87 1.60279
20 0.975 333.69 1.67823
24000 20 8000
1 0.025 313.08 1.84953
2 0.075 322.69 1.75898
3 0.125 342.02 1.66484
4 0.175 364.41 1.58799
5 0.225 390.21 1.42313
6 0.275 407.48 1.36668
7 0.325 421.23 1.32315
8 0.375 433.8 1.24529
9 0.425 452.64 1.16974
10 0.475 469.09 1.10316
11 0.525 479.85 1.04586
12 0.575 466.74 1.13969
13 0.625 451.65 1.19762
14 0.675 430.04 1.24855
15 0.725 417.52 1.33226
16 0.775 405.4 1.37658
17 0.825 376.39 1.48678
18 0.875 368.15 1.57185
19 0.925 354.02 1.62554
20 0.975 333.59 1.75438
25000 20 8000
1 0.025 313.64 1.80634
2 0.075 324.05 1.71933
3 0.125 343.33 1.61348
4 0.175 362.61 1.54935
5 0.225 378.09 1.47967
6 0.275 414.43 1.36693
7 0.325 426.07 1.3185
8 0.375 437.75 1.25993
9 0.425 452.65 1.19176
10 0.475 471.59 1.12704
11 0.525 481.75 1.07074
12 0.575 470.23 1.11916
13 0.625 454.89 1.19264
14 0.675 435.34 1.25786
15 0.725 414.67 1.34818
16 0.775 394.04 1.42967
17 0.825 380.43 1.48823
18 0.875 370.28 1.55129
19 0.925 351.59 1.6162
20 0.975 322.57 1.74431
26000 20 8000
1 0.025 308.4 1.84473
2 0.075 325.48 1.75272
3 0.125 352.47 1.62619
4 0.175 365.99 1.57145
5 0.225 391.95 1.43775
6 0.275 404.62 1.39548
7 0.325 425.2 1.28537
8 0.375 437.1 1.25101
9 0.425 451.1 1.18608
10 0.475 462.46 1.12661
11 0.525 472.15 1.07181
12 0.575 465.21 1.09039
13 0.625 453.32 1.17063
14 0.675 441.42 1.24972
15 0.725 426.69 1.2962
16 0.775 399.87 1.42231
17 0.825 381.98 1.5052
18 0.875 366.7 1.55646
19 0.925 348.6 1.65512
20 0.975 319.29 1.77302
27000 20 8000
1 0.025 316.89 1.80217
2 0.075 321.88 1.73459
3 0.125 345.89 1.63671
4 0.175 370.85 1.56358
5 0.225 385.03 1.51454
6 0.275 401.79 1.42111
7 0.325 427.07 1.29632
8 0.375 439.6 1.25415
9 0.425 460.68 1.16856
10 0.475 462.89 1.14373
11 0.525 478.74 1.08743
12 0.575 471.03 1.13158
13 0.625 452.29 1.17223
14 0.675 441.54 1.23337
15 0.725 419.13 1.30045
16 0.775 407.11 1.37184
17 0.825 373 1.4903
18 0.875 355.4 1.56997
19 0.925 343.35 1.66031
20 0.975 325.84 1.74216
28000 20 8000
1 0.025 321.78 1.81229
2 0.075 324.74 1.73651
3 0.125 331.52 1.62837
4 0.175 356.73 1.54322
5 0.225 381.28 1.46292
6 0.275 400.35 1.38037
7 0.325 422 1.30288
8 0.375 439.97 1.23362
9 0.425 452.44 1.17339
10 0.475 468.53 1.12737
11 0.525 475.51 1.07379
12 0.575 470.06 1.09663
13 0.625 453.01 1.17756
14 0.675 436.94 1.25999
15 0.725 420.08 1.31529
16 0.775 413.04 1.37827
17 0.825 383.82 1.50271
18 0.875 357.85 1.61812
19 0.925 346.4 1.67671
20 0.975 343.95 1.7339
29000 20 8000
1 0.025 313.72 1.80333
2 0.075 328.8 1.75365
3 0.125 344.37 1.64779
4 0.175 359.25 1.56358
5 0.225 376.97 1.51736
6 0.275 405.82 1.3928
7 0.325 421.01 1.32598
8 0.375 438.17 1.23454
9 0.425 455.16 1.17988
10 0.475 470.1 1.14151
11 0.525 480 1.09958
12 0.575 472.1 1.13539
13 0.625 452.87 1.20604
14 0.675 436.55 1.23619
15 0.725 422.97 1.30866
16 0.775 405.26 1.35304
17 0.825 376 1.46238
18 0.875 367.09 1.55544
19 0.925 346.8 1.60587
20 0.975 326.99 1.68799
30000 20 8000
1 0.025 321.81 1.8543
2 0.075 331.9 1.76078
3 0.125 347.94 1.66013
4 0.175 370.04 1.53588
5 0.225 381.88 1.46325
6 0.275 397.98 1.3956
7 0.325 413.72 1.3412
8 0.375 436.53 1.25271
9 0.425 452.86 1.17908
10 0.475 463.63 1.10328
11 0.525 478.33 1.05537
12 0.575 468.82 1.10703
13 0.625 448.79 1.18322
14 0.675 433.19 1.23552
15 0.725 420.5 1.27843
16 0.775 399.18 1.38097
17 0.825 385.02 1.44688
18 0.875 364.12 1.5531
19 0.925 345.41 1.63361
20 0.975 338.35 1.75235
31000 20 8000
1 0.025 310.01 1.8588
2 0.075 325.87 1.67311
3 0.125 346.92 1.60961
4 0.175 359.66 1.57966
5 0.225 383.88 1.48165
6 0.275 402.2 1.39282
7 0.325 414.39 1.33445
8 0.375 435.04 1.25763
9 0.425 454.62 1.19412
10 0.475 463.43 1.14685
11 0.525 480.51 1.09661
12 0.575 472.75 1.11081
13 0.625 462.84 1.1498
14 0.675 443.69 1.22548
15 0.725 419.56 1.35579
16 0.775 402.69 1.42284
17 0.825 384.78 1.4833
18 0.875 356.57 1.56438
19 0.925 351.55 1.60594
20 0.975 329.04 1.69805