From ed627579f6b78147e952767b65bc95375f8ca21b Mon Sep 17 00:00:00 2001 From: jtclemm Date: Thu, 30 Jan 2025 16:44:11 -0700 Subject: [PATCH] Adding bpm/plastic demo, fixing mistake in poisson example --- examples/bpm/plasticity/in.bpm.plasticity | 96 + .../log.19Nov2024.bpm.plasticity.g++.1 | 1573 +++++++++++++++++ .../bpm/poissons_ratio/in.bpm.poissons_ratio | 2 +- 3 files changed, 1670 insertions(+), 1 deletion(-) create mode 100644 examples/bpm/plasticity/in.bpm.plasticity create mode 100644 examples/bpm/plasticity/log.19Nov2024.bpm.plasticity.g++.1 diff --git a/examples/bpm/plasticity/in.bpm.plasticity b/examples/bpm/plasticity/in.bpm.plasticity new file mode 100644 index 0000000000..47e82455ed --- /dev/null +++ b/examples/bpm/plasticity/in.bpm.plasticity @@ -0,0 +1,96 @@ +# Load then unload four cantilevers with different plastic yield strains +# model from DOI: 10.1016/j.powtec.2024.120563 + +units lj +dimension 3 +boundary p s s +atom_style bond +special_bonds lj 0.0 1.0 1.0 coul 0.0 1.0 1.0 +newton on off +comm_modify vel yes cutoff 2.6 +lattice fcc 1.41 +region box block 0 20 -1 41 -5 5 units box +create_box 4 box bond/types 4 extra/bond/per/atom 20 extra/special/per/atom 50 +mass * 1.0 + +# Create initial disordered geometry somewhat near jamming + +region cantilever block 3.0 12.0 0 40 -2.5 2.5 units box +region cantileverw block 2.0 13.0 -1 41 -3.5 3.5 units box +create_atoms 1 region cantilever + +velocity all create 0.1 345910 + +pair_style bpm/spring +pair_coeff * * 1.0 1.0 1.0 + +fix 1 all nve +fix wtemp all wall/region cantileverw harmonic 1.0 1.0 1.0 + +thermo_style custom step ke pe pxx pyy pzz +thermo 100 + +timestep 0.1 +run 20000 +unfix wtemp + +# Replicate cylinder and add bonds + +replicate 4 1 1 + +region r1 block 0.0 20.0 EDGE EDGE EDGE EDGE side in units box +region r2 block 20.0 40.0 EDGE EDGE EDGE EDGE side in units box +region r3 block 40.0 60.0 EDGE EDGE EDGE EDGE side in units box +region r4 block 60.0 80.0 EDGE EDGE EDGE EDGE side in units box + +group c1 region r1 +group c2 region r2 +group c3 region r3 +group c4 region r4 + +set group c2 type 2 +set group c3 type 3 +set group c4 type 4 + +velocity all set 0.0 0.0 0.0 +neighbor 1.0 bin + +create_bonds many c1 c1 1 0.0 1.5 +create_bonds many c2 c2 2 0.0 1.5 +create_bonds many c3 c3 3 0.0 1.5 +create_bonds many c4 c4 4 0.0 1.5 + +neighbor 0.3 bin +special_bonds lj 0.0 1.0 1.0 coul 1.0 1.0 1.0 + +bond_style bpm/spring/plastic break no smooth no +bond_coeff 1 1.0 0 1.0 1.0 +bond_coeff 2 1.0 0 1.0 0.03 +bond_coeff 3 1.0 0 1.0 0.02 +bond_coeff 4 1.0 0 1.0 0.01 + +# apply load + +region anchor block EDGE EDGE EDGE 2.0 EDGE EDGE units box +region load block EDGE EDGE 38.0 EDGE EDGE EDGE units box +group anchor region anchor +group load region load + +fix 2 anchor setforce 0.0 0.0 0.0 +fix 3 load addforce 0.0 0.0 -0.0005 +fix 4 all viscous 0.001 + +compute zmin1 c1 reduce min z +compute zmin2 c2 reduce min z +compute zmin3 c3 reduce min z +compute zmin4 c4 reduce min z + +thermo_style custom step ke c_zmin1 c_zmin2 c_zmin3 c_zmin4 +#dump 1 all custom 200 atomDump id type x y z + +run 50000 + +# remove load + +unfix 3 +run 50000 diff --git a/examples/bpm/plasticity/log.19Nov2024.bpm.plasticity.g++.1 b/examples/bpm/plasticity/log.19Nov2024.bpm.plasticity.g++.1 new file mode 100644 index 0000000000..45fd829e70 --- /dev/null +++ b/examples/bpm/plasticity/log.19Nov2024.bpm.plasticity.g++.1 @@ -0,0 +1,1573 @@ +LAMMPS (19 Nov 2024 - Development - patch_19Nov2024-1028-g800a5f6310-modified) +# Load then unload four cantilevers with different plastic yield strains +# model from DOI: 10.1016/j.powtec.2024.120563 + +units lj +dimension 3 +boundary p s s +atom_style bond +special_bonds lj 0.0 1.0 1.0 coul 0.0 1.0 1.0 +newton on off +comm_modify vel yes cutoff 2.6 +lattice fcc 1.41 +Lattice spacing in x,y,z = 1.4156209 1.4156209 1.4156209 +region box block 0 20 -1 41 -5 5 units box +create_box 4 box bond/types 4 extra/bond/per/atom 20 extra/special/per/atom 50 +Created orthogonal box = (0 -1 -5) to (20 41 5) + 1 by 1 by 1 MPI processor grid +mass * 1.0 + +# Create initial disordered geometry somewhat near jamming + +region cantilever block 3.0 12.0 0 40 -2.5 2.5 units box +region cantileverw block 2.0 13.0 -1 41 -3.5 3.5 units box +create_atoms 1 region cantilever +Created 2394 atoms + using lattice units in orthogonal box = (0 -1.0042 -5.001) to (20 41.0042 5.001) + create_atoms CPU = 0.001 seconds + +velocity all create 0.1 345910 + +pair_style bpm/spring +pair_coeff * * 1.0 1.0 1.0 + +fix 1 all nve +fix wtemp all wall/region cantileverw harmonic 1.0 1.0 1.0 + +thermo_style custom step ke pe pxx pyy pzz +thermo 100 + +timestep 0.1 +run 20000 +Generated 0 of 6 mixed pair_coeff terms from geometric mixing rule +Neighbor list info ... + update: every = 1 steps, delay = 0 steps, check = yes + max neighbors/atom: 2000, page size: 100000 + master list distance cutoff = 1.3 + ghost atom cutoff = 2.6 + binsize = 0.65, bins = 31 65 16 + 1 neighbor lists, perpetual/occasional/extra = 1 0 0 + (1) pair bpm/spring, perpetual + attributes: half, newton on + pair build: half/bin/newton + stencil: half/bin/3d + bin: standard +Per MPI rank memory allocation (min/avg/max) = 9.752 | 9.752 | 9.752 Mbytes + Step KinEng PotEng Pxx Pyy Pzz + 0 0.14993734 0 0.071514461 0.070121562 0.071454556 + 100 0.0066788839 0.0029982274 0.016637041 0.016868542 0.015823802 + 200 0.0019715508 0.00071490141 0.007843495 0.0076223765 0.0067662021 + 300 0.00087217257 0.00026043228 0.003985349 0.004259066 0.0040664397 + 400 0.00050307484 0.00011882009 0.0024899485 0.0025393633 0.0024477801 + 500 0.00033241763 6.5955939e-05 0.0016895447 0.0018048761 0.0017308071 + 600 0.00023663671 4.4828066e-05 0.0013515511 0.001353668 0.0013208904 + 700 0.0001837694 2.7185024e-05 0.00094852029 0.0010286312 0.00099370739 + 800 0.00014655165 1.8873405e-05 0.00075979673 0.0008230931 0.00075345331 + 900 0.00012145909 1.4344012e-05 0.00060702604 0.00068683252 0.00061213385 + 1000 0.00010239299 1.2153676e-05 0.00055078741 0.00060838425 0.00054159268 + 1100 8.8075366e-05 9.9954299e-06 0.00045915185 0.00057179945 0.00047014896 + 1200 7.8336822e-05 7.4598572e-06 0.00038286385 0.000423282 0.00040077777 + 1300 6.8617848e-05 7.4071518e-06 0.00038962894 0.00045042653 0.00042176853 + 1400 6.1816066e-05 5.7757188e-06 0.00035626235 0.00037452754 0.00033148835 + 1500 5.5158763e-05 5.3010672e-06 0.00032659745 0.00035059437 0.0003143055 + 1600 5.0680523e-05 4.4203877e-06 0.00027420802 0.00032194903 0.00030193132 + 1700 4.6020086e-05 4.0228214e-06 0.00026909527 0.00029969818 0.00027492066 + 1800 4.1968387e-05 3.9751303e-06 0.00027184181 0.00028203292 0.00026747334 + 1900 3.8741903e-05 3.5539155e-06 0.00024183418 0.00026550737 0.00025356968 + 2000 3.6204768e-05 2.9462894e-06 0.00022277001 0.00025408288 0.0002123446 + 2100 3.3266509e-05 2.8028832e-06 0.00021917117 0.00024069099 0.00020369452 + 2200 3.102664e-05 2.6085678e-06 0.0001882006 0.00023685724 0.00019134947 + 2300 2.903812e-05 2.3143895e-06 0.0001749797 0.00021704715 0.00019021214 + 2400 2.7141792e-05 2.1848123e-06 0.00018014226 0.00020933228 0.00017930569 + 2500 2.5612882e-05 1.8545406e-06 0.00017807794 0.00016882959 0.0001750063 + 2600 2.4068378e-05 1.7104525e-06 0.00015853683 0.000178264 0.00015848288 + 2700 2.2573196e-05 1.7771307e-06 0.00015577573 0.00017747772 0.00016654472 + 2800 2.1422175e-05 1.4847142e-06 0.00014322637 0.00013897114 0.00014490629 + 2900 2.0433268e-05 1.3514585e-06 0.00012745407 0.00014678277 0.000135221 + 3000 1.9204132e-05 1.3849721e-06 0.00012986068 0.00015084177 0.00014142185 + 3100 1.8212704e-05 1.2748732e-06 0.00013317833 0.00014368027 0.00013451894 + 3200 1.7326097e-05 1.3213286e-06 0.00012777819 0.00014793494 0.00012886678 + 3300 1.6615969e-05 1.1158927e-06 0.00011223033 0.00012323335 0.0001210878 + 3400 1.5855041e-05 1.1229014e-06 0.00011468711 0.00013243589 0.00012051018 + 3500 1.532374e-05 8.888678e-07 0.00010638549 0.00011430069 0.00010089698 + 3600 1.467535e-05 8.7654654e-07 8.8094885e-05 0.00011642682 0.00010671988 + 3700 1.4062105e-05 8.1807269e-07 9.4265348e-05 0.0001161673 9.3404926e-05 + 3800 1.3552196e-05 8.3487647e-07 8.6030873e-05 0.00011280194 0.0001010521 + 3900 1.3022591e-05 7.9989391e-07 0.00010119027 0.00010258667 9.1434393e-05 + 4000 1.2525445e-05 8.1520281e-07 9.5211467e-05 0.00010785495 9.97298e-05 + 4100 1.2181734e-05 6.4410557e-07 9.4373715e-05 8.9279858e-05 7.9893316e-05 + 4200 1.1704251e-05 7.0450874e-07 8.0872404e-05 9.4334475e-05 8.8302993e-05 + 4300 1.1280286e-05 7.1408139e-07 9.1119222e-05 8.5870034e-05 8.8445587e-05 + 4400 1.1016643e-05 6.0873467e-07 7.9024773e-05 8.1967988e-05 8.5977338e-05 + 4500 1.0668349e-05 5.8516931e-07 7.8294095e-05 8.0043527e-05 8.4715888e-05 + 4600 1.0331041e-05 5.7376843e-07 7.5679063e-05 8.2398328e-05 7.4278168e-05 + 4700 1.0004816e-05 5.4597537e-07 7.2607537e-05 7.7956907e-05 6.9095147e-05 + 4800 9.7594591e-06 5.2681078e-07 7.6737892e-05 6.79129e-05 7.2013864e-05 + 4900 9.4372215e-06 5.108678e-07 7.1926492e-05 7.6336467e-05 7.2404427e-05 + 5000 9.2301398e-06 4.9943356e-07 6.9272822e-05 7.4437007e-05 6.8591392e-05 + 5100 8.9087719e-06 5.0137676e-07 7.1953364e-05 7.0112306e-05 7.1894235e-05 + 5200 8.6733384e-06 4.6530133e-07 5.9160431e-05 7.2957759e-05 6.7731265e-05 + 5300 8.5237067e-06 4.0166784e-07 5.6976938e-05 6.4780167e-05 6.2493655e-05 + 5400 8.2585748e-06 3.9144457e-07 5.6287251e-05 6.5281415e-05 6.4678701e-05 + 5500 8.1089159e-06 3.5903718e-07 6.2738749e-05 5.8042636e-05 5.1954549e-05 + 5600 7.8605737e-06 3.9176106e-07 6.1198736e-05 6.5959291e-05 5.8081374e-05 + 5700 7.6788462e-06 3.8505614e-07 5.5766462e-05 6.4100303e-05 6.0012403e-05 + 5800 7.5163424e-06 3.4531571e-07 5.326349e-05 6.1833688e-05 4.8218455e-05 + 5900 7.3552892e-06 3.1726542e-07 5.2593444e-05 5.3015093e-05 5.2366381e-05 + 6000 7.1550139e-06 3.188599e-07 4.7713401e-05 5.6322276e-05 5.2424598e-05 + 6100 7.0361757e-06 3.0641052e-07 4.7286992e-05 5.8028248e-05 5.2063931e-05 + 6200 6.8801441e-06 2.7974487e-07 5.3580486e-05 4.6047586e-05 4.7327222e-05 + 6300 6.6802143e-06 3.3499395e-07 5.8400669e-05 5.5627535e-05 4.7978034e-05 + 6400 6.5789434e-06 2.8602311e-07 5.0198228e-05 5.5844532e-05 4.877662e-05 + 6500 6.4436475e-06 2.6008325e-07 4.7602086e-05 4.4486849e-05 4.6404674e-05 + 6600 6.2843629e-06 2.7832981e-07 4.7074622e-05 4.6606369e-05 4.9936757e-05 + 6700 6.1703075e-06 2.8089174e-07 4.7083973e-05 4.6575052e-05 5.0006622e-05 + 6800 6.0794707e-06 2.5065496e-07 4.3818817e-05 4.8409259e-05 3.8917908e-05 + 6900 5.9095362e-06 2.8634305e-07 5.2986415e-05 4.8027488e-05 5.1770643e-05 + 7000 5.8364704e-06 2.4891212e-07 4.3004439e-05 4.6236407e-05 4.203024e-05 + 7100 5.676111e-06 2.8000071e-07 4.2647162e-05 5.4717075e-05 4.5375773e-05 + 7200 5.6485498e-06 2.0372933e-07 3.7420101e-05 4.5276234e-05 3.716119e-05 + 7300 5.4828422e-06 2.3962736e-07 4.4493988e-05 4.3446503e-05 4.673755e-05 + 7400 5.3980604e-06 2.3856758e-07 4.0222008e-05 4.6533781e-05 4.3072176e-05 + 7500 5.3006008e-06 2.2358346e-07 4.0025366e-05 3.6897725e-05 4.4834626e-05 + 7600 5.197958e-06 2.3660321e-07 4.0784719e-05 4.4961031e-05 4.2087955e-05 + 7700 5.1806799e-06 1.7718602e-07 3.2411759e-05 3.9110371e-05 3.8802708e-05 + 7800 5.0962887e-06 1.7412195e-07 3.6552483e-05 3.5361276e-05 3.2904976e-05 + 7900 4.9703998e-06 2.1337014e-07 4.3936409e-05 4.0886722e-05 3.9732939e-05 + 8000 4.9053377e-06 1.7953438e-07 3.7608294e-05 3.4731494e-05 3.6081071e-05 + 8100 4.8171492e-06 1.9697727e-07 3.3678657e-05 4.1808664e-05 4.1007093e-05 + 8200 4.7721132e-06 1.7544811e-07 3.4665105e-05 3.2777884e-05 3.774731e-05 + 8300 4.708456e-06 1.6203898e-07 3.7062724e-05 3.4881585e-05 2.9315416e-05 + 8400 4.6397062e-06 1.3777361e-07 3.1993062e-05 2.840598e-05 2.8943249e-05 + 8500 4.5983716e-06 1.3333705e-07 3.2032331e-05 3.1515185e-05 2.7257332e-05 + 8600 4.5175708e-06 1.4915416e-07 3.0728382e-05 3.0468802e-05 3.1961748e-05 + 8700 4.4347789e-06 1.5492856e-07 3.0961313e-05 3.4774356e-05 2.957729e-05 + 8800 4.3438643e-06 1.9416152e-07 3.5496445e-05 3.650361e-05 3.7686819e-05 + 8900 4.2916882e-06 1.7972527e-07 3.4999135e-05 3.4301174e-05 3.369003e-05 + 9000 4.2420397e-06 1.653212e-07 3.2002391e-05 3.5994842e-05 3.6333896e-05 + 9100 4.156708e-06 1.7541257e-07 3.1781091e-05 3.3713908e-05 3.6414424e-05 + 9200 4.1206942e-06 1.6684543e-07 3.1148086e-05 3.4009984e-05 3.6543294e-05 + 9300 4.0739027e-06 1.5661214e-07 3.2100841e-05 3.1002639e-05 3.2377644e-05 + 9400 4.0749842e-06 1.1325533e-07 2.7132088e-05 2.6729932e-05 2.6505394e-05 + 9500 4.0182431e-06 1.1916661e-07 2.8267149e-05 2.7495849e-05 2.5348393e-05 + 9600 3.9394743e-06 1.4897486e-07 3.1303396e-05 2.8986697e-05 3.3906299e-05 + 9700 3.9202642e-06 1.1265772e-07 2.6904853e-05 2.9419644e-05 2.3388387e-05 + 9800 3.8563029e-06 1.4201267e-07 2.9742419e-05 3.0333264e-05 2.8597551e-05 + 9900 3.8407551e-06 1.0925008e-07 2.4813089e-05 2.518408e-05 2.5132946e-05 + 10000 3.8031732e-06 1.0429945e-07 2.4587006e-05 2.4160366e-05 2.2475372e-05 + 10100 3.7403802e-06 1.188634e-07 2.621425e-05 2.747381e-05 2.7423478e-05 + 10200 3.6933569e-06 1.177636e-07 2.3506114e-05 3.0373875e-05 2.7330288e-05 + 10300 3.6359757e-06 1.3540565e-07 2.4720627e-05 3.1264106e-05 3.0431809e-05 + 10400 3.6158282e-06 1.1027073e-07 2.6627036e-05 2.5041131e-05 2.3909283e-05 + 10500 3.5857824e-06 1.0049365e-07 2.5306577e-05 2.7251152e-05 2.1635787e-05 + 10600 3.5260536e-06 1.1579533e-07 2.581535e-05 2.7130923e-05 2.4935895e-05 + 10700 3.5024547e-06 1.0954559e-07 2.5027035e-05 2.8379233e-05 2.2447535e-05 + 10800 3.4853309e-06 9.4720082e-08 2.2724388e-05 2.0056563e-05 2.3309701e-05 + 10900 3.4396529e-06 1.0829632e-07 2.5958339e-05 2.6491256e-05 2.1617336e-05 + 11000 3.3987023e-06 9.7092727e-08 2.2470175e-05 2.2846816e-05 2.091066e-05 + 11100 3.3606128e-06 9.8398405e-08 2.2648156e-05 2.4731296e-05 2.0418175e-05 + 11200 3.3396876e-06 7.6470057e-08 1.9873816e-05 1.9797248e-05 1.9115439e-05 + 11300 3.3041431e-06 9.3155692e-08 2.3746612e-05 2.3286585e-05 1.93667e-05 + 11400 3.2797866e-06 9.0584145e-08 2.1717142e-05 2.5320734e-05 2.118408e-05 + 11500 3.2454003e-06 8.5422835e-08 2.2434639e-05 1.920118e-05 2.0617104e-05 + 11600 3.1819911e-06 1.1244129e-07 1.791276e-05 2.6915583e-05 2.355823e-05 + 11700 3.1635572e-06 1.0031236e-07 2.2122734e-05 2.3651576e-05 2.2707601e-05 + 11800 3.1420932e-06 8.9733842e-08 2.0691568e-05 2.1346189e-05 2.3543388e-05 + 11900 3.1066705e-06 9.1501672e-08 2.0236573e-05 2.2932414e-05 1.9990374e-05 + 12000 3.0745943e-06 8.9277079e-08 2.0150526e-05 2.2352883e-05 1.8398403e-05 + 12100 3.0733981e-06 7.4616245e-08 1.9385797e-05 1.7454795e-05 1.9264248e-05 + 12200 3.0257929e-06 9.2174543e-08 1.869441e-05 2.4202442e-05 2.2015e-05 + 12300 3.0234068e-06 6.5090899e-08 1.767132e-05 1.9462735e-05 1.6725567e-05 + 12400 3.0027086e-06 6.2193553e-08 1.4880016e-05 2.0662066e-05 1.8614087e-05 + 12500 2.9319141e-06 1.0097547e-07 2.1488942e-05 2.1615501e-05 2.148531e-05 + 12600 2.9289151e-06 7.6302482e-08 1.7692348e-05 1.8642681e-05 1.7425701e-05 + 12700 2.90702e-06 7.6215717e-08 1.9509559e-05 1.8415327e-05 2.2776532e-05 + 12800 2.8758381e-06 7.9105433e-08 1.8521166e-05 2.0454146e-05 2.0370962e-05 + 12900 2.8317602e-06 9.1385197e-08 2.0335336e-05 2.4975614e-05 2.0240701e-05 + 13000 2.8410057e-06 6.0044423e-08 1.629434e-05 1.823554e-05 1.511721e-05 + 13100 2.782039e-06 8.4410675e-08 2.0919093e-05 2.0809882e-05 1.8657917e-05 + 13200 2.791875e-06 5.337143e-08 1.5825524e-05 1.9356399e-05 1.2911019e-05 + 13300 2.7487173e-06 7.6547916e-08 1.962803e-05 1.8981737e-05 2.0142921e-05 + 13400 2.7299962e-06 6.7372367e-08 1.7153528e-05 1.8103989e-05 1.7732609e-05 + 13500 2.7057441e-06 7.2429583e-08 1.9019906e-05 1.9691015e-05 1.9578145e-05 + 13600 2.6778573e-06 8.0573212e-08 1.8433433e-05 2.0729038e-05 2.0267435e-05 + 13700 2.6826077e-06 4.4511734e-08 1.3084576e-05 1.5000823e-05 1.3750559e-05 + 13800 2.647424e-06 5.9574367e-08 1.4631775e-05 1.4306273e-05 1.7157491e-05 + 13900 2.6110162e-06 7.4881186e-08 1.7818332e-05 2.148169e-05 1.7463563e-05 + 14000 2.6014529e-06 6.3294079e-08 1.5635773e-05 1.7333417e-05 1.6871959e-05 + 14100 2.5707329e-06 7.3706564e-08 1.691795e-05 2.2276988e-05 1.7435609e-05 + 14200 2.5693783e-06 5.5754478e-08 1.6045272e-05 1.5266604e-05 1.6872656e-05 + 14300 2.543378e-06 5.6841288e-08 1.7342484e-05 1.6861913e-05 1.5333195e-05 + 14400 2.5180009e-06 6.7438588e-08 1.7569821e-05 1.8124358e-05 1.8789645e-05 + 14500 2.505371e-06 5.8372098e-08 1.5370256e-05 1.7944439e-05 1.3727925e-05 + 14600 2.50052e-06 5.2960357e-08 1.2874059e-05 1.6976938e-05 1.5523433e-05 + 14700 2.4677314e-06 5.6098859e-08 1.5234684e-05 1.4810739e-05 1.6137648e-05 + 14800 2.4431342e-06 6.9542809e-08 1.7307233e-05 1.675592e-05 1.8785122e-05 + 14900 2.4217996e-06 6.5392872e-08 1.8214833e-05 1.7872976e-05 1.6079527e-05 + 15000 2.4223306e-06 4.8254257e-08 1.6243518e-05 1.3430592e-05 1.4053815e-05 + 15100 2.3971545e-06 5.8012806e-08 1.5024579e-05 1.5105347e-05 1.3560674e-05 + 15200 2.3734628e-06 5.4310001e-08 1.3063191e-05 1.5414145e-05 1.5458277e-05 + 15300 2.3523286e-06 6.1640634e-08 1.7876627e-05 1.6249813e-05 1.687417e-05 + 15400 2.3489693e-06 4.7278336e-08 1.2011731e-05 1.2230688e-05 1.678822e-05 + 15500 2.3284108e-06 5.2715812e-08 1.3176372e-05 1.6503583e-05 1.3955123e-05 + 15600 2.2947375e-06 6.7465869e-08 1.5326712e-05 1.9881502e-05 1.6487331e-05 + 15700 2.3019796e-06 4.8005103e-08 1.4590446e-05 1.4276946e-05 1.4055723e-05 + 15800 2.2826912e-06 4.9641879e-08 1.3481867e-05 1.6276418e-05 1.3419255e-05 + 15900 2.2583072e-06 5.4614844e-08 1.7382114e-05 1.604672e-05 1.4741544e-05 + 16000 2.2430797e-06 4.8420384e-08 1.2843624e-05 1.2618834e-05 1.2491054e-05 + 16100 2.2264678e-06 4.796181e-08 1.4438989e-05 1.5552177e-05 1.2048946e-05 + 16200 2.2123007e-06 5.2907699e-08 1.2747032e-05 1.4431874e-05 1.6152188e-05 + 16300 2.1942666e-06 5.1538008e-08 1.3243122e-05 1.4793006e-05 1.4855916e-05 + 16400 2.1800741e-06 5.0280071e-08 1.1549921e-05 1.5836438e-05 1.5087605e-05 + 16500 2.1718081e-06 4.2144732e-08 1.4603756e-05 1.3084635e-05 1.2983558e-05 + 16600 2.15145e-06 4.6326815e-08 1.4904114e-05 1.4500088e-05 1.2283839e-05 + 16700 2.1563978e-06 3.2862475e-08 1.1841441e-05 1.2463519e-05 1.0662253e-05 + 16800 2.1128692e-06 5.462978e-08 1.3455978e-05 1.6216741e-05 1.4636363e-05 + 16900 2.1178726e-06 3.6500222e-08 1.1666564e-05 1.2789709e-05 9.6393609e-06 + 17000 2.1058751e-06 3.6428095e-08 1.3109685e-05 1.304962e-05 1.0397847e-05 + 17100 2.0817867e-06 5.1894295e-08 1.5527039e-05 1.442429e-05 1.1662417e-05 + 17200 2.0697616e-06 4.8400908e-08 1.2051247e-05 1.4854095e-05 1.4277423e-05 + 17300 2.0496506e-06 5.1765629e-08 1.5484172e-05 1.6516861e-05 1.3787806e-05 + 17400 2.0318945e-06 5.0178134e-08 1.2610435e-05 1.5926848e-05 1.3146463e-05 + 17500 2.0312975e-06 4.2775082e-08 1.1413756e-05 1.2303892e-05 1.332512e-05 + 17600 2.0286562e-06 3.3246228e-08 9.9623166e-06 1.1190246e-05 1.0713751e-05 + 17700 2.0057117e-06 4.1465626e-08 1.28564e-05 1.2354674e-05 1.1732572e-05 + 17800 1.9775891e-06 5.3253383e-08 1.3320474e-05 1.6084721e-05 1.1697579e-05 + 17900 1.9763707e-06 4.0037754e-08 1.0787927e-05 1.0905895e-05 1.2524822e-05 + 18000 1.9659819e-06 4.2213997e-08 1.0786016e-05 1.356023e-05 1.3052018e-05 + 18100 1.962808e-06 3.1950711e-08 1.0829314e-05 1.246617e-05 1.1079602e-05 + 18200 1.9360469e-06 4.3810974e-08 1.382493e-05 1.245849e-05 1.2288243e-05 + 18300 1.9407685e-06 3.0042901e-08 7.3716082e-06 1.109871e-05 1.226921e-05 + 18400 1.910359e-06 4.7042217e-08 1.3644803e-05 1.4884969e-05 1.1956804e-05 + 18500 1.8946497e-06 5.0743345e-08 1.1974771e-05 1.4385508e-05 1.7610709e-05 + 18600 1.8962987e-06 3.9757677e-08 1.0582884e-05 1.2673246e-05 1.3298389e-05 + 18700 1.8742478e-06 4.4718119e-08 1.5344401e-05 1.2546575e-05 1.1764246e-05 + 18800 1.8650878e-06 3.7704721e-08 1.07802e-05 1.2917662e-05 9.5373738e-06 + 18900 1.8685713e-06 3.6504734e-08 1.0815769e-05 1.3268273e-05 9.2228416e-06 + 19000 1.8510631e-06 3.5157873e-08 1.2534306e-05 1.2128038e-05 1.0355377e-05 + 19100 1.844547e-06 3.5807341e-08 1.1011196e-05 1.2288462e-05 1.042938e-05 + 19200 1.8237482e-06 4.4374495e-08 1.0976301e-05 1.0525434e-05 1.3496339e-05 + 19300 1.809658e-06 4.663351e-08 1.0427019e-05 1.2761301e-05 1.2443509e-05 + 19400 1.8047595e-06 3.7730816e-08 1.2549441e-05 9.4564734e-06 1.159109e-05 + 19500 1.7912254e-06 3.9882052e-08 1.3940679e-05 1.220604e-05 1.1210472e-05 + 19600 1.789224e-06 3.7125586e-08 1.0231473e-05 1.2070753e-05 1.1047121e-05 + 19700 1.7845362e-06 3.456485e-08 1.12335e-05 1.042505e-05 9.0484867e-06 + 19800 1.7683141e-06 3.8353201e-08 1.29417e-05 1.0764969e-05 1.116828e-05 + 19900 1.7664071e-06 3.2153664e-08 9.4515931e-06 1.1051242e-05 1.0924767e-05 + 20000 1.7618892e-06 2.2689236e-08 8.4011955e-06 7.5721128e-06 8.9521609e-06 +Loop time of 2.54444 on 1 procs for 20000 steps with 2394 atoms + +Performance: 67912724.488 tau/day, 7860.269 timesteps/s, 18.817 Matom-step/s +99.7% CPU use with 1 MPI tasks x no OpenMP threads + +MPI task timing breakdown: +Section | min time | avg time | max time |%varavg| %total +--------------------------------------------------------------- +Pair | 1.2052 | 1.2052 | 1.2052 | 0.0 | 47.37 +Bond | 0.00069514 | 0.00069514 | 0.00069514 | 0.0 | 0.03 +Neigh | 0.18155 | 0.18155 | 0.18155 | 0.0 | 7.14 +Comm | 0.0076275 | 0.0076275 | 0.0076275 | 0.0 | 0.30 +Output | 0.0043813 | 0.0043813 | 0.0043813 | 0.0 | 0.17 +Modify | 1.0429 | 1.0429 | 1.0429 | 0.0 | 40.99 +Other | | 0.1021 | | | 4.01 + +Nlocal: 2394 ave 2394 max 2394 min +Histogram: 1 0 0 0 0 0 0 0 0 0 +Nghost: 0 ave 0 max 0 min +Histogram: 1 0 0 0 0 0 0 0 0 0 +Neighs: 11170 ave 11170 max 11170 min +Histogram: 1 0 0 0 0 0 0 0 0 0 + +Total # of neighbors = 11170 +Ave neighs/atom = 4.6658312 +Ave special neighs/atom = 0 +Neighbor list builds = 199 +Dangerous builds = 0 +unfix wtemp + +# Replicate cylinder and add bonds + +replicate 4 1 1 +Replication is creating a 4x1x1 = 4 times larger system... + orthogonal box = (0 -0.0019118212 -2.5018545) to (80 40.003414 2.5021045) + 1 by 1 by 1 MPI processor grid + 9576 atoms +Finding 1-2 1-3 1-4 neighbors ... + special bond factors lj: 0 1 1 + special bond factors coul: 0 1 1 + 0 = max # of 1-2 neighbors + 101 = max # of special neighbors + special bonds CPU = 0.000 seconds + replicate CPU = 0.002 seconds + +region r1 block 0.0 20.0 EDGE EDGE EDGE EDGE side in units box +region r2 block 20.0 40.0 EDGE EDGE EDGE EDGE side in units box +region r3 block 40.0 60.0 EDGE EDGE EDGE EDGE side in units box +region r4 block 60.0 80.0 EDGE EDGE EDGE EDGE side in units box + +group c1 region r1 +2394 atoms in group c1 +group c2 region r2 +2394 atoms in group c2 +group c3 region r3 +2394 atoms in group c3 +group c4 region r4 +2394 atoms in group c4 + +set group c2 type 2 +Setting atom values ... + 2394 settings made for type +set group c3 type 3 +Setting atom values ... + 2394 settings made for type +set group c4 type 4 +Setting atom values ... + 2394 settings made for type + +velocity all set 0.0 0.0 0.0 +neighbor 1.0 bin + +create_bonds many c1 c1 1 0.0 1.5 +Generated 0 of 6 mixed pair_coeff terms from geometric mixing rule +Neighbor list info ... + update: every = 1 steps, delay = 0 steps, check = yes + max neighbors/atom: 2000, page size: 100000 + master list distance cutoff = 2 + ghost atom cutoff = 2.6 + binsize = 1, bins = 80 41 6 + 2 neighbor lists, perpetual/occasional/extra = 1 1 0 + (1) command create_bonds, occasional + attributes: full, newton on + pair build: full/bin + stencil: full/bin/3d + bin: standard + (2) pair bpm/spring, perpetual + attributes: half, newton on + pair build: half/bin/newton + stencil: half/bin/3d + bin: standard +Added 13760 bonds, new total = 13760 +Finding 1-2 1-3 1-4 neighbors ... + special bond factors lj: 0 1 1 + special bond factors coul: 0 1 1 + 17 = max # of 1-2 neighbors + 101 = max # of special neighbors + special bonds CPU = 0.001 seconds +create_bonds many c2 c2 2 0.0 1.5 +Generated 0 of 6 mixed pair_coeff terms from geometric mixing rule +WARNING: Bonds are defined but no bond style is set (../force.cpp:197) +WARNING: Likewise 1-2 special neighbor interactions != 1.0 (../force.cpp:199) +Neighbor list info ... + update: every = 1 steps, delay = 0 steps, check = yes + max neighbors/atom: 2000, page size: 100000 + master list distance cutoff = 2 + ghost atom cutoff = 2.6 + binsize = 1, bins = 80 41 6 + 2 neighbor lists, perpetual/occasional/extra = 1 1 0 + (1) command create_bonds, occasional + attributes: full, newton on + pair build: full/bin + stencil: full/bin/3d + bin: standard + (2) pair bpm/spring, perpetual + attributes: half, newton on + pair build: half/bin/newton + stencil: half/bin/3d + bin: standard +Added 13760 bonds, new total = 27520 +Finding 1-2 1-3 1-4 neighbors ... + special bond factors lj: 0 1 1 + special bond factors coul: 0 1 1 + 17 = max # of 1-2 neighbors + 101 = max # of special neighbors + special bonds CPU = 0.002 seconds +create_bonds many c3 c3 3 0.0 1.5 +Generated 0 of 6 mixed pair_coeff terms from geometric mixing rule +WARNING: Bonds are defined but no bond style is set (../force.cpp:197) +WARNING: Likewise 1-2 special neighbor interactions != 1.0 (../force.cpp:199) +Neighbor list info ... + update: every = 1 steps, delay = 0 steps, check = yes + max neighbors/atom: 2000, page size: 100000 + master list distance cutoff = 2 + ghost atom cutoff = 2.6 + binsize = 1, bins = 80 41 6 + 2 neighbor lists, perpetual/occasional/extra = 1 1 0 + (1) command create_bonds, occasional + attributes: full, newton on + pair build: full/bin + stencil: full/bin/3d + bin: standard + (2) pair bpm/spring, perpetual + attributes: half, newton on + pair build: half/bin/newton + stencil: half/bin/3d + bin: standard +Added 13760 bonds, new total = 41280 +Finding 1-2 1-3 1-4 neighbors ... + special bond factors lj: 0 1 1 + special bond factors coul: 0 1 1 + 17 = max # of 1-2 neighbors + 101 = max # of special neighbors + special bonds CPU = 0.003 seconds +create_bonds many c4 c4 4 0.0 1.5 +Generated 0 of 6 mixed pair_coeff terms from geometric mixing rule +WARNING: Bonds are defined but no bond style is set (../force.cpp:197) +WARNING: Likewise 1-2 special neighbor interactions != 1.0 (../force.cpp:199) +Neighbor list info ... + update: every = 1 steps, delay = 0 steps, check = yes + max neighbors/atom: 2000, page size: 100000 + master list distance cutoff = 2 + ghost atom cutoff = 2.6 + binsize = 1, bins = 80 41 6 + 2 neighbor lists, perpetual/occasional/extra = 1 1 0 + (1) command create_bonds, occasional + attributes: full, newton on + pair build: full/bin + stencil: full/bin/3d + bin: standard + (2) pair bpm/spring, perpetual + attributes: half, newton on + pair build: half/bin/newton + stencil: half/bin/3d + bin: standard +Added 13760 bonds, new total = 55040 +Finding 1-2 1-3 1-4 neighbors ... + special bond factors lj: 0 1 1 + special bond factors coul: 0 1 1 + 17 = max # of 1-2 neighbors + 101 = max # of special neighbors + special bonds CPU = 0.003 seconds + +neighbor 0.3 bin +special_bonds lj 0.0 1.0 1.0 coul 1.0 1.0 1.0 + +bond_style bpm/spring/plastic break no smooth no +bond_coeff 1 1.0 0 1.0 1.0 +bond_coeff 2 1.0 0 1.0 0.03 +bond_coeff 3 1.0 0 1.0 0.02 +bond_coeff 4 1.0 0 1.0 0.01 + +# apply load + +region anchor block EDGE EDGE EDGE 2.0 EDGE EDGE units box +region load block EDGE EDGE 38.0 EDGE EDGE EDGE units box +group anchor region anchor +512 atoms in group anchor +group load region load +592 atoms in group load + +fix 2 anchor setforce 0.0 0.0 0.0 +fix 3 load addforce 0.0 0.0 -0.0005 +fix 4 all viscous 0.001 + +compute zmin1 c1 reduce min z +compute zmin2 c2 reduce min z +compute zmin3 c3 reduce min z +compute zmin4 c4 reduce min z + +thermo_style custom step ke c_zmin1 c_zmin2 c_zmin3 c_zmin4 +#dump 1 all custom 200 atomDump id type x y z + +run 50000 + +CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE + +Your simulation uses code contributions which should be cited: + +- BPM bond style: doi:10.1039/D3SM01373A + +@Article{Clemmer2024, + author = {Clemmer, Joel T. and Monti, Joseph M. and Lechman, Jeremy B.}, + title = {A soft departure from jamming: the compaction of deformable + granular matter under high pressures}, + journal = {Soft Matter}, + year = 2024, + volume = 20, + number = 8, + pages = {1702--1718} +} + +CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE + +Generated 0 of 6 mixed pair_coeff terms from geometric mixing rule +Neighbor list info ... + update: every = 1 steps, delay = 0 steps, check = yes + max neighbors/atom: 2000, page size: 100000 + master list distance cutoff = 1.3 + ghost atom cutoff = 2.6 + binsize = 0.65, bins = 124 62 8 + 1 neighbor lists, perpetual/occasional/extra = 1 0 0 + (1) pair bpm/spring, perpetual + attributes: half, newton on + pair build: half/bin/newton + stencil: half/bin/3d + bin: standard +Per MPI rank memory allocation (min/avg/max) = 20.02 | 20.02 | 20.02 Mbytes + Step KinEng c_zmin1 c_zmin2 c_zmin3 c_zmin4 + 20000 0 -2.5009426 -2.5009426 -2.5009426 -2.5009426 + 20100 4.0956567e-07 -2.5160602 -2.5160602 -2.5160602 -2.5160602 + 20200 1.3143434e-06 -2.555086 -2.555086 -2.555086 -2.555086 + 20300 2.5328497e-06 -2.6114806 -2.6114806 -2.6114806 -2.6114806 + 20400 4.0044252e-06 -2.6807339 -2.6807339 -2.6807339 -2.6807339 + 20500 5.683827e-06 -2.7603965 -2.7603965 -2.7603965 -2.7603965 + 20600 7.5588865e-06 -2.8493381 -2.8493381 -2.8493381 -2.8493381 + 20700 9.5460405e-06 -2.9466232 -2.9466232 -2.9466232 -2.9466232 + 20800 1.1434976e-05 -3.0516294 -3.0516294 -3.0516294 -3.0516294 + 20900 1.3298452e-05 -3.1649086 -3.1649086 -3.1649086 -3.1649086 + 21000 1.5334087e-05 -3.2853774 -3.2853774 -3.2853774 -3.2853774 + 21100 1.7644543e-05 -3.4124904 -3.4124904 -3.4124904 -3.4124904 + 21200 2.0249305e-05 -3.5471314 -3.5471314 -3.5471314 -3.5471314 + 21300 2.3096201e-05 -3.6870368 -3.6870368 -3.6870368 -3.6870367 + 21400 2.6137865e-05 -3.8263018 -3.8263018 -3.8263018 -3.8262991 + 21500 2.94271e-05 -3.9605484 -3.9605484 -3.9605484 -3.9605778 + 21600 3.2932776e-05 -4.0900526 -4.0900526 -4.0900526 -4.0902275 + 21700 3.6482197e-05 -4.2183658 -4.2183658 -4.2183658 -4.2187208 + 21800 3.9901933e-05 -4.3498077 -4.3498077 -4.3498077 -4.3503142 + 21900 4.3096586e-05 -4.4878953 -4.4878953 -4.4878953 -4.4885197 + 22000 4.6065281e-05 -4.6348263 -4.6348263 -4.6348263 -4.6355379 + 22100 4.8854721e-05 -4.7916676 -4.7916676 -4.7916676 -4.7924466 + 22200 5.1467377e-05 -4.958865 -4.958865 -4.958865 -4.9597307 + 22300 5.3899117e-05 -5.1366375 -5.1366375 -5.1366375 -5.1376323 + 22400 5.6210492e-05 -5.3246709 -5.3246709 -5.3246709 -5.3257786 + 22500 5.8495105e-05 -5.5216133 -5.5216133 -5.5216133 -5.5227109 + 22600 6.0812028e-05 -5.7258178 -5.7258178 -5.7258178 -5.7265308 + 22700 6.3143749e-05 -5.9368597 -5.9368597 -5.9368595 -5.9365948 + 22800 6.5410125e-05 -6.1554904 -6.1554904 -6.155484 -6.1532615 + 22900 6.7542278e-05 -6.3835712 -6.3835712 -6.3835479 -6.3779477 + 23000 6.9507997e-05 -6.6216992 -6.6216992 -6.6216707 -6.6104365 + 23100 7.1266388e-05 -6.8682926 -6.8682926 -6.8682125 -6.8496968 + 23200 7.2755292e-05 -7.1203614 -7.1203614 -7.1199808 -7.0940328 + 23300 7.3926113e-05 -7.3743208 -7.3743208 -7.3732718 -7.3415245 + 23400 7.477841e-05 -7.6268907 -7.6268907 -7.624851 -7.5907534 + 23500 7.5367785e-05 -7.8755433 -7.8755433 -7.8722996 -7.840737 + 23600 7.5767316e-05 -8.118477 -8.118477 -8.1140037 -8.0907012 + 23700 7.602765e-05 -8.3545004 -8.3545004 -8.3489055 -8.3399761 + 23800 7.6179748e-05 -8.5831449 -8.5831448 -8.5766169 -8.588118 + 23900 7.6254392e-05 -8.8047933 -8.8047901 -8.7975178 -8.8348184 + 24000 7.6282088e-05 -9.0203998 -9.02038 -9.0124888 -9.0799295 + 24100 7.6273769e-05 -9.2308282 -9.2307693 -9.2224638 -9.3235303 + 24200 7.6210421e-05 -9.4363163 -9.4361182 -9.4281809 -9.5657774 + 24300 7.6054129e-05 -9.6364815 -9.6359364 -9.6298499 -9.8069317 + 24400 7.5770911e-05 -9.8308075 -9.8296133 -9.8273813 -10.047332 + 24500 7.5346158e-05 -10.019212 -10.016965 -10.020848 -10.287277 + 24600 7.4773894e-05 -10.202327 -10.198639 -10.21069 -10.526948 + 24700 7.4050263e-05 -10.381365 -10.375876 -10.397763 -10.766375 + 24800 7.3180356e-05 -10.557735 -10.550055 -10.58316 -11.005505 + 24900 7.2185526e-05 -10.732679 -10.722522 -10.767885 -11.244347 + 25000 7.1106011e-05 -10.907087 -10.894308 -10.952726 -11.483049 + 25100 6.9985212e-05 -11.081419 -11.066111 -11.138362 -11.721873 + 25200 6.8853108e-05 -11.255636 -11.238316 -11.325516 -11.961107 + 25300 6.7725304e-05 -11.429152 -11.410944 -11.514757 -12.201026 + 25400 6.660963e-05 -11.600933 -11.583606 -11.706214 -12.441881 + 25500 6.5511741e-05 -11.769787 -11.755527 -11.899618 -12.683904 + 25600 6.4431312e-05 -11.93471 -11.925846 -12.09448 -12.927272 + 25700 6.3358881e-05 -12.095051 -12.094001 -12.290238 -13.17205 + 25800 6.2279783e-05 -12.250417 -12.259724 -12.486358 -13.41819 + 25900 6.1182597e-05 -12.400395 -12.422931 -12.682434 -13.665588 + 26000 6.0064731e-05 -12.544284 -12.583582 -12.878121 -13.914096 + 26100 5.8931468e-05 -12.681003 -12.74145 -13.072993 -14.163547 + 26200 5.7791958e-05 -12.809231 -12.895994 -13.266591 -14.413818 + 26300 5.6657846e-05 -12.92764 -13.04641 -13.458453 -14.664819 + 26400 5.5543753e-05 -13.035138 -13.191802 -13.648055 -14.916483 + 26500 5.4465863e-05 -13.13102 -13.331331 -13.834847 -15.168771 + 26600 5.3437869e-05 -13.215073 -13.464305 -14.018282 -15.421628 + 26700 5.2467582e-05 -13.287603 -13.590164 -14.197897 -15.674928 + 26800 5.1557173e-05 -13.34938 -13.708555 -14.373328 -15.928477 + 26900 5.0705883e-05 -13.40147 -13.819378 -14.544261 -16.182031 + 27000 4.9911134e-05 -13.444985 -13.922817 -14.710472 -16.435349 + 27100 4.9168609e-05 -13.480839 -14.019362 -14.871829 -16.688246 + 27200 4.8472969e-05 -13.50963 -14.109672 -15.028327 -16.940585 + 27300 4.7820336e-05 -13.53168 -14.194403 -15.180117 -17.192235 + 27400 4.7210393e-05 -13.547187 -14.274111 -15.327496 -17.443055 + 27500 4.6646131e-05 -13.556406 -14.349198 -15.470825 -17.692916 + 27600 4.6131822e-05 -13.55974 -14.419959 -15.61046 -17.941734 + 27700 4.5671463e-05 -13.557748 -14.486644 -15.746709 -18.189438 + 27800 4.526841e-05 -13.55106 -14.549496 -15.879832 -18.435965 + 27900 4.4925027e-05 -13.540261 -14.608781 -16.01005 -18.681274 + 28000 4.4641149e-05 -13.525772 -14.664791 -16.137576 -18.925336 + 28100 4.4413443e-05 -13.507739 -14.717881 -16.262673 -19.168131 + 28200 4.4236093e-05 -13.485997 -14.768473 -16.385611 -19.409639 + 28300 4.4103038e-05 -13.460087 -14.816971 -16.506607 -19.649817 + 28400 4.4009677e-05 -13.429366 -14.863676 -16.625875 -19.88859 + 28500 4.395268e-05 -13.393145 -14.908709 -16.743629 -20.125895 + 28600 4.3929199e-05 -13.35081 -14.952007 -16.860112 -20.361701 + 28700 4.3931357e-05 -13.301876 -14.99336 -16.97559 -20.596003 + 28800 4.3958849e-05 -13.245964 -15.03248 -17.090288 -20.828842 + 28900 4.4017767e-05 -13.182747 -15.069077 -17.204323 -21.060354 + 29000 4.4112107e-05 -13.111905 -15.102891 -17.317673 -21.290863 + 29100 4.4238616e-05 -13.033134 -15.133697 -17.430203 -21.520613 + 29200 4.437609e-05 -12.946196 -15.161312 -17.541762 -21.749803 + 29300 4.4521265e-05 -12.851002 -15.185556 -17.652227 -21.97852 + 29400 4.4663842e-05 -12.747688 -15.206178 -17.761476 -22.206898 + 29500 4.4799305e-05 -12.63667 -15.222853 -17.869356 -22.435207 + 29600 4.4928852e-05 -12.518653 -15.235188 -17.975691 -22.66335 + 29700 4.5048862e-05 -12.394588 -15.242741 -18.080302 -22.89138 + 29800 4.5151469e-05 -12.265571 -15.245083 -18.183038 -23.119178 + 29900 4.5238694e-05 -12.132715 -15.241915 -18.283755 -23.346422 + 30000 4.5310976e-05 -11.997019 -15.23315 -18.382297 -23.572939 + 30100 4.5367895e-05 -11.859285 -15.218941 -18.478497 -23.798746 + 30200 4.5409758e-05 -11.720105 -15.199629 -18.572189 -24.023589 + 30300 4.5437914e-05 -11.579914 -15.175689 -18.663215 -24.246982 + 30400 4.5454781e-05 -11.439079 -15.147661 -18.751468 -24.468408 + 30500 4.546054e-05 -11.297967 -15.116074 -18.8369 -24.687405 + 30600 4.5448468e-05 -11.156981 -15.08141 -18.919516 -24.903683 + 30700 4.5410939e-05 -11.016532 -15.044084 -18.999368 -25.117042 + 30800 4.5339487e-05 -10.876969 -15.004442 -19.076532 -25.32748 + 30900 4.5224936e-05 -10.7385 -14.962771 -19.151114 -25.53528 + 31000 4.5064834e-05 -10.601114 -14.919327 -19.223245 -25.740748 + 31100 4.4858022e-05 -10.464556 -14.874359 -19.293077 -25.944182 + 31200 4.4603577e-05 -10.328349 -14.82811 -19.360777 -26.145888 + 31300 4.4302138e-05 -10.191878 -14.780786 -19.426527 -26.34594 + 31400 4.3955895e-05 -10.054501 -14.732515 -19.490517 -26.544193 + 31500 4.356891e-05 -9.9156707 -14.683313 -19.552928 -26.740381 + 31600 4.3145785e-05 -9.7750191 -14.633087 -19.613944 -26.934353 + 31700 4.2689819e-05 -9.6323955 -14.581655 -19.67376 -27.126104 + 31800 4.2204394e-05 -9.4878602 -14.528796 -19.73257 -27.31555 + 31900 4.169173e-05 -9.3416467 -14.474281 -19.790543 -27.502505 + 32000 4.1153109e-05 -9.1941157 -14.417899 -19.847789 -27.686714 + 32100 4.0588966e-05 -9.0457218 -14.359461 -19.904349 -27.86792 + 32200 3.9998911e-05 -8.8969982 -14.298815 -19.960257 -28.04599 + 32300 3.9382527e-05 -8.7485603 -14.235852 -20.01557 -28.220844 + 32400 3.8739886e-05 -8.6011166 -14.170496 -20.070338 -28.392428 + 32500 3.807179e-05 -8.455473 -14.102695 -20.124582 -28.560692 + 32600 3.7380597e-05 -8.3125174 -14.032425 -20.178286 -28.725584 + 32700 3.6669507e-05 -8.173176 -13.959714 -20.231415 -28.887079 + 32800 3.5942269e-05 -8.0383455 -13.884682 -20.283914 -29.045243 + 32900 3.5203089e-05 -7.9088151 -13.807571 -20.335713 -29.200239 + 33000 3.4455796e-05 -7.7851998 -13.728756 -20.386736 -29.352283 + 33100 3.3703315e-05 -7.6679015 -13.648728 -20.436897 -29.501552 + 33200 3.2947524e-05 -7.5571093 -13.568052 -20.486094 -29.648164 + 33300 3.2188895e-05 -7.4528303 -13.487314 -20.534215 -29.792195 + 33400 3.1427363e-05 -7.3549392 -13.407076 -20.581164 -29.933685 + 33500 3.0663873e-05 -7.2632256 -13.327826 -20.626861 -30.072653 + 33600 2.9898968e-05 -7.1774267 -13.249955 -20.671239 -30.20916 + 33700 2.9133275e-05 -7.0972388 -13.173745 -20.714241 -30.343255 + 33800 2.8368812e-05 -7.0223118 -13.099383 -20.755828 -30.47491 + 33900 2.7609003e-05 -6.9522386 -13.026985 -20.79601 -30.604079 + 34000 2.6857675e-05 -6.8865525 -12.956623 -20.834836 -30.730759 + 34100 2.6118839e-05 -6.8247453 -12.888336 -20.87237 -30.854988 + 34200 2.5396622e-05 -6.7663094 -12.822128 -20.908707 -30.97674 + 34300 2.469467e-05 -6.710799 -12.757955 -20.943955 -31.095936 + 34400 2.4015785e-05 -6.6578961 -12.695711 -20.978205 -31.212492 + 34500 2.3351522e-05 -6.6074644 -12.635225 -21.011508 -31.326423 + 34600 2.2693793e-05 -6.5595753 -12.576275 -21.043901 -31.437755 + 34700 2.2042845e-05 -6.5145006 -12.518607 -21.075419 -31.546556 + 34800 2.1400448e-05 -6.4726739 -12.461972 -21.106109 -31.654648 + 34900 2.076709e-05 -6.4346333 -12.406153 -21.136024 -31.762833 + 35000 2.0143611e-05 -6.4009581 -12.350999 -21.165239 -31.872069 + 35100 1.9532155e-05 -6.3722162 -12.296436 -21.193882 -31.979576 + 35200 1.8935937e-05 -6.3489286 -12.242478 -21.222151 -32.084957 + 35300 1.8358185e-05 -6.331553 -12.18921 -21.250135 -32.188072 + 35400 1.7801738e-05 -6.3204815 -12.136783 -21.277903 -32.289118 + 35500 1.7266114e-05 -6.3160416 -12.085391 -21.305528 -32.38804 + 35600 1.675254e-05 -6.318492 -12.035274 -21.333023 -32.484516 + 35700 1.6263012e-05 -6.3280071 -11.986708 -21.360341 -32.578294 + 35800 1.5799035e-05 -6.3446531 -11.940009 -21.387371 -32.669409 + 35900 1.5360163e-05 -6.368361 -11.895533 -21.413946 -32.757946 + 36000 1.4944282e-05 -6.3989087 -11.853657 -21.439897 -32.843827 + 36100 1.454895e-05 -6.4359198 -11.814763 -21.465082 -32.926809 + 36200 1.4172256e-05 -6.4788836 -11.779207 -21.489395 -33.00665 + 36300 1.3812971e-05 -6.5271927 -11.747285 -21.512758 -33.083362 + 36400 1.3469641e-05 -6.5801919 -11.719208 -21.535122 -33.157216 + 36500 1.3140865e-05 -6.6372266 -11.695081 -21.556456 -33.228519 + 36600 1.282636e-05 -6.6976818 -11.674904 -21.576758 -33.297421 + 36700 1.2526815e-05 -6.7610085 -11.65858 -21.596045 -33.363994 + 36800 1.2242856e-05 -6.8267339 -11.645938 -21.614336 -33.428392 + 36900 1.1974517e-05 -6.8944642 -11.636757 -21.631625 -33.490797 + 37000 1.1721571e-05 -6.9638828 -11.630788 -21.647883 -33.551316 + 37100 1.1483256e-05 -7.0347523 -11.62777 -21.66307 -33.610018 + 37200 1.1258337e-05 -7.1069211 -11.627437 -21.67715 -33.666966 + 37300 1.1045047e-05 -7.180334 -11.629521 -21.690113 -33.722209 + 37400 1.0841233e-05 -7.2550403 -11.633755 -21.701985 -33.775826 + 37500 1.0644622e-05 -7.3311938 -11.639882 -21.712812 -33.828108 + 37600 1.0453073e-05 -7.4090393 -11.647665 -21.722656 -33.879698 + 37700 1.026458e-05 -7.4888842 -11.656902 -21.731577 -33.92968 + 37800 1.0078232e-05 -7.5710581 -11.667446 -21.739637 -33.978107 + 37900 9.8937342e-06 -7.6558676 -11.679216 -21.746914 -34.025034 + 38000 9.7112276e-06 -7.7435547 -11.692204 -21.753491 -34.070539 + 38100 9.531055e-06 -7.834266 -11.706475 -21.759429 -34.114628 + 38200 9.3534226e-06 -7.9280367 -11.722153 -21.764777 -34.15724 + 38300 9.1781939e-06 -8.02479 -11.739403 -21.769582 -34.19834 + 38400 9.0048139e-06 -8.1243473 -11.758418 -21.773899 -34.237904 + 38500 8.8323047e-06 -8.2264447 -11.779393 -21.777796 -34.275926 + 38600 8.6593378e-06 -8.3307479 -11.802514 -21.781335 -34.312426 + 38700 8.4842616e-06 -8.4368643 -11.827949 -21.784562 -34.34743 + 38800 8.3052026e-06 -8.544351 -11.855839 -21.787506 -34.381302 + 38900 8.1209916e-06 -8.6527202 -11.886292 -21.790172 -34.414268 + 39000 7.9309937e-06 -8.7614486 -11.919374 -21.79254 -34.445847 + 39100 7.7350519e-06 -8.8699912 -11.955104 -21.79457 -34.476114 + 39200 7.5334096e-06 -8.9778021 -11.99344 -21.796214 -34.505084 + 39300 7.3266446e-06 -9.0843611 -12.034275 -21.797414 -34.532715 + 39400 7.1155426e-06 -9.1892007 -12.077433 -21.79811 -34.558975 + 39500 6.9009068e-06 -9.2919305 -12.122675 -21.798242 -34.583889 + 39600 6.683376e-06 -9.3922532 -12.169711 -21.797759 -34.607507 + 39700 6.4633394e-06 -9.4899717 -12.218216 -21.796601 -34.62988 + 39800 6.2409324e-06 -9.5849851 -12.26785 -21.794705 -34.65104 + 39900 6.0160834e-06 -9.6772783 -12.318281 -21.792015 -34.671019 + 40000 5.788596e-06 -9.7669066 -12.369199 -21.788497 -34.689868 + 40100 5.5582747e-06 -9.8539797 -12.420331 -21.784138 -34.707653 + 40200 5.325066e-06 -9.9386467 -12.47144 -21.778948 -34.724427 + 40300 5.0891883e-06 -10.021084 -12.522336 -21.772958 -34.740229 + 40400 4.8512027e-06 -10.10148 -12.572868 -21.766214 -34.755086 + 40500 4.6120028e-06 -10.180028 -12.622928 -21.758781 -34.769037 + 40600 4.3727307e-06 -10.256904 -12.67245 -21.750735 -34.782131 + 40700 4.1346564e-06 -10.332254 -12.721408 -21.742158 -34.794412 + 40800 3.8990492e-06 -10.406178 -12.769818 -21.733127 -34.80591 + 40900 3.6670602e-06 -10.47871 -12.81773 -21.723713 -34.816648 + 41000 3.4396317e-06 -10.54981 -12.865232 -21.713989 -34.826651 + 41100 3.2174552e-06 -10.619362 -12.912429 -21.704028 -34.835948 + 41200 3.0009848e-06 -10.687175 -12.95944 -21.693906 -34.844567 + 41300 2.7905083e-06 -10.753 -13.006379 -21.683694 -34.852523 + 41400 2.5862488e-06 -10.816542 -13.053344 -21.673446 -34.859823 + 41500 2.3884725e-06 -10.877482 -13.100405 -21.6632 -34.866469 + 41600 2.1975779e-06 -10.935493 -13.147596 -21.652968 -34.87247 + 41700 2.0141477e-06 -10.990252 -13.194912 -21.642745 -34.877833 + 41800 1.8389545e-06 -11.041457 -13.242307 -21.632508 -34.882566 + 41900 1.6729135e-06 -11.088833 -13.289693 -21.622218 -34.886671 + 42000 1.5169969e-06 -11.13214 -13.336945 -21.611829 -34.890156 + 42100 1.3721296e-06 -11.171183 -13.383902 -21.60129 -34.893022 + 42200 1.2390898e-06 -11.205818 -13.430374 -21.590558 -34.895275 + 42300 1.1184368e-06 -11.235964 -13.476143 -21.579594 -34.896914 + 42400 1.0104789e-06 -11.261605 -13.520974 -21.568371 -34.897933 + 42500 9.1528123e-07 -11.282793 -13.564626 -21.556865 -34.898324 + 42600 8.3271128e-07 -11.299646 -13.606862 -21.545064 -34.898075 + 42700 7.6250806e-07 -11.312343 -13.647464 -21.532961 -34.89718 + 42800 7.0436072e-07 -11.321105 -13.686241 -21.520567 -34.89563 + 42900 6.5798437e-07 -11.326185 -13.723044 -21.50791 -34.893406 + 43000 6.2317463e-07 -11.327851 -13.757765 -21.495036 -34.890484 + 43100 5.99827e-07 -11.326369 -13.790344 -21.482008 -34.886839 + 43200 5.8791483e-07 -11.32199 -13.820764 -21.4689 -34.88247 + 43300 5.8743308e-07 -11.314943 -13.849045 -21.455795 -34.877393 + 43400 5.9832137e-07 -11.305423 -13.875238 -21.442777 -34.871632 + 43500 6.2038589e-07 -11.293589 -13.89942 -21.429927 -34.865214 + 43600 6.5324191e-07 -11.279557 -13.921686 -21.417316 -34.858184 + 43700 6.962899e-07 -11.263395 -13.942143 -21.405008 -34.850595 + 43800 7.48728e-07 -11.245126 -13.960907 -21.393055 -34.842505 + 43900 8.095946e-07 -11.224719 -13.978097 -21.381504 -34.833968 + 44000 8.7783193e-07 -11.2021 -13.99383 -21.370389 -34.825024 + 44100 9.5235889e-07 -11.177157 -14.008212 -21.359738 -34.815711 + 44200 1.0321412e-06 -11.149748 -14.021335 -21.349566 -34.806057 + 44300 1.1162454e-06 -11.119718 -14.033267 -21.339877 -34.796076 + 44400 1.2038703e-06 -11.086917 -14.044049 -21.330656 -34.785763 + 44500 1.2943411e-06 -11.051211 -14.05369 -21.321878 -34.775106 + 44600 1.3870678e-06 -11.012496 -14.062168 -21.313504 -34.764089 + 44700 1.4814963e-06 -10.970711 -14.069429 -21.305491 -34.752697 + 44800 1.5770468e-06 -10.925843 -14.075397 -21.297788 -34.740916 + 44900 1.6730595e-06 -10.877929 -14.079972 -21.290349 -34.72874 + 45000 1.7687641e-06 -10.827058 -14.083043 -21.283132 -34.716176 + 45100 1.8632791e-06 -10.773368 -14.084492 -21.276104 -34.703246 + 45200 1.9556417e-06 -10.717045 -14.084203 -21.269242 -34.689988 + 45300 2.0448667e-06 -10.658317 -14.082067 -21.262535 -34.676447 + 45400 2.1300137e-06 -10.59745 -14.077992 -21.255984 -34.662668 + 45500 2.210249e-06 -10.534738 -14.071906 -21.249602 -34.648696 + 45600 2.284911e-06 -10.470497 -14.063769 -21.243414 -34.634575 + 45700 2.3535289e-06 -10.40505 -14.053569 -21.237457 -34.620345 + 45800 2.4157755e-06 -10.338716 -14.041333 -21.231779 -34.60604 + 45900 2.4714568e-06 -10.271796 -14.02712 -21.226434 -34.591699 + 46000 2.5204501e-06 -10.204564 -14.011022 -21.221479 -34.57735 + 46100 2.5626655e-06 -10.137255 -13.993158 -21.216974 -34.563002 + 46200 2.598002e-06 -10.070062 -13.973667 -21.212978 -34.548664 + 46300 2.6263075e-06 -10.003128 -13.952702 -21.209538 -34.534347 + 46400 2.6473967e-06 -9.9365471 -13.930421 -21.206696 -34.520073 + 46500 2.6610253e-06 -9.8703701 -13.906984 -21.204481 -34.505864 + 46600 2.6669427e-06 -9.8046037 -13.882544 -21.20291 -34.491745 + 46700 2.6649472e-06 -9.7392186 -13.857245 -21.201985 -34.477755 + 46800 2.6549354e-06 -9.6741563 -13.831219 -21.201699 -34.463941 + 46900 2.6369413e-06 -9.6093378 -13.80458 -21.20203 -34.450356 + 47000 2.6111538e-06 -9.5446744 -13.777424 -21.202947 -34.437043 + 47100 2.5779103e-06 -9.4800794 -13.749827 -21.20441 -34.424039 + 47200 2.5376682e-06 -9.4154798 -13.721839 -21.206368 -34.411379 + 47300 2.4909592e-06 -9.3508289 -13.693486 -21.208768 -34.399098 + 47400 2.4383406e-06 -9.2861156 -13.664771 -21.211552 -34.387218 + 47500 2.3803377e-06 -9.2213727 -13.635675 -21.214662 -34.375757 + 47600 2.3174227e-06 -9.1566799 -13.60616 -21.218044 -34.364719 + 47700 2.2500069e-06 -9.0921645 -13.57618 -21.221648 -34.354101 + 47800 2.1784399e-06 -9.0279971 -13.54568 -21.225434 -34.343896 + 47900 2.1030418e-06 -8.9643856 -13.51461 -21.229373 -34.334101 + 48000 2.0241445e-06 -8.9015668 -13.48293 -21.233446 -34.324715 + 48100 1.9421273e-06 -8.8397967 -13.450613 -21.237647 -34.315737 + 48200 1.8574451e-06 -8.7793411 -13.417654 -21.241982 -34.307181 + 48300 1.7706403e-06 -8.720466 -13.384069 -21.246467 -34.299069 + 48400 1.6823382e-06 -8.6634283 -13.349902 -21.251127 -34.291431 + 48500 1.593223e-06 -8.6084673 -13.315221 -21.255993 -34.284298 + 48600 1.5039994e-06 -8.5557957 -13.280119 -21.261099 -34.277697 + 48700 1.4153491e-06 -8.5055928 -13.24471 -21.266481 -34.271649 + 48800 1.3278916e-06 -8.4579971 -13.209128 -21.272173 -34.266175 + 48900 1.2421559e-06 -8.4131019 -13.173518 -21.278201 -34.261292 + 49000 1.1585716e-06 -8.3709529 -13.138035 -21.284589 -34.257347 + 49100 1.0774783e-06 -8.331548 -13.102833 -21.291349 -34.254075 + 49200 9.991497e-07 -8.2948411 -13.068063 -21.298484 -34.251377 + 49300 9.2382685e-07 -8.2607479 -13.033864 -21.305986 -34.24924 + 49400 8.5175121e-07 -8.2291539 -13.000361 -21.313836 -34.247649 + 49500 7.8318971e-07 -8.1999246 -12.967661 -21.322004 -34.246593 + 49600 7.1844633e-07 -8.1729156 -12.935851 -21.330453 -34.246057 + 49700 6.5785849e-07 -8.1479832 -12.904995 -21.339136 -34.246033 + 49800 6.01779e-07 -8.1249945 -12.875137 -21.348002 -34.246515 + 49900 5.5054733e-07 -8.1038357 -12.8463 -21.356994 -34.2475 + 50000 5.0445677e-07 -8.0844203 -12.818485 -21.366056 -34.24899 + 50100 4.6372437e-07 -8.066694 -12.791677 -21.375131 -34.250986 + 50200 4.2847076e-07 -8.0506387 -12.765845 -21.384166 -34.253491 + 50300 3.9871392e-07 -8.036274 -12.740947 -21.393111 -34.256503 + 50400 3.7437835e-07 -8.0236554 -12.716933 -21.401923 -34.260017 + 50500 3.553174e-07 -8.0128711 -12.693748 -21.410567 -34.264019 + 50600 3.413435e-07 -8.0040356 -12.671343 -21.419017 -34.268492 + 50700 3.3225923e-07 -7.9972817 -12.649675 -21.427258 -34.273416 + 50800 3.2788239e-07 -7.9927513 -12.628714 -21.435284 -34.278907 + 50900 3.280597e-07 -7.990585 -12.608448 -21.443102 -34.285043 + 51000 3.326667e-07 -7.9909133 -12.588883 -21.450724 -34.291557 + 51100 3.4159474e-07 -7.9938474 -12.570045 -21.458173 -34.298424 + 51200 3.5472924e-07 -7.9994724 -12.551979 -21.465475 -34.305625 + 51300 3.7192462e-07 -8.0078416 -12.534748 -21.472659 -34.313143 + 51400 3.9298227e-07 -8.0189729 -12.518427 -21.479754 -34.320966 + 51500 4.176367e-07 -8.0328468 -12.503104 -21.486787 -34.329084 + 51600 4.455534e-07 -8.0494058 -12.488869 -21.493779 -34.337484 + 51700 4.7633879e-07 -8.0685554 -12.475814 -21.500746 -34.346158 + 51800 5.0956009e-07 -8.0901672 -12.464026 -21.507696 -34.355093 + 51900 5.4477085e-07 -8.1140827 -12.453585 -21.514629 -34.364276 + 52000 5.815338e-07 -8.140119 -12.444556 -21.521535 -34.373691 + 52100 6.1944112e-07 -8.1680767 -12.436991 -21.528395 -34.38332 + 52200 6.5813507e-07 -8.1977473 -12.430921 -21.535182 -34.393142 + 52300 6.9730824e-07 -8.228923 -12.426359 -21.541863 -34.403134 + 52400 7.366909e-07 -8.2614049 -12.423298 -21.548397 -34.41327 + 52500 7.7603538e-07 -8.2950118 -12.42171 -21.554741 -34.423513 + 52600 8.1510322e-07 -8.3295868 -12.421551 -21.560851 -34.433831 + 52700 8.536414e-07 -8.3650023 -12.422763 -21.566681 -34.444194 + 52800 8.9136981e-07 -8.4011628 -12.425273 -21.572191 -34.454578 + 52900 9.2797815e-07 -8.4380059 -12.429002 -21.577344 -34.464966 + 53000 9.6313291e-07 -8.4755004 -12.433865 -21.58211 -34.475334 + 53100 9.9649361e-07 -8.5136428 -12.439775 -21.586464 -34.485662 + 53200 1.0277337e-06 -8.5524529 -12.446649 -21.590392 -34.495934 + 53300 1.0565618e-06 -8.5919672 -12.45441 -21.593887 -34.506139 + 53400 1.0827383e-06 -8.6322321 -12.462988 -21.59695 -34.516258 + 53500 1.1060844e-06 -8.6732967 -12.472326 -21.599594 -34.526272 + 53600 1.1264823e-06 -8.7152053 -12.482383 -21.601835 -34.536161 + 53700 1.1438662e-06 -8.7579904 -12.49313 -21.6037 -34.545904 + 53800 1.1582004e-06 -8.8016669 -12.504554 -21.605217 -34.555481 + 53900 1.1694765e-06 -8.846227 -12.516657 -21.606421 -34.564872 + 54000 1.1776909e-06 -8.8916365 -12.529452 -21.607345 -34.574062 + 54100 1.1828304e-06 -8.9378335 -12.542962 -21.608025 -34.583039 + 54200 1.1848722e-06 -8.9847277 -12.557216 -21.608491 -34.591791 + 54300 1.1837845e-06 -9.0322019 -12.572243 -21.608769 -34.600309 + 54400 1.1795365e-06 -9.0801152 -12.588069 -21.608883 -34.608583 + 54500 1.1721131e-06 -9.1283068 -12.604716 -21.608845 -34.6166 + 54600 1.1615269e-06 -9.1766005 -12.622195 -21.608664 -34.624352 + 54700 1.147836e-06 -9.2248111 -12.640506 -21.608339 -34.631834 + 54800 1.1311509e-06 -9.2727499 -12.659636 -21.607863 -34.639044 + 54900 1.1116306e-06 -9.3202309 -12.679556 -21.607224 -34.645983 + 55000 1.0894686e-06 -9.3670775 -12.700224 -21.606401 -34.652644 + 55100 1.0648885e-06 -9.413128 -12.721582 -21.605372 -34.659018 + 55200 1.0381271e-06 -9.4582409 -12.743563 -21.604112 -34.665094 + 55300 1.0094209e-06 -9.5022994 -12.766086 -21.602595 -34.670866 + 55400 9.789906e-07 -9.545214 -12.789064 -21.600794 -34.676323 + 55500 9.4703646e-07 -9.5869241 -12.812403 -21.598688 -34.681447 + 55600 9.137395e-07 -9.6273972 -12.836007 -21.596259 -34.686219 + 55700 8.7926831e-07 -9.666627 -12.859783 -21.593493 -34.690631 + 55800 8.4378974e-07 -9.7046295 -12.883639 -21.590385 -34.694678 + 55900 8.0748027e-07 -9.7414384 -12.907492 -21.586937 -34.698356 + 56000 7.705351e-07 -9.7770988 -12.931267 -21.583156 -34.701659 + 56100 7.3317279e-07 -9.811662 -12.954898 -21.579059 -34.704589 + 56200 6.9563428e-07 -9.8451789 -12.978333 -21.574669 -34.707151 + 56300 6.5817669e-07 -9.8776945 -13.001532 -21.570015 -34.709352 + 56400 6.2106264e-07 -9.9092437 -13.024468 -21.56513 -34.711197 + 56500 5.8454726e-07 -9.9398467 -13.047123 -21.560052 -34.712689 + 56600 5.4886558e-07 -9.9695068 -13.069493 -21.55482 -34.713835 + 56700 5.1422298e-07 -9.9982083 -13.091579 -21.549474 -34.714641 + 56800 4.8079053e-07 -10.025916 -13.11339 -21.54405 -34.715114 + 56900 4.4870591e-07 -10.052575 -13.134939 -21.538584 -34.715255 + 57000 4.1807927e-07 -10.078115 -13.156238 -21.533107 -34.715066 + 57100 3.8900229e-07 -10.102449 -13.177298 -21.527644 -34.714547 + 57200 3.6155842e-07 -10.125481 -13.198124 -21.522212 -34.713698 + 57300 3.3583144e-07 -10.147105 -13.218714 -21.516824 -34.712518 + 57400 3.1191058e-07 -10.167217 -13.239058 -21.511485 -34.711007 + 57500 2.8989098e-07 -10.185712 -13.259135 -21.506192 -34.709168 + 57600 2.6986966e-07 -10.202495 -13.278915 -21.500938 -34.707008 + 57700 2.5193814e-07 -10.217483 -13.298357 -21.495711 -34.704539 + 57800 2.3617333e-07 -10.230608 -13.317413 -21.490495 -34.701772 + 57900 2.2262857e-07 -10.241822 -13.336026 -21.485272 -34.69872 + 58000 2.1132698e-07 -10.2511 -13.354136 -21.480025 -34.695395 + 58100 2.022583e-07 -10.258438 -13.371679 -21.474736 -34.691808 + 58200 1.9538022e-07 -10.263855 -13.388591 -21.469393 -34.687968 + 58300 1.9062358e-07 -10.267394 -13.40481 -21.463984 -34.683885 + 58400 1.8790049e-07 -10.269116 -13.420277 -21.458504 -34.679564 + 58500 1.8711356e-07 -10.269097 -13.434942 -21.452956 -34.675013 + 58600 1.8816432e-07 -10.267428 -13.448761 -21.447344 -34.670241 + 58700 1.9095913e-07 -10.264207 -13.461699 -21.441683 -34.665255 + 58800 1.9541341e-07 -10.259532 -13.473735 -21.43599 -34.660062 + 58900 2.0144991e-07 -10.253502 -13.484854 -21.430291 -34.65467 + 59000 2.0899126e-07 -10.246208 -13.495055 -21.424611 -34.649091 + 59100 2.1795238e-07 -10.237732 -13.504345 -21.418984 -34.643339 + 59200 2.2823416e-07 -10.22814 -13.512741 -21.41344 -34.637428 + 59300 2.397203e-07 -10.217487 -13.520266 -21.408014 -34.631375 + 59400 2.5227554e-07 -10.205808 -13.52695 -21.402738 -34.625192 + 59500 2.6574729e-07 -10.193124 -13.532826 -21.39764 -34.61889 + 59600 2.7996905e-07 -10.179441 -13.537927 -21.392747 -34.612484 + 59700 2.9476817e-07 -10.164751 -13.542288 -21.388082 -34.605988 + 59800 3.099747e-07 -10.149035 -13.545938 -21.383659 -34.599414 + 59900 3.2542651e-07 -10.132269 -13.548905 -21.37949 -34.592773 + 60000 3.4097512e-07 -10.114421 -13.551209 -21.375579 -34.586077 + 60100 3.5648824e-07 -10.095461 -13.552862 -21.371926 -34.579337 + 60200 3.7184901e-07 -10.075365 -13.553873 -21.368522 -34.572563 + 60300 3.8695243e-07 -10.054112 -13.554239 -21.365356 -34.56576 + 60400 4.0170024e-07 -10.031696 -13.553952 -21.362413 -34.558941 + 60500 4.1599766e-07 -10.008124 -13.552998 -21.359675 -34.55212 + 60600 4.2974819e-07 -9.9834177 -13.551359 -21.35712 -34.545314 + 60700 4.428495e-07 -9.9576181 -13.549014 -21.354728 -34.538539 + 60800 4.5519287e-07 -9.9307832 -13.545942 -21.352482 -34.531814 + 60900 4.6666819e-07 -9.9029881 -13.542121 -21.350363 -34.525157 + 61000 4.7716916e-07 -9.8743233 -13.537537 -21.34836 -34.518582 + 61100 4.8659927e-07 -9.8448926 -13.532175 -21.346462 -34.5121 + 61200 4.948774e-07 -9.8148092 -13.526032 -21.344666 -34.505718 + 61300 5.0194163e-07 -9.7841929 -13.519108 -21.342973 -34.499446 + 61400 5.0775067e-07 -9.7531657 -13.511414 -21.341389 -34.493295 + 61500 5.1228302e-07 -9.7218479 -13.50297 -21.339924 -34.487275 + 61600 5.1553383e-07 -9.6903541 -13.493801 -21.338592 -34.481392 + 61700 5.1750997e-07 -9.6587902 -13.483943 -21.337412 -34.475654 + 61800 5.1822469e-07 -9.6272499 -13.473438 -21.336402 -34.470069 + 61900 5.1769313e-07 -9.5958132 -13.462333 -21.335584 -34.464643 + 62000 5.159296e-07 -9.5645442 -13.450681 -21.334978 -34.459382 + 62100 5.1294706e-07 -9.5334913 -13.438535 -21.334602 -34.454293 + 62200 5.0875877e-07 -9.5026877 -13.42595 -21.334472 -34.449386 + 62300 5.0338157e-07 -9.472152 -13.41298 -21.334601 -34.44467 + 62400 4.9684005e-07 -9.4418912 -13.399677 -21.334997 -34.440158 + 62500 4.8917046e-07 -9.4119026 -13.386087 -21.335664 -34.435861 + 62600 4.8042341e-07 -9.3821771 -13.372253 -21.336599 -34.431789 + 62700 4.7066452e-07 -9.3527029 -13.358211 -21.337796 -34.427948 + 62800 4.5997305e-07 -9.3234682 -13.343989 -21.339244 -34.424347 + 62900 4.4843861e-07 -9.2944652 -13.329609 -21.340926 -34.42099 + 63000 4.3615693e-07 -9.2656926 -13.315086 -21.342823 -34.417879 + 63100 4.23225e-07 -9.2371584 -13.300432 -21.344913 -34.415017 + 63200 4.0973696e-07 -9.2088819 -13.285649 -21.347172 -34.412406 + 63300 3.9578127e-07 -9.1808945 -13.270739 -21.349576 -34.410046 + 63400 3.8143977e-07 -9.1532408 -13.255702 -21.3521 -34.407939 + 63500 3.6678858e-07 -9.1259778 -13.240536 -21.354722 -34.406089 + 63600 3.5190033e-07 -9.099174 -13.225242 -21.357422 -34.404497 + 63700 3.3684723e-07 -9.0729073 -13.209823 -21.360183 -34.403165 + 63800 3.2170403e-07 -9.0472626 -13.194287 -21.362992 -34.402096 + 63900 3.0655022e-07 -9.022329 -13.178647 -21.365842 -34.40129 + 64000 2.9147055e-07 -8.9981963 -13.162922 -21.368728 -34.400747 + 64100 2.7655408e-07 -8.9749518 -13.147137 -21.371649 -34.400466 + 64200 2.6189153e-07 -8.952677 -13.131325 -21.374609 -34.400445 + 64300 2.4757175e-07 -8.9314448 -13.115523 -21.377614 -34.400682 + 64400 2.3367792e-07 -8.9113167 -13.099774 -21.380675 -34.401175 + 64500 2.2028422e-07 -8.8923414 -13.084123 -21.3838 -34.401918 + 64600 2.0745373e-07 -8.8745529 -13.06862 -21.387001 -34.402906 + 64700 1.9523782e-07 -8.8579706 -13.053315 -21.390288 -34.404131 + 64800 1.8367711e-07 -8.842599 -13.038256 -21.393669 -34.405585 + 64900 1.7280379e-07 -8.8284291 -13.023493 -21.397152 -34.407258 + 65000 1.6264447e-07 -8.8154396 -13.009068 -21.400739 -34.409144 + 65100 1.5322292e-07 -8.8035997 -12.995021 -21.404432 -34.411233 + 65200 1.4456224e-07 -8.7928711 -12.981387 -21.408227 -34.41352 + 65300 1.3668573e-07 -8.7832113 -12.968193 -21.412116 -34.416 + 65400 1.2961657e-07 -8.7745765 -12.95546 -21.416089 -34.418668 + 65500 1.2337608e-07 -8.7669245 -12.943203 -21.420131 -34.421519 + 65600 1.1798139e-07 -8.7602174 -12.93143 -21.424225 -34.424549 + 65700 1.1344272e-07 -8.7544238 -12.920141 -21.428351 -34.427751 + 65800 1.0976115e-07 -8.7495209 -12.909333 -21.432488 -34.431119 + 65900 1.0692732e-07 -8.7454949 -12.899 -21.436614 -34.434647 + 66000 1.0492138e-07 -8.7423423 -12.88913 -21.440707 -34.438326 + 66100 1.0371413e-07 -8.740069 -12.879713 -21.444746 -34.442147 + 66200 1.032693e-07 -8.73869 -12.870736 -21.448712 -34.446103 + 66300 1.0354625e-07 -8.7382278 -12.862191 -21.452587 -34.450184 + 66400 1.0450274e-07 -8.7387102 -12.854069 -21.45636 -34.45438 + 66500 1.06097e-07 -8.7401685 -12.846368 -21.460018 -34.45868 + 66600 1.0828893e-07 -8.7426347 -12.839088 -21.463556 -34.463074 + 66700 1.1104016e-07 -8.746139 -12.832234 -21.466971 -34.467552 + 66800 1.1431301e-07 -8.7507069 -12.825819 -21.470264 -34.472103 + 66900 1.1806887e-07 -8.7563574 -12.819856 -21.473438 -34.476717 + 67000 1.222663e-07 -8.7631004 -12.814365 -21.476499 -34.481385 + 67100 1.2685945e-07 -8.7709356 -12.809367 -21.479457 -34.486097 + 67200 1.3179721e-07 -8.7798511 -12.804885 -21.482319 -34.490843 + 67300 1.3702336e-07 -8.7898235 -12.800945 -21.485097 -34.495616 + 67400 1.4247769e-07 -8.8008173 -12.797568 -21.487799 -34.500406 + 67500 1.4809796e-07 -8.8127864 -12.794775 -21.490433 -34.505203 + 67600 1.5382219e-07 -8.825675 -12.792586 -21.493006 -34.509995 + 67700 1.5959101e-07 -8.8394198 -12.791013 -21.49552 -34.514774 + 67800 1.653494e-07 -8.8539516 -12.790064 -21.497979 -34.519527 + 67900 1.7104771e-07 -8.8691981 -12.789743 -21.500379 -34.524247 + 68000 1.7664163e-07 -8.8850862 -12.790046 -21.502716 -34.528922 + 68100 1.8209146e-07 -8.9015444 -12.790964 -21.504983 -34.533546 + 68200 1.8736063e-07 -8.9185053 -12.792483 -21.507171 -34.53811 + 68300 1.9241418e-07 -8.9359072 -12.794581 -21.509267 -34.542607 + 68400 1.9721738e-07 -8.953696 -12.797235 -21.511259 -34.547032 + 68500 2.0173492e-07 -8.9718259 -12.800416 -21.513133 -34.551378 + 68600 2.0593036e-07 -8.99026 -12.804095 -21.514874 -34.55564 + 68700 2.0976636e-07 -9.0089704 -12.808238 -21.516471 -34.559812 + 68800 2.1320809e-07 -9.0279373 -12.812815 -21.51791 -34.563891 + 68900 2.1622499e-07 -9.0471484 -12.817796 -21.519181 -34.567872 + 69000 2.1879138e-07 -9.066597 -12.823153 -21.520277 -34.571752 + 69100 2.2088772e-07 -9.0862802 -12.82886 -21.521192 -34.575526 + 69200 2.225012e-07 -9.1061973 -12.834898 -21.521923 -34.579185 + 69300 2.236256e-07 -9.1263478 -12.84125 -21.522472 -34.582726 + 69400 2.2426041e-07 -9.1467288 -12.847903 -21.522842 -34.586141 + 69500 2.2440945e-07 -9.167334 -12.85485 -21.523038 -34.589424 + 69600 2.2407935e-07 -9.1881513 -12.862086 -21.52307 -34.592569 + 69700 2.2327811e-07 -9.2091624 -12.869609 -21.522947 -34.595571 + 69800 2.2201421e-07 -9.230341 -12.877422 -21.522682 -34.598423 + 69900 2.2029621e-07 -9.2516534 -12.885525 -21.522287 -34.601124 + 70000 2.1813311e-07 -9.2730579 -12.893919 -21.521774 -34.603669 +Loop time of 68.4982 on 1 procs for 50000 steps with 9576 atoms + +Performance: 6306730.826 tau/day, 729.946 timesteps/s, 6.990 Matom-step/s +99.7% CPU use with 1 MPI tasks x no OpenMP threads + +MPI task timing breakdown: +Section | min time | avg time | max time |%varavg| %total +--------------------------------------------------------------- +Pair | 6.3143 | 6.3143 | 6.3143 | 0.0 | 9.22 +Bond | 50.992 | 50.992 | 50.992 | 0.0 | 74.44 +Neigh | 1.9091 | 1.9091 | 1.9091 | 0.0 | 2.79 +Comm | 0.12117 | 0.12117 | 0.12117 | 0.0 | 0.18 +Output | 0.043458 | 0.043458 | 0.043458 | 0.0 | 0.06 +Modify | 7.4924 | 7.4924 | 7.4924 | 0.0 | 10.94 +Other | | 1.626 | | | 2.37 + +Nlocal: 9576 ave 9576 max 9576 min +Histogram: 1 0 0 0 0 0 0 0 0 0 +Nghost: 0 ave 0 max 0 min +Histogram: 1 0 0 0 0 0 0 0 0 0 +Neighs: 44768 ave 44768 max 44768 min +Histogram: 1 0 0 0 0 0 0 0 0 0 + +Total # of neighbors = 44768 +Ave neighs/atom = 4.6750209 +Ave special neighs/atom = 11.495405 +Neighbor list builds = 365 +Dangerous builds = 0 + +# remove load + +unfix 3 +run 50000 +Generated 0 of 6 mixed pair_coeff terms from geometric mixing rule +Per MPI rank memory allocation (min/avg/max) = 21.82 | 21.82 | 21.82 Mbytes + Step KinEng c_zmin1 c_zmin2 c_zmin3 c_zmin4 + 70000 2.1813311e-07 -9.2730579 -12.893919 -21.521774 -34.603669 + 70100 2.7416894e-07 -9.2778729 -12.8872 -21.508977 -34.598374 + 70200 6.6335855e-07 -9.2611204 -12.861214 -21.482732 -34.591084 + 70300 1.2627672e-06 -9.2298939 -12.821921 -21.445965 -34.581892 + 70400 1.9905726e-06 -9.1871577 -12.772186 -21.400939 -34.570948 + 70500 2.7955954e-06 -9.1344585 -12.713505 -21.348764 -34.55824 + 70600 3.7193627e-06 -9.072659 -12.646828 -21.290459 -34.544303 + 70700 4.7902891e-06 -9.0025065 -12.573171 -21.227795 -34.531964 + 70800 5.9255152e-06 -8.9247979 -12.493414 -21.161896 -34.524412 + 70900 7.1463035e-06 -8.8402444 -12.408235 -21.092862 -34.520546 + 71000 8.4987987e-06 -8.7500719 -12.318272 -21.02089 -34.5176 + 71100 9.9679997e-06 -8.6529242 -12.222656 -20.944645 -34.515152 + 71200 1.1556155e-05 -8.5474358 -12.119754 -20.862458 -34.513186 + 71300 1.3274773e-05 -8.4349355 -12.010461 -20.774236 -34.510169 + 71400 1.5140577e-05 -8.3200716 -11.89855 -20.680776 -34.503754 + 71500 1.7213748e-05 -8.2066613 -11.787147 -20.583327 -34.492199 + 71600 1.9474236e-05 -8.094834 -11.676235 -20.482879 -34.475507 + 71700 2.1804764e-05 -7.9820249 -11.563587 -20.379511 -34.454603 + 71800 2.4087422e-05 -7.8650707 -11.44673 -20.272786 -34.429938 + 71900 2.6250109e-05 -7.7414839 -11.323935 -20.162217 -34.40189 + 72000 2.8282526e-05 -7.6098137 -11.194489 -20.0479 -34.371761 + 72100 3.0215375e-05 -7.4694383 -11.058334 -19.930379 -34.341565 + 72200 3.2059121e-05 -7.3199949 -10.915298 -19.809722 -34.312965 + 72300 3.3810773e-05 -7.1609696 -10.764694 -19.685224 -34.286642 + 72400 3.5491637e-05 -6.9927121 -10.605812 -19.555971 -34.262596 + 72500 3.7130799e-05 -6.8151209 -10.438646 -19.42141 -34.240774 + 72600 3.8743785e-05 -6.6289755 -10.263622 -19.28124 -34.220996 + 72700 4.0324092e-05 -6.4343589 -10.080529 -19.134916 -34.202367 + 72800 4.1847068e-05 -6.2300723 -9.8881888 -18.981521 -34.183319 + 72900 4.329335e-05 -6.0146378 -9.6853669 -18.820111 -34.162309 + 73000 4.4646752e-05 -5.7877259 -9.4719896 -18.650386 -34.138336 + 73100 4.5868964e-05 -5.5508283 -9.2496659 -18.473222 -34.110917 + 73200 4.6904622e-05 -5.3069204 -9.0213313 -18.290458 -34.079809 + 73300 4.771078e-05 -5.0595313 -8.7903532 -18.104406 -34.045051 + 73400 4.8277744e-05 -4.8126505 -8.559649 -17.917282 -34.00728 + 73500 4.8629391e-05 -4.568759 -8.3312922 -17.730865 -33.967722 + 73600 4.8799887e-05 -4.3288349 -8.1066643 -17.54648 -33.92781 + 73700 4.881458e-05 -4.0939477 -7.8867039 -17.365046 -33.888768 + 73800 4.869771e-05 -3.8644512 -7.6718523 -17.187029 -33.851448 + 73900 4.8480856e-05 -3.639927 -7.4618911 -17.012349 -33.81647 + 74000 4.8192265e-05 -3.4195074 -7.2561496 -16.84047 -33.784115 + 74100 4.7838296e-05 -3.2025614 -7.0541038 -16.670784 -33.75417 + 74200 4.7395899e-05 -2.9892342 -6.855871 -16.503048 -33.725929 + 74300 4.6826781e-05 -2.7803886 -6.662165 -16.337535 -33.698396 + 74400 4.6102731e-05 -2.5929199 -6.4737509 -16.174814 -33.670563 + 74500 4.5212041e-05 -2.5004006 -6.2908018 -16.015347 -33.641579 + 74600 4.415534e-05 -2.5004006 -6.1125837 -15.859213 -33.610792 + 74700 4.2941694e-05 -2.5004006 -5.9376331 -15.706122 -33.577758 + 74800 4.159009e-05 -2.5004006 -5.764227 -15.555559 -33.542369 + 74900 4.0130778e-05 -2.5004006 -5.5908516 -15.407018 -33.505 + 75000 3.8597973e-05 -2.5004006 -5.4164924 -15.259954 -33.46642 + 75100 3.7020341e-05 -2.5004006 -5.2406689 -15.113731 -33.427524 + 75200 3.5415765e-05 -2.5004006 -5.0632786 -14.967729 -33.389216 + 75300 3.3797341e-05 -2.5004006 -4.8848086 -14.821644 -33.352335 + 75400 3.2177997e-05 -2.5004006 -4.7062922 -14.675562 -33.317619 + 75500 3.0568725e-05 -2.5004006 -4.5290262 -14.529776 -33.285516 + 75600 2.8969586e-05 -2.5004006 -4.3543706 -14.384706 -33.25606 + 75700 2.7366145e-05 -2.5004006 -4.1842559 -14.240768 -33.22894 + 75800 2.5737424e-05 -2.5004006 -4.0228204 -14.098359 -33.203651 + 75900 2.4070025e-05 -2.5004006 -3.8657517 -13.957859 -33.179564 + 76000 2.2367906e-05 -2.5004006 -3.7134777 -13.819847 -33.155959 + 76100 2.0650719e-05 -2.5004006 -3.5667989 -13.685129 -33.132123 + 76200 1.8946116e-05 -2.5004006 -3.4267708 -13.554616 -33.107467 + 76300 1.7282439e-05 -2.5004006 -3.3101325 -13.429211 -33.081682 + 76400 1.568469e-05 -2.5004006 -3.2452779 -13.30972 -33.054767 + 76500 1.4171528e-05 -2.5004006 -3.1836568 -13.196845 -33.027017 + 76600 1.2752052e-05 -2.5004006 -3.1258382 -13.091151 -32.998932 + 76700 1.1425314e-05 -2.5004006 -3.0790348 -12.992978 -32.971118 + 76800 1.0184013e-05 -2.5004006 -3.0322214 -12.902374 -32.94421 + 76900 9.0204418e-06 -2.5004006 -2.9951747 -12.81913 -32.918857 + 77000 7.9291293e-06 -2.5004006 -2.9586958 -12.74286 -32.895624 + 77100 6.9053423e-06 -2.5004006 -2.9223822 -12.673103 -32.87484 + 77200 5.943679e-06 -2.5004006 -2.886878 -12.609336 -32.856548 + 77300 5.0389494e-06 -2.5004006 -2.8530124 -12.551002 -32.840573 + 77400 4.1906274e-06 -2.5004006 -2.8216962 -12.49764 -32.826569 + 77500 3.4056552e-06 -2.5004006 -2.7938113 -12.448923 -32.814065 + 77600 2.6964734e-06 -2.5004006 -2.7700595 -12.404564 -32.802575 + 77700 2.0773116e-06 -2.5004006 -2.7509054 -12.364239 -32.791623 + 77800 1.5603947e-06 -2.5004006 -2.7389168 -12.327578 -32.780779 + 77900 1.1536435e-06 -2.5004006 -2.7324424 -12.294218 -32.769736 + 78000 8.5913247e-07 -2.5004006 -2.7298568 -12.263793 -32.758356 + 78100 6.7197913e-07 -2.5004006 -2.7309094 -12.235918 -32.746643 + 78200 5.8074058e-07 -2.5004006 -2.7350041 -12.210245 -32.734744 + 78300 5.7014189e-07 -2.5004006 -2.7418314 -12.186576 -32.722966 + 78400 6.2534729e-07 -2.5004006 -2.7512369 -12.164896 -32.711731 + 78500 7.3549757e-07 -2.5004006 -2.7628012 -12.145326 -32.701569 + 78600 8.9487844e-07 -2.5004006 -2.777618 -12.128128 -32.692881 + 78700 1.1020774e-06 -2.5004006 -2.7960631 -12.113691 -32.685976 + 78800 1.3584927e-06 -2.5004006 -2.8155025 -12.102461 -32.681053 + 78900 1.6672136e-06 -2.5004006 -2.8352023 -12.094917 -32.678147 + 79000 2.0319882e-06 -2.5004006 -2.8544683 -12.091542 -32.677131 + 79100 2.4558006e-06 -2.5004006 -2.8727545 -12.092777 -32.677735 + 79200 2.9393301e-06 -2.5004006 -2.8897361 -12.099054 -32.679571 + 79300 3.4800305e-06 -2.5004006 -2.9053135 -12.110801 -32.68218 + 79400 4.0718633e-06 -2.5004006 -2.9195724 -12.128376 -32.685153 + 79500 4.7054146e-06 -2.5004006 -2.9327397 -12.152035 -32.688217 + 79600 5.3685887e-06 -2.5004006 -2.9451562 -12.181882 -32.691255 + 79700 6.0479312e-06 -2.5004006 -2.9581594 -12.217843 -32.694301 + 79800 6.7304839e-06 -2.5004006 -2.980518 -12.259659 -32.697456 + 79900 7.4061484e-06 -2.5004006 -3.014133 -12.306891 -32.700882 + 80000 8.0693275e-06 -2.5004006 -3.054981 -12.35896 -32.704898 + 80100 8.7189192e-06 -2.5004006 -3.1333394 -12.415209 -32.709901 + 80200 9.3569252e-06 -2.5004006 -3.2342664 -12.474964 -32.716211 + 80300 9.986466e-06 -2.5004006 -3.336316 -12.537602 -32.724009 + 80400 1.0609962e-05 -2.5004006 -3.4389447 -12.602587 -32.733338 + 80500 1.1227761e-05 -2.5004006 -3.5418262 -12.669484 -32.74413 + 80600 1.1837222e-05 -2.5004006 -3.6448003 -12.737933 -32.756209 + 80700 1.2432416e-05 -2.5004006 -3.7478202 -12.807611 -32.769286 + 80800 1.3004796e-05 -2.5004006 -3.850915 -12.878214 -32.782992 + 80900 1.354463e-05 -2.5004006 -3.9542357 -12.949464 -32.796975 + 81000 1.4042611e-05 -2.5004006 -4.0598928 -13.021131 -32.810953 + 81100 1.4491109e-05 -2.5004006 -4.1659576 -13.093069 -32.82474 + 81200 1.4884855e-05 -2.5004006 -4.2727509 -13.165235 -32.838249 + 81300 1.522131e-05 -2.5004006 -4.3806836 -13.237697 -32.85151 + 81400 1.5500821e-05 -2.5004006 -4.490242 -13.310629 -32.864686 + 81500 1.5726329e-05 -2.5004006 -4.601939 -13.384292 -32.878162 + 81600 1.5902495e-05 -2.5004006 -4.7162335 -13.459008 -32.89211 + 81700 1.6034321e-05 -2.5004006 -4.8334403 -13.535111 -32.906782 + 81800 1.6125714e-05 -2.5004006 -4.9536647 -13.612907 -32.922376 + 81900 1.6178445e-05 -2.5004006 -5.0767806 -13.692637 -32.938998 + 82000 1.61917e-05 -2.5004006 -5.2024517 -13.774461 -32.956646 + 82100 1.6162276e-05 -2.5004006 -5.3301785 -13.858461 -32.975202 + 82200 1.6085459e-05 -2.5004006 -5.4593529 -13.944631 -32.99444 + 82300 1.5956479e-05 -2.5010432 -5.5893034 -14.032877 -33.014068 + 82400 1.5772147e-05 -2.5860931 -5.7193287 -14.123002 -33.03378 + 82500 1.5532042e-05 -2.7427755 -5.8487168 -14.214698 -33.053294 + 82600 1.5238792e-05 -2.8967382 -5.9767558 -14.307552 -33.072397 + 82700 1.4897488e-05 -3.0471521 -6.1027397 -14.401066 -33.090962 + 82800 1.4514647e-05 -3.1932042 -6.2259798 -14.494679 -33.108972 + 82900 1.4097212e-05 -3.3341629 -6.3458284 -14.587807 -33.126511 + 83000 1.3651746e-05 -3.4694575 -6.4617181 -14.679885 -33.143746 + 83100 1.3183838e-05 -3.5987454 -6.5732089 -14.770398 -33.160883 + 83200 1.2697643e-05 -3.7219456 -6.6800295 -14.858922 -33.178136 + 83300 1.2195641e-05 -3.8392287 -6.7820952 -14.945141 -33.195681 + 83400 1.1678698e-05 -3.9509667 -6.8794988 -15.028842 -33.213634 + 83500 1.1146457e-05 -4.0576572 -6.9724748 -15.109904 -33.232024 + 83600 1.0598043e-05 -4.1598417 -7.0613513 -15.188277 -33.250789 + 83700 1.0032989e-05 -4.2580386 -7.1464999 -15.263958 -33.269783 + 83800 9.4522503e-06 -4.3527026 -7.228295 -15.336983 -33.2888 + 83900 8.8589236e-06 -4.4442133 -7.3070864 -15.407417 -33.307608 + 84000 8.2583636e-06 -4.5328857 -7.3831864 -15.475355 -33.325977 + 84100 7.6575676e-06 -4.6189859 -7.4568677 -15.540922 -33.343713 + 84200 7.0640356e-06 -4.7027382 -7.5283677 -15.604273 -33.360691 + 84300 6.4845715e-06 -4.7843111 -7.5978881 -15.66559 -33.376874 + 84400 5.9243908e-06 -4.8637864 -7.6655868 -15.725079 -33.392315 + 84500 5.3867541e-06 -4.9411178 -7.7315557 -15.782953 -33.407131 + 84600 4.8731791e-06 -5.0160978 -7.795793 -15.839411 -33.421468 + 84700 4.3838933e-06 -5.0883447 -7.8581775 -15.894627 -33.43551 + 84800 3.9181935e-06 -5.1573188 -7.9184559 -15.948719 -33.449442 + 84900 3.475115e-06 -5.2223663 -7.9762491 -16.001729 -33.463396 + 85000 3.0540118e-06 -5.2827818 -8.0310766 -16.053625 -33.477416 + 85100 2.6550467e-06 -5.337876 -8.0823966 -16.104293 -33.491481 + 85200 2.279518e-06 -5.3870346 -8.1296519 -16.153515 -33.505511 + 85300 1.9299553e-06 -5.4297582 -8.172316 -16.201007 -33.519371 + 85400 1.6099081e-06 -5.4657787 -8.2099325 -16.246444 -33.532875 + 85500 1.3234198e-06 -5.4947923 -8.2421424 -16.289491 -33.545818 + 85600 1.0742562e-06 -5.5166676 -8.2686986 -16.32981 -33.558019 + 85700 8.6509493e-07 -5.5313922 -8.2894677 -16.367077 -33.569358 + 85800 6.9695078e-07 -5.5390521 -8.3044253 -16.401002 -33.579778 + 85900 5.6903348e-07 -5.5398336 -8.3136472 -16.431345 -33.589288 + 86000 4.7904183e-07 -5.5340277 -8.3173016 -16.457935 -33.597956 + 86100 4.2376266e-07 -5.5220333 -8.315641 -16.480674 -33.60591 + 86200 3.9977185e-07 -5.5043487 -8.3089927 -16.499543 -33.613305 + 86300 4.0406617e-07 -5.4815499 -8.2977464 -16.514591 -33.620299 + 86400 4.3451e-07 -5.4542534 -8.2823352 -16.525934 -33.627026 + 86500 4.9002551e-07 -5.423071 -8.2632126 -16.533738 -33.633571 + 86600 5.7050242e-07 -5.3887073 -8.2408268 -16.538213 -33.639963 + 86700 6.7647032e-07 -5.3514757 -8.2155949 -16.539592 -33.646172 + 86800 8.0863336e-07 -5.3116716 -8.1878807 -16.538109 -33.652113 + 86900 9.6738244e-07 -5.2695331 -8.1579793 -16.533984 -33.657659 + 87000 1.1523769e-06 -5.2251923 -8.126109 -16.527423 -33.662666 + 87100 1.3622566e-06 -5.178691 -8.0924113 -16.518618 -33.666997 + 87200 1.5945312e-06 -5.1302188 -8.0569553 -16.507752 -33.670543 + 87300 1.8456841e-06 -5.0795185 -8.0197439 -16.495005 -33.673238 + 87400 2.1114927e-06 -5.0267242 -7.9807189 -16.480544 -33.67507 + 87500 2.3875087e-06 -4.9713658 -7.939764 -16.464519 -33.676082 + 87600 2.6695964e-06 -4.9130837 -7.8967073 -16.447048 -33.676364 + 87700 2.9543991e-06 -4.8515463 -7.8513277 -16.428208 -33.676036 + 87800 3.2396175e-06 -4.7864078 -7.8033655 -16.408021 -33.675235 + 87900 3.5240252e-06 -4.7173428 -7.7525426 -16.386444 -33.674086 + 88000 3.8072131e-06 -4.6440852 -7.6985888 -16.363371 -33.672686 + 88100 4.0891272e-06 -4.5664667 -7.6412733 -16.338636 -33.67109 + 88200 4.3695427e-06 -4.4844452 -7.5804352 -16.312035 -33.669306 + 88300 4.6476362e-06 -4.3981188 -7.5160094 -16.283359 -33.667291 + 88400 4.9217735e-06 -4.3077234 -7.4480414 -16.252411 -33.664962 + 88500 5.18954e-06 -4.2136145 -7.3766897 -16.219033 -33.662216 + 88600 5.4479572e-06 -4.1162415 -7.3022169 -16.183113 -33.658944 + 88700 5.6937964e-06 -4.0161173 -7.2249694 -16.144601 -33.655057 + 88800 5.9239171e-06 -3.9137928 -7.1453545 -16.103504 -33.650491 + 88900 6.1355882e-06 -3.8098358 -7.063816 -16.059893 -33.645227 + 89000 6.3267599e-06 -3.7048168 -6.9808139 -16.013895 -33.639287 + 89100 6.4962463e-06 -3.5992968 -6.896808 -15.965687 -33.632738 + 89200 6.6437691e-06 -3.4938158 -6.8122455 -15.915486 -33.625677 + 89300 6.7698261e-06 -3.3888758 -6.7275506 -15.863545 -33.618215 + 89400 6.8753929e-06 -3.2849207 -6.6431133 -15.810137 -33.610469 + 89500 6.9615167e-06 -3.1823134 -6.5592776 -15.755552 -33.602537 + 89600 7.0289135e-06 -3.0813151 -6.4763283 -15.700075 -33.594488 + 89700 7.0776907e-06 -2.9820715 -6.3944791 -15.643978 -33.586357 + 89800 7.1072828e-06 -2.8846104 -6.3138641 -15.587509 -33.578134 + 89900 7.1166145e-06 -2.789937 -6.2345345 -15.530884 -33.569776 + 90000 7.1044239e-06 -2.6970821 -6.1564634 -15.474287 -33.561217 + 90100 7.0696316e-06 -2.6053965 -6.0795571 -15.417871 -33.552378 + 90200 7.0116529e-06 -2.5146927 -6.0036712 -15.361757 -33.543183 + 90300 6.9305698e-06 -2.5004006 -5.9286302 -15.306035 -33.533579 + 90400 6.827186e-06 -2.5004006 -5.854246 -15.250766 -33.523543 + 90500 6.7029652e-06 -2.5004006 -5.780336 -15.195975 -33.513091 + 90600 6.5598756e-06 -2.5004006 -5.7067378 -15.141655 -33.502268 + 90700 6.4001617e-06 -2.5004006 -5.6333219 -15.087762 -33.491142 + 90800 6.2260766e-06 -2.5004006 -5.5600032 -15.03422 -33.479809 + 90900 6.0396236e-06 -2.5004006 -5.4867516 -14.980925 -33.468369 + 91000 5.8423603e-06 -2.5004006 -5.4136004 -14.927757 -33.456911 + 91100 5.6353156e-06 -2.5004006 -5.3406533 -14.874596 -33.445499 + 91200 5.4190491e-06 -2.5004006 -5.2680875 -14.82134 -33.434174 + 91300 5.1938435e-06 -2.5004006 -5.1961503 -14.76792 -33.422945 + 91400 4.9599708e-06 -2.5004006 -5.1251511 -14.714309 -33.411794 + 91500 4.7179704e-06 -2.5004006 -5.0554463 -14.660539 -33.400678 + 91600 4.4688521e-06 -2.5004006 -4.9874193 -14.606693 -33.389541 + 91700 4.2141788e-06 -2.5004006 -4.9214594 -14.552907 -33.378326 + 91800 3.9560121e-06 -2.5004006 -4.8579398 -14.499361 -33.366992 + 91900 3.6967289e-06 -2.5004006 -4.7971993 -14.44627 -33.355518 + 92000 3.4387974e-06 -2.5004006 -4.7395281 -14.393868 -33.343905 + 92100 3.1845304e-06 -2.5004006 -4.685158 -14.3424 -33.332185 + 92200 2.9358778e-06 -2.5004006 -4.6342572 -14.292112 -33.320407 + 92300 2.6943026e-06 -2.5004006 -4.5869277 -14.243236 -33.308638 + 92400 2.4607576e-06 -2.5004006 -4.5432063 -14.195989 -33.296968 + 92500 2.2357671e-06 -2.5004006 -4.5030653 -14.150561 -33.285483 + 92600 2.0195819e-06 -2.5004006 -4.4664168 -14.107111 -33.27426 + 92700 1.8123641e-06 -2.5004006 -4.4331183 -14.065765 -33.263347 + 92800 1.6143761e-06 -2.5004006 -4.4029809 -14.026607 -33.252764 + 92900 1.4261208e-06 -2.5004006 -4.3757815 -13.989684 -33.242497 + 93000 1.2484102e-06 -2.5004006 -4.3512773 -13.955005 -33.232517 + 93100 1.0823488e-06 -2.5004006 -4.3292232 -13.922543 -33.222778 + 93200 9.2923589e-07 -2.5004006 -4.3093889 -13.892241 -33.213232 + 93300 7.9040001e-07 -2.5004006 -4.291575 -13.864013 -33.203834 + 93400 6.6699612e-07 -2.5004006 -4.275626 -13.837756 -33.19455 + 93500 5.5983172e-07 -2.5004006 -4.2614395 -13.813353 -33.185368 + 93600 4.6924843e-07 -2.5004006 -4.2489704 -13.790679 -33.176299 + 93700 3.9509573e-07 -2.5004006 -4.2382306 -13.769612 -33.167372 + 93800 3.3680487e-07 -2.5004006 -4.2292859 -13.750034 -33.158633 + 93900 2.9354212e-07 -2.5004006 -4.2222491 -13.731846 -33.150141 + 94000 2.6439722e-07 -2.5004006 -4.2172717 -13.714969 -33.141959 + 94100 2.4855516e-07 -2.5004006 -4.2145349 -13.699353 -33.134141 + 94200 2.454113e-07 -2.5004006 -4.2142377 -13.684984 -33.126728 + 94300 2.5461072e-07 -2.5004006 -4.2165856 -13.671886 -33.119744 + 94400 2.7601438e-07 -2.5004006 -4.2217774 -13.660118 -33.113191 + 94500 3.0960954e-07 -2.5004006 -4.2299918 -13.649775 -33.107051 + 94600 3.5538783e-07 -2.5004006 -4.2413741 -13.640981 -33.101291 + 94700 4.1321548e-07 -2.5004006 -4.2560242 -13.633873 -33.09587 + 94800 4.8271828e-07 -2.5004006 -4.2739875 -13.628598 -33.090747 + 94900 5.6320367e-07 -2.5004006 -4.2952488 -13.625294 -33.085887 + 95000 6.5363882e-07 -2.5004006 -4.319731 -13.624087 -33.081268 + 95100 7.5269568e-07 -2.5004006 -4.3472981 -13.625073 -33.076891 + 95200 8.5885868e-07 -2.5004006 -4.3777617 -13.628322 -33.072773 + 95300 9.705736e-07 -2.5004006 -4.410891 -13.63387 -33.06895 + 95400 1.0864014e-06 -2.5004006 -4.4464236 -13.641717 -33.065471 + 95500 1.205137e-06 -2.5004006 -4.4840785 -13.651832 -33.062385 + 95600 1.3258634e-06 -2.5004006 -4.5235672 -13.664151 -33.059738 + 95700 1.4479293e-06 -2.5004006 -4.5646058 -13.678584 -33.057562 + 95800 1.5708652e-06 -2.5004006 -4.6069251 -13.695014 -33.055872 + 95900 1.6942639e-06 -2.5004006 -4.6502806 -13.713305 -33.054666 + 96000 1.8176604e-06 -2.5004006 -4.694461 -13.733304 -33.053923 + 96100 1.9404391e-06 -2.5004006 -4.7392952 -13.75485 -33.053605 + 96200 2.0617855e-06 -2.5004006 -4.7846582 -13.777776 -33.053669 + 96300 2.1806885e-06 -2.5004006 -4.8304729 -13.80192 -33.054066 + 96400 2.2959893e-06 -2.5004006 -4.8767108 -13.827127 -33.054751 + 96500 2.4064669e-06 -2.5004006 -4.9233877 -13.853258 -33.055688 + 96600 2.5109437e-06 -2.5004006 -4.9705574 -13.880193 -33.056855 + 96700 2.6083922e-06 -2.5004006 -5.0183018 -13.907834 -33.058247 + 96800 2.6980208e-06 -2.5004006 -5.0667204 -13.93611 -33.059876 + 96900 2.7793185e-06 -2.5004006 -5.1159176 -13.964974 -33.061765 + 97000 2.8520473e-06 -2.5004006 -5.1659916 -13.994405 -33.063947 + 97100 2.9161831e-06 -2.5004006 -5.2170233 -14.024406 -33.066455 + 97200 2.9718191e-06 -2.5004006 -5.2690671 -14.055001 -33.069316 + 97300 3.0190574e-06 -2.5004006 -5.3221438 -14.086232 -33.072549 + 97400 3.0579167e-06 -2.5004006 -5.3762351 -14.118151 -33.076158 + 97500 3.0882798e-06 -2.5004006 -5.43128 -14.150815 -33.080131 + 97600 3.1098902e-06 -2.5004006 -5.4871738 -14.184277 -33.084443 + 97700 3.1223974e-06 -2.5004006 -5.5437694 -14.218578 -33.089059 + 97800 3.1254339e-06 -2.5004006 -5.6008804 -14.253742 -33.093938 + 97900 3.1187057e-06 -2.5004006 -5.6582882 -14.289767 -33.099039 + 98000 3.1020739e-06 -2.5004006 -5.7157494 -14.326627 -33.104326 + 98100 3.0756099e-06 -2.5064869 -5.7730067 -14.364262 -33.109772 + 98200 3.0396159e-06 -2.5732628 -5.8297997 -14.402591 -33.115365 + 98300 2.9946057e-06 -2.6388327 -5.8858759 -14.441505 -33.121101 + 98400 2.9412519e-06 -2.7033969 -5.9410015 -14.480879 -33.126992 + 98500 2.8803103e-06 -2.7677637 -5.9949693 -14.52057 -33.133054 + 98600 2.8125364e-06 -2.8302739 -6.0476057 -14.560432 -33.139306 + 98700 2.7386116e-06 -2.8908094 -6.0987749 -14.600311 -33.145763 + 98800 2.6590957e-06 -2.9493014 -6.148381 -14.640058 -33.152435 + 98900 2.5744147e-06 -3.0057286 -6.1963679 -14.679528 -33.15932 + 99000 2.484888e-06 -3.060113 -6.2427178 -14.718588 -33.166404 + 99100 2.3907859e-06 -3.1125141 -6.287448 -14.757116 -33.173663 + 99200 2.292403e-06 -3.1630222 -6.3306067 -14.795007 -33.181063 + 99300 2.1901277e-06 -3.2117487 -6.3722665 -14.832176 -33.188564 + 99400 2.0844912e-06 -3.2588163 -6.4125175 -14.868555 -33.196127 + 99500 1.9761812e-06 -3.3043473 -6.4514592 -14.9041 -33.203712 + 99600 1.8660216e-06 -3.3484527 -6.4891914 -14.938785 -33.211292 + 99700 1.7549199e-06 -3.3912217 -6.5258053 -14.972606 -33.218846 + 99800 1.6437984e-06 -3.4327129 -6.5613749 -15.005572 -33.226367 + 99900 1.5335215e-06 -3.4729488 -6.5959502 -15.037711 -33.233858 + 100000 1.4248393e-06 -3.5119118 -6.629551 -15.069057 -33.241331 + 100100 1.3183578e-06 -3.5495445 -6.6621645 -15.099651 -33.248801 + 100200 1.2145333e-06 -3.585752 -6.6937442 -15.129537 -33.256287 + 100300 1.1136966e-06 -3.6204064 -6.7242112 -15.158752 -33.263802 + 100400 1.0160984e-06 -3.6533534 -6.753458 -15.187327 -33.271351 + 100500 9.219661e-07 -3.6844197 -6.7813535 -15.215281 -33.278931 + 100600 8.3155685e-07 -3.7134215 -6.8077502 -15.242614 -33.28653 + 100700 7.4519693e-07 -3.7401734 -6.832491 -15.26931 -33.294128 + 100800 6.6329422e-07 -3.7644978 -6.8554181 -15.295335 -33.301698 + 100900 5.8632697e-07 -3.7862335 -6.8763804 -15.320633 -33.309208 + 101000 5.14808e-07 -3.8052446 -6.8952419 -15.345134 -33.316627 + 101100 4.492345e-07 -3.8214278 -6.9118885 -15.368754 -33.323928 + 101200 3.9002454e-07 -3.8347177 -6.9262338 -15.391402 -33.331093 + 101300 3.3747189e-07 -3.8450908 -6.9382236 -15.412982 -33.33811 + 101400 2.9171997e-07 -3.8525658 -6.9478381 -15.433399 -33.344971 + 101500 2.5276094e-07 -3.8572029 -6.9550922 -15.452566 -33.351676 + 101600 2.2046129e-07 -3.8590989 -6.9600343 -15.470404 -33.358233 + 101700 1.9460679e-07 -3.8583819 -6.9627423 -15.486848 -33.364652 + 101800 1.7495364e-07 -3.855204 -6.9633197 -15.501849 -33.370943 + 101900 1.612756e-07 -3.8497328 -6.9618893 -15.515373 -33.377114 + 102000 1.5339526e-07 -3.8421433 -6.9585869 -15.527402 -33.383166 + 102100 1.511964e-07 -3.8326098 -6.9535549 -15.537937 -33.389093 + 102200 1.5461595e-07 -3.8212973 -6.9469354 -15.546992 -33.39488 + 102300 1.6361336e-07 -3.8083553 -6.938864 -15.554596 -33.400511 + 102400 1.7813048e-07 -3.793911 -6.9294636 -15.56079 -33.405962 + 102500 1.9805311e-07 -3.7780644 -6.9188395 -15.565629 -33.411208 + 102600 2.2317591e-07 -3.7608852 -6.907075 -15.569172 -33.416227 + 102700 2.5318267e-07 -3.7424111 -6.8942283 -15.571486 -33.421007 + 102800 2.8764868e-07 -3.7226486 -6.880331 -15.57264 -33.425539 + 102900 3.2606274e-07 -3.7015752 -6.865388 -15.572701 -33.429824 + 103000 3.6786316e-07 -3.6791446 -6.8493799 -15.571735 -33.43387 + 103100 4.1248375e-07 -3.6552928 -6.8322662 -15.569796 -33.437686 + 103200 4.5939555e-07 -3.6299456 -6.8139904 -15.56693 -33.441284 + 103300 5.0813645e-07 -3.603027 -6.7944864 -15.56317 -33.444677 + 103400 5.5832257e-07 -3.5744669 -6.7736857 -15.558538 -33.447876 + 103500 6.0963997e-07 -3.5442088 -6.7515235 -15.553037 -33.450889 + 103600 6.6181983e-07 -3.5122161 -6.7279464 -15.546662 -33.453718 + 103700 7.1460491e-07 -3.4784776 -6.7029176 -15.539392 -33.456356 + 103800 7.6771598e-07 -3.4430109 -6.676422 -15.531201 -33.458794 + 103900 8.2082628e-07 -3.4058648 -6.6484691 -15.522055 -33.461019 + 104000 8.7354951e-07 -3.3671198 -6.6190958 -15.511918 -33.463017 + 104100 9.2544355e-07 -3.326887 -6.5883657 -15.500756 -33.464774 + 104200 9.7602904e-07 -3.2853057 -6.5563691 -15.48854 -33.466277 + 104300 1.0248188e-06 -3.2425391 -6.5232197 -15.475253 -33.467518 + 104400 1.0713521e-06 -3.1987689 -6.4890516 -15.460886 -33.468495 + 104500 1.1152268e-06 -3.1541887 -6.4540141 -15.445447 -33.469209 + 104600 1.1561221e-06 -3.1089966 -6.4182665 -15.428956 -33.46967 + 104700 1.1938078e-06 -3.063388 -6.3819724 -15.411449 -33.46989 + 104800 1.2281376e-06 -3.0175481 -6.3452929 -15.392974 -33.469883 + 104900 1.2590281e-06 -2.9716456 -6.308382 -15.373594 -33.469663 + 105000 1.2864298e-06 -2.9258274 -6.2713803 -15.353379 -33.469243 + 105100 1.3102962e-06 -2.8802147 -6.2344114 -15.332406 -33.468631 + 105200 1.3305604e-06 -2.8349009 -6.1975781 -15.310761 -33.467829 + 105300 1.3471233e-06 -2.7899508 -6.1609606 -15.288529 -33.466835 + 105400 1.3598556e-06 -2.7454017 -6.1246153 -15.265797 -33.465645 + 105500 1.3686136e-06 -2.7012659 -6.0885754 -15.242649 -33.464254 + 105600 1.3732626e-06 -2.6585488 -6.0528524 -15.219167 -33.462654 + 105700 1.3737039e-06 -2.6162001 -6.0174388 -15.195423 -33.460842 + 105800 1.3698978e-06 -2.5741793 -5.9823122 -15.171485 -33.458818 + 105900 1.3618787e-06 -2.5324483 -5.9474396 -15.147407 -33.456583 + 106000 1.3497592e-06 -2.5004006 -5.9127828 -15.123235 -33.454147 + 106100 1.3337214e-06 -2.5004006 -5.878304 -15.099001 -33.451519 + 106200 1.313999e-06 -2.5004006 -5.8439715 -15.074727 -33.448713 + 106300 1.2908518e-06 -2.5004006 -5.8097644 -15.050423 -33.445745 + 106400 1.2645404e-06 -2.5004006 -5.7756772 -15.026093 -33.442631 + 106500 1.2353047e-06 -2.5004006 -5.741723 -15.001731 -33.439384 + 106600 1.2033522e-06 -2.5004006 -5.7079352 -14.977327 -33.436015 + 106700 1.1688577e-06 -2.5004006 -5.6743684 -14.952873 -33.432531 + 106800 1.1319736e-06 -2.5004006 -5.641097 -14.928358 -33.428935 + 106900 1.0928486e-06 -2.5004006 -5.6082138 -14.903781 -33.425226 + 107000 1.0516487e-06 -2.5004006 -5.575826 -14.879144 -33.421403 + 107100 1.0085757e-06 -2.5004006 -5.5440517 -14.854461 -33.417461 + 107200 9.6387785e-07 -2.5004006 -5.5130153 -14.829755 -33.413396 + 107300 9.1785266e-07 -2.5004006 -5.4828423 -14.805063 -33.409206 + 107400 8.7083734e-07 -2.5004006 -5.4536544 -14.780429 -33.404894 + 107500 8.2319132e-07 -2.5004006 -5.425565 -14.755909 -33.400466 + 107600 7.7527543e-07 -2.5004006 -5.3986747 -14.731569 -33.395933 + 107700 7.2743149e-07 -2.5004006 -5.3730678 -14.707478 -33.391307 + 107800 6.7996411e-07 -2.5004006 -5.3488093 -14.683712 -33.386603 + 107900 6.3313122e-07 -2.5004006 -5.3259432 -14.660344 -33.381839 + 108000 5.8714599e-07 -2.5004006 -5.3044918 -14.63745 -33.37703 + 108100 5.4218381e-07 -2.5004006 -5.2844564 -14.615099 -33.37219 + 108200 4.9839695e-07 -2.5004006 -5.2658187 -14.593356 -33.367329 + 108300 4.559321e-07 -2.5004006 -5.2485432 -14.572279 -33.362457 + 108400 4.1494542e-07 -2.5004006 -5.2325809 -14.551913 -33.357576 + 108500 3.7561221e-07 -2.5004006 -5.2178732 -14.5323 -33.352686 + 108600 3.3812879e-07 -2.5004006 -5.2043565 -14.513465 -33.347788 + 108700 3.0270652e-07 -2.5004006 -5.1919663 -14.495427 -33.342881 + 108800 2.6955871e-07 -2.5004006 -5.180642 -14.478192 -33.337969 + 108900 2.3888338e-07 -2.5004006 -5.1703304 -14.461759 -33.333055 + 109000 2.1084663e-07 -2.5004006 -5.1609896 -14.446117 -33.328143 + 109100 1.8556964e-07 -2.5004006 -5.1525916 -14.431248 -33.323242 + 109200 1.6312328e-07 -2.5004006 -5.1451235 -14.417132 -33.318362 + 109300 1.4352767e-07 -2.5004006 -5.1385886 -14.403742 -33.313514 + 109400 1.267618e-07 -2.5004006 -5.133006 -14.391055 -33.308712 + 109500 1.1277879e-07 -2.5004006 -5.1284086 -14.379049 -33.303968 + 109600 1.0151997e-07 -2.5004006 -5.1248416 -14.367705 -33.299294 + 109700 9.2927869e-08 -2.5004006 -5.1223582 -14.357012 -33.294697 + 109800 8.6955973e-08 -2.5004006 -5.1210168 -14.346966 -33.290188 + 109900 8.3569025e-08 -2.5004006 -5.1208763 -14.337572 -33.285767 + 110000 8.2740535e-08 -2.5004006 -5.1219918 -14.328843 -33.281432 + 110100 8.4443675e-08 -2.5004006 -5.1244115 -14.3208 -33.277178 + 110200 8.8639129e-08 -2.5004006 -5.1281723 -14.313471 -33.273001 + 110300 9.5262978e-08 -2.5004006 -5.1332973 -14.306888 -33.268903 + 110400 1.0421742e-07 -2.5004006 -5.1397937 -14.301088 -33.264886 + 110500 1.1536651e-07 -2.5004006 -5.1476515 -14.296108 -33.26095 + 110600 1.2853831e-07 -2.5004006 -5.1568431 -14.291983 -33.257099 + 110700 1.4353299e-07 -2.5004006 -5.1673241 -14.288746 -33.25334 + 110800 1.6013514e-07 -2.5004006 -5.1790343 -14.286423 -33.249685 + 110900 1.7812767e-07 -2.5004006 -5.1919001 -14.285035 -33.246146 + 111000 1.9730403e-07 -2.5004006 -5.2058375 -14.284593 -33.242731 + 111100 2.174765e-07 -2.5004006 -5.220755 -14.285098 -33.239452 + 111200 2.3847885e-07 -2.5004006 -5.2365577 -14.286543 -33.236316 + 111300 2.6016334e-07 -2.5004006 -5.2531507 -14.288914 -33.233329 + 111400 2.8239334e-07 -2.5004006 -5.2704429 -14.292184 -33.230496 + 111500 3.0503358e-07 -2.5004006 -5.2883503 -14.29632 -33.227817 + 111600 3.2794066e-07 -2.5004006 -5.3067986 -14.301285 -33.225289 + 111700 3.5095611e-07 -2.5004006 -5.3257256 -14.307032 -33.222909 + 111800 3.739037e-07 -2.5004006 -5.3450822 -14.313516 -33.220672 + 111900 3.9659154e-07 -2.5004006 -5.3648334 -14.320686 -33.218575 + 112000 4.1881831e-07 -2.5004006 -5.3849576 -14.328494 -33.216615 + 112100 4.4038335e-07 -2.5004006 -5.4054455 -14.336891 -33.214793 + 112200 4.6109763e-07 -2.5004006 -5.4262983 -14.345835 -33.213108 + 112300 4.8079364e-07 -2.5004006 -5.4475254 -14.355286 -33.211565 + 112400 4.9933188e-07 -2.5004006 -5.4691411 -14.365213 -33.210167 + 112500 5.1660295e-07 -2.5004006 -5.4911615 -14.375588 -33.20892 + 112600 5.3252481e-07 -2.5004006 -5.5136014 -14.386395 -33.207829 + 112700 5.4703619e-07 -2.5004006 -5.5364708 -14.397621 -33.206897 + 112800 5.6008784e-07 -2.5004006 -5.5597723 -14.409261 -33.206127 + 112900 5.7163385e-07 -2.5004006 -5.5834987 -14.421317 -33.205518 + 113000 5.8162508e-07 -2.5004006 -5.6076309 -14.43379 -33.205068 + 113100 5.900063e-07 -2.5004006 -5.6321376 -14.446688 -33.204772 + 113200 5.9671738e-07 -2.5004006 -5.6569743 -14.460016 -33.204625 + 113300 6.0169828e-07 -2.5004006 -5.6820843 -14.473778 -33.20462 + 113400 6.0489638e-07 -2.5004006 -5.7073998 -14.487975 -33.204751 + 113500 6.0627464e-07 -2.5004006 -5.7328438 -14.502602 -33.205013 + 113600 6.0581853e-07 -2.5004006 -5.7583325 -14.517649 -33.205403 + 113700 6.0354028e-07 -2.5004006 -5.783778 -14.533099 -33.20592 + 113800 5.9947954e-07 -2.5004006 -5.8090911 -14.548926 -33.206563 + 113900 5.9370019e-07 -2.5159493 -5.8341843 -14.565099 -33.207336 + 114000 5.8628419e-07 -2.5451136 -5.8589749 -14.58158 -33.208242 + 114100 5.773237e-07 -2.5737325 -5.883387 -14.598325 -33.209282 + 114200 5.6691316e-07 -2.6017422 -5.9073542 -14.615285 -33.21046 + 114300 5.5514312e-07 -2.6290924 -5.9308213 -14.632407 -33.211774 + 114400 5.4209684e-07 -2.6557466 -5.9537448 -14.649638 -33.213224 + 114500 5.2785052e-07 -2.681682 -5.9760937 -14.666924 -33.214806 + 114600 5.1247659e-07 -2.7070928 -5.9978496 -14.684213 -33.216515 + 114700 4.9604917e-07 -2.7320576 -6.0190051 -14.701453 -33.218343 + 114800 4.7865027e-07 -2.7562929 -6.0395625 -14.718599 -33.220283 + 114900 4.6037496e-07 -2.779819 -6.0595325 -14.735608 -33.222326 + 115000 4.4133424e-07 -2.8026588 -6.0789306 -14.752447 -33.224465 + 115100 4.2165484e-07 -2.8248348 -6.0977757 -14.769085 -33.226692 + 115200 4.0147597e-07 -2.8463656 -6.1160866 -14.7855 -33.229 + 115300 3.8094391e-07 -2.8672636 -6.1338801 -14.801677 -33.231385 + 115400 3.6020502e-07 -2.8875328 -6.1511682 -14.817604 -33.233844 + 115500 3.3939905e-07 -2.907167 -6.1679562 -14.833278 -33.236375 + 115600 3.1865396e-07 -2.9261493 -6.1842417 -14.848697 -33.238976 + 115700 2.9808318e-07 -2.9444517 -6.2000131 -14.863865 -33.241648 + 115800 2.7778573e-07 -2.9620356 -6.2152499 -14.878786 -33.244388 + 115900 2.5784886e-07 -2.9788534 -6.2299226 -14.893463 -33.247194 + 116000 2.3835251e-07 -2.99485 -6.243994 -14.9079 -33.250064 + 116100 2.1937431e-07 -3.0099654 -6.2574204 -14.922095 -33.252991 + 116200 2.0099393e-07 -3.0241375 -6.2701532 -14.936046 -33.25597 + 116300 1.8329564e-07 -3.0373046 -6.282142 -14.949743 -33.258994 + 116400 1.6636846e-07 -3.0494088 -6.2933358 -14.963172 -33.262056 + 116500 1.5030376e-07 -3.0603984 -6.3036862 -14.976314 -33.265149 + 116600 1.3519094e-07 -3.0702307 -6.3131493 -14.989147 -33.268265 + 116700 1.2111212e-07 -3.0788734 -6.3216877 -15.00164 -33.271401 + 116800 1.081369e-07 -3.0863067 -6.3292725 -15.013763 -33.274552 + 116900 9.6318621e-08 -3.0925235 -6.3358839 -15.025483 -33.277715 + 117000 8.569256e-08 -3.09753 -6.3415129 -15.036763 -33.280888 + 117100 7.6276644e-08 -3.1013449 -6.3461607 -15.047569 -33.284068 + 117200 6.8074253e-08 -3.1039982 -6.3498389 -15.057869 -33.287255 + 117300 6.1078459e-08 -3.1055299 -6.3525685 -15.06763 -33.290446 + 117400 5.52767e-08 -3.1059874 -6.354379 -15.076828 -33.293639 + 117500 5.0654836e-08 -3.1054234 -6.3553064 -15.085439 -33.29683 + 117600 4.7199554e-08 -3.103893 -6.3553917 -15.093447 -33.300015 + 117700 4.4898805e-08 -3.1014511 -6.3546782 -15.100841 -33.303188 + 117800 4.3741386e-08 -3.0981502 -6.3532103 -15.107617 -33.306343 + 117900 4.3713431e-08 -3.0940376 -6.3510303 -15.113773 -33.309473 + 118000 4.4794472e-08 -3.0891536 -6.3481775 -15.119318 -33.312571 + 118100 4.6953371e-08 -3.0835306 -6.3446858 -15.124261 -33.315631 + 118200 5.0146958e-08 -3.0771919 -6.3405829 -15.128618 -33.318647 + 118300 5.4318621e-08 -3.0701515 -6.3358895 -15.132406 -33.321616 + 118400 5.9399274e-08 -3.0624147 -6.3306186 -15.135646 -33.324535 + 118500 6.5310223e-08 -3.0539791 -6.3247763 -15.138359 -33.327399 + 118600 7.1967003e-08 -3.044836 -6.3183621 -15.140565 -33.330205 + 118700 7.9283741e-08 -3.034972 -6.3113698 -15.142283 -33.332952 + 118800 8.7177355e-08 -3.0243717 -6.3037897 -15.143528 -33.335639 + 118900 9.5569716e-08 -3.0130198 -6.2956093 -15.144315 -33.338267 + 119000 1.0438844e-07 -3.000903 -6.2868162 -15.144651 -33.340835 + 119100 1.1356601e-07 -2.9880129 -6.277399 -15.144543 -33.34334 + 119200 1.2303752e-07 -2.9743476 -6.26735 -15.14399 -33.345778 + 119300 1.327378e-07 -2.9599133 -6.2566661 -15.142991 -33.348147 + 119400 1.4259868e-07 -2.9447257 -6.2453506 -15.141538 -33.350444 + 119500 1.5254705e-07 -2.9288102 -6.2334142 -15.139624 -33.352663 + 119600 1.6250448e-07 -2.9122025 -6.2208752 -15.137238 -33.3548 + 119700 1.7238807e-07 -2.8949478 -6.2077601 -15.13437 -33.356853 + 119800 1.8211288e-07 -2.8771001 -6.1941032 -15.131011 -33.35882 + 119900 1.9159516e-07 -2.8587205 -6.1799459 -15.127154 -33.3607 + 120000 2.0075596e-07 -2.8398756 -6.1653354 -15.122793 -33.362492 +Loop time of 66.7366 on 1 procs for 50000 steps with 9576 atoms + +Performance: 6473210.162 tau/day, 749.214 timesteps/s, 7.174 Matom-step/s +99.7% CPU use with 1 MPI tasks x no OpenMP threads + +MPI task timing breakdown: +Section | min time | avg time | max time |%varavg| %total +--------------------------------------------------------------- +Pair | 6.3674 | 6.3674 | 6.3674 | 0.0 | 9.54 +Bond | 50.92 | 50.92 | 50.92 | 0.0 | 76.30 +Neigh | 1.1547 | 1.1547 | 1.1547 | 0.0 | 1.73 +Comm | 0.11662 | 0.11662 | 0.11662 | 0.0 | 0.17 +Output | 0.041399 | 0.041399 | 0.041399 | 0.0 | 0.06 +Modify | 6.3682 | 6.3682 | 6.3682 | 0.0 | 9.54 +Other | | 1.768 | | | 2.65 + +Nlocal: 9576 ave 9576 max 9576 min +Histogram: 1 0 0 0 0 0 0 0 0 0 +Nghost: 0 ave 0 max 0 min +Histogram: 1 0 0 0 0 0 0 0 0 0 +Neighs: 44738 ave 44738 max 44738 min +Histogram: 1 0 0 0 0 0 0 0 0 0 + +Total # of neighbors = 44738 +Ave neighs/atom = 4.6718881 +Ave special neighs/atom = 11.495405 +Neighbor list builds = 219 +Dangerous builds = 0 +Total wall time: 0:02:17 diff --git a/examples/bpm/poissons_ratio/in.bpm.poissons_ratio b/examples/bpm/poissons_ratio/in.bpm.poissons_ratio index 4f31291a6a..71052ad417 100644 --- a/examples/bpm/poissons_ratio/in.bpm.poissons_ratio +++ b/examples/bpm/poissons_ratio/in.bpm.poissons_ratio @@ -1,4 +1,4 @@ -# Squish three pillars with different Poisson's ratios +# Squish four pillars with different Poisson's ratios # going from nearly incompressible to auxetic # from left to right, nu ~ 0.5, 0.25, 0, <0 # calibration from DOI: 10.1039/d3sm01373a