From 0129ff9696c792016eaa4a7501420b436ac31a81 Mon Sep 17 00:00:00 2001 From: sjplimp Date: Tue, 7 Jun 2016 18:17:05 +0000 Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@15149 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- examples/KAPPA/log.ehex.7Jun16 | 225 +++++++++++++++ examples/KAPPA/log.heat.7Jun16 | 225 +++++++++++++++ examples/KAPPA/profile.heat.7Jun16 | 423 +++++++++++++++++++++++++++++ 3 files changed, 873 insertions(+) create mode 100644 examples/KAPPA/log.ehex.7Jun16 create mode 100644 examples/KAPPA/log.heat.7Jun16 create mode 100644 examples/KAPPA/profile.heat.7Jun16 diff --git a/examples/KAPPA/log.ehex.7Jun16 b/examples/KAPPA/log.ehex.7Jun16 new file mode 100644 index 0000000000..64d6340d18 --- /dev/null +++ b/examples/KAPPA/log.ehex.7Jun16 @@ -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 diff --git a/examples/KAPPA/log.heat.7Jun16 b/examples/KAPPA/log.heat.7Jun16 new file mode 100644 index 0000000000..c9a035ddc4 --- /dev/null +++ b/examples/KAPPA/log.heat.7Jun16 @@ -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 diff --git a/examples/KAPPA/profile.heat.7Jun16 b/examples/KAPPA/profile.heat.7Jun16 new file mode 100644 index 0000000000..32ddbec13d --- /dev/null +++ b/examples/KAPPA/profile.heat.7Jun16 @@ -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