Added the 2D crack propagation program. Works on 1 partition, to be tested with
authorMilind Bhandarkar <milind@cs.uiuc.edu>
Thu, 24 Aug 2000 20:05:33 +0000 (20:05 +0000)
committerMilind Bhandarkar <milind@cs.uiuc.edu>
Thu, 24 Aug 2000 20:05:33 +0000 (20:05 +0000)
multiple partitions.

examples/fem/crack2D/Makefile [new file with mode: 0644]
examples/fem/crack2D/cohesive.inp [new file with mode: 0644]
examples/fem/crack2D/crack.h [new file with mode: 0644]
examples/fem/crack2D/crck_bar.inp [new file with mode: 0644]
examples/fem/crack2D/driver.C [new file with mode: 0644]
examples/fem/crack2D/getmesh.C [new file with mode: 0644]
examples/fem/crack2D/lst_NL.C [new file with mode: 0644]
examples/fem/crack2D/lst_coh2.C [new file with mode: 0644]
examples/fem/crack2D/readfile.C [new file with mode: 0644]
examples/fem/crack2D/volume.C [new file with mode: 0644]

diff --git a/examples/fem/crack2D/Makefile b/examples/fem/crack2D/Makefile
new file mode 100644 (file)
index 0000000..d348871
--- /dev/null
@@ -0,0 +1,28 @@
+CHARMC=../../../../bin/charmc $(OPTS)
+all: getmesh pgm
+
+getmesh: getmesh.C
+       $(CHARMC) -seq -o getmesh getmesh.C -language c++
+
+pgm: driver.o readfile.o volume.o lst_NL.o lst_coh2.o
+       $(CHARMC) -o pgm readfile.o volume.o lst_NL.o lst_coh2.o driver.o \
+         -language fem 
+
+readfile.o: readfile.C crack.h
+       $(CHARMC) -c readfile.C
+
+driver.o: driver.C crack.h
+       $(CHARMC) -c driver.C
+
+volume.o: volume.C crack.h
+       $(CHARMC) -c volume.C
+
+lst_NL.o: lst_NL.C crack.h
+       $(CHARMC) -c lst_NL.C
+
+lst_coh2.o: lst_coh2.C crack.h
+       $(CHARMC) -c lst_coh2.C
+
+clean:
+       rm -rf *.o pgm *.sts *.log *.bak *~ SunWS_cache ii_files ti_files
+       rm -f conv-host getmesh *.mesh meshdata.*
diff --git a/examples/fem/crack2D/cohesive.inp b/examples/fem/crack2D/cohesive.inp
new file mode 100644 (file)
index 0000000..5bc76b1
--- /dev/null
@@ -0,0 +1,43 @@
+prefix of output files
+crck_bar
+Execution: ntime,nsteps,ntsint_full,ntsint_energy,restart(0=no,1=yes)
+        2000 20 500 15003 0
+nplane(0=stress,1=strain,2=axi) ncoh(0=exp,1=lin) lin(0=nonlin,1=linear)
+        0 1 0
+number of volumetric material sets : numat_vol
+        1
+each material set : e1,e2,g12,g23,nu12,nu23,rho,alpha1,alpha2,theta(deg)
+       3.24e9 3.24e9 1.20e9 1.20e9 0.35 0.35 1190 0.0 0.0 0.0
+number of cohesive material sets : numat_coh
+        1
+for each material set : deltan, deltat, sigmax, taumax, mu, sinit
+       2.416e-6 2.416e-6 3.24e8 3.24e8 0.0 0.9
+impact data : 1 or 0, vo_imp (m/s), mass_imp (kg/m),  radius_imp (m)
+        0  0.0  1.0  0.0
+Young's modulus and Poisson's ratio of impactor: E_imp, xnu_imp
+        1.0 0.3
+thermal load
+        0.0
+boundary conditions (0 for imp. velo, 1 for imp. force)
+1  0 0 0.0 0.0
+2  0 0 0.0 0.0
+3  0 0 1.0 1.0
+4  0 0 1.0 1.0
+5  0 0 1.0 1.0
+6  1 0 0.0 0.0
+7  1 0 0.0 20.0
+8  1 0 0.0 0.0
+9  1 0 0.0 0.0
+10 1 0 0.0 0.0
+11 0 1 0.0 0.0
+12 0 1 0.0 0.0
+13 0 1 0.0 0.0
+14 0 1 0.0 0.0
+15 0 1 0.0 0.0
+16 1 1 0.0 0.0
+17 1 1 0.0 0.0
+18 1 1 0.0 0.0
+19 1 1 0.0 0.0
+20 1 1 0.0 0.0
+
+
diff --git a/examples/fem/crack2D/crack.h b/examples/fem/crack2D/crack.h
new file mode 100644 (file)
index 0000000..711e397
--- /dev/null
@@ -0,0 +1,135 @@
+#include <string>
+#include <stdlib.h>
+#include "fem.h"
+
+//Coord describes a property with X and Y components
+struct Coord
+{
+  double x;
+  double y;
+};
+
+// Declaration for class VolElement (Volumetric Elements) 
+// These are 6-vertex triangles (each side has an extra vertex)
+struct Vol
+{
+  double s11l[3], s22l[3], s12l[3];//Stress coefficients
+};
+
+
+// Declaration for class CohElement (Cohesive Elements) 
+//  These are 6-vertex rectangles (long sides have an extra vertex)
+struct Coh
+{
+  Coord T[3];//Tractions at each sample point
+  double sidel[3];//[0]->length of element; 
+                //[1],[2] give cosine and sine of orientation
+  double Sthresh[3];//The threshold, and damage for this edge
+};
+
+// Declaration for class Element (generic elements)
+struct Element
+{
+  int material;          // matlst, matclst, matc
+  union {
+    Vol v;
+    Coh c;
+  };
+};
+
+  
+// Declaration for class Material
+struct Material
+{
+  double c[4], ci[4];          // Elastic stiffness constants & inverses
+};
+
+struct VolMaterial : public Material
+{
+  double e1, e2;               // Young's moduli 
+  double g12, g23;             // Shear moduli 
+  double xnu12, xnu23;         // Poisson's ratios 
+  double rho;                  // density 
+  double alpha1, alpha2;       // thermal expansion coefficients 
+  double theta;                // principal material direction 
+};
+
+struct CohMaterial : public Material
+{
+  double deltan;           // normal characteristic length
+  double deltat;           // tangent characteristic length
+  double sigmax;           // max. normal stress
+  double taumax;           // max. shearing stress
+  double mu;               // coefficient of friction
+  double Sinit;            // initial Sthreshold
+};
+  
+struct Node
+{
+  Coord Rin;             //Internal force
+  Coord Rco;             //Cohesive traction load
+  Coord xM;              //Mass at this node (xM.x==xM.y, too)
+  Coord vel;
+  Coord accel;
+  Coord disp;
+  Coord pos;
+  Coord r;
+  unsigned char id1, id2;
+  unsigned char isbnd;
+};
+
+//Global constants
+const int    numBoundMax = 1123;   // Maximum number of boundary nodes
+const double   g1          = -0.774596669241483;
+const double   g3          = 0.774596669241483;
+const double   w1          = 0.555555555555555;
+const double   w2          = 0.888888888888888;
+const double   w3          = 0.555555555555555;
+const double   pi          = 3.14159265358979;
+
+struct GlobalData {
+  int numNP;               //number of nodal points (numnp)
+  int numLST;              //number of LST elements (numlst)
+  int numCLST;             //number of LST cohesive elements (numclst)
+  int numBound;            //number of boundary nodes w/loads
+  int nTime;               //total number of time steps
+  int steps;               //ratio of delta and Courant cond
+  double delta;            //timestep
+  double delta2;           //delta squared times 2????
+  int *ts_proportion;      //time step for given constant
+  double *proportion;      //load proportionality constant
+
+  int lin;                 //elastic formulation, 0=nl, 1=linear
+  int ncoh;                //type of cohesive law, 0=exp., 1=linear
+  int nplane;              //type of plane analysis, 0=stress, 1=strain
+  int tsintFull;          //output interval in time steps
+  int tsintEnergy;        //partial output interval
+  int restart;
+  int imp;                //Is there impact?
+  double voImp;         //Velocity of impactor
+  double dImp;             //Displacement
+  double aImp;             //Acceleration
+  double vImp;             //Velocity
+  double xnuImp,eImp,eTop,radiusImp;   //Impact parameters
+  double indent;           //indentation of impactor
+  int nImp;            //node hit
+  double fImp;             //contact force
+  double massImp;          //mass
+
+  int numMatVol;           //number of volumetric materials (numat_vol)
+  int numMatCoh;           //number of cohesive materials (numat_coh)
+  int numProp;             //number of proportionality constants
+  VolMaterial *volm;
+  CohMaterial *cohm;
+
+  int nn, ne, npere;
+  int *nnums, *enums, *conn;
+  Node *nodes;
+  Element *elements;
+  int scoh, ecoh, svol, evol;
+};
+
+extern void readFile(GlobalData *gd);
+extern void vol_elem(GlobalData *gd);
+extern void lst_NL(GlobalData *gd);
+extern void lst_coh2(GlobalData *gd);
diff --git a/examples/fem/crack2D/crck_bar.inp b/examples/fem/crack2D/crck_bar.inp
new file mode 100644 (file)
index 0000000..30c1d3b
--- /dev/null
@@ -0,0 +1,3347 @@
+ 2196     366
+    2.91979169935698332E-8
+    3
+    1     0
+    20     1
+    2000     1
+    2196
+    1     2.70476109685850318E-4     1.47470775643085319E-4
+    2     2.62862460053854216E-4     2.58243951636348947E-4
+    3     1.80546315170733955E-4     1.83727399972074352E-4
+    4     2.66669284869852294E-4     2.02857363639717133E-4
+    5     2.21704387612294072E-4     2.2098567580421165E-4
+    6     2.2551121242829215E-4     1.65599087807579836E-4
+    7     9.21159286470410715E-4     6.74999999999999927E-4
+    8     8.35489473019320351E-4     7.27605773008054794E-4
+    9     8.14703805291646039E-4     6.47042209754021894E-4
+    10     8.78324379744865533E-4     7.01302886504027306E-4
+    11     8.25096639155483195E-4     6.87323991381038344E-4
+    12     8.67931545881028377E-4     6.61021104877010856E-4
+    13     2.70476109685850318E-4     1.47470775643085319E-4
+    14     3.64813480721653041E-4     1.25057137488653904E-4
+    15     3.3457565114347566E-4     2.07524696466039707E-4
+    16     3.17644795203751652E-4     1.36263956565869611E-4
+    17     3.49694565932564378E-4     1.66290916977346792E-4
+    18     3.02525880414662989E-4     1.774977360545625E-4
+    19     1.1927993626285668E-4     9.2101752411100806E-4
+    20     1.80463733243880608E-4     1.01635139957644981E-3
+    21     8.06327547355035015E-5     1.013121025387205E-3
+    22     1.4987183475336863E-4     9.6868446184372888E-4
+    23     1.30548243989692055E-4     1.0147362124818273E-3
+    24     9.99563454991800906E-5     9.67069274749106584E-4
+    25     4.25259790828855515E-4     1.93997913774442133E-4
+    26     3.662109375E-4     2.98836517333984323E-4
+    27     3.3457565114347566E-4     2.07524696466039707E-4
+    28     3.95735364164427757E-4     2.46417215554213228E-4
+    29     3.5039329432173783E-4     2.53180606900012001E-4
+    30     3.79917720986165588E-4     2.00761305120240906E-4
+    31     8.70539573309993481E-4     3.74999999999999954E-4
+    32     8.08703962690366035E-4     4.17299929575055159E-4
+    33     7.77749415513451876E-4     3.13113370342386432E-4
+    34     8.39621768000179758E-4     3.96149964787527583E-4
+    35     7.93226689101908956E-4     3.65206649958720823E-4
+    36     8.24144494411722679E-4     3.44056685171193193E-4
+    37     1.74429316483144684E-4     8.22070312499999993E-4
+    38     8.72146582415723147E-5     8.12372124900119985E-4
+    39     1.36565743695333484E-4     7.39812343895404297E-4
+    40     1.30821987362358499E-4     8.17221218700059935E-4
+    41     1.11890200968452899E-4     7.76092234397762195E-4
+    42     1.55497530089239097E-4     7.80941328197702145E-4
+    43     2.50000000000000005E-4     5.99999999999999947E-4
+    44     3.79565818807853228E-4     5.99785262317693206E-4
+    45     3.14898739406715305E-4     6.69780733668936504E-4
+    46     3.14782909403926617E-4     5.99892631158846631E-4
+    47     3.47232279107284267E-4     6.34782997993314909E-4
+    48     2.82449369703357655E-4     6.34890366834468226E-4
+    49     1.220703125E-4     1.19999999999999989E-3
+    50     1.27766480316576791E-4     1.1007034017909631E-3
+    51     1.8310546875E-4     1.19999999999999989E-3
+    52     1.24918396408288396E-4     1.15035170089548139E-3
+    53     1.55435974533288396E-4     1.15035170089548139E-3
+    54     1.52587890625E-4     1.19999999999999989E-3
+    55     7.24825013803433236E-4     8.92866282184746245E-4
+    56     6.94054151555351649E-4     7.95336883715192389E-4
+    57     7.56891934562431879E-4     7.83576045089881936E-4
+    58     7.09439582679392442E-4     8.44101582949969317E-4
+    59     7.25473043058891764E-4     7.89456464402537163E-4
+    60     7.40858474182932557E-4     8.3822116363731409E-4
+    61     8.92816374246439338E-4     5.66266606957709193E-4
+    62     7.90472897978866684E-4     5.45453687327126634E-4
+    63     8.58425917146639784E-4     4.73341467989030294E-4
+    64     8.41644636112653011E-4     5.55860147142417859E-4
+    65     8.2444940756275318E-4     5.09397577658078464E-4
+    66     8.75621145696539507E-4     5.19804037473369689E-4
+    67     0     8.22070312499999993E-4
+    68     6.75304712410491789E-5     8.77236181658321186E-4
+    69     0     8.83105468749999993E-4
+    70     3.37652356205245894E-5     8.49653247079160644E-4
+    71     3.37652356205245894E-5     8.80170825204160644E-4
+    72     0     8.52587890624999993E-4
+    73     5.19966569338154057E-4     8.58840618576917807E-4
+    74     4.51818848350113347E-4     9.0027268151548134E-4
+    75     4.10887142695853748E-4     8.117939619078767E-4
+    76     4.85892708844133702E-4     8.79556650046199573E-4
+    77     4.31352995522983574E-4     8.56033321711678966E-4
+    78     4.65426856017003929E-4     8.35317290242397308E-4
+    79     6.27929687500000016E-4     9.68445999237773649E-4
+    80     5.67809578584186053E-4     1.01842229960029352E-3
+    81     5.42759630809160629E-4     9.44364496429779174E-4
+    82     5.97869633042093034E-4     9.9343414941903364E-4
+    83     5.55284604696673341E-4     9.81393398015036403E-4
+    84     5.85344659154580322E-4     9.56405247833776412E-4
+    85     4.8087086624308442E-4     1.00763880523435656E-3
+    86     5.63964843750000013E-4     1.0965073564807807E-3
+    87     4.64295327735872877E-4     1.106418874786976E-3
+    88     5.22417854996542162E-4     1.05207308085756863E-3
+    89     5.14130085742936418E-4     1.10146311563387835E-3
+    90     4.72583096989478675E-4     1.05702884001066628E-3
+    91     7.24825013803433236E-4     8.92866282184746245E-4
+    92     8.37129054430499212E-4     8.81217371382065837E-4
+    93     8.1517205990673709E-4     9.40785388914013596E-4
+    94     7.8097703411696617E-4     8.87041826783406041E-4
+    95     8.26150557168618205E-4     9.11001380148039716E-4
+    96     7.69998536855085163E-4     9.1682583554937992E-4
+    97     5.59110560621724098E-4     5.93289128784703487E-4
+    98     4.45119450916024088E-4     5.73219517420557926E-4
+    99     5.77009580214111854E-4     5.23102021510448161E-4
+    100     5.0211500576887412E-4     5.83254323102630652E-4
+    101     5.11064515565067944E-4     5.48160769465503043E-4
+    102     5.6806007041791803E-4     5.58195575147575878E-4
+    103     1.80546315170733955E-4     1.83727399972074352E-4
+    104     9.02731575853669911E-5     2.19793387486037181E-4
+    105     1.01956056377857214E-4     1.1802584118703223E-4
+    106     1.3540973637805046E-4     2.01760393729055753E-4
+    107     9.61146069816121027E-5     1.68909614336534706E-4
+    108     1.41251185774295585E-4     1.50876620579553278E-4
+    109     2.39680731901657833E-4     5.14383951086584473E-4
+    110     1.87500000000000004E-4     5.7499999999999999E-4
+    111     1.72756436172360502E-4     4.70585859703467847E-4
+    112     2.13590365950828932E-4     5.44691975543292232E-4
+    113     1.80128218086180239E-4     5.22792929851733919E-4
+    114     2.06218584037009167E-4     4.9248490539502616E-4
+    115     7.42631330235908388E-4     4.52617545905931022E-4
+    116     6.86868266859239985E-4     5.40099367412093744E-4
+    117     6.55870398270680236E-4     4.58827279636259491E-4
+    118     7.14749798547574186E-4     4.96358456659012383E-4
+    119     6.71369332564960165E-4     4.99463323524176645E-4
+    120     6.99250864253294366E-4     4.55722412771095283E-4
+    121     7.66063425991097419E-4     2.25618978565681171E-4
+    122     6.73095146243751124E-4     1.78431349377634624E-4
+    123     7.40024751298205222E-4     1.41990090695903195E-4
+    124     7.19579286117424271E-4     2.02025163971657911E-4
+    125     7.06559948770978173E-4     1.6021072003676891E-4
+    126     7.53044088644651321E-4     1.83804534630792197E-4
+    127     9.37074205321842733E-4     6.1279086004113368E-4
+    128     1.00000000000000002E-3     5.99999999999999947E-4
+    129     1.00000000000000002E-3     6.75000000000000036E-4
+    130     9.68537102660921431E-4     6.06395430020566814E-4
+    131     1.00000000000000002E-3     6.37500000000000046E-4
+    132     9.68537102660921431E-4     6.43895430020566912E-4
+    133     6.91894531250000018E-4     1.19999999999999989E-3
+    134     7.63546512545414027E-4     1.12019473601693395E-3
+    135     7.55859375000000021E-4     1.19999999999999989E-3
+    136     7.27720521897706968E-4     1.16009736800846692E-3
+    137     7.59702943772707024E-4     1.16009736800846692E-3
+    138     7.23876953124999965E-4     1.19999999999999989E-3
+    139     5.63964843750000013E-4     1.0965073564807807E-3
+    140     5.67809578584186053E-4     1.01842229960029352E-3
+    141     6.3613763360561681E-4     1.03946949707743879E-3
+    142     5.65887211167093087E-4     1.05746482804053711E-3
+    143     6.01973606094901486E-4     1.02894589833886626E-3
+    144     6.00051238677808412E-4     1.06798842677910963E-3
+    145     6.28726770153248599E-4     8.69566754220071274E-4
+    146     6.94054151555351649E-4     7.95336883715192389E-4
+    147     7.24825013803433236E-4     8.92866282184746245E-4
+    148     6.61390460854300124E-4     8.32451818967631886E-4
+    149     7.09439582679392442E-4     8.44101582949969317E-4
+    150     6.76775891978340918E-4     8.81216518202408705E-4
+    151     2.16128514871350092E-4     6.96571287178375043E-4
+    152     1.25000000000000003E-4     6.4999999999999997E-4
+    153     1.87500000000000004E-4     6.24999999999999905E-4
+    154     1.70564257435675047E-4     6.73285643589187507E-4
+    155     1.56250000000000003E-4     6.37499999999999937E-4
+    156     2.01814257435675048E-4     6.60785643589187474E-4
+    157     2.97380987249783093E-4     1.12592184443437615E-3
+    158     3.88700445385212667E-4     1.12415429991009354E-3
+    159     3.72070312500000005E-4     1.19999999999999989E-3
+    160     3.43040716317497907E-4     1.12503807217223485E-3
+    161     3.80385378942606363E-4     1.16207714995504672E-3
+    162     3.34725649874891549E-4     1.16296092221718802E-3
+    163     1.74429316483144684E-4     8.22070312499999993E-4
+    164     2.29766755665683817E-4     7.70779974744912363E-4
+    165     2.56894004484586392E-4     9.02486475935557219E-4
+    166     2.02098036074414237E-4     7.96425143622456178E-4
+    167     2.43330380075135091E-4     8.36633225340234791E-4
+    168     2.15661660483865538E-4     8.62278394217778606E-4
+    169     1.27766480316576791E-4     1.1007034017909631E-3
+    170     8.06327547355035015E-5     1.013121025387205E-3
+    171     1.80463733243880608E-4     1.01635139957644981E-3
+    172     1.04199617526040146E-4     1.05691221358908405E-3
+    173     1.30548243989692055E-4     1.0147362124818273E-3
+    174     1.541151067802287E-4     1.05852740068370635E-3
+    175     1.73262011544884361E-4     2.84564040120359029E-4
+    176     2.62862460053854216E-4     2.58243951636348947E-4
+    177     2.31915182548602798E-4     3.185630410561948E-4
+    178     2.18062235799369275E-4     2.71403995878353988E-4
+    179     2.47388821301228494E-4     2.88403496346271873E-4
+    180     2.0258859704674358E-4     3.01563540588276914E-4
+    181     3.78329331029246172E-4     4.03040117695809171E-4
+    182     3.15575840731535889E-4     4.21078468404699436E-4
+    183     2.80075730378807853E-4     3.66277694634743636E-4
+    184     3.46952585880391058E-4     4.12059293050254303E-4
+    185     2.97825785555171844E-4     3.93678081519721536E-4
+    186     3.29202530704027013E-4     3.84658906165276403E-4
+    187     1.35354037059463752E-4     6.29907264959934823E-5
+    188     0     6.39648437500000026E-5
+    189     1.220703125E-4     0
+    190     6.76770185297318762E-5     6.34777851229967425E-5
+    191     6.103515625E-5     3.19824218750000013E-5
+    192     1.28712174779731876E-4     3.14953632479967412E-5
+    193     1.95469003684387412E-4     8.60228725173496005E-5
+    194     2.70476109685850318E-4     1.47470775643085319E-4
+    195     1.80546315170733955E-4     1.83727399972074352E-4
+    196     2.32972556685118865E-4     1.1674682408021746E-4
+    197     2.2551121242829215E-4     1.65599087807579836E-4
+    198     1.88007659427560697E-4     1.34875136244711976E-4
+    199     4.11905532296878403E-4     9.69320888710480284E-4
+    200     4.14780023525604651E-4     1.04732934218306818E-3
+    201     3.33510000366525681E-4     9.78187713756657664E-4
+    202     4.13342777911241527E-4     1.00832511544677412E-3
+    203     3.74145011946065166E-4     1.01275852796986303E-3
+    204     3.72707766331702042E-4     9.73754301233568974E-4
+    205     5.0000000000000001E-4     1.19999999999999989E-3
+    206     3.88700445385212667E-4     1.12415429991009354E-3
+    207     4.64295327735872877E-4     1.106418874786976E-3
+    208     4.44350222692606366E-4     1.16207714995504672E-3
+    209     4.26497886560542772E-4     1.11528658734853488E-3
+    210     4.82147663867936471E-4     1.15320943739348784E-3
+    211     4.11905532296878403E-4     9.69320888710480284E-4
+    212     4.51818848350113347E-4     9.0027268151548134E-4
+    213     4.8087086624308442E-4     1.00763880523435656E-3
+    214     4.31862190323495875E-4     9.34796785112980812E-4
+    215     4.66344857296598883E-4     9.5395574337491895E-4
+    216     4.46388199269981412E-4     9.88479846972418422E-4
+    217     5.63964843750000013E-4     1.0965073564807807E-3
+    218     5.63964843750000013E-4     1.19999999999999989E-3
+    219     5.0000000000000001E-4     1.19999999999999989E-3
+    220     5.63964843750000013E-4     1.14825367824039019E-3
+    221     5.31982421875000066E-4     1.19999999999999989E-3
+    222     5.31982421875000066E-4     1.14825367824039019E-3
+    223     3.14898739406715305E-4     6.69780733668936504E-4
+    224     3.98747049720822426E-4     7.15064925451776041E-4
+    225     3.3342684329849511E-4     7.35743010166048958E-4
+    226     3.56822894563768893E-4     6.92422829560356218E-4
+    227     3.66086946509658795E-4     7.25403967808912553E-4
+    228     3.24162791352605208E-4     7.02761871917492731E-4
+    229     6.28726770153248599E-4     8.69566754220071274E-4
+    230     5.79668274240937665E-4     8.10244013239797283E-4
+    231     6.3251725796544824E-4     7.69458270396488286E-4
+    232     6.04197522197093186E-4     8.39905383729934278E-4
+    233     6.06092766103192952E-4     7.89851141818142784E-4
+    234     6.30622014059348474E-4     8.1951251230827978E-4
+    235     3.79565818807853228E-4     5.99785262317693206E-4
+    236     3.75421084408168141E-4     5.26840780869595538E-4
+    237     4.45119450916024088E-4     5.73219517420557926E-4
+    238     3.77493451608010657E-4     5.63313021593644372E-4
+    239     4.10270267662096115E-4     5.50030149145076678E-4
+    240     4.12342634861938631E-4     5.86502389869125512E-4
+    241     5.60394849123839033E-4     7.35715772264671007E-4
+    242     4.75850298603600714E-4     7.54514336975707829E-4
+    243     4.62057889734597506E-4     6.86398137230062556E-4
+    244     5.18122573863719874E-4     7.45115054620189418E-4
+    245     4.6895409416909911E-4     7.20456237102885247E-4
+    246     5.11226369429218269E-4     7.11056954747366727E-4
+    247     7.20535236875600961E-4     6.29731480771863312E-4
+    248     7.52655853090133443E-4     7.19787645578877193E-4
+    249     6.60634466318554738E-4     7.08911474490989688E-4
+    250     7.36595544982867256E-4     6.74759563175370307E-4
+    251     7.06645159704344036E-4     7.1434956003493344E-4
+    252     6.90584851597077849E-4     6.69321477631426446E-4
+    253     9.00310858258578334E-4     8.8742365740778872E-4
+    254     8.1517205990673709E-4     9.40785388914013596E-4
+    255     8.37129054430499212E-4     8.81217371382065837E-4
+    256     8.57741459082657712E-4     9.14104523160901104E-4
+    257     8.26150557168618205E-4     9.11001380148039716E-4
+    258     8.68719956344538718E-4     8.84320514394927332E-4
+    259     4.45119450916024088E-4     5.73219517420557926E-4
+    260     3.75421084408168141E-4     5.26840780869595538E-4
+    261     4.44306653512529842E-4     5.02492175855262643E-4
+    262     4.10270267662096115E-4     5.50030149145076678E-4
+    263     4.09863868960348991E-4     5.1466647836242909E-4
+    264     4.44713052214276965E-4     5.3785584663791023E-4
+    265     5.469449348011982E-4     3.56411724829354306E-4
+    266     5.18154572622838414E-4     4.64528827497301226E-4
+    267     4.57065334871855982E-4     4.13720111675068352E-4
+    268     5.32549753712018362E-4     4.10470276163327793E-4
+    269     4.87609953747347225E-4     4.39124469586184762E-4
+    270     5.02005134836527064E-4     3.85065918252211329E-4
+    271     7.90472897978866684E-4     5.45453687327126634E-4
+    272     6.86868266859239985E-4     5.40099367412093744E-4
+    273     7.42631330235908388E-4     4.52617545905931022E-4
+    274     7.38670582419053388E-4     5.42776527369610243E-4
+    275     7.14749798547574186E-4     4.96358456659012383E-4
+    276     7.6655211410738759E-4     4.99035616616528774E-4
+    277     8.37129054430499212E-4     8.81217371382065837E-4
+    278     8.9886962494534965E-4     7.81508352249731973E-4
+    279     9.00310858258578334E-4     8.8742365740778872E-4
+    280     8.67999339687924431E-4     8.31362861815898905E-4
+    281     8.99590241601963938E-4     8.34466004828760292E-4
+    282     8.68719956344538718E-4     8.84320514394927332E-4
+    283     6.27455298609753205E-4     5.60579639099527925E-4
+    284     6.86868266859239985E-4     5.40099367412093744E-4
+    285     6.25000000000000013E-4     6.23375449734813098E-4
+    286     6.57161782734496649E-4     5.50339503255810889E-4
+    287     6.55934133429619945E-4     5.81737408573453421E-4
+    288     6.26227649304876609E-4     5.91977544417170458E-4
+    289     8.70539573309993481E-4     3.74999999999999954E-4
+    290     9.33976498687064506E-4     3.35267605171179269E-4
+    291     1.00000000000000002E-3     4.49999999999999988E-4
+    292     9.02258035998528994E-4     3.55133802585589611E-4
+    293     9.66988249343532264E-4     3.92633802585589655E-4
+    294     9.35269786654996751E-4     4.12499999999999998E-4
+    295     9.22298701458379552E-4     1.09448242187499991E-3
+    296     9.12446330413049594E-4     9.87167427587815364E-4
+    297     1.00000000000000002E-3     1.04999999999999994E-3
+    298     9.17372515935714573E-4     1.04082492473140764E-3
+    299     9.56223165206524861E-4     1.01858371379390765E-3
+    300     9.61149350729189732E-4     1.07224121093749993E-3
+    301     7.63546512545414027E-4     1.12019473601693395E-3
+    302     6.91894531250000018E-4     1.19999999999999989E-3
+    303     6.91894531250000018E-4     1.08422299961888672E-3
+    304     7.27720521897706968E-4     1.16009736800846692E-3
+    305     6.91894531250000018E-4     1.14211149980944331E-3
+    306     7.27720521897706968E-4     1.10220886781791034E-3
+    307     7.29985750852659397E-4     9.92266244781275222E-4
+    308     6.8729190829219008E-4     9.44648911063777994E-4
+    309     7.24825013803433236E-4     8.92866282184746245E-4
+    310     7.08638829572424793E-4     9.68457577922526608E-4
+    311     7.06058461047811658E-4     9.1875759662426212E-4
+    312     7.27405382328046371E-4     9.42566263483010734E-4
+    313     5.0000000000000001E-4     1.19999999999999989E-3
+    314     4.64295327735872877E-4     1.106418874786976E-3
+    315     5.63964843750000013E-4     1.0965073564807807E-3
+    316     4.82147663867936471E-4     1.15320943739348784E-3
+    317     5.14130085742936418E-4     1.10146311563387835E-3
+    318     5.31982421875000066E-4     1.14825367824039019E-3
+    319     7.70557602221913452E-4     1.04998903577999537E-3
+    320     8.33723092676693115E-4     1.08142335682909598E-3
+    321     7.63546512545414027E-4     1.12019473601693395E-3
+    322     8.02140347449303284E-4     1.06570619630454557E-3
+    323     7.98634802611053571E-4     1.10080904642301507E-3
+    324     7.67052057383663739E-4     1.08509188589846466E-3
+    325     8.42471384914615607E-4     1.00872691912019851E-3
+    326     9.22298701458379552E-4     1.09448242187499991E-3
+    327     8.33723092676693115E-4     1.08142335682909598E-3
+    328     8.82385043186497634E-4     1.0516046704975991E-3
+    329     8.78010897067536388E-4     1.08795288935204784E-3
+    330     8.38097238795654361E-4     1.04507513797464724E-3
+    331     8.33723092676693115E-4     1.08142335682909598E-3
+    332     8.16894531250000021E-4     1.19999999999999989E-3
+    333     7.63546512545414027E-4     1.12019473601693395E-3
+    334     8.25308811963346568E-4     1.14071167841454804E-3
+    335     7.90220521897707024E-4     1.16009736800846692E-3
+    336     7.98634802611053571E-4     1.10080904642301507E-3
+    337     4.8087086624308442E-4     1.00763880523435656E-3
+    338     5.42759630809160629E-4     9.44364496429779174E-4
+    339     5.67809578584186053E-4     1.01842229960029352E-3
+    340     5.11815248526122524E-4     9.76001650832067922E-4
+    341     5.55284604696673341E-4     9.81393398015036403E-4
+    342     5.24340222413635236E-4     1.01303055241732504E-3
+    343     7.70557602221913452E-4     1.04998903577999537E-3
+    344     6.91894531250000018E-4     1.08422299961888672E-3
+    345     7.29985750852659397E-4     9.92266244781275222E-4
+    346     7.31226066735956789E-4     1.06710601769944105E-3
+    347     7.10940141051329653E-4     1.03824462220008097E-3
+    348     7.50271676537286425E-4     1.0211276402806353E-3
+    349     7.55859375000000021E-4     1.19999999999999989E-3
+    350     7.63546512545414027E-4     1.12019473601693395E-3
+    351     8.16894531250000021E-4     1.19999999999999989E-3
+    352     7.59702943772707024E-4     1.16009736800846692E-3
+    353     7.90220521897707024E-4     1.16009736800846692E-3
+    354     7.86376953125000021E-4     1.19999999999999989E-3
+    355     1.00000000000000002E-3     1.13896484375000011E-3
+    356     9.22298701458379552E-4     1.09448242187499991E-3
+    357     1.00000000000000002E-3     1.04999999999999994E-3
+    358     9.61149350729189732E-4     1.1167236328124999E-3
+    359     9.61149350729189732E-4     1.07224121093749993E-3
+    360     1.00000000000000002E-3     1.09448242187500013E-3
+    361     9.22298701458379552E-4     1.09448242187499991E-3
+    362     8.42471384914615607E-4     1.00872691912019851E-3
+    363     9.12446330413049594E-4     9.87167427587815364E-4
+    364     8.82385043186497634E-4     1.0516046704975991E-3
+    365     8.77458857663832546E-4     9.97947173354006827E-4
+    366     9.17372515935714573E-4     1.04082492473140764E-3
+    367     1.00000000000000002E-3     1.13896484375000011E-3
+    368     1.00000000000000002E-3     1.19999999999999989E-3
+    369     9.38964843750000021E-4     1.19999999999999989E-3
+    370     1.00000000000000002E-3     1.16948242187500011E-3
+    371     9.69482421875000021E-4     1.19999999999999989E-3
+    372     9.69482421875000021E-4     1.16948242187500011E-3
+    373     9.12446330413049594E-4     9.87167427587815364E-4
+    374     8.42471384914615607E-4     1.00872691912019851E-3
+    375     8.1517205990673709E-4     9.40785388914013596E-4
+    376     8.77458857663832546E-4     9.97947173354006827E-4
+    377     8.28821722410676403E-4     9.74756154017106051E-4
+    378     8.63809195159893342E-4     9.6397640825091448E-4
+    379     9.00310858258578334E-4     8.8742365740778872E-4
+    380     9.12446330413049594E-4     9.87167427587815364E-4
+    381     8.1517205990673709E-4     9.40785388914013596E-4
+    382     9.06378594335813964E-4     9.37295542497801988E-4
+    383     8.63809195159893342E-4     9.6397640825091448E-4
+    384     8.57741459082657712E-4     9.14104523160901104E-4
+    385     9.12446330413049594E-4     9.87167427587815364E-4
+    386     1.00000000000000002E-3     9.0E-4
+    387     1.00000000000000002E-3     9.74999999999999955E-4
+    388     9.56223165206524861E-4     9.4358371379390767E-4
+    389     1.00000000000000002E-3     9.37499999999999965E-4
+    390     9.56223165206524861E-4     9.8108371379390766E-4
+    391     7.24825013803433236E-4     8.92866282184746245E-4
+    392     8.03413482021300262E-4     8.27424077205904536E-4
+    393     8.37129054430499212E-4     8.81217371382065837E-4
+    394     7.64119247912366803E-4     8.60145179695325337E-4
+    395     8.20271268225899737E-4     8.54320724293985241E-4
+    396     7.8097703411696617E-4     8.87041826783406041E-4
+    397     7.52655853090133443E-4     7.19787645578877193E-4
+    398     7.56891934562431879E-4     7.83576045089881936E-4
+    399     6.94054151555351649E-4     7.95336883715192389E-4
+    400     7.54773893826282661E-4     7.51681845334379564E-4
+    401     7.25473043058891764E-4     7.89456464402537163E-4
+    402     7.23355002322742546E-4     7.57562264647034791E-4
+    403     1.00000000000000002E-3     8.24999999999999995E-4
+    404     9.36226394521655E-4     8.3396748156078118E-4
+    405     1.00000000000000002E-3     7.50000000000000016E-4
+    406     9.68113197260827511E-4     8.29483740780390588E-4
+    407     9.68113197260827511E-4     7.91983740780390598E-4
+    408     1.00000000000000002E-3     7.87500000000000006E-4
+    409     8.9886962494534965E-4     7.81508352249731973E-4
+    410     9.36226394521655E-4     8.3396748156078118E-4
+    411     9.00310858258578334E-4     8.8742365740778872E-4
+    412     9.17548009733502325E-4     8.07737916905256577E-4
+    413     9.18268626390116721E-4     8.60695569484284896E-4
+    414     8.99590241601963938E-4     8.34466004828760292E-4
+    415     9.00310858258578334E-4     8.8742365740778872E-4
+    416     9.36226394521655E-4     8.3396748156078118E-4
+    417     1.00000000000000002E-3     9.0E-4
+    418     9.18268626390116721E-4     8.60695569484284896E-4
+    419     9.68113197260827511E-4     8.66983740780390578E-4
+    420     9.50155429129289231E-4     8.93711828703894293E-4
+    421     8.03413482021300262E-4     8.27424077205904536E-4
+    422     7.56891934562431879E-4     7.83576045089881936E-4
+    423     8.35489473019320351E-4     7.27605773008054794E-4
+    424     7.80152708291866016E-4     8.05500061147893182E-4
+    425     7.9619070379087606E-4     7.55590909048968365E-4
+    426     8.19451477520310306E-4     7.77514925106979719E-4
+    427     9.37074205321842733E-4     6.1279086004113368E-4
+    428     1.00000000000000002E-3     6.75000000000000036E-4
+    429     9.21159286470410715E-4     6.74999999999999927E-4
+    430     9.68537102660921431E-4     6.43895430020566912E-4
+    431     9.60579643235205368E-4     6.75000000000000036E-4
+    432     9.29116745896126778E-4     6.43895430020566804E-4
+    433     9.21159286470410715E-4     6.74999999999999927E-4
+    434     9.39685899131729899E-4     7.34463702578968722E-4
+    435     8.35489473019320351E-4     7.27605773008054794E-4
+    436     9.30422592801070361E-4     7.04731851289484325E-4
+    437     8.87587686075525179E-4     7.31034737793511812E-4
+    438     8.78324379744865533E-4     7.01302886504027306E-4
+    439     9.37074205321842733E-4     6.1279086004113368E-4
+    440     8.14703805291646039E-4     6.47042209754021894E-4
+    441     8.92816374246439338E-4     5.66266606957709193E-4
+    442     8.75889005306744331E-4     6.29916534897577733E-4
+    443     8.53760089769042742E-4     6.06654408355865489E-4
+    444     9.14945289784141035E-4     5.89528733499421437E-4
+    445     8.14703805291646039E-4     6.47042209754021894E-4
+    446     7.90472897978866684E-4     5.45453687327126634E-4
+    447     8.92816374246439338E-4     5.66266606957709193E-4
+    448     8.02588351635256415E-4     5.96247948540574264E-4
+    449     8.41644636112653011E-4     5.55860147142417859E-4
+    450     8.53760089769042742E-4     6.06654408355865489E-4
+    451     8.08703962690366035E-4     4.17299929575055159E-4
+    452     7.47347240742430642E-4     3.78837751313385344E-4
+    453     7.77749415513451876E-4     3.13113370342386432E-4
+    454     7.78025601716398393E-4     3.98068840444220251E-4
+    455     7.62548328127941313E-4     3.45975560827885861E-4
+    456     7.93226689101908956E-4     3.65206649958720823E-4
+    457     9.33471456834415085E-4     4.87499999999999978E-4
+    458     1.00000000000000002E-3     4.49999999999999988E-4
+    459     1.00000000000000002E-3     5.24999999999999968E-4
+    460     9.66735728417207553E-4     4.68749999999999983E-4
+    461     1.00000000000000002E-3     4.87499999999999978E-4
+    462     9.66735728417207553E-4     5.06249999999999973E-4
+    463     9.37074205321842733E-4     6.1279086004113368E-4
+    464     1.00000000000000002E-3     5.24999999999999968E-4
+    465     1.00000000000000002E-3     5.99999999999999947E-4
+    466     9.68537102660921431E-4     5.68895430020566824E-4
+    467     1.00000000000000002E-3     5.62499999999999957E-4
+    468     9.68537102660921431E-4     6.06395430020566814E-4
+    469     6.15127063129517924E-4     1.28961775828680109E-4
+    470     5.52246093750000003E-4     0
+    471     6.16210937500000005E-4     0
+    472     5.83686578439758909E-4     6.44808879143400545E-5
+    473     5.8422851562499995E-4     0
+    474     6.15669000314758965E-4     6.44808879143400545E-5
+    475     9.30681494601841297E-4     1.85796537326297451E-4
+    476     9.35392068859210385E-4     6.524280030581833E-5
+    477     1.00000000000000002E-3     1.49999999999999987E-4
+    478     9.33036781730525841E-4     1.2551966881605789E-4
+    479     9.67696034429605203E-4     1.07621400152909158E-4
+    480     9.65340747300920659E-4     1.67898268663148732E-4
+    481     1.00000000000000002E-3     2.99999999999999974E-4
+    482     9.33976498687064506E-4     3.35267605171179269E-4
+    483     9.16432772236257481E-4     2.62499999999999984E-4
+    484     9.66988249343532264E-4     3.17633802585589621E-4
+    485     9.2520463546166094E-4     2.98883802585589626E-4
+    486     9.58216386118128805E-4     2.81249999999999979E-4
+    487     8.79865520605121044E-4     9.96839116576813268E-5
+    488     8.13728732664855108E-4     1.33343473922578902E-4
+    489     7.69147152428217916E-4     7.17795642973629926E-5
+    490     8.4679712663498813E-4     1.16513692790130114E-4
+    491     7.91437942546536566E-4     1.02561519109970947E-4
+    492     8.2450633651666948E-4     8.57317379775221597E-5
+    493     6.86868266859239985E-4     5.40099367412093744E-4
+    494     7.20535236875600961E-4     6.29731480771863312E-4
+    495     6.25000000000000013E-4     6.23375449734813098E-4
+    496     7.03701751867420473E-4     5.84915424091978528E-4
+    497     6.72767618437800433E-4     6.26553465253338205E-4
+    498     6.55934133429619945E-4     5.81737408573453421E-4
+    499     8.58425917146639784E-4     4.73341467989030294E-4
+    500     8.08703962690366035E-4     4.17299929575055159E-4
+    501     8.70539573309993481E-4     3.74999999999999954E-4
+    502     8.3356493991850291E-4     4.45320698782042726E-4
+    503     8.39621768000179758E-4     3.96149964787527583E-4
+    504     8.64482745228316632E-4     4.24170733994515096E-4
+    505     7.20535236875600961E-4     6.29731480771863312E-4
+    506     7.90472897978866684E-4     5.45453687327126634E-4
+    507     8.14703805291646039E-4     6.47042209754021894E-4
+    508     7.55504067427233876E-4     5.87592584049495027E-4
+    509     8.02588351635256415E-4     5.96247948540574264E-4
+    510     7.676195210836235E-4     6.38386845262942657E-4
+    511     5.77009580214111854E-4     5.23102021510448161E-4
+    512     5.82521566988404701E-4     4.12072931624541923E-4
+    513     6.55870398270680236E-4     4.58827279636259491E-4
+    514     5.79765573601258278E-4     4.67587476567495042E-4
+    515     6.19195982629542523E-4     4.3545010563040068E-4
+    516     6.16439989242395991E-4     4.90964650573353853E-4
+    517     4.15683819919475818E-4     3.45889508016053047E-4
+    518     4.57065334871855982E-4     4.13720111675068352E-4
+    519     3.78329331029246172E-4     4.03040117695809171E-4
+    520     4.363745773956659E-4     3.798048098455607E-4
+    521     4.17697332950551104E-4     4.08380114685438734E-4
+    522     3.97006575474360968E-4     3.74464812855931136E-4
+    523     6.149060889867273E-4     3.58911302456159137E-4
+    524     5.8651429788997224E-4     3.03514350561540362E-4
+    525     6.48721061501020204E-4     3.0580228455003082E-4
+    526     6.00710193438349716E-4     3.3121282650884975E-4
+    527     6.17617679695496276E-4     3.04658317555785564E-4
+    528     6.31813575243873752E-4     3.32356793503094951E-4
+    529     5.77009580214111854E-4     5.23102021510448161E-4
+    530     6.27455298609753205E-4     5.60579639099527925E-4
+    531     5.59110560621724098E-4     5.93289128784703487E-4
+    532     6.02232439411932476E-4     5.41840830304987989E-4
+    533     5.93282929615738652E-4     5.76934383942115706E-4
+    534     5.6806007041791803E-4     5.58195575147575878E-4
+    535     5.18154572622838414E-4     4.64528827497301226E-4
+    536     5.82521566988404701E-4     4.12072931624541923E-4
+    537     5.77009580214111854E-4     5.23102021510448161E-4
+    538     5.50338069805621612E-4     4.38300879560921575E-4
+    539     5.79765573601258278E-4     4.67587476567495042E-4
+    540     5.4758207641847508E-4     4.93815424503874639E-4
+    541     4.9377437948036596E-4     6.24558736453494859E-4
+    542     5.59110560621724098E-4     5.93289128784703487E-4
+    543     5.60394849123839033E-4     7.35715772264671007E-4
+    544     5.26442470051044975E-4     6.08923932619099173E-4
+    545     5.5975270487278162E-4     6.64502450524687193E-4
+    546     5.27084614302102497E-4     6.80137254359082987E-4
+    547     3.79565818807853228E-4     5.99785262317693206E-4
+    548     4.45119450916024088E-4     5.73219517420557926E-4
+    549     4.62057889734597506E-4     6.86398137230062556E-4
+    550     4.12342634861938631E-4     5.86502389869125512E-4
+    551     4.5358867032531077E-4     6.29808827325310187E-4
+    552     4.20811854271225367E-4     6.43091699773877881E-4
+    553     5.60394849123839033E-4     7.35715772264671007E-4
+    554     6.25000000000000013E-4     6.23375449734813098E-4
+    555     6.60634466318554738E-4     7.08911474490989688E-4
+    556     5.92697424561919577E-4     6.79545610999742053E-4
+    557     6.4281723315927743E-4     6.66143462112901447E-4
+    558     6.10514657721196886E-4     7.22313623377830293E-4
+    559     5.19966569338154057E-4     8.58840618576917807E-4
+    560     4.10887142695853748E-4     8.117939619078767E-4
+    561     4.75850298603600714E-4     7.54514336975707829E-4
+    562     4.65426856017003929E-4     8.35317290242397308E-4
+    563     4.43368720649727258E-4     7.83154149441792264E-4
+    564     4.97908433970877386E-4     8.06677477776312872E-4
+    565     7.42631330235908388E-4     4.52617545905931022E-4
+    566     8.58425917146639784E-4     4.73341467989030294E-4
+    567     7.90472897978866684E-4     5.45453687327126634E-4
+    568     8.00528623691274086E-4     4.62979506947480658E-4
+    569     8.2444940756275318E-4     5.09397577658078464E-4
+    570     7.6655211410738759E-4     4.99035616616528774E-4
+    571     9.33471456834415085E-4     4.87499999999999978E-4
+    572     8.92816374246439338E-4     5.66266606957709193E-4
+    573     8.58425917146639784E-4     4.73341467989030294E-4
+    574     9.13143915540427266E-4     5.26883303478854585E-4
+    575     8.75621145696539507E-4     5.19804037473369689E-4
+    576     8.95948686990527435E-4     4.80420733994515136E-4
+    577     7.90472897978866684E-4     5.45453687327126634E-4
+    578     7.20535236875600961E-4     6.29731480771863312E-4
+    579     6.86868266859239985E-4     5.40099367412093744E-4
+    580     7.55504067427233876E-4     5.87592584049495027E-4
+    581     7.03701751867420473E-4     5.84915424091978528E-4
+    582     7.38670582419053388E-4     5.42776527369610243E-4
+    583     7.52655853090133443E-4     7.19787645578877193E-4
+    584     7.20535236875600961E-4     6.29731480771863312E-4
+    585     8.14703805291646039E-4     6.47042209754021894E-4
+    586     7.36595544982867256E-4     6.74759563175370307E-4
+    587     7.676195210836235E-4     6.38386845262942657E-4
+    588     7.83679829190889795E-4     6.83414927666449543E-4
+    589     7.52655853090133443E-4     7.19787645578877193E-4
+    590     8.14703805291646039E-4     6.47042209754021894E-4
+    591     8.35489473019320351E-4     7.27605773008054794E-4
+    592     7.83679829190889795E-4     6.83414927666449543E-4
+    593     8.25096639155483195E-4     6.87323991381038344E-4
+    594     7.94072663054726843E-4     7.23696709293465993E-4
+    595     6.60634466318554738E-4     7.08911474490989688E-4
+    596     7.52655853090133443E-4     7.19787645578877193E-4
+    597     6.94054151555351649E-4     7.95336883715192389E-4
+    598     7.06645159704344036E-4     7.1434956003493344E-4
+    599     7.23355002322742546E-4     7.57562264647034791E-4
+    600     6.77344308936953139E-4     7.52124179103091039E-4
+    601     6.60634466318554738E-4     7.08911474490989688E-4
+    602     6.94054151555351649E-4     7.95336883715192389E-4
+    603     6.3251725796544824E-4     7.69458270396488286E-4
+    604     6.77344308936953139E-4     7.52124179103091039E-4
+    605     6.63285704760399998E-4     7.82397577055840392E-4
+    606     6.46575862142001489E-4     7.39184872443738933E-4
+    607     7.20535236875600961E-4     6.29731480771863312E-4
+    608     6.60634466318554738E-4     7.08911474490989688E-4
+    609     6.25000000000000013E-4     6.23375449734813098E-4
+    610     6.90584851597077849E-4     6.69321477631426446E-4
+    611     6.4281723315927743E-4     6.66143462112901447E-4
+    612     6.72767618437800433E-4     6.26553465253338205E-4
+    613     7.29985750852659397E-4     9.92266244781275222E-4
+    614     7.24825013803433236E-4     8.92866282184746245E-4
+    615     8.1517205990673709E-4     9.40785388914013596E-4
+    616     7.27405382328046371E-4     9.42566263483010734E-4
+    617     7.69998536855085163E-4     9.1682583554937992E-4
+    618     7.72578905379698189E-4     9.66525816847644409E-4
+    619     7.24825013803433236E-4     8.92866282184746245E-4
+    620     6.8729190829219008E-4     9.44648911063777994E-4
+    621     6.28726770153248599E-4     8.69566754220071274E-4
+    622     7.06058461047811658E-4     9.1875759662426212E-4
+    623     6.5800933922271934E-4     9.07107832641924688E-4
+    624     6.76775891978340918E-4     8.81216518202408705E-4
+    625     8.42471384914615607E-4     1.00872691912019851E-3
+    626     7.70557602221913452E-4     1.04998903577999537E-3
+    627     7.29985750852659397E-4     9.92266244781275222E-4
+    628     8.0651449356826453E-4     1.02935797745009683E-3
+    629     7.50271676537286425E-4     1.0211276402806353E-3
+    630     7.86228567883637502E-4     1.00049658195073676E-3
+    631     7.29985750852659397E-4     9.92266244781275222E-4
+    632     6.91894531250000018E-4     1.08422299961888672E-3
+    633     6.3613763360561681E-4     1.03946949707743879E-3
+    634     7.10940141051329653E-4     1.03824462220008097E-3
+    635     6.64016082427808414E-4     1.06184624834816275E-3
+    636     6.83061692229138158E-4     1.015867870929357E-3
+    637     5.42759630809160629E-4     9.44364496429779174E-4
+    638     6.28726770153248599E-4     8.69566754220071274E-4
+    639     6.27929687500000016E-4     9.68445999237773649E-4
+    640     5.85743200481204668E-4     9.06965625324925224E-4
+    641     6.28328228826624362E-4     9.19006376728922461E-4
+    642     5.85344659154580322E-4     9.56405247833776412E-4
+    643     5.79668274240937665E-4     8.10244013239797283E-4
+    644     6.28726770153248599E-4     8.69566754220071274E-4
+    645     5.19966569338154057E-4     8.58840618576917807E-4
+    646     6.04197522197093186E-4     8.39905383729934278E-4
+    647     5.74346669745701274E-4     8.64203686398494486E-4
+    648     5.49817421789545861E-4     8.34542315908357599E-4
+    649     3.14898739406715305E-4     6.69780733668936504E-4
+    650     3.79565818807853228E-4     5.99785262317693206E-4
+    651     3.98747049720822426E-4     7.15064925451776041E-4
+    652     3.47232279107284267E-4     6.34782997993314909E-4
+    653     3.891564342643378E-4     6.57425093884734624E-4
+    654     3.56822894563768893E-4     6.92422829560356218E-4
+    655     3.75421084408168141E-4     5.26840780869595538E-4
+    656     2.50000000000000005E-4     5.99999999999999947E-4
+    657     3.14594753959753046E-4     4.86365680855438756E-4
+    658     3.12710542204084046E-4     5.63420390434797797E-4
+    659     2.82297376979876499E-4     5.43182840427719298E-4
+    660     3.45007919183960593E-4     5.06603230862517147E-4
+    661     7.66063425991097419E-4     2.25618978565681171E-4
+    662     8.13728732664855108E-4     1.33343473922578902E-4
+    663     8.53348088789513369E-4     1.96091760522196289E-4
+    664     7.89896079327976263E-4     1.79481226244130037E-4
+    665     8.33538410727184239E-4     1.64717617222387609E-4
+    666     8.09705757390305394E-4     2.10855369543938717E-4
+    667     1.73262011544884361E-4     2.84564040120359029E-4
+    668     2.31915182548602798E-4     3.185630410561948E-4
+    669     1.71301572768688224E-4     3.7792968750000001E-4
+    670     2.0258859704674358E-4     3.01563540588276914E-4
+    671     2.01608377658645498E-4     3.48246364278097405E-4
+    672     1.72281792156786279E-4     3.3124686381017952E-4
+    673     6.4572363977386294E-5     3.7792968750000001E-4
+    674     1.71301572768688224E-4     3.7792968750000001E-4
+    675     8.56507863843441119E-5     4.3896484375000001E-4
+    676     1.17936968373037266E-4     3.7792968750000001E-4
+    677     1.28476179576516181E-4     4.0844726562500001E-4
+    678     7.51115751808651962E-5     4.0844726562500001E-4
+    679     8.56507863843441119E-5     4.3896484375000001E-4
+    680     1.25000000000000003E-4     5.49999999999999925E-4
+    681     6.25000000000000013E-5     5.24999999999999968E-4
+    682     1.05325393192172057E-4     4.94482421874999968E-4
+    683     9.3750000000000002E-5     5.375E-4
+    684     7.40753931921720566E-5     4.81982421874999989E-4
+    685     0     5.0000000000000001E-4
+    686     0     4.3896484375000001E-4
+    687     8.56507863843441119E-5     4.3896484375000001E-4
+    688     0     4.6948242187500001E-4
+    689     4.28253931921720559E-5     4.3896484375000001E-4
+    690     4.28253931921720559E-5     4.6948242187500001E-4
+    691     8.03413482021300262E-4     8.27424077205904536E-4
+    692     8.35489473019320351E-4     7.27605773008054794E-4
+    693     8.9886962494534965E-4     7.81508352249731973E-4
+    694     8.19451477520310306E-4     7.77514925106979719E-4
+    695     8.67179548982335E-4     7.54557062628893383E-4
+    696     8.51141553483324956E-4     8.044662147278182E-4
+    697     0     2.5585937500000001E-4
+    698     9.02731575853669911E-5     2.19793387486037181E-4
+    699     8.56507863843441119E-5     3.1689453125000001E-4
+    700     4.51365787926834956E-5     2.37826381243018609E-4
+    701     8.79619719848555515E-5     2.68343959368018609E-4
+    702     4.28253931921720559E-5     2.8637695312500001E-4
+    703     6.4572363977386294E-5     3.7792968750000001E-4
+    704     0     3.1689453125000001E-4
+    705     8.56507863843441119E-5     3.1689453125000001E-4
+    706     3.2286181988693147E-5     3.4741210937500001E-4
+    707     4.28253931921720559E-5     3.1689453125000001E-4
+    708     7.51115751808651962E-5     3.4741210937500001E-4
+    709     0     1.91894531250000008E-4
+    710     1.01956056377857214E-4     1.1802584118703223E-4
+    711     9.02731575853669911E-5     2.19793387486037181E-4
+    712     5.09780281889286071E-5     1.54960186218516132E-4
+    713     9.61146069816121027E-5     1.68909614336534706E-4
+    714     4.51365787926834956E-5     2.05843959368018608E-4
+    715     1.73262011544884361E-4     2.84564040120359029E-4
+    716     9.02731575853669911E-5     2.19793387486037181E-4
+    717     1.80546315170733955E-4     1.83727399972074352E-4
+    718     1.31767584565125663E-4     2.52178713803198119E-4
+    719     1.3540973637805046E-4     2.01760393729055753E-4
+    720     1.76904163357809158E-4     2.34145720046216691E-4
+    721     0     1.27929687500000005E-4
+    722     1.01956056377857214E-4     1.1802584118703223E-4
+    723     0     1.91894531250000008E-4
+    724     5.09780281889286071E-5     1.22977764343516131E-4
+    725     5.09780281889286071E-5     1.54960186218516132E-4
+    726     0     1.59912109375000007E-4
+    727     1.80546315170733955E-4     1.83727399972074352E-4
+    728     1.01956056377857214E-4     1.1802584118703223E-4
+    729     1.95469003684387412E-4     8.60228725173496005E-5
+    730     1.41251185774295585E-4     1.50876620579553278E-4
+    731     1.487125300311223E-4     1.02024356852190915E-4
+    732     1.88007659427560697E-4     1.34875136244711976E-4
+    733     0     6.39648437500000026E-5
+    734     0     0
+    735     1.220703125E-4     0
+    736     0     3.19824218750000013E-5
+    737     6.103515625E-5     0
+    738     6.103515625E-5     3.19824218750000013E-5
+    739     1.35354037059463752E-4     6.29907264959934823E-5
+    740     2.44140625E-4     0
+    741     1.95469003684387412E-4     8.60228725173496005E-5
+    742     1.89747331029731876E-4     3.14953632479967412E-5
+    743     2.19804814342193706E-4     4.30114362586748002E-5
+    744     1.65411520371925582E-4     7.45067995066715414E-5
+    745     1.35354037059463752E-4     6.29907264959934823E-5
+    746     1.95469003684387412E-4     8.60228725173496005E-5
+    747     1.01956056377857214E-4     1.1802584118703223E-4
+    748     1.65411520371925582E-4     7.45067995066715414E-5
+    749     1.487125300311223E-4     1.02024356852190915E-4
+    750     1.18655046718660483E-4     9.05082838415128562E-5
+    751     3.64813480721653041E-4     1.25057137488653904E-4
+    752     2.70476109685850318E-4     1.47470775643085319E-4
+    753     3.00937473641497586E-4     6.59440565847406544E-5
+    754     3.17644795203751652E-4     1.36263956565869611E-4
+    755     2.85706791663673952E-4     1.06707416113912987E-4
+    756     3.32875477181575341E-4     9.55005970366972791E-5
+    757     4.2724609375E-4     0
+    758     4.8828125E-4     0
+    759     4.20397371381042544E-4     7.46936636915224194E-5
+    760     4.57763671875E-4     0
+    761     4.54339310690521272E-4     3.73468318457612097E-5
+    762     4.23821732565521272E-4     3.73468318457612097E-5
+    763     7.66063425991097419E-4     2.25618978565681171E-4
+    764     8.53348088789513369E-4     1.96091760522196289E-4
+    765     8.40965569411703819E-4     2.71022418711743549E-4
+    766     8.09705757390305394E-4     2.10855369543938717E-4
+    767     8.47156829100608648E-4     2.33557089616969919E-4
+    768     8.03514497701400673E-4     2.48320698638712374E-4
+    769     5.1540545604394721E-4     2.10741452330906448E-4
+    770     4.25259790828855515E-4     1.93997913774442133E-4
+    771     4.88018586665576635E-4     1.0715009043865536E-4
+    772     4.70332623436401363E-4     2.02369683052674304E-4
+    773     4.56639188747216102E-4     1.50574002106548739E-4
+    774     5.01712021354761895E-4     1.58945771384780911E-4
+    775     5.8651429788997224E-4     3.03514350561540362E-4
+    776     6.149060889867273E-4     3.58911302456159137E-4
+    777     5.469449348011982E-4     3.56411724829354306E-4
+    778     6.00710193438349716E-4     3.3121282650884975E-4
+    779     5.80925511893962696E-4     3.57661513642756695E-4
+    780     5.6672961634558522E-4     3.29963037695447307E-4
+    781     3.662109375E-4     0
+    782     4.20397371381042544E-4     7.46936636915224194E-5
+    783     3.00937473641497586E-4     6.59440565847406544E-5
+    784     3.93304154440521272E-4     3.73468318457612097E-5
+    785     3.60667422511270065E-4     7.03188601381315369E-5
+    786     3.33574205570748793E-4     3.29720282923703272E-5
+    787     2.70476109685850318E-4     1.47470775643085319E-4
+    788     3.3457565114347566E-4     2.07524696466039707E-4
+    789     2.62862460053854216E-4     2.58243951636348947E-4
+    790     3.02525880414662989E-4     1.774977360545625E-4
+    791     2.98719055598664911E-4     2.32884324051194341E-4
+    792     2.66669284869852294E-4     2.02857363639717133E-4
+    793     2.31915182548602798E-4     3.185630410561948E-4
+    794     3.03283108388862287E-4     3.07191821602677501E-4
+    795     2.80075730378807853E-4     3.66277694634743636E-4
+    796     2.67599145468732529E-4     3.12877431329436177E-4
+    797     2.91679419383835043E-4     3.36734758118710596E-4
+    798     2.55995456463705339E-4     3.42420367845469218E-4
+    799     1.73262011544884361E-4     2.84564040120359029E-4
+    800     1.71301572768688224E-4     3.7792968750000001E-4
+    801     8.56507863843441119E-5     3.1689453125000001E-4
+    802     1.72281792156786279E-4     3.3124686381017952E-4
+    803     1.28476179576516181E-4     3.4741210937500001E-4
+    804     1.29456398964614236E-4     3.0072928568517952E-4
+    805     2.80075730378807853E-4     3.66277694634743636E-4
+    806     3.03283108388862287E-4     3.07191821602677501E-4
+    807     3.78329331029246172E-4     4.03040117695809171E-4
+    808     2.91679419383835043E-4     3.36734758118710596E-4
+    809     3.40806219709054257E-4     3.55115969649243308E-4
+    810     3.29202530704027013E-4     3.84658906165276403E-4
+    811     1.71301572768688224E-4     3.7792968750000001E-4
+    812     2.31175962939123751E-4     4.23329063260392325E-4
+    813     1.72756436172360502E-4     4.70585859703467847E-4
+    814     2.01238767853905988E-4     4.00629375380196168E-4
+    815     2.0196619955574214E-4     4.46957461481930086E-4
+    816     1.72029004470524349E-4     4.24257773601733929E-4
+    817     3.78329331029246172E-4     4.03040117695809171E-4
+    818     3.662109375E-4     2.98836517333984323E-4
+    819     4.15683819919475818E-4     3.45889508016053047E-4
+    820     3.72270134264623086E-4     3.5093831751489672E-4
+    821     3.90947378709737936E-4     3.22363012675018685E-4
+    822     3.97006575474360968E-4     3.74464812855931136E-4
+    823     3.75421084408168141E-4     5.26840780869595538E-4
+    824     3.78329331029246172E-4     4.03040117695809171E-4
+    825     4.44306653512529842E-4     5.02492175855262643E-4
+    826     3.76875207718707183E-4     4.64940449282702354E-4
+    827     4.1131799227088798E-4     4.52766146775535907E-4
+    828     4.09863868960348991E-4     5.1466647836242909E-4
+    829     4.58849660034789083E-4     2.64192897481815869E-4
+    830     3.662109375E-4     2.98836517333984323E-4
+    831     4.25259790828855515E-4     1.93997913774442133E-4
+    832     4.12530298767394541E-4     2.81514707407900123E-4
+    833     3.95735364164427757E-4     2.46417215554213228E-4
+    834     4.42054725431822299E-4     2.29095405628129001E-4
+    835     5.1540545604394721E-4     2.10741452330906448E-4
+    836     6.03238036253151279E-4     2.17057585313148337E-4
+    837     5.8651429788997224E-4     3.03514350561540362E-4
+    838     5.5932174614854919E-4     2.13899518822027393E-4
+    839     5.94876167071561705E-4     2.60285967937344336E-4
+    840     5.50959876966959725E-4     2.57127901446223392E-4
+    841     4.25259790828855515E-4     1.93997913774442133E-4
+    842     3.3457565114347566E-4     2.07524696466039707E-4
+    843     3.64813480721653041E-4     1.25057137488653904E-4
+    844     3.79917720986165588E-4     2.00761305120240906E-4
+    845     3.49694565932564378E-4     1.66290916977346792E-4
+    846     3.95036635775254251E-4     1.59527525631548018E-4
+    847     3.03283108388862287E-4     3.07191821602677501E-4
+    848     3.3457565114347566E-4     2.07524696466039707E-4
+    849     3.662109375E-4     2.98836517333984323E-4
+    850     3.18929379766168947E-4     2.5735825903435859E-4
+    851     3.5039329432173783E-4     2.53180606900012001E-4
+    852     3.34747022944431117E-4     3.03014169468330912E-4
+    853     3.03283108388862287E-4     3.07191821602677501E-4
+    854     2.31915182548602798E-4     3.185630410561948E-4
+    855     2.62862460053854216E-4     2.58243951636348947E-4
+    856     2.67599145468732529E-4     3.12877431329436177E-4
+    857     2.47388821301228494E-4     2.88403496346271873E-4
+    858     2.83072784221358252E-4     2.82717886619513197E-4
+    859     4.25259790828855515E-4     1.93997913774442133E-4
+    860     5.1540545604394721E-4     2.10741452330906448E-4
+    861     4.58849660034789083E-4     2.64192897481815869E-4
+    862     4.70332623436401363E-4     2.02369683052674304E-4
+    863     4.87127558039368146E-4     2.37467174906361145E-4
+    864     4.42054725431822299E-4     2.29095405628129001E-4
+    865     4.20397371381042544E-4     7.46936636915224194E-5
+    866     4.25259790828855515E-4     1.93997913774442133E-4
+    867     3.64813480721653041E-4     1.25057137488653904E-4
+    868     4.22828581104949029E-4     1.34345788732982262E-4
+    869     3.95036635775254251E-4     1.59527525631548018E-4
+    870     3.9260542605134782E-4     9.98754005900881616E-5
+    871     4.88018586665576635E-4     1.0715009043865536E-4
+    872     5.46592219681165661E-4     1.47080610942937849E-4
+    873     5.1540545604394721E-4     2.10741452330906448E-4
+    874     5.17305403173371121E-4     1.27115350690796598E-4
+    875     5.30998837862556436E-4     1.78911031636922162E-4
+    876     5.01712021354761895E-4     1.58945771384780911E-4
+    877     4.88018586665576635E-4     1.0715009043865536E-4
+    878     4.8828125E-4     0
+    879     5.52246093750000003E-4     0
+    880     4.8814991833278829E-4     5.35750452193276799E-5
+    881     5.20263671875000056E-4     0
+    882     5.20132340207788346E-4     5.35750452193276799E-5
+    883     6.80175781250000008E-4     6.50230427205121634E-5
+    884     6.16210937500000005E-4     0
+    885     7.4414062500000001E-4     0
+    886     6.48193359375000061E-4     3.25115213602560817E-5
+    887     6.80175781250000008E-4     0
+    888     7.12158203124999955E-4     3.25115213602560817E-5
+    889     6.03238036253151279E-4     2.17057585313148337E-4
+    890     6.83188905983988506E-4     2.53967124884082096E-4
+    891     6.48721061501020204E-4     3.0580228455003082E-4
+    892     6.43213471118569838E-4     2.3551235509861523E-4
+    893     6.65954983742504409E-4     2.79884704717056458E-4
+    894     6.25979548877085741E-4     2.61429934931589592E-4
+    895     3.00937473641497586E-4     6.59440565847406544E-5
+    896     1.95469003684387412E-4     8.60228725173496005E-5
+    897     2.44140625E-4     0
+    898     2.48203238662942499E-4     7.59834645510451274E-5
+    899     2.19804814342193706E-4     4.30114362586748002E-5
+    900     2.72539049320748793E-4     3.29720282923703272E-5
+    901     1.95469003684387412E-4     8.60228725173496005E-5
+    902     3.00937473641497586E-4     6.59440565847406544E-5
+    903     2.70476109685850318E-4     1.47470775643085319E-4
+    904     2.48203238662942499E-4     7.59834645510451274E-5
+    905     2.85706791663673952E-4     1.06707416113912987E-4
+    906     2.32972556685118865E-4     1.1674682408021746E-4
+    907     3.662109375E-4     0
+    908     3.00937473641497586E-4     6.59440565847406544E-5
+    909     2.44140625E-4     0
+    910     3.33574205570748793E-4     3.29720282923703272E-5
+    911     2.72539049320748793E-4     3.29720282923703272E-5
+    912     3.0517578125E-4     0
+    913     9.02731575853669911E-5     2.19793387486037181E-4
+    914     1.73262011544884361E-4     2.84564040120359029E-4
+    915     8.56507863843441119E-5     3.1689453125000001E-4
+    916     1.31767584565125663E-4     2.52178713803198119E-4
+    917     1.29456398964614236E-4     3.0072928568517952E-4
+    918     8.79619719848555515E-5     2.68343959368018609E-4
+    919     1.80546315170733955E-4     1.83727399972074352E-4
+    920     2.62862460053854216E-4     2.58243951636348947E-4
+    921     1.73262011544884361E-4     2.84564040120359029E-4
+    922     2.21704387612294072E-4     2.2098567580421165E-4
+    923     2.18062235799369275E-4     2.71403995878353988E-4
+    924     1.76904163357809158E-4     2.34145720046216691E-4
+    925     5.8651429788997224E-4     3.03514350561540362E-4
+    926     4.83502295899330232E-4     3.38002535036003411E-4
+    927     4.58849660034789083E-4     2.64192897481815869E-4
+    928     5.35008296894651236E-4     3.20758442798771914E-4
+    929     4.71175977967059657E-4     3.01097716258909667E-4
+    930     5.22681978962380661E-4     2.83853624021678116E-4
+    931     6.48721061501020204E-4     3.0580228455003082E-4
+    932     6.83188905983988506E-4     2.53967124884082096E-4
+    933     7.14223074903462172E-4     3.09513800853735191E-4
+    934     6.65954983742504409E-4     2.79884704717056458E-4
+    935     6.98705990443725285E-4     2.81740462868908617E-4
+    936     6.81472068202241188E-4     3.07658042701882978E-4
+    937     6.48721061501020204E-4     3.0580228455003082E-4
+    938     7.14223074903462172E-4     3.09513800853735191E-4
+    939     6.77112852597775263E-4     3.61199236444649594E-4
+    940     6.81472068202241188E-4     3.07658042701882978E-4
+    941     6.95667963750618663E-4     3.35356518649192366E-4
+    942     6.62916957049397787E-4     3.33500760497340207E-4
+    943     6.55870398270680236E-4     4.58827279636259491E-4
+    944     6.77112852597775263E-4     3.61199236444649594E-4
+    945     7.42631330235908388E-4     4.52617545905931022E-4
+    946     6.66491625434227804E-4     4.10013258040454543E-4
+    947     7.09872091416841825E-4     4.06908391175290281E-4
+    948     6.99250864253294366E-4     4.55722412771095283E-4
+    949     6.83188905983988506E-4     2.53967124884082096E-4
+    950     6.73095146243751124E-4     1.78431349377634624E-4
+    951     7.66063425991097419E-4     2.25618978565681171E-4
+    952     6.78142026113869815E-4     2.1619923713085836E-4
+    953     7.19579286117424271E-4     2.02025163971657911E-4
+    954     7.24626165987542962E-4     2.3979305172488162E-4
+    955     9.16432772236257481E-4     2.62499999999999984E-4
+    956     8.40965569411703819E-4     2.71022418711743549E-4
+    957     8.53348088789513369E-4     1.96091760522196289E-4
+    958     8.7869917082398065E-4     2.66761209355871739E-4
+    959     8.47156829100608648E-4     2.33557089616969919E-4
+    960     8.84890430512885479E-4     2.29295880261098136E-4
+    961     6.16210937500000005E-4     0
+    962     6.80175781250000008E-4     6.50230427205121634E-5
+    963     6.15127063129517924E-4     1.28961775828680109E-4
+    964     6.48193359375000061E-4     3.25115213602560817E-5
+    965     6.4765142218975902E-4     9.69924092745961294E-5
+    966     6.15669000314758965E-4     6.44808879143400545E-5
+    967     8.08105468750000013E-4     0
+    968     8.72070312500000016E-4     0
+    969     8.79865520605121044E-4     9.96839116576813268E-5
+    970     8.4008789062499996E-4     0
+    971     8.75967916552560476E-4     4.98419558288406634E-5
+    972     8.43985494677560529E-4     4.98419558288406634E-5
+    973     1.72756436172360502E-4     4.70585859703467847E-4
+    974     8.56507863843441119E-5     4.3896484375000001E-4
+    975     1.71301572768688224E-4     3.7792968750000001E-4
+    976     1.29203611278352307E-4     4.54775351726733929E-4
+    977     1.28476179576516181E-4     4.0844726562500001E-4
+    978     1.72029004470524349E-4     4.24257773601733929E-4
+    979     8.79865520605121044E-4     9.96839116576813268E-5
+    980     9.35392068859210385E-4     6.524280030581833E-5
+    981     9.30681494601841297E-4     1.85796537326297451E-4
+    982     9.07628794732165661E-4     8.24633559817498284E-5
+    983     9.33036781730525841E-4     1.2551966881605789E-4
+    984     9.05273507603481225E-4     1.42740224491989375E-4
+    985     1.00000000000000002E-3     2.24999999999999994E-4
+    986     9.16432772236257481E-4     2.62499999999999984E-4
+    987     9.30681494601841297E-4     1.85796537326297451E-4
+    988     9.58216386118128805E-4     2.43749999999999989E-4
+    989     9.23557133419049443E-4     2.24148268663148717E-4
+    990     9.65340747300920659E-4     2.05398268663148722E-4
+    991     8.56507863843441119E-5     4.3896484375000001E-4
+    992     0     4.3896484375000001E-4
+    993     6.4572363977386294E-5     3.7792968750000001E-4
+    994     4.28253931921720559E-5     4.3896484375000001E-4
+    995     3.2286181988693147E-5     4.0844726562500001E-4
+    996     7.51115751808651962E-5     4.0844726562500001E-4
+    997     1.71301572768688224E-4     3.7792968750000001E-4
+    998     2.31915182548602798E-4     3.185630410561948E-4
+    999     2.31175962939123751E-4     4.23329063260392325E-4
+    1000     2.01608377658645498E-4     3.48246364278097405E-4
+    1001     2.31545572743863289E-4     3.70946052158293562E-4
+    1002     2.01238767853905988E-4     4.00629375380196168E-4
+    1003     8.40965569411703819E-4     2.71022418711743549E-4
+    1004     9.16432772236257481E-4     2.62499999999999984E-4
+    1005     9.33976498687064506E-4     3.35267605171179269E-4
+    1006     8.7869917082398065E-4     2.66761209355871739E-4
+    1007     9.2520463546166094E-4     2.98883802585589626E-4
+    1008     8.87471034049384109E-4     3.03145011941461382E-4
+    1009     5.8651429788997224E-4     3.03514350561540362E-4
+    1010     4.58849660034789083E-4     2.64192897481815869E-4
+    1011     5.1540545604394721E-4     2.10741452330906448E-4
+    1012     5.22681978962380661E-4     2.83853624021678116E-4
+    1013     4.87127558039368146E-4     2.37467174906361145E-4
+    1014     5.50959876966959725E-4     2.57127901446223392E-4
+    1015     4.58849660034789083E-4     2.64192897481815869E-4
+    1016     4.83502295899330232E-4     3.38002535036003411E-4
+    1017     4.15683819919475818E-4     3.45889508016053047E-4
+    1018     4.71175977967059657E-4     3.01097716258909667E-4
+    1019     4.49593057909403052E-4     3.41946021526028229E-4
+    1020     4.37266739977132477E-4     3.05041202748934431E-4
+    1021     6.55870398270680236E-4     4.58827279636259491E-4
+    1022     6.86868266859239985E-4     5.40099367412093744E-4
+    1023     6.27455298609753205E-4     5.60579639099527925E-4
+    1024     6.71369332564960165E-4     4.99463323524176645E-4
+    1025     6.57161782734496649E-4     5.50339503255810889E-4
+    1026     6.41662848440216721E-4     5.09703459367893681E-4
+    1027     6.55870398270680236E-4     4.58827279636259491E-4
+    1028     5.82521566988404701E-4     4.12072931624541923E-4
+    1029     6.149060889867273E-4     3.58911302456159137E-4
+    1030     6.19195982629542523E-4     4.3545010563040068E-4
+    1031     5.98713827987565946E-4     3.8549211704035053E-4
+    1032     6.35388243628703768E-4     4.08869291046209341E-4
+    1033     9.16432772236257481E-4     2.62499999999999984E-4
+    1034     8.53348088789513369E-4     1.96091760522196289E-4
+    1035     9.30681494601841297E-4     1.85796537326297451E-4
+    1036     8.84890430512885479E-4     2.29295880261098136E-4
+    1037     8.92014791695677333E-4     1.9094414892424687E-4
+    1038     9.23557133419049443E-4     2.24148268663148717E-4
+    1039     1.00000000000000002E-3     2.24999999999999994E-4
+    1040     1.00000000000000002E-3     2.99999999999999974E-4
+    1041     9.16432772236257481E-4     2.62499999999999984E-4
+    1042     1.00000000000000002E-3     2.62499999999999984E-4
+    1043     9.58216386118128805E-4     2.81249999999999979E-4
+    1044     9.58216386118128805E-4     2.43749999999999989E-4
+    1045     8.40965569411703819E-4     2.71022418711743549E-4
+    1046     8.70539573309993481E-4     3.74999999999999954E-4
+    1047     7.77749415513451876E-4     3.13113370342386432E-4
+    1048     8.55752571360848596E-4     3.23011209355871724E-4
+    1049     8.24144494411722679E-4     3.44056685171193193E-4
+    1050     8.09357492462577793E-4     2.92067894527064963E-4
+    1051     8.79865520605121044E-4     9.96839116576813268E-5
+    1052     9.30681494601841297E-4     1.85796537326297451E-4
+    1053     8.53348088789513369E-4     1.96091760522196289E-4
+    1054     9.05273507603481225E-4     1.42740224491989375E-4
+    1055     8.92014791695677333E-4     1.9094414892424687E-4
+    1056     8.66606804697317261E-4     1.47887836089938794E-4
+    1057     1.00000000000000002E-3     2.24999999999999994E-4
+    1058     9.30681494601841297E-4     1.85796537326297451E-4
+    1059     1.00000000000000002E-3     1.49999999999999987E-4
+    1060     9.65340747300920659E-4     2.05398268663148722E-4
+    1061     9.65340747300920659E-4     1.67898268663148732E-4
+    1062     1.00000000000000002E-3     1.87500000000000004E-4
+    1063     8.13728732664855108E-4     1.33343473922578902E-4
+    1064     8.79865520605121044E-4     9.96839116576813268E-5
+    1065     8.53348088789513369E-4     1.96091760522196289E-4
+    1066     8.4679712663498813E-4     1.16513692790130114E-4
+    1067     8.66606804697317261E-4     1.47887836089938794E-4
+    1068     8.33538410727184239E-4     1.64717617222387609E-4
+    1069     1.00000000000000002E-3     7.49999999999999934E-5
+    1070     9.35392068859210385E-4     6.524280030581833E-5
+    1071     1.00000000000000002E-3     0
+    1072     9.67696034429605203E-4     7.01214001529091549E-5
+    1073     9.67696034429605203E-4     3.2621400152909165E-5
+    1074     1.00000000000000002E-3     3.74999999999999967E-5
+    1075     2.39680731901657833E-4     5.14383951086584473E-4
+    1076     2.50000000000000005E-4     5.99999999999999947E-4
+    1077     1.87500000000000004E-4     5.7499999999999999E-4
+    1078     2.44840365950828906E-4     5.57191975543292264E-4
+    1079     2.18750000000000005E-4     5.87499999999999915E-4
+    1080     2.13590365950828932E-4     5.44691975543292232E-4
+    1081     7.69147152428217916E-4     7.17795642973629926E-5
+    1082     8.08105468750000013E-4     0
+    1083     8.79865520605121044E-4     9.96839116576813268E-5
+    1084     7.88626310589108964E-4     3.58897821486814963E-5
+    1085     8.43985494677560529E-4     4.98419558288406634E-5
+    1086     8.2450633651666948E-4     8.57317379775221597E-5
+    1087     1.87500000000000004E-4     5.7499999999999999E-4
+    1088     1.25000000000000003E-4     5.49999999999999925E-4
+    1089     1.72756436172360502E-4     4.70585859703467847E-4
+    1090     1.56250000000000003E-4     5.62499999999999957E-4
+    1091     1.48878218086180266E-4     5.10292929851733886E-4
+    1092     1.80128218086180239E-4     5.22792929851733919E-4
+    1093     7.14223074903462172E-4     3.09513800853735191E-4
+    1094     7.66063425991097419E-4     2.25618978565681171E-4
+    1095     7.77749415513451876E-4     3.13113370342386432E-4
+    1096     7.40143250447279849E-4     2.67566389709708195E-4
+    1097     7.71906420752274647E-4     2.69366174454033788E-4
+    1098     7.4598624520845697E-4     3.11313585598060839E-4
+    1099     2.97380987249783093E-4     1.12592184443437615E-3
+    1100     2.25989861581488047E-4     1.1034297389022796E-3
+    1101     2.63550467795857588E-4     1.0283522649845746E-3
+    1102     2.61685424415635556E-4     1.11467579166832798E-3
+    1103     2.44770164688672804E-4     1.06589100194342699E-3
+    1104     2.80465727522820341E-4     1.07713705470947538E-3
+    1105     4.14780023525604651E-4     1.04732934218306818E-3
+    1106     4.8087086624308442E-4     1.00763880523435656E-3
+    1107     4.64295327735872877E-4     1.106418874786976E-3
+    1108     4.47825444884344536E-4     1.02748407370871226E-3
+    1109     4.72583096989478675E-4     1.05702884001066628E-3
+    1110     4.39537675630738791E-4     1.07687410848502198E-3
+    1111     5.0000000000000001E-4     1.19999999999999989E-3
+    1112     3.72070312500000005E-4     1.19999999999999989E-3
+    1113     3.88700445385212667E-4     1.12415429991009354E-3
+    1114     4.36035156250000008E-4     1.19999999999999989E-3
+    1115     3.80385378942606363E-4     1.16207714995504672E-3
+    1116     4.44350222692606366E-4     1.16207714995504672E-3
+    1117     1.00000000000000002E-3     4.49999999999999988E-4
+    1118     9.33976498687064506E-4     3.35267605171179269E-4
+    1119     1.00000000000000002E-3     2.99999999999999974E-4
+    1120     9.66988249343532264E-4     3.92633802585589655E-4
+    1121     9.66988249343532264E-4     3.17633802585589621E-4
+    1122     1.00000000000000002E-3     3.75000000000000008E-4
+    1123     8.40965569411703819E-4     2.71022418711743549E-4
+    1124     9.33976498687064506E-4     3.35267605171179269E-4
+    1125     8.70539573309993481E-4     3.74999999999999954E-4
+    1126     8.87471034049384109E-4     3.03145011941461382E-4
+    1127     9.02258035998528994E-4     3.55133802585589611E-4
+    1128     8.55752571360848596E-4     3.23011209355871724E-4
+    1129     9.33471456834415085E-4     4.87499999999999978E-4
+    1130     8.70539573309993481E-4     3.74999999999999954E-4
+    1131     1.00000000000000002E-3     4.49999999999999988E-4
+    1132     9.02005515072204283E-4     4.31249999999999993E-4
+    1133     9.35269786654996751E-4     4.12499999999999998E-4
+    1134     9.66735728417207553E-4     4.68749999999999983E-4
+    1135     8.92816374246439338E-4     5.66266606957709193E-4
+    1136     9.33471456834415085E-4     4.87499999999999978E-4
+    1137     1.00000000000000002E-3     5.24999999999999968E-4
+    1138     9.13143915540427266E-4     5.26883303478854585E-4
+    1139     9.66735728417207553E-4     5.06249999999999973E-4
+    1140     9.46408187123219625E-4     5.4563330347885458E-4
+    1141     8.70539573309993481E-4     3.74999999999999954E-4
+    1142     9.33471456834415085E-4     4.87499999999999978E-4
+    1143     8.58425917146639784E-4     4.73341467989030294E-4
+    1144     9.02005515072204283E-4     4.31249999999999993E-4
+    1145     8.95948686990527435E-4     4.80420733994515136E-4
+    1146     8.64482745228316632E-4     4.24170733994515096E-4
+    1147     4.10887142695853748E-4     8.117939619078767E-4
+    1148     3.00868592099089952E-4     7.96027940566129363E-4
+    1149     3.3342684329849511E-4     7.35743010166048958E-4
+    1150     3.55877867397471877E-4     8.03910951237003031E-4
+    1151     3.17147717698792531E-4     7.6588547536608916E-4
+    1152     3.72156992997174402E-4     7.73768486036962829E-4
+    1153     2.29766755665683817E-4     7.70779974744912363E-4
+    1154     3.00868592099089952E-4     7.96027940566129363E-4
+    1155     2.56894004484586392E-4     9.02486475935557219E-4
+    1156     2.65317673882386898E-4     7.83403957655520863E-4
+    1157     2.78881298291838199E-4     8.49257208250843291E-4
+    1158     2.43330380075135091E-4     8.36633225340234791E-4
+    1159     2.56894004484586392E-4     9.02486475935557219E-4
+    1160     3.00868592099089952E-4     7.96027940566129363E-4
+    1161     3.14063844356556141E-4     8.63789983373621688E-4
+    1162     2.78881298291838199E-4     8.49257208250843291E-4
+    1163     3.07466218227823019E-4     8.2990896196987558E-4
+    1164     2.85478924420571266E-4     8.83138229654589508E-4
+    1165     1.36565743695333484E-4     7.39812343895404297E-4
+    1166     1.25000000000000003E-4     6.4999999999999997E-4
+    1167     2.16128514871350092E-4     6.96571287178375043E-4
+    1168     1.3078287184766673E-4     6.9490617194770208E-4
+    1169     1.70564257435675047E-4     6.73285643589187507E-4
+    1170     1.76347129283341774E-4     7.18191815536889724E-4
+    1171     2.16128514871350092E-4     6.96571287178375043E-4
+    1172     3.14898739406715305E-4     6.69780733668936504E-4
+    1173     3.3342684329849511E-4     7.35743010166048958E-4
+    1174     2.65513627139032698E-4     6.83176010423655828E-4
+    1175     3.24162791352605208E-4     7.02761871917492731E-4
+    1176     2.74777679084922601E-4     7.16157148672211946E-4
+    1177     3.14898739406715305E-4     6.69780733668936504E-4
+    1178     2.16128514871350092E-4     6.96571287178375043E-4
+    1179     2.50000000000000005E-4     5.99999999999999947E-4
+    1180     2.65513627139032698E-4     6.83176010423655828E-4
+    1181     2.33064257435675048E-4     6.48285643589187549E-4
+    1182     2.82449369703357655E-4     6.34890366834468226E-4
+    1183     2.16128514871350092E-4     6.96571287178375043E-4
+    1184     2.29766755665683817E-4     7.70779974744912363E-4
+    1185     1.36565743695333484E-4     7.39812343895404297E-4
+    1186     2.22947635268516968E-4     7.33675630961643649E-4
+    1187     1.8316624968050865E-4     7.5529615932015833E-4
+    1188     1.76347129283341774E-4     7.18191815536889724E-4
+    1189     2.63550467795857588E-4     1.0283522649845746E-3
+    1190     3.33510000366525681E-4     9.78187713756657664E-4
+    1191     3.39525065135284183E-4     1.0640634227390719E-3
+    1192     2.98530234081191635E-4     1.00326998937061602E-3
+    1193     3.36517532750904959E-4     1.02112556824786489E-3
+    1194     3.01537766465570858E-4     1.04620784386182325E-3
+    1195     2.29999438806578629E-4     9.67017715356922923E-4
+    1196     1.1927993626285668E-4     9.2101752411100806E-4
+    1197     1.85578989176374238E-4     8.9312739269816724E-4
+    1198     1.74639687534717668E-4     9.44017619733965438E-4
+    1199     1.52429462719615459E-4     9.07072458404587596E-4
+    1200     2.07789213991476447E-4     9.30072554027545082E-4
+    1201     1.74429316483144684E-4     8.22070312499999993E-4
+    1202     1.1927993626285668E-4     9.2101752411100806E-4
+    1203     6.75304712410491789E-5     8.77236181658321186E-4
+    1204     1.46854626373000695E-4     8.71543918305504027E-4
+    1205     9.34052037519529361E-5     8.99126852884664677E-4
+    1206     1.20979893862096924E-4     8.49653247079160644E-4
+    1207     1.25000000000000003E-4     6.4999999999999997E-4
+    1208     1.36565743695333484E-4     7.39812343895404297E-4
+    1209     6.25000000000000013E-5     6.75000000000000036E-4
+    1210     1.3078287184766673E-4     6.9490617194770208E-4
+    1211     9.95328718476667424E-5     7.07406171947702221E-4
+    1212     9.3750000000000002E-5     6.62500000000000003E-4
+    1213     1.85578989176374238E-4     8.9312739269816724E-4
+    1214     1.1927993626285668E-4     9.2101752411100806E-4
+    1215     1.74429316483144684E-4     8.22070312499999993E-4
+    1216     1.52429462719615459E-4     9.07072458404587596E-4
+    1217     1.46854626373000695E-4     8.71543918305504027E-4
+    1218     1.80004152829759474E-4     8.57598852599083562E-4
+    1219     3.79565818807853228E-4     5.99785262317693206E-4
+    1220     4.62057889734597506E-4     6.86398137230062556E-4
+    1221     3.98747049720822426E-4     7.15064925451776041E-4
+    1222     4.20811854271225367E-4     6.43091699773877881E-4
+    1223     4.30402469727709939E-4     7.00731531340919298E-4
+    1224     3.891564342643378E-4     6.57425093884734624E-4
+    1225     8.72146582415723147E-5     8.12372124900119985E-4
+    1226     0     7.61035156249999993E-4
+    1227     6.7318724869494569E-5     7.44219603001817923E-4
+    1228     4.36073291207861573E-5     7.86703640575059935E-4
+    1229     3.36593624347472845E-5     7.52627379625908904E-4
+    1230     7.72666915555334486E-5     7.78295863950968954E-4
+    1231     2.25989861581488047E-4     1.1034297389022796E-3
+    1232     1.8310546875E-4     1.19999999999999989E-3
+    1233     1.27766480316576791E-4     1.1007034017909631E-3
+    1234     2.0454766516574401E-4     1.15171486945113985E-3
+    1235     1.55435974533288396E-4     1.15035170089548139E-3
+    1236     1.76878170949032432E-4     1.10206657034662135E-3
+    1237     0     8.22070312499999993E-4
+    1238     8.72146582415723147E-5     8.12372124900119985E-4
+    1239     6.75304712410491789E-5     8.77236181658321186E-4
+    1240     4.36073291207861573E-5     8.17221218700059935E-4
+    1241     7.737256474131074E-5     8.44804153279220585E-4
+    1242     3.37652356205245894E-5     8.49653247079160644E-4
+    1243     1.27766480316576791E-4     1.1007034017909631E-3
+    1244     6.77802435850866764E-5     1.13246159963255083E-3
+    1245     8.06327547355035015E-5     1.013121025387205E-3
+    1246     9.77733619508317338E-5     1.11658250071175686E-3
+    1247     7.42064991602950889E-5     1.07279131250987781E-3
+    1248     1.04199617526040146E-4     1.05691221358908405E-3
+    1249     1.1927993626285668E-4     9.2101752411100806E-4
+    1250     8.06327547355035015E-5     1.013121025387205E-3
+    1251     6.33482805365148801E-5     9.51708311133460578E-4
+    1252     9.99563454991800906E-5     9.67069274749106584E-4
+    1253     7.19905176360091908E-5     9.82414668260332734E-4
+    1254     9.131410839968578E-5     9.36362917622234319E-4
+    1255     0     1.19999999999999989E-3
+    1256     0     1.13603515625E-3
+    1257     6.103515625E-5     1.19999999999999989E-3
+    1258     0     1.16801757812500006E-3
+    1259     3.0517578125E-5     1.16801757812500006E-3
+    1260     3.0517578125E-5     1.19999999999999989E-3
+    1261     8.06327547355035015E-5     1.013121025387205E-3
+    1262     0     1.07207031249999989E-3
+    1263     0     1.00810546875E-3
+    1264     4.03163773677517507E-5     1.04259566894360244E-3
+    1265     0     1.04008789062499983E-3
+    1266     4.03163773677517507E-5     1.01061324706860261E-3
+    1267     1.27766480316576791E-4     1.1007034017909631E-3
+    1268     1.80463733243880608E-4     1.01635139957644981E-3
+    1269     2.25989861581488047E-4     1.1034297389022796E-3
+    1270     1.541151067802287E-4     1.05852740068370635E-3
+    1271     2.03226797412684314E-4     1.05989056923936481E-3
+    1272     1.76878170949032432E-4     1.10206657034662135E-3
+    1273     6.77802435850866764E-5     1.13246159963255083E-3
+    1274     1.220703125E-4     1.19999999999999989E-3
+    1275     6.103515625E-5     1.19999999999999989E-3
+    1276     9.49252780425433382E-5     1.16623079981627536E-3
+    1277     9.1552734375E-5     1.19999999999999989E-3
+    1278     6.44076999175433382E-5     1.16623079981627536E-3
+    1279     2.63550467795857588E-4     1.0283522649845746E-3
+    1280     2.25989861581488047E-4     1.1034297389022796E-3
+    1281     1.80463733243880608E-4     1.01635139957644981E-3
+    1282     2.44770164688672804E-4     1.06589100194342699E-3
+    1283     2.03226797412684314E-4     1.05989056923936481E-3
+    1284     2.22007100519869098E-4     1.02235183228051221E-3
+    1285     3.39525065135284183E-4     1.0640634227390719E-3
+    1286     2.97380987249783093E-4     1.12592184443437615E-3
+    1287     2.63550467795857588E-4     1.0283522649845746E-3
+    1288     3.18453026192533611E-4     1.09499263358672403E-3
+    1289     2.80465727522820341E-4     1.07713705470947538E-3
+    1290     3.01537766465570858E-4     1.04620784386182325E-3
+    1291     1.80463733243880608E-4     1.01635139957644981E-3
+    1292     2.29999438806578629E-4     9.67017715356922923E-4
+    1293     2.63550467795857588E-4     1.0283522649845746E-3
+    1294     2.05231586025229619E-4     9.91684557466686366E-4
+    1295     2.46774953301218082E-4     9.97684990170748763E-4
+    1296     2.22007100519869098E-4     1.02235183228051221E-3
+    1297     2.56894004484586392E-4     9.02486475935557219E-4
+    1298     3.14063844356556141E-4     8.63789983373621688E-4
+    1299     3.33510000366525681E-4     9.78187713756657664E-4
+    1300     2.85478924420571266E-4     8.83138229654589508E-4
+    1301     3.23786922361540938E-4     9.20988848565139676E-4
+    1302     2.9520200242555601E-4     9.40337094846107387E-4
+    1303     3.33510000366525681E-4     9.78187713756657664E-4
+    1304     2.63550467795857588E-4     1.0283522649845746E-3
+    1305     2.29999438806578629E-4     9.67017715356922923E-4
+    1306     2.98530234081191635E-4     1.00326998937061602E-3
+    1307     2.46774953301218082E-4     9.97684990170748763E-4
+    1308     2.81754719586552183E-4     9.72602714556790294E-4
+    1309     4.11905532296878403E-4     9.69320888710480284E-4
+    1310     3.33510000366525681E-4     9.78187713756657664E-4
+    1311     3.72070312500000005E-4     9.01220703124999949E-4
+    1312     3.72707766331702042E-4     9.73754301233568974E-4
+    1313     3.52790156433262843E-4     9.39704208440828807E-4
+    1314     3.91987922398439204E-4     9.35270795917740117E-4
+    1315     4.51818848350113347E-4     9.0027268151548134E-4
+    1316     4.11905532296878403E-4     9.69320888710480284E-4
+    1317     3.72070312500000005E-4     9.01220703124999949E-4
+    1318     4.31862190323495875E-4     9.34796785112980812E-4
+    1319     3.91987922398439204E-4     9.35270795917740117E-4
+    1320     4.11944580425056676E-4     9.00746692320240645E-4
+    1321     4.11905532296878403E-4     9.69320888710480284E-4
+    1322     4.8087086624308442E-4     1.00763880523435656E-3
+    1323     4.14780023525604651E-4     1.04732934218306818E-3
+    1324     4.46388199269981412E-4     9.88479846972418422E-4
+    1325     4.47825444884344536E-4     1.02748407370871226E-3
+    1326     4.13342777911241527E-4     1.00832511544677412E-3
+    1327     4.10887142695853748E-4     8.117939619078767E-4
+    1328     4.51818848350113347E-4     9.0027268151548134E-4
+    1329     3.72070312500000005E-4     9.01220703124999949E-4
+    1330     4.31352995522983574E-4     8.56033321711678966E-4
+    1331     4.11944580425056676E-4     9.00746692320240645E-4
+    1332     3.91478727597926904E-4     8.5650733251643827E-4
+    1333     4.51818848350113347E-4     9.0027268151548134E-4
+    1334     5.19966569338154057E-4     8.58840618576917807E-4
+    1335     5.42759630809160629E-4     9.44364496429779174E-4
+    1336     4.85892708844133702E-4     8.79556650046199573E-4
+    1337     5.31363100073657343E-4     9.01602557503348436E-4
+    1338     4.97289239579636988E-4     9.22318588972630311E-4
+    1339     4.10887142695853748E-4     8.117939619078767E-4
+    1340     3.72070312500000005E-4     9.01220703124999949E-4
+    1341     3.14063844356556141E-4     8.63789983373621688E-4
+    1342     3.91478727597926904E-4     8.5650733251643827E-4
+    1343     3.43067078428278046E-4     8.82505343249310819E-4
+    1344     3.62475493526204944E-4     8.37791972640749248E-4
+    1345     4.10887142695853748E-4     8.117939619078767E-4
+    1346     3.98747049720822426E-4     7.15064925451776041E-4
+    1347     4.75850298603600714E-4     7.54514336975707829E-4
+    1348     4.04817096208338087E-4     7.63429443679826316E-4
+    1349     4.37298674162211543E-4     7.3478963121374188E-4
+    1350     4.43368720649727258E-4     7.83154149441792264E-4
+    1351     4.62057889734597506E-4     6.86398137230062556E-4
+    1352     4.75850298603600714E-4     7.54514336975707829E-4
+    1353     3.98747049720822426E-4     7.15064925451776041E-4
+    1354     4.6895409416909911E-4     7.20456237102885247E-4
+    1355     4.37298674162211543E-4     7.3478963121374188E-4
+    1356     4.30402469727709939E-4     7.00731531340919298E-4
+    1357     4.75850298603600714E-4     7.54514336975707829E-4
+    1358     5.60394849123839033E-4     7.35715772264671007E-4
+    1359     5.79668274240937665E-4     8.10244013239797283E-4
+    1360     5.18122573863719874E-4     7.45115054620189418E-4
+    1361     5.70031561682388349E-4     7.72979892752234145E-4
+    1362     5.27759286422269189E-4     7.82379175107752556E-4
+    1363     0     7.0E-4
+    1364     6.25000000000000013E-5     6.75000000000000036E-4
+    1365     6.7318724869494569E-5     7.44219603001817923E-4
+    1366     3.12500000000000007E-5     6.8749999999999996E-4
+    1367     6.49093624347472919E-5     7.09609801500908979E-4
+    1368     3.36593624347472845E-5     7.22109801500908904E-4
+    1369     2.25989861581488047E-4     1.1034297389022796E-3
+    1370     2.97380987249783093E-4     1.12592184443437615E-3
+    1371     2.44140625E-4     1.19999999999999989E-3
+    1372     2.61685424415635556E-4     1.11467579166832798E-3
+    1373     2.70760806124891546E-4     1.16296092221718802E-3
+    1374     2.3506524329074401E-4     1.15171486945113985E-3
+    1375     1.1927993626285668E-4     9.2101752411100806E-4
+    1376     2.29999438806578629E-4     9.67017715356922923E-4
+    1377     1.80463733243880608E-4     1.01635139957644981E-3
+    1378     1.74639687534717668E-4     9.44017619733965438E-4
+    1379     2.05231586025229619E-4     9.91684557466686366E-4
+    1380     1.4987183475336863E-4     9.6868446184372888E-4
+    1381     3.33510000366525681E-4     9.78187713756657664E-4
+    1382     2.29999438806578629E-4     9.67017715356922923E-4
+    1383     2.56894004484586392E-4     9.02486475935557219E-4
+    1384     2.81754719586552183E-4     9.72602714556790294E-4
+    1385     2.43446721645582511E-4     9.34752095646240126E-4
+    1386     2.9520200242555601E-4     9.40337094846107387E-4
+    1387     5.46592219681165661E-4     1.47080610942937849E-4
+    1388     6.03238036253151279E-4     2.17057585313148337E-4
+    1389     5.1540545604394721E-4     2.10741452330906448E-4
+    1390     5.74915127967158416E-4     1.82069098128043107E-4
+    1391     5.5932174614854919E-4     2.13899518822027393E-4
+    1392     5.30998837862556436E-4     1.78911031636922162E-4
+    1393     6.73095146243751124E-4     1.78431349377634624E-4
+    1394     6.03238036253151279E-4     2.17057585313148337E-4
+    1395     6.15127063129517924E-4     1.28961775828680109E-4
+    1396     6.38166591248451147E-4     1.97744467345391467E-4
+    1397     6.09182549691334547E-4     1.73009680570914223E-4
+    1398     6.44111104686634524E-4     1.53696562603157353E-4
+    1399     3.15575840731535889E-4     4.21078468404699436E-4
+    1400     2.31175962939123751E-4     4.23329063260392325E-4
+    1401     2.80075730378807853E-4     3.66277694634743636E-4
+    1402     2.73375901835329793E-4     4.22203765832545881E-4
+    1403     2.55625846658965802E-4     3.94803378947567981E-4
+    1404     2.97825785555171844E-4     3.93678081519721536E-4
+    1405     2.39680731901657833E-4     5.14383951086584473E-4
+    1406     2.31175962939123751E-4     4.23329063260392325E-4
+    1407     3.14594753959753046E-4     4.86365680855438756E-4
+    1408     2.35428347420390779E-4     4.68856507173488399E-4
+    1409     2.72885358449438372E-4     4.54847372057915541E-4
+    1410     2.77137742930705453E-4     5.00374815971011615E-4
+    1411     6.77802435850866764E-5     1.13246159963255083E-3
+    1412     6.103515625E-5     1.19999999999999989E-3
+    1413     0     1.13603515625E-3
+    1414     6.44076999175433382E-5     1.16623079981627536E-3
+    1415     3.0517578125E-5     1.16801757812500006E-3
+    1416     3.38901217925433382E-5     1.13424837794127552E-3
+    1417     6.77802435850866764E-5     1.13246159963255083E-3
+    1418     0     1.07207031249999989E-3
+    1419     8.06327547355035015E-5     1.013121025387205E-3
+    1420     3.38901217925433382E-5     1.10226595606627547E-3
+    1421     4.03163773677517507E-5     1.04259566894360244E-3
+    1422     7.42064991602950889E-5     1.07279131250987781E-3
+    1423     6.7318724869494569E-5     7.44219603001817923E-4
+    1424     6.25000000000000013E-5     6.75000000000000036E-4
+    1425     1.36565743695333484E-4     7.39812343895404297E-4
+    1426     6.49093624347472919E-5     7.09609801500908979E-4
+    1427     9.95328718476667424E-5     7.07406171947702221E-4
+    1428     1.0194223428241402E-4     7.42015973448611164E-4
+    1429     8.72146582415723147E-5     8.12372124900119985E-4
+    1430     0     8.22070312499999993E-4
+    1431     0     7.61035156249999993E-4
+    1432     4.36073291207861573E-5     8.17221218700059935E-4
+    1433     0     7.91552734374999993E-4
+    1434     4.36073291207861573E-5     7.86703640575059935E-4
+    1435     5.60394849123839033E-4     7.35715772264671007E-4
+    1436     4.62057889734597506E-4     6.86398137230062556E-4
+    1437     4.9377437948036596E-4     6.24558736453494859E-4
+    1438     5.11226369429218269E-4     7.11056954747366727E-4
+    1439     4.77916134607481733E-4     6.55478436841778708E-4
+    1440     5.27084614302102497E-4     6.80137254359082987E-4
+    1441     4.51818848350113347E-4     9.0027268151548134E-4
+    1442     5.42759630809160629E-4     9.44364496429779174E-4
+    1443     4.8087086624308442E-4     1.00763880523435656E-3
+    1444     4.97289239579636988E-4     9.22318588972630311E-4
+    1445     5.11815248526122524E-4     9.76001650832067922E-4
+    1446     4.66344857296598883E-4     9.5395574337491895E-4
+    1447     6.28726770153248599E-4     8.69566754220071274E-4
+    1448     5.42759630809160629E-4     9.44364496429779174E-4
+    1449     5.19966569338154057E-4     8.58840618576917807E-4
+    1450     5.85743200481204668E-4     9.06965625324925224E-4
+    1451     5.31363100073657343E-4     9.01602557503348436E-4
+    1452     5.74346669745701274E-4     8.64203686398494486E-4
+    1453     6.3613763360561681E-4     1.03946949707743879E-3
+    1454     5.67809578584186053E-4     1.01842229960029352E-3
+    1455     6.27929687500000016E-4     9.68445999237773649E-4
+    1456     6.01973606094901486E-4     1.02894589833886626E-3
+    1457     5.97869633042093034E-4     9.9343414941903364E-4
+    1458     6.32033660552808359E-4     1.00395774815760616E-3
+    1459     4.8087086624308442E-4     1.00763880523435656E-3
+    1460     5.67809578584186053E-4     1.01842229960029352E-3
+    1461     5.63964843750000013E-4     1.0965073564807807E-3
+    1462     5.24340222413635236E-4     1.01303055241732504E-3
+    1463     5.65887211167093087E-4     1.05746482804053711E-3
+    1464     5.22417854996542162E-4     1.05207308085756863E-3
+    1465     7.29985750852659397E-4     9.92266244781275222E-4
+    1466     6.3613763360561681E-4     1.03946949707743879E-3
+    1467     6.27929687500000016E-4     9.68445999237773649E-4
+    1468     6.83061692229138158E-4     1.015867870929357E-3
+    1469     6.32033660552808359E-4     1.00395774815760616E-3
+    1470     6.78957719176329706E-4     9.80356122009524381E-4
+    1471     6.3613763360561681E-4     1.03946949707743879E-3
+    1472     6.91894531250000018E-4     1.08422299961888672E-3
+    1473     6.31385186434097476E-4     1.12635085509955684E-3
+    1474     6.64016082427808414E-4     1.06184624834816275E-3
+    1475     6.61639858842048801E-4     1.10528692735922178E-3
+    1476     6.33761410019857089E-4     1.08291017608849781E-3
+    1477     6.27929687500000016E-4     1.19999999999999989E-3
+    1478     5.63964843750000013E-4     1.19999999999999989E-3
+    1479     6.31385186434097476E-4     1.12635085509955684E-3
+    1480     5.9594726562499996E-4     1.19999999999999989E-3
+    1481     5.9767501509204869E-4     1.16317542754977837E-3
+    1482     6.29657436967048746E-4     1.16317542754977837E-3
+    1483     6.91894531250000018E-4     1.19999999999999989E-3
+    1484     6.27929687500000016E-4     1.19999999999999989E-3
+    1485     6.31385186434097476E-4     1.12635085509955684E-3
+    1486     6.59912109375000071E-4     1.19999999999999989E-3
+    1487     6.29657436967048746E-4     1.16317542754977837E-3
+    1488     6.61639858842048801E-4     1.16317542754977837E-3
+    1489     8.33723092676693115E-4     1.08142335682909598E-3
+    1490     8.77929687500000021E-4     1.19999999999999989E-3
+    1491     8.16894531250000021E-4     1.19999999999999989E-3
+    1492     8.55826390088346568E-4     1.14071167841454804E-3
+    1493     8.47412109375000021E-4     1.19999999999999989E-3
+    1494     8.25308811963346568E-4     1.14071167841454804E-3
+    1495     6.91894531250000018E-4     1.08422299961888672E-3
+    1496     7.70557602221913452E-4     1.04998903577999537E-3
+    1497     7.63546512545414027E-4     1.12019473601693395E-3
+    1498     7.31226066735956789E-4     1.06710601769944105E-3
+    1499     7.67052057383663739E-4     1.08509188589846466E-3
+    1500     7.27720521897706968E-4     1.10220886781791034E-3
+    1501     8.42471384914615607E-4     1.00872691912019851E-3
+    1502     7.29985750852659397E-4     9.92266244781275222E-4
+    1503     8.1517205990673709E-4     9.40785388914013596E-4
+    1504     7.86228567883637502E-4     1.00049658195073676E-3
+    1505     7.72578905379698189E-4     9.66525816847644409E-4
+    1506     8.28821722410676403E-4     9.74756154017106051E-4
+    1507     8.77929687500000021E-4     1.19999999999999989E-3
+    1508     8.33723092676693115E-4     1.08142335682909598E-3
+    1509     9.22298701458379552E-4     1.09448242187499991E-3
+    1510     8.55826390088346568E-4     1.14071167841454804E-3
+    1511     8.78010897067536388E-4     1.08795288935204784E-3
+    1512     9.00114194479189732E-4     1.1472412109374999E-3
+    1513     6.31385186434097476E-4     1.12635085509955684E-3
+    1514     5.63964843750000013E-4     1.19999999999999989E-3
+    1515     5.63964843750000013E-4     1.0965073564807807E-3
+    1516     5.9767501509204869E-4     1.16317542754977837E-3
+    1517     5.63964843750000013E-4     1.14825367824039019E-3
+    1518     5.9767501509204869E-4     1.11142910579016888E-3
+    1519     9.22298701458379552E-4     1.09448242187499991E-3
+    1520     9.38964843750000021E-4     1.19999999999999989E-3
+    1521     8.77929687500000021E-4     1.19999999999999989E-3
+    1522     9.30631772604189732E-4     1.1472412109374999E-3
+    1523     9.08447265625000021E-4     1.19999999999999989E-3
+    1524     9.00114194479189732E-4     1.1472412109374999E-3
+    1525     9.22298701458379552E-4     1.09448242187499991E-3
+    1526     1.00000000000000002E-3     1.13896484375000011E-3
+    1527     9.38964843750000021E-4     1.19999999999999989E-3
+    1528     9.61149350729189732E-4     1.1167236328124999E-3
+    1529     9.69482421875000021E-4     1.16948242187500011E-3
+    1530     9.30631772604189732E-4     1.1472412109374999E-3
+    1531     7.70557602221913452E-4     1.04998903577999537E-3
+    1532     8.42471384914615607E-4     1.00872691912019851E-3
+    1533     8.33723092676693115E-4     1.08142335682909598E-3
+    1534     8.0651449356826453E-4     1.02935797745009683E-3
+    1535     8.38097238795654361E-4     1.04507513797464724E-3
+    1536     8.02140347449303284E-4     1.06570619630454557E-3
+    1537     1.00000000000000002E-3     9.0E-4
+    1538     9.12446330413049594E-4     9.87167427587815364E-4
+    1539     9.00310858258578334E-4     8.8742365740778872E-4
+    1540     9.56223165206524861E-4     9.4358371379390767E-4
+    1541     9.06378594335813964E-4     9.37295542497801988E-4
+    1542     9.50155429129289231E-4     8.93711828703894293E-4
+    1543     1.00000000000000002E-3     9.74999999999999955E-4
+    1544     1.00000000000000002E-3     1.04999999999999994E-3
+    1545     9.12446330413049594E-4     9.87167427587815364E-4
+    1546     1.00000000000000002E-3     1.01249999999999995E-3
+    1547     9.56223165206524861E-4     1.01858371379390765E-3
+    1548     9.56223165206524861E-4     9.8108371379390766E-4
+    1549     8.9886962494534965E-4     7.81508352249731973E-4
+    1550     8.37129054430499212E-4     8.81217371382065837E-4
+    1551     8.03413482021300262E-4     8.27424077205904536E-4
+    1552     8.67999339687924431E-4     8.31362861815898905E-4
+    1553     8.20271268225899737E-4     8.54320724293985241E-4
+    1554     8.51141553483324956E-4     8.044662147278182E-4
+    1555     9.36226394521655E-4     8.3396748156078118E-4
+    1556     9.39685899131729899E-4     7.34463702578968722E-4
+    1557     1.00000000000000002E-3     7.50000000000000016E-4
+    1558     9.37956146826692504E-4     7.84215592069874897E-4
+    1559     9.69842949565865014E-4     7.42231851289484315E-4
+    1560     9.68113197260827511E-4     7.91983740780390598E-4
+    1561     8.35489473019320351E-4     7.27605773008054794E-4
+    1562     9.39685899131729899E-4     7.34463702578968722E-4
+    1563     8.9886962494534965E-4     7.81508352249731973E-4
+    1564     8.87587686075525179E-4     7.31034737793511812E-4
+    1565     9.1927776203853972E-4     7.57986027414350293E-4
+    1566     8.67179548982335E-4     7.54557062628893383E-4
+    1567     9.39685899131729899E-4     7.34463702578968722E-4
+    1568     1.00000000000000002E-3     6.75000000000000036E-4
+    1569     1.00000000000000002E-3     7.50000000000000016E-4
+    1570     9.69842949565865014E-4     7.04731851289484325E-4
+    1571     1.00000000000000002E-3     7.12500000000000026E-4
+    1572     9.69842949565865014E-4     7.42231851289484315E-4
+    1573     1.00000000000000002E-3     6.75000000000000036E-4
+    1574     9.39685899131729899E-4     7.34463702578968722E-4
+    1575     9.21159286470410715E-4     6.74999999999999927E-4
+    1576     9.69842949565865014E-4     7.04731851289484325E-4
+    1577     9.30422592801070361E-4     7.04731851289484325E-4
+    1578     9.60579643235205368E-4     6.75000000000000036E-4
+    1579     8.14703805291646039E-4     6.47042209754021894E-4
+    1580     9.37074205321842733E-4     6.1279086004113368E-4
+    1581     9.21159286470410715E-4     6.74999999999999927E-4
+    1582     8.75889005306744331E-4     6.29916534897577733E-4
+    1583     9.29116745896126778E-4     6.43895430020566804E-4
+    1584     8.67931545881028377E-4     6.61021104877010856E-4
+    1585     1.00000000000000002E-3     5.24999999999999968E-4
+    1586     9.37074205321842733E-4     6.1279086004113368E-4
+    1587     8.92816374246439338E-4     5.66266606957709193E-4
+    1588     9.68537102660921431E-4     5.68895430020566824E-4
+    1589     9.14945289784141035E-4     5.89528733499421437E-4
+    1590     9.46408187123219625E-4     5.4563330347885458E-4
+    1591     5.59110560621724098E-4     5.93289128784703487E-4
+    1592     6.27455298609753205E-4     5.60579639099527925E-4
+    1593     6.25000000000000013E-4     6.23375449734813098E-4
+    1594     5.93282929615738652E-4     5.76934383942115706E-4
+    1595     6.26227649304876609E-4     5.91977544417170458E-4
+    1596     5.92055280310862055E-4     6.08332289259758347E-4
+    1597     6.55870398270680236E-4     4.58827279636259491E-4
+    1598     6.27455298609753205E-4     5.60579639099527925E-4
+    1599     5.77009580214111854E-4     5.23102021510448161E-4
+    1600     6.41662848440216721E-4     5.09703459367893681E-4
+    1601     6.02232439411932476E-4     5.41840830304987989E-4
+    1602     6.16439989242395991E-4     4.90964650573353853E-4
+    1603     5.77009580214111854E-4     5.23102021510448161E-4
+    1604     4.44306653512529842E-4     5.02492175855262643E-4
+    1605     5.18154572622838414E-4     4.64528827497301226E-4
+    1606     5.10658116863320875E-4     5.12797098682855456E-4
+    1607     4.81230613067684101E-4     4.83510501676281935E-4
+    1608     5.4758207641847508E-4     4.93815424503874639E-4
+    1609     4.44306653512529842E-4     5.02492175855262643E-4
+    1610     3.78329331029246172E-4     4.03040117695809171E-4
+    1611     4.57065334871855982E-4     4.13720111675068352E-4
+    1612     4.1131799227088798E-4     4.52766146775535907E-4
+    1613     4.17697332950551104E-4     4.08380114685438734E-4
+    1614     4.50685994192192912E-4     4.5810614376516547E-4
+    1615     5.60394849123839033E-4     7.35715772264671007E-4
+    1616     5.59110560621724098E-4     5.93289128784703487E-4
+    1617     6.25000000000000013E-4     6.23375449734813098E-4
+    1618     5.5975270487278162E-4     6.64502450524687193E-4
+    1619     5.92055280310862055E-4     6.08332289259758347E-4
+    1620     5.92697424561919577E-4     6.79545610999742053E-4
+    1621     5.59110560621724098E-4     5.93289128784703487E-4
+    1622     4.9377437948036596E-4     6.24558736453494859E-4
+    1623     4.45119450916024088E-4     5.73219517420557926E-4
+    1624     5.26442470051044975E-4     6.08923932619099173E-4
+    1625     4.69446915198194997E-4     5.98889126937026447E-4
+    1626     5.0211500576887412E-4     5.83254323102630652E-4
+    1627     8.58425917146639784E-4     4.73341467989030294E-4
+    1628     7.42631330235908388E-4     4.52617545905931022E-4
+    1629     8.08703962690366035E-4     4.17299929575055159E-4
+    1630     8.00528623691274086E-4     4.62979506947480658E-4
+    1631     7.75667646463137211E-4     4.3495873774049309E-4
+    1632     8.3356493991850291E-4     4.45320698782042726E-4
+    1633     7.42631330235908388E-4     4.52617545905931022E-4
+    1634     6.77112852597775263E-4     3.61199236444649594E-4
+    1635     7.47347240742430642E-4     3.78837751313385344E-4
+    1636     7.09872091416841825E-4     4.06908391175290281E-4
+    1637     7.12230046670103006E-4     3.70018493879017496E-4
+    1638     7.44989285489169569E-4     4.15727648609658183E-4
+    1639     8.35489473019320351E-4     7.27605773008054794E-4
+    1640     7.56891934562431879E-4     7.83576045089881936E-4
+    1641     7.52655853090133443E-4     7.19787645578877193E-4
+    1642     7.9619070379087606E-4     7.55590909048968365E-4
+    1643     7.54773893826282661E-4     7.51681845334379564E-4
+    1644     7.94072663054726843E-4     7.23696709293465993E-4
+    1645     7.24825013803433236E-4     8.92866282184746245E-4
+    1646     7.56891934562431879E-4     7.83576045089881936E-4
+    1647     8.03413482021300262E-4     8.27424077205904536E-4
+    1648     7.40858474182932557E-4     8.3822116363731409E-4
+    1649     7.80152708291866016E-4     8.05500061147893182E-4
+    1650     7.64119247912366803E-4     8.60145179695325337E-4
+    1651     7.29985750852659397E-4     9.92266244781275222E-4
+    1652     6.27929687500000016E-4     9.68445999237773649E-4
+    1653     6.8729190829219008E-4     9.44648911063777994E-4
+    1654     6.78957719176329706E-4     9.80356122009524381E-4
+    1655     6.57610797896094994E-4     9.56547455150775876E-4
+    1656     7.08638829572424793E-4     9.68457577922526608E-4
+    1657     6.28726770153248599E-4     8.69566754220071274E-4
+    1658     6.8729190829219008E-4     9.44648911063777994E-4
+    1659     6.27929687500000016E-4     9.68445999237773649E-4
+    1660     6.5800933922271934E-4     9.07107832641924688E-4
+    1661     6.57610797896094994E-4     9.56547455150775876E-4
+    1662     6.28328228826624362E-4     9.19006376728922461E-4
+    1663     4.75850298603600714E-4     7.54514336975707829E-4
+    1664     5.79668274240937665E-4     8.10244013239797283E-4
+    1665     5.19966569338154057E-4     8.58840618576917807E-4
+    1666     5.27759286422269189E-4     7.82379175107752556E-4
+    1667     5.49817421789545861E-4     8.34542315908357599E-4
+    1668     4.97908433970877386E-4     8.06677477776312872E-4
+    1669     5.60394849123839033E-4     7.35715772264671007E-4
+    1670     6.60634466318554738E-4     7.08911474490989688E-4
+    1671     6.3251725796544824E-4     7.69458270396488286E-4
+    1672     6.10514657721196886E-4     7.22313623377830293E-4
+    1673     6.46575862142001489E-4     7.39184872443738933E-4
+    1674     5.96456053544643637E-4     7.52587021330579646E-4
+    1675     4.45119450916024088E-4     5.73219517420557926E-4
+    1676     4.44306653512529842E-4     5.02492175855262643E-4
+    1677     5.77009580214111854E-4     5.23102021510448161E-4
+    1678     4.44713052214276965E-4     5.3785584663791023E-4
+    1679     5.10658116863320875E-4     5.12797098682855456E-4
+    1680     5.11064515565067944E-4     5.48160769465503043E-4
+    1681     4.62057889734597506E-4     6.86398137230062556E-4
+    1682     4.45119450916024088E-4     5.73219517420557926E-4
+    1683     4.9377437948036596E-4     6.24558736453494859E-4
+    1684     4.5358867032531077E-4     6.29808827325310187E-4
+    1685     4.69446915198194997E-4     5.98889126937026447E-4
+    1686     4.77916134607481733E-4     6.55478436841778708E-4
+    1687     2.50000000000000005E-4     5.99999999999999947E-4
+    1688     3.75421084408168141E-4     5.26840780869595538E-4
+    1689     3.79565818807853228E-4     5.99785262317693206E-4
+    1690     3.12710542204084046E-4     5.63420390434797797E-4
+    1691     3.77493451608010657E-4     5.63313021593644372E-4
+    1692     3.14782909403926617E-4     5.99892631158846631E-4
+    1693     3.78329331029246172E-4     4.03040117695809171E-4
+    1694     3.75421084408168141E-4     5.26840780869595538E-4
+    1695     3.14594753959753046E-4     4.86365680855438756E-4
+    1696     3.76875207718707183E-4     4.64940449282702354E-4
+    1697     3.45007919183960593E-4     5.06603230862517147E-4
+    1698     3.46462042494499636E-4     4.44702899275623963E-4
+    1699     9.02731575853669911E-5     2.19793387486037181E-4
+    1700     0     2.5585937500000001E-4
+    1701     0     1.91894531250000008E-4
+    1702     4.51365787926834956E-5     2.37826381243018609E-4
+    1703     0     2.23876953125000009E-4
+    1704     4.51365787926834956E-5     2.05843959368018608E-4
+    1705     0     2.5585937500000001E-4
+    1706     8.56507863843441119E-5     3.1689453125000001E-4
+    1707     0     3.1689453125000001E-4
+    1708     4.28253931921720559E-5     2.8637695312500001E-4
+    1709     4.28253931921720559E-5     3.1689453125000001E-4
+    1710     0     2.8637695312500001E-4
+    1711     1.01956056377857214E-4     1.1802584118703223E-4
+    1712     0     1.27929687500000005E-4
+    1713     0     6.39648437500000026E-5
+    1714     5.09780281889286071E-5     1.22977764343516131E-4
+    1715     0     9.59472656250000039E-5
+    1716     5.09780281889286071E-5     9.09953424685161163E-5
+    1717     0     6.39648437500000026E-5
+    1718     1.35354037059463752E-4     6.29907264959934823E-5
+    1719     1.01956056377857214E-4     1.1802584118703223E-4
+    1720     6.76770185297318762E-5     6.34777851229967425E-5
+    1721     1.18655046718660483E-4     9.05082838415128562E-5
+    1722     5.09780281889286071E-5     9.09953424685161163E-5
+    1723     2.44140625E-4     0
+    1724     1.35354037059463752E-4     6.29907264959934823E-5
+    1725     1.220703125E-4     0
+    1726     1.89747331029731876E-4     3.14953632479967412E-5
+    1727     1.28712174779731876E-4     3.14953632479967412E-5
+    1728     1.8310546875E-4     0
+    1729     3.3457565114347566E-4     2.07524696466039707E-4
+    1730     3.03283108388862287E-4     3.07191821602677501E-4
+    1731     2.62862460053854216E-4     2.58243951636348947E-4
+    1732     3.18929379766168947E-4     2.5735825903435859E-4
+    1733     2.83072784221358252E-4     2.82717886619513197E-4
+    1734     2.98719055598664911E-4     2.32884324051194341E-4
+    1735     3.78329331029246172E-4     4.03040117695809171E-4
+    1736     3.03283108388862287E-4     3.07191821602677501E-4
+    1737     3.662109375E-4     2.98836517333984323E-4
+    1738     3.40806219709054257E-4     3.55115969649243308E-4
+    1739     3.34747022944431117E-4     3.03014169468330912E-4
+    1740     3.72270134264623086E-4     3.5093831751489672E-4
+    1741     3.78329331029246172E-4     4.03040117695809171E-4
+    1742     3.14594753959753046E-4     4.86365680855438756E-4
+    1743     3.15575840731535889E-4     4.21078468404699436E-4
+    1744     3.46462042494499636E-4     4.44702899275623963E-4
+    1745     3.15085297345644468E-4     4.53722074630069096E-4
+    1746     3.46952585880391058E-4     4.12059293050254303E-4
+    1747     2.31175962939123751E-4     4.23329063260392325E-4
+    1748     3.15575840731535889E-4     4.21078468404699436E-4
+    1749     3.14594753959753046E-4     4.86365680855438756E-4
+    1750     2.73375901835329793E-4     4.22203765832545881E-4
+    1751     3.15085297345644468E-4     4.53722074630069096E-4
+    1752     2.72885358449438372E-4     4.54847372057915541E-4
+    1753     4.58849660034789083E-4     2.64192897481815869E-4
+    1754     4.15683819919475818E-4     3.45889508016053047E-4
+    1755     3.662109375E-4     2.98836517333984323E-4
+    1756     4.37266739977132477E-4     3.05041202748934431E-4
+    1757     3.90947378709737936E-4     3.22363012675018685E-4
+    1758     4.12530298767394541E-4     2.81514707407900123E-4
+    1759     5.469449348011982E-4     3.56411724829354306E-4
+    1760     4.57065334871855982E-4     4.13720111675068352E-4
+    1761     4.83502295899330232E-4     3.38002535036003411E-4
+    1762     5.02005134836527064E-4     3.85065918252211329E-4
+    1763     4.7028381538559308E-4     3.75861323355535881E-4
+    1764     5.15223615350264216E-4     3.47207129932678858E-4
+    1765     4.83502295899330232E-4     3.38002535036003411E-4
+    1766     4.57065334871855982E-4     4.13720111675068352E-4
+    1767     4.15683819919475818E-4     3.45889508016053047E-4
+    1768     4.7028381538559308E-4     3.75861323355535881E-4
+    1769     4.363745773956659E-4     3.798048098455607E-4
+    1770     4.49593057909403052E-4     3.41946021526028229E-4
+    1771     4.44306653512529842E-4     5.02492175855262643E-4
+    1772     4.57065334871855982E-4     4.13720111675068352E-4
+    1773     5.18154572622838414E-4     4.64528827497301226E-4
+    1774     4.50685994192192912E-4     4.5810614376516547E-4
+    1775     4.87609953747347225E-4     4.39124469586184762E-4
+    1776     4.81230613067684101E-4     4.83510501676281935E-4
+    1777     5.52246093750000003E-4     0
+    1778     5.46592219681165661E-4     1.47080610942937849E-4
+    1779     4.88018586665576635E-4     1.0715009043865536E-4
+    1780     5.49419156715582778E-4     7.35403054714689247E-5
+    1781     5.17305403173371121E-4     1.27115350690796598E-4
+    1782     5.20132340207788346E-4     5.35750452193276799E-5
+    1783     5.52246093750000003E-4     0
+    1784     6.15127063129517924E-4     1.28961775828680109E-4
+    1785     5.46592219681165661E-4     1.47080610942937849E-4
+    1786     5.83686578439758909E-4     6.44808879143400545E-5
+    1787     5.80859641405341793E-4     1.38021193385808993E-4
+    1788     5.49419156715582778E-4     7.35403054714689247E-5
+    1789     6.03238036253151279E-4     2.17057585313148337E-4
+    1790     5.46592219681165661E-4     1.47080610942937849E-4
+    1791     6.15127063129517924E-4     1.28961775828680109E-4
+    1792     5.74915127967158416E-4     1.82069098128043107E-4
+    1793     5.80859641405341793E-4     1.38021193385808993E-4
+    1794     6.09182549691334547E-4     1.73009680570914223E-4
+    1795     7.69147152428217916E-4     7.17795642973629926E-5
+    1796     6.80175781250000008E-4     6.50230427205121634E-5
+    1797     7.4414062500000001E-4     0
+    1798     7.24661466839108962E-4     6.84013035089375712E-5
+    1799     7.12158203124999955E-4     3.25115213602560817E-5
+    1800     7.56643888714109017E-4     3.58897821486814963E-5
+    1801     6.73095146243751124E-4     1.78431349377634624E-4
+    1802     6.80175781250000008E-4     6.50230427205121634E-5
+    1803     7.40024751298205222E-4     1.41990090695903195E-4
+    1804     6.7663546374687562E-4     1.21727196049073401E-4
+    1805     7.10100266274102561E-4     1.03506566708207686E-4
+    1806     7.06559948770978173E-4     1.6021072003676891E-4
+    1807     4.20397371381042544E-4     7.46936636915224194E-5
+    1808     3.662109375E-4     0
+    1809     4.2724609375E-4     0
+    1810     3.93304154440521272E-4     3.73468318457612097E-5
+    1811     3.96728515625E-4     0
+    1812     4.23821732565521272E-4     3.73468318457612097E-5
+    1813     4.8828125E-4     0
+    1814     4.88018586665576635E-4     1.0715009043865536E-4
+    1815     4.20397371381042544E-4     7.46936636915224194E-5
+    1816     4.8814991833278829E-4     5.35750452193276799E-5
+    1817     4.54207979023309562E-4     9.09218770650888964E-5
+    1818     4.54339310690521272E-4     3.73468318457612097E-5
+    1819     2.31915182548602798E-4     3.185630410561948E-4
+    1820     2.80075730378807853E-4     3.66277694634743636E-4
+    1821     2.31175962939123751E-4     4.23329063260392325E-4
+    1822     2.55995456463705339E-4     3.42420367845469218E-4
+    1823     2.55625846658965802E-4     3.94803378947567981E-4
+    1824     2.31545572743863289E-4     3.70946052158293562E-4
+    1825     5.469449348011982E-4     3.56411724829354306E-4
+    1826     5.82521566988404701E-4     4.12072931624541923E-4
+    1827     5.18154572622838414E-4     4.64528827497301226E-4
+    1828     5.64733250894801451E-4     3.84242328226948142E-4
+    1829     5.50338069805621612E-4     4.38300879560921575E-4
+    1830     5.32549753712018362E-4     4.10470276163327793E-4
+    1831     4.83502295899330232E-4     3.38002535036003411E-4
+    1832     5.8651429788997224E-4     3.03514350561540362E-4
+    1833     5.469449348011982E-4     3.56411724829354306E-4
+    1834     5.35008296894651236E-4     3.20758442798771914E-4
+    1835     5.6672961634558522E-4     3.29963037695447307E-4
+    1836     5.15223615350264216E-4     3.47207129932678858E-4
+    1837     6.55870398270680236E-4     4.58827279636259491E-4
+    1838     6.149060889867273E-4     3.58911302456159137E-4
+    1839     6.77112852597775263E-4     3.61199236444649594E-4
+    1840     6.35388243628703768E-4     4.08869291046209341E-4
+    1841     6.46009470792251227E-4     3.60055269450404339E-4
+    1842     6.66491625434227804E-4     4.10013258040454543E-4
+    1843     5.469449348011982E-4     3.56411724829354306E-4
+    1844     6.149060889867273E-4     3.58911302456159137E-4
+    1845     5.82521566988404701E-4     4.12072931624541923E-4
+    1846     5.80925511893962696E-4     3.57661513642756695E-4
+    1847     5.98713827987565946E-4     3.8549211704035053E-4
+    1848     5.64733250894801451E-4     3.84242328226948142E-4
+    1849     6.149060889867273E-4     3.58911302456159137E-4
+    1850     6.48721061501020204E-4     3.0580228455003082E-4
+    1851     6.77112852597775263E-4     3.61199236444649594E-4
+    1852     6.31813575243873752E-4     3.32356793503094951E-4
+    1853     6.62916957049397787E-4     3.33500760497340207E-4
+    1854     6.46009470792251227E-4     3.60055269450404339E-4
+    1855     6.03238036253151279E-4     2.17057585313148337E-4
+    1856     6.48721061501020204E-4     3.0580228455003082E-4
+    1857     5.8651429788997224E-4     3.03514350561540362E-4
+    1858     6.25979548877085741E-4     2.61429934931589592E-4
+    1859     6.17617679695496276E-4     3.04658317555785564E-4
+    1860     5.94876167071561705E-4     2.60285967937344336E-4
+    1861     7.66063425991097419E-4     2.25618978565681171E-4
+    1862     7.14223074903462172E-4     3.09513800853735191E-4
+    1863     6.83188905983988506E-4     2.53967124884082096E-4
+    1864     7.40143250447279849E-4     2.67566389709708195E-4
+    1865     6.98705990443725285E-4     2.81740462868908617E-4
+    1866     7.24626165987542962E-4     2.3979305172488162E-4
+    1867     7.47347240742430642E-4     3.78837751313385344E-4
+    1868     7.14223074903462172E-4     3.09513800853735191E-4
+    1869     7.77749415513451876E-4     3.13113370342386432E-4
+    1870     7.30785157822946407E-4     3.44175776083560268E-4
+    1871     7.4598624520845697E-4     3.11313585598060839E-4
+    1872     7.62548328127941313E-4     3.45975560827885861E-4
+    1873     6.03238036253151279E-4     2.17057585313148337E-4
+    1874     6.73095146243751124E-4     1.78431349377634624E-4
+    1875     6.83188905983988506E-4     2.53967124884082096E-4
+    1876     6.38166591248451147E-4     1.97744467345391467E-4
+    1877     6.78142026113869815E-4     2.1619923713085836E-4
+    1878     6.43213471118569838E-4     2.3551235509861523E-4
+    1879     6.80175781250000008E-4     6.50230427205121634E-5
+    1880     6.73095146243751124E-4     1.78431349377634624E-4
+    1881     6.15127063129517924E-4     1.28961775828680109E-4
+    1882     6.7663546374687562E-4     1.21727196049073401E-4
+    1883     6.44111104686634524E-4     1.53696562603157353E-4
+    1884     6.4765142218975902E-4     9.69924092745961294E-5
+    1885     0     3.1689453125000001E-4
+    1886     6.4572363977386294E-5     3.7792968750000001E-4
+    1887     0     3.7792968750000001E-4
+    1888     3.2286181988693147E-5     3.4741210937500001E-4
+    1889     3.2286181988693147E-5     3.7792968750000001E-4
+    1890     0     3.4741210937500001E-4
+    1891     6.4572363977386294E-5     3.7792968750000001E-4
+    1892     0     4.3896484375000001E-4
+    1893     0     3.7792968750000001E-4
+    1894     3.2286181988693147E-5     4.0844726562500001E-4
+    1895     0     4.0844726562500001E-4
+    1896     3.2286181988693147E-5     3.7792968750000001E-4
+    1897     1.71301572768688224E-4     3.7792968750000001E-4
+    1898     6.4572363977386294E-5     3.7792968750000001E-4
+    1899     8.56507863843441119E-5     3.1689453125000001E-4
+    1900     1.17936968373037266E-4     3.7792968750000001E-4
+    1901     7.51115751808651962E-5     3.4741210937500001E-4
+    1902     1.28476179576516181E-4     3.4741210937500001E-4
+    1903     1.25000000000000003E-4     5.49999999999999925E-4
+    1904     8.56507863843441119E-5     4.3896484375000001E-4
+    1905     1.72756436172360502E-4     4.70585859703467847E-4
+    1906     1.05325393192172057E-4     4.94482421874999968E-4
+    1907     1.29203611278352307E-4     4.54775351726733929E-4
+    1908     1.48878218086180266E-4     5.10292929851733886E-4
+    1909     6.25000000000000013E-5     5.24999999999999968E-4
+    1910     0     5.0000000000000001E-4
+    1911     8.56507863843441119E-5     4.3896484375000001E-4
+    1912     3.12500000000000007E-5     5.12499999999999935E-4
+    1913     4.28253931921720559E-5     4.6948242187500001E-4
+    1914     7.40753931921720566E-5     4.81982421874999989E-4
+    1915     8.13728732664855108E-4     1.33343473922578902E-4
+    1916     7.40024751298205222E-4     1.41990090695903195E-4
+    1917     7.69147152428217916E-4     7.17795642973629926E-5
+    1918     7.76876741981530165E-4     1.37666782309241035E-4
+    1919     7.54585951863211623E-4     1.06884827496633094E-4
+    1920     7.91437942546536566E-4     1.02561519109970947E-4
+    1921     6.80175781250000008E-4     6.50230427205121634E-5
+    1922     7.69147152428217916E-4     7.17795642973629926E-5
+    1923     7.40024751298205222E-4     1.41990090695903195E-4
+    1924     7.24661466839108962E-4     6.84013035089375712E-5
+    1925     7.54585951863211623E-4     1.06884827496633094E-4
+    1926     7.10100266274102561E-4     1.03506566708207686E-4
+    1927     8.08105468750000013E-4     0
+    1928     7.69147152428217916E-4     7.17795642973629926E-5
+    1929     7.4414062500000001E-4     0
+    1930     7.88626310589108964E-4     3.58897821486814963E-5
+    1931     7.56643888714109017E-4     3.58897821486814963E-5
+    1932     7.76123046875000066E-4     0
+    1933     9.35392068859210385E-4     6.524280030581833E-5
+    1934     8.72070312500000016E-4     0
+    1935     9.38964843750000021E-4     0
+    1936     9.037311906796052E-4     3.2621400152909165E-5
+    1937     9.05517578125000018E-4     0
+    1938     9.37178456304605203E-4     3.2621400152909165E-5
+    1939     9.38964843750000021E-4     0
+    1940     1.00000000000000002E-3     0
+    1941     9.35392068859210385E-4     6.524280030581833E-5
+    1942     9.69482421875000021E-4     0
+    1943     9.67696034429605203E-4     3.2621400152909165E-5
+    1944     9.37178456304605203E-4     3.2621400152909165E-5
+    1945     2.31175962939123751E-4     4.23329063260392325E-4
+    1946     2.39680731901657833E-4     5.14383951086584473E-4
+    1947     1.72756436172360502E-4     4.70585859703467847E-4
+    1948     2.35428347420390779E-4     4.68856507173488399E-4
+    1949     2.06218584037009167E-4     4.9248490539502616E-4
+    1950     2.0196619955574214E-4     4.46957461481930086E-4
+    1951     2.50000000000000005E-4     5.99999999999999947E-4
+    1952     2.39680731901657833E-4     5.14383951086584473E-4
+    1953     3.14594753959753046E-4     4.86365680855438756E-4
+    1954     2.44840365950828906E-4     5.57191975543292264E-4
+    1955     2.77137742930705453E-4     5.00374815971011615E-4
+    1956     2.82297376979876499E-4     5.43182840427719298E-4
+    1957     7.66063425991097419E-4     2.25618978565681171E-4
+    1958     7.40024751298205222E-4     1.41990090695903195E-4
+    1959     8.13728732664855108E-4     1.33343473922578902E-4
+    1960     7.53044088644651321E-4     1.83804534630792197E-4
+    1961     7.76876741981530165E-4     1.37666782309241035E-4
+    1962     7.89896079327976263E-4     1.79481226244130037E-4
+    1963     7.77749415513451876E-4     3.13113370342386432E-4
+    1964     7.66063425991097419E-4     2.25618978565681171E-4
+    1965     8.40965569411703819E-4     2.71022418711743549E-4
+    1966     7.71906420752274647E-4     2.69366174454033788E-4
+    1967     8.03514497701400673E-4     2.48320698638712374E-4
+    1968     8.09357492462577793E-4     2.92067894527064963E-4
+    1969     3.88700445385212667E-4     1.12415429991009354E-3
+    1970     4.14780023525604651E-4     1.04732934218306818E-3
+    1971     4.64295327735872877E-4     1.106418874786976E-3
+    1972     4.01740234455408686E-4     1.08574182104658086E-3
+    1973     4.39537675630738791E-4     1.07687410848502198E-3
+    1974     4.26497886560542772E-4     1.11528658734853488E-3
+    1975     4.14780023525604651E-4     1.04732934218306818E-3
+    1976     3.39525065135284183E-4     1.0640634227390719E-3
+    1977     3.33510000366525681E-4     9.78187713756657664E-4
+    1978     3.7715254433044439E-4     1.05569638246107004E-3
+    1979     3.36517532750904959E-4     1.02112556824786489E-3
+    1980     3.74145011946065166E-4     1.01275852796986303E-3
+    1981     2.97380987249783093E-4     1.12592184443437615E-3
+    1982     3.72070312500000005E-4     1.19999999999999989E-3
+    1983     3.08105468750000003E-4     1.19999999999999989E-3
+    1984     3.34725649874891549E-4     1.16296092221718802E-3
+    1985     3.40087890625000004E-4     1.19999999999999989E-3
+    1986     3.02743227999891548E-4     1.16296092221718802E-3
+    1987     4.14780023525604651E-4     1.04732934218306818E-3
+    1988     3.88700445385212667E-4     1.12415429991009354E-3
+    1989     3.39525065135284183E-4     1.0640634227390719E-3
+    1990     4.01740234455408686E-4     1.08574182104658086E-3
+    1991     3.64112755260248425E-4     1.09410886132458272E-3
+    1992     3.7715254433044439E-4     1.05569638246107004E-3
+    1993     4.10887142695853748E-4     8.117939619078767E-4
+    1994     3.3342684329849511E-4     7.35743010166048958E-4
+    1995     3.98747049720822426E-4     7.15064925451776041E-4
+    1996     3.72156992997174402E-4     7.73768486036962829E-4
+    1997     3.66086946509658795E-4     7.25403967808912553E-4
+    1998     4.04817096208338087E-4     7.63429443679826316E-4
+    1999     3.3342684329849511E-4     7.35743010166048958E-4
+    2000     3.00868592099089952E-4     7.96027940566129363E-4
+    2001     2.29766755665683817E-4     7.70779974744912363E-4
+    2002     3.17147717698792531E-4     7.6588547536608916E-4
+    2003     2.65317673882386898E-4     7.83403957655520863E-4
+    2004     2.81596799482089477E-4     7.53261492455480661E-4
+    2005     1.36565743695333484E-4     7.39812343895404297E-4
+    2006     2.29766755665683817E-4     7.70779974744912363E-4
+    2007     1.74429316483144684E-4     8.22070312499999993E-4
+    2008     1.8316624968050865E-4     7.5529615932015833E-4
+    2009     2.02098036074414237E-4     7.96425143622456178E-4
+    2010     1.55497530089239097E-4     7.80941328197702145E-4
+    2011     1.87500000000000004E-4     6.24999999999999905E-4
+    2012     2.50000000000000005E-4     5.99999999999999947E-4
+    2013     2.16128514871350092E-4     6.96571287178375043E-4
+    2014     2.18750000000000005E-4     6.1249999999999998E-4
+    2015     2.33064257435675048E-4     6.48285643589187549E-4
+    2016     2.01814257435675048E-4     6.60785643589187474E-4
+    2017     3.3342684329849511E-4     7.35743010166048958E-4
+    2018     2.29766755665683817E-4     7.70779974744912363E-4
+    2019     2.16128514871350092E-4     6.96571287178375043E-4
+    2020     2.81596799482089477E-4     7.53261492455480661E-4
+    2021     2.22947635268516968E-4     7.33675630961643649E-4
+    2022     2.74777679084922601E-4     7.16157148672211946E-4
+    2023     4.10887142695853748E-4     8.117939619078767E-4
+    2024     3.14063844356556141E-4     8.63789983373621688E-4
+    2025     3.00868592099089952E-4     7.96027940566129363E-4
+    2026     3.62475493526204944E-4     8.37791972640749248E-4
+    2027     3.07466218227823019E-4     8.2990896196987558E-4
+    2028     3.55877867397471877E-4     8.03910951237003031E-4
+    2029     3.33510000366525681E-4     9.78187713756657664E-4
+    2030     3.14063844356556141E-4     8.63789983373621688E-4
+    2031     3.72070312500000005E-4     9.01220703124999949E-4
+    2032     3.23786922361540938E-4     9.20988848565139676E-4
+    2033     3.43067078428278046E-4     8.82505343249310819E-4
+    2034     3.52790156433262843E-4     9.39704208440828807E-4
+    2035     1.74429316483144684E-4     8.22070312499999993E-4
+    2036     2.56894004484586392E-4     9.02486475935557219E-4
+    2037     1.85578989176374238E-4     8.9312739269816724E-4
+    2038     2.15661660483865538E-4     8.62278394217778606E-4
+    2039     2.21236496830480329E-4     8.97806934316862284E-4
+    2040     1.80004152829759474E-4     8.57598852599083562E-4
+    2041     2.29999438806578629E-4     9.67017715356922923E-4
+    2042     1.85578989176374238E-4     8.9312739269816724E-4
+    2043     2.56894004484586392E-4     9.02486475935557219E-4
+    2044     2.07789213991476447E-4     9.30072554027545082E-4
+    2045     2.21236496830480329E-4     8.97806934316862284E-4
+    2046     2.43446721645582511E-4     9.34752095646240126E-4
+    2047     0     8.83105468749999993E-4
+    2048     6.75304712410491789E-5     8.77236181658321186E-4
+    2049     6.33482805365148801E-5     9.51708311133460578E-4
+    2050     3.37652356205245894E-5     8.80170825204160644E-4
+    2051     6.54393758887820227E-5     9.14472246395890828E-4
+    2052     3.16741402682574401E-5     9.17406889941730285E-4
+    2053     0     1.00810546875E-3
+    2054     6.33482805365148801E-5     9.51708311133460578E-4
+    2055     8.06327547355035015E-5     1.013121025387205E-3
+    2056     3.16741402682574401E-5     9.79906889941730341E-4
+    2057     7.19905176360091908E-5     9.82414668260332734E-4
+    2058     4.03163773677517507E-5     1.01061324706860261E-3
+    2059     6.33482805365148801E-5     9.51708311133460578E-4
+    2060     0     1.00810546875E-3
+    2061     0     9.44140624999999993E-4
+    2062     3.16741402682574401E-5     9.79906889941730341E-4
+    2063     0     9.7612304687499994E-4
+    2064     3.16741402682574401E-5     9.47924468066730285E-4
+    2065     0     8.83105468749999993E-4
+    2066     6.33482805365148801E-5     9.51708311133460578E-4
+    2067     0     9.44140624999999993E-4
+    2068     3.16741402682574401E-5     9.17406889941730285E-4
+    2069     3.16741402682574401E-5     9.47924468066730285E-4
+    2070     0     9.13623046874999993E-4
+    2071     2.25989861581488047E-4     1.1034297389022796E-3
+    2072     2.44140625E-4     1.19999999999999989E-3
+    2073     1.8310546875E-4     1.19999999999999989E-3
+    2074     2.3506524329074401E-4     1.15171486945113985E-3
+    2075     2.13623046875E-4     1.19999999999999989E-3
+    2076     2.0454766516574401E-4     1.15171486945113985E-3
+    2077     3.88700445385212667E-4     1.12415429991009354E-3
+    2078     2.97380987249783093E-4     1.12592184443437615E-3
+    2079     3.39525065135284183E-4     1.0640634227390719E-3
+    2080     3.43040716317497907E-4     1.12503807217223485E-3
+    2081     3.18453026192533611E-4     1.09499263358672403E-3
+    2082     3.64112755260248425E-4     1.09410886132458272E-3
+    2083     2.44140625E-4     1.19999999999999989E-3
+    2084     2.97380987249783093E-4     1.12592184443437615E-3
+    2085     3.08105468750000003E-4     1.19999999999999989E-3
+    2086     2.70760806124891546E-4     1.16296092221718802E-3
+    2087     3.02743227999891548E-4     1.16296092221718802E-3
+    2088     2.76123046875000001E-4     1.19999999999999989E-3
+    2089     0     1.07207031249999989E-3
+    2090     6.77802435850866764E-5     1.13246159963255083E-3
+    2091     0     1.13603515625E-3
+    2092     3.38901217925433382E-5     1.10226595606627547E-3
+    2093     3.38901217925433382E-5     1.13424837794127552E-3
+    2094     0     1.10405273437499995E-3
+    2095     1.220703125E-4     1.19999999999999989E-3
+    2096     6.77802435850866764E-5     1.13246159963255083E-3
+    2097     1.27766480316576791E-4     1.1007034017909631E-3
+    2098     9.49252780425433382E-5     1.16623079981627536E-3
+    2099     9.77733619508317338E-5     1.11658250071175686E-3
+    2100     1.24918396408288396E-4     1.15035170089548139E-3
+    2101     8.72146582415723147E-5     8.12372124900119985E-4
+    2102     6.7318724869494569E-5     7.44219603001817923E-4
+    2103     1.36565743695333484E-4     7.39812343895404297E-4
+    2104     7.72666915555334486E-5     7.78295863950968954E-4
+    2105     1.0194223428241402E-4     7.42015973448611164E-4
+    2106     1.11890200968452899E-4     7.76092234397762195E-4
+    2107     0     7.0E-4
+    2108     6.7318724869494569E-5     7.44219603001817923E-4
+    2109     0     7.61035156249999993E-4
+    2110     3.36593624347472845E-5     7.22109801500908904E-4
+    2111     3.36593624347472845E-5     7.52627379625908904E-4
+    2112     0     7.30517578124999993E-4
+    2113     6.3613763360561681E-4     1.03946949707743879E-3
+    2114     6.31385186434097476E-4     1.12635085509955684E-3
+    2115     5.63964843750000013E-4     1.0965073564807807E-3
+    2116     6.33761410019857089E-4     1.08291017608849781E-3
+    2117     5.9767501509204869E-4     1.11142910579016888E-3
+    2118     6.00051238677808412E-4     1.06798842677910963E-3
+    2119     6.91894531250000018E-4     1.19999999999999989E-3
+    2120     6.31385186434097476E-4     1.12635085509955684E-3
+    2121     6.91894531250000018E-4     1.08422299961888672E-3
+    2122     6.61639858842048801E-4     1.16317542754977837E-3
+    2123     6.61639858842048801E-4     1.10528692735922178E-3
+    2124     6.91894531250000018E-4     1.14211149980944331E-3
+    2125     9.39685899131729899E-4     7.34463702578968722E-4
+    2126     9.36226394521655E-4     8.3396748156078118E-4
+    2127     8.9886962494534965E-4     7.81508352249731973E-4
+    2128     9.37956146826692504E-4     7.84215592069874897E-4
+    2129     9.17548009733502325E-4     8.07737916905256577E-4
+    2130     9.1927776203853972E-4     7.57986027414350293E-4
+    2131     1.00000000000000002E-3     9.0E-4
+    2132     9.36226394521655E-4     8.3396748156078118E-4
+    2133     1.00000000000000002E-3     8.24999999999999995E-4
+    2134     9.68113197260827511E-4     8.66983740780390578E-4
+    2135     9.68113197260827511E-4     8.29483740780390588E-4
+    2136     1.00000000000000002E-3     8.62499999999999985E-4
+    2137     7.42631330235908388E-4     4.52617545905931022E-4
+    2138     7.47347240742430642E-4     3.78837751313385344E-4
+    2139     8.08703962690366035E-4     4.17299929575055159E-4
+    2140     7.44989285489169569E-4     4.15727648609658183E-4
+    2141     7.78025601716398393E-4     3.98068840444220251E-4
+    2142     7.75667646463137211E-4     4.3495873774049309E-4
+    2143     7.14223074903462172E-4     3.09513800853735191E-4
+    2144     7.47347240742430642E-4     3.78837751313385344E-4
+    2145     6.77112852597775263E-4     3.61199236444649594E-4
+    2146     7.30785157822946407E-4     3.44175776083560268E-4
+    2147     7.12230046670103006E-4     3.70018493879017496E-4
+    2148     6.95667963750618663E-4     3.35356518649192366E-4
+    2149     6.28726770153248599E-4     8.69566754220071274E-4
+    2150     6.3251725796544824E-4     7.69458270396488286E-4
+    2151     6.94054151555351649E-4     7.95336883715192389E-4
+    2152     6.30622014059348474E-4     8.1951251230827978E-4
+    2153     6.63285704760399998E-4     7.82397577055840392E-4
+    2154     6.61390460854300124E-4     8.32451818967631886E-4
+    2155     5.60394849123839033E-4     7.35715772264671007E-4
+    2156     6.3251725796544824E-4     7.69458270396488286E-4
+    2157     5.79668274240937665E-4     8.10244013239797283E-4
+    2158     5.96456053544643637E-4     7.52587021330579646E-4
+    2159     6.06092766103192952E-4     7.89851141818142784E-4
+    2160     5.70031561682388349E-4     7.72979892752234145E-4
+    2161     4.25259790828855515E-4     1.93997913774442133E-4
+    2162     4.20397371381042544E-4     7.46936636915224194E-5
+    2163     4.88018586665576635E-4     1.0715009043865536E-4
+    2164     4.22828581104949029E-4     1.34345788732982262E-4
+    2165     4.54207979023309562E-4     9.09218770650888964E-5
+    2166     4.56639188747216102E-4     1.50574002106548739E-4
+    2167     4.20397371381042544E-4     7.46936636915224194E-5
+    2168     3.64813480721653041E-4     1.25057137488653904E-4
+    2169     3.00937473641497586E-4     6.59440565847406544E-5
+    2170     3.9260542605134782E-4     9.98754005900881616E-5
+    2171     3.32875477181575341E-4     9.55005970366972791E-5
+    2172     3.60667422511270065E-4     7.03188601381315369E-5
+    2173     1.00000000000000002E-3     1.49999999999999987E-4
+    2174     9.35392068859210385E-4     6.524280030581833E-5
+    2175     1.00000000000000002E-3     7.49999999999999934E-5
+    2176     9.67696034429605203E-4     1.07621400152909158E-4
+    2177     9.67696034429605203E-4     7.01214001529091549E-5
+    2178     1.00000000000000002E-3     1.12499999999999997E-4
+    2179     8.72070312500000016E-4     0
+    2180     9.35392068859210385E-4     6.524280030581833E-5
+    2181     8.79865520605121044E-4     9.96839116576813268E-5
+    2182     9.037311906796052E-4     3.2621400152909165E-5
+    2183     9.07628794732165661E-4     8.24633559817498284E-5
+    2184     8.75967916552560476E-4     4.98419558288406634E-5
+    2185     1.74429316483144684E-4     8.22070312499999993E-4
+    2186     6.75304712410491789E-5     8.77236181658321186E-4
+    2187     8.72146582415723147E-5     8.12372124900119985E-4
+    2188     1.20979893862096924E-4     8.49653247079160644E-4
+    2189     7.737256474131074E-5     8.44804153279220585E-4
+    2190     1.30821987362358499E-4     8.17221218700059935E-4
+    2191     6.33482805365148801E-5     9.51708311133460578E-4
+    2192     6.75304712410491789E-5     8.77236181658321186E-4
+    2193     1.1927993626285668E-4     9.2101752411100806E-4
+    2194     6.54393758887820227E-5     9.14472246395890828E-4
+    2195     9.34052037519529361E-5     8.99126852884664677E-4
+    2196     9.131410839968578E-5     9.36362917622234319E-4
+    254
+    43     1     1     0     0
+    49     1     0     0     20.
+    51     1     0     0     20.
+    54     1     0     0     20.
+    67     1     1     0     0
+    69     1     1     0     0
+    72     1     1     0     0
+    110     1     1     0     0
+    128     1     1     0     0
+    129     1     1     0     0
+    131     1     1     0     0
+    133     1     0     0     20.
+    135     1     0     0     20.
+    138     1     0     0     20.
+    152     1     1     0     0
+    153     1     1     0     0
+    155     1     1     0     0
+    159     1     0     0     20.
+    188     1     1     0     0
+    189     1     0     0     0
+    205     1     0     0     20.
+    218     1     0     0     20.
+    219     1     0     0     20.
+    221     1     0     0     20.
+    291     1     1     0     0
+    297     1     1     0     0
+    302     1     0     0     20.
+    313     1     0     0     20.
+    332     1     0     0     20.
+    349     1     0     0     20.
+    351     1     0     0     20.
+    354     1     0     0     20.
+    355     1     1     0     0
+    357     1     1     0     0
+    360     1     1     0     0
+    367     1     1     0     0
+    368     1     0     0     20.
+    369     1     0     0     20.
+    371     1     0     0     20.
+    370     1     1     0     0
+    386     1     1     0     0
+    387     1     1     0     0
+    389     1     1     0     0
+    403     1     1     0     0
+    405     1     1     0     0
+    408     1     1     0     0
+    417     1     1     0     0
+    428     1     1     0     0
+    458     1     1     0     0
+    459     1     1     0     0
+    461     1     1     0     0
+    464     1     1     0     0
+    465     1     1     0     0
+    467     1     1     0     0
+    470     1     0     0     0
+    471     1     0     0     0
+    473     1     0     0     0
+    477     1     1     0     0
+    481     1     1     0     0
+    656     1     1     0     0
+    680     1     1     0     0
+    681     1     1     0     0
+    683     1     1     0     0
+    685     1     1     0     0
+    686     1     1     0     0
+    688     1     1     0     0
+    697     1     1     0     0
+    704     1     1     0     0
+    709     1     1     0     0
+    721     1     1     0     0
+    723     1     1     0     0
+    726     1     1     0     0
+    733     1     1     0     0
+    734     1     0     0     0
+    735     1     0     0     0
+    737     1     0     0     0
+    736     1     1     0     0
+    740     1     0     0     0
+    757     1     0     0     0
+    758     1     0     0     0
+    760     1     0     0     0
+    781     1     0     0     0
+    878     1     0     0     0
+    879     1     0     0     0
+    881     1     0     0     0
+    884     1     0     0     0
+    885     1     0     0     0
+    887     1     0     0     0
+    897     1     0     0     0
+    907     1     0     0     0
+    909     1     0     0     0
+    912     1     0     0     0
+    961     1     0     0     0
+    967     1     0     0     0
+    968     1     0     0     0
+    970     1     0     0     0
+    985     1     1     0     0
+    992     1     1     0     0
+    1039     1     1     0     0
+    1040     1     1     0     0
+    1042     1     1     0     0
+    1057     1     1     0     0
+    1059     1     1     0     0
+    1062     1     1     0     0
+    1069     1     1     0     0
+    1071     1     0     0     0
+    1074     1     1     0     0
+    1076     1     1     0     0
+    1077     1     1     0     0
+    1079     1     1     0     0
+    1082     1     0     0     0
+    1087     1     1     0     0
+    1088     1     1     0     0
+    1090     1     1     0     0
+    1111     1     0     0     20.
+    1112     1     0     0     20.
+    1114     1     0     0     20.
+    1117     1     1     0     0
+    1119     1     1     0     0
+    1122     1     1     0     0
+    1131     1     1     0     0
+    1137     1     1     0     0
+    1166     1     1     0     0
+    1179     1     1     0     0
+    1207     1     1     0     0
+    1209     1     1     0     0
+    1212     1     1     0     0
+    1226     1     1     0     0
+    1232     1     0     0     20.
+    1237     1     1     0     0
+    1255     1     0     0     20.
+    1256     1     1     0     0
+    1257     1     0     0     20.
+    1260     1     0     0     20.
+    1258     1     1     0     0
+    1262     1     1     0     0
+    1263     1     1     0     0
+    1265     1     1     0     0
+    1274     1     0     0     20.
+    1275     1     0     0     20.
+    1277     1     0     0     20.
+    1363     1     1     0     0
+    1364     1     1     0     0
+    1366     1     1     0     0
+    1371     1     0     0     20.
+    1412     1     0     0     20.
+    1413     1     1     0     0
+    1418     1     1     0     0
+    1424     1     1     0     0
+    1430     1     1     0     0
+    1431     1     1     0     0
+    1433     1     1     0     0
+    1477     1     0     0     20.
+    1478     1     0     0     20.
+    1480     1     0     0     20.
+    1483     1     0     0     20.
+    1484     1     0     0     20.
+    1486     1     0     0     20.
+    1490     1     0     0     20.
+    1491     1     0     0     20.
+    1493     1     0     0     20.
+    1507     1     0     0     20.
+    1514     1     0     0     20.
+    1520     1     0     0     20.
+    1521     1     0     0     20.
+    1523     1     0     0     20.
+    1526     1     1     0     0
+    1527     1     0     0     20.
+    1537     1     1     0     0
+    1543     1     1     0     0
+    1544     1     1     0     0
+    1546     1     1     0     0
+    1557     1     1     0     0
+    1568     1     1     0     0
+    1569     1     1     0     0
+    1571     1     1     0     0
+    1573     1     1     0     0
+    1585     1     1     0     0
+    1687     1     1     0     0
+    1700     1     1     0     0
+    1701     1     1     0     0
+    1703     1     1     0     0
+    1705     1     1     0     0
+    1707     1     1     0     0
+    1710     1     1     0     0
+    1712     1     1     0     0
+    1713     1     1     0     0
+    1715     1     1     0     0
+    1717     1     1     0     0
+    1723     1     0     0     0
+    1725     1     0     0     0
+    1728     1     0     0     0
+    1777     1     0     0     0
+    1783     1     0     0     0
+    1797     1     0     0     0
+    1808     1     0     0     0
+    1809     1     0     0     0
+    1811     1     0     0     0
+    1813     1     0     0     0
+    1885     1     1     0     0
+    1887     1     1     0     0
+    1890     1     1     0     0
+    1892     1     1     0     0
+    1893     1     1     0     0
+    1895     1     1     0     0
+    1903     1     1     0     0
+    1909     1     1     0     0
+    1910     1     1     0     0
+    1912     1     1     0     0
+    1927     1     0     0     0
+    1929     1     0     0     0
+    1932     1     0     0     0
+    1934     1     0     0     0
+    1935     1     0     0     0
+    1937     1     0     0     0
+    1939     1     0     0     0
+    1940     1     0     0     0
+    1942     1     0     0     0
+    1951     1     1     0     0
+    1982     1     0     0     20.
+    1983     1     0     0     20.
+    1985     1     0     0     20.
+    2011     1     1     0     0
+    2012     1     1     0     0
+    2014     1     1     0     0
+    2047     1     1     0     0
+    2053     1     1     0     0
+    2060     1     1     0     0
+    2061     1     1     0     0
+    2063     1     1     0     0
+    2065     1     1     0     0
+    2067     1     1     0     0
+    2070     1     1     0     0
+    2072     1     0     0     20.
+    2073     1     0     0     20.
+    2075     1     0     0     20.
+    2083     1     0     0     20.
+    2085     1     0     0     20.
+    2088     1     0     0     20.
+    2089     1     1     0     0
+    2091     1     1     0     0
+    2094     1     1     0     0
+    2095     1     0     0     20.
+    2107     1     1     0     0
+    2109     1     1     0     0
+    2112     1     1     0     0
+    2119     1     0     0     20.
+    2131     1     1     0     0
+    2133     1     1     0     0
+    2136     1     1     0     0
+    2173     1     1     0     0
+    2175     1     1     0     0
+    2178     1     1     0     0
+    2179     1     0     0     0
+    0     516     0     0     0
+    1     2     3     919     920     5     922     -1
+    1     1     3     195     194     6     197     -1
+    1     1     2     789     787     4     792     -1
+    1     8     9     590     591     11     593     -1
+    1     7     9     1579     1581     12     1584     -1
+    1     7     8     435     433     10     438     -1
+    1     14     15     842     843     17     845     -1
+    1     13     15     788     787     18     790     -1
+    1     13     14     751     752     16     754     -1
+    1     20     21     170     171     23     173     -1
+    1     19     21     1250     1249     24     1252     -1
+    1     19     20     1377     1375     22     1380     -1
+    1     26     27     848     849     29     851     -1
+    1     25     27     842     841     30     844     -1
+    1     25     26     830     831     28     833     -1
+    1     32     33     453     451     35     456     -1
+    1     31     33     1047     1046     36     1049     -1
+    1     31     32     500     501     34     503     -1
+    1     38     39     2103     2101     41     2106     -1
+    1     37     39     2005     2007     42     2010     -1
+    1     37     38     2187     2185     40     2190     -1
+    1     44     45     649     650     47     652     -1
+    1     43     45     1177     1179     48     1182     -1
+    1     43     44     1689     1687     46     1692     -1
+    1     50     51     1232     1233     53     1235     -1
+    1     49     50     2097     2095     52     2100     -1
+    1     56     57     398     399     59     401     -1
+    1     55     57     1646     1645     60     1648     -1
+    1     55     56     146     147     58     149     -1
+    1     62     63     566     567     65     569     -1
+    1     61     63     573     572     66     575     -1
+    1     61     62     446     447     64     449     -1
+    1     68     69     2047     2048     71     2050     -1
+    1     67     68     1239     1237     70     1242     -1
+    1     74     75     1327     1328     77     1330     -1
+    1     73     75     560     559     78     562     -1
+    1     73     74     1333     1334     76     1336     -1
+    1     80     81     338     339     83     341     -1
+    1     79     81     637     639     84     642     -1
+    1     79     80     1454     1455     82     1457     -1
+    1     86     87     314     315     89     317     -1
+    1     85     87     1107     1106     90     1109     -1
+    1     85     86     1461     1459     88     1464     -1
+    1     92     93     254     255     95     257     -1
+    1     91     93     615     614     96     617     -1
+    1     91     92     393     391     94     396     -1
+    1     98     99     1677     1675     101     1680     -1
+    1     97     99     529     531     102     534     -1
+    1     97     98     1623     1621     100     1626     -1
+    1     104     105     710     711     107     713     -1
+    1     103     105     728     727     108     730     -1
+    1     103     104     716     717     106     719     -1
+    1     110     111     1089     1087     113     1092     -1
+    1     109     111     1947     1946     114     1949     -1
+    1     109     110     1077     1075     112     1080     -1
+    1     116     117     1021     1022     119     1024     -1
+    1     115     117     943     945     120     948     -1
+    1     115     116     272     273     118     275     -1
+    1     122     123     1803     1801     125     1806     -1
+    1     121     123     1958     1957     126     1960     -1
+    1     121     122     950     951     124     953     -1
+    1     127     129     428     427     132     430     -1
+    1     127     128     465     463     130     468     -1
+    1     134     135     349     350     137     352     -1
+    1     133     134     301     302     136     304     -1
+    1     140     141     1453     1454     143     1456     -1
+    1     139     141     2113     2115     144     2118     -1
+    1     139     140     1460     1461     142     1463     -1
+    1     145     147     619     621     150     624     -1
+    1     145     146     2151     2149     148     2154     -1
+    1     151     153     2011     2013     156     2016     -1
+    1     151     152     1166     1167     154     1169     -1
+    1     158     159     1112     1113     161     1115     -1
+    1     157     159     1982     1981     162     1984     -1
+    1     157     158     2077     2078     160     2080     -1
+    1     164     165     1155     1153     167     1158     -1
+    1     163     165     2036     2035     168     2038     -1
+    1     163     164     2006     2007     166     2009     -1
+    1     169     171     1268     1267     174     1270     -1
+    1     169     170     1245     1243     172     1248     -1
+    1     176     177     854     855     179     857     -1
+    1     175     177     668     667     180     670     -1
+    1     175     176     920     921     178     923     -1
+    1     182     183     1401     1399     185     1404     -1
+    1     181     183     805     807     186     810     -1
+    1     181     182     1743     1741     184     1746     -1
+    1     188     189     735     733     191     738     -1
+    1     187     189     1725     1724     192     1727     -1
+    1     187     188     1717     1718     190     1720     -1
+    1     193     195     727     729     198     732     -1
+    1     193     194     903     901     196     906     -1
+    1     200     201     1977     1975     203     1980     -1
+    1     199     201     1310     1309     204     1312     -1
+    1     199     200     1323     1321     202     1326     -1
+    1     206     207     1971     1969     209     1974     -1
+    1     205     207     314     313     210     316     -1
+    1     205     206     1113     1111     208     1116     -1
+    1     212     213     1443     1441     215     1446     -1
+    1     211     213     1322     1321     216     1324     -1
+    1     211     212     1315     1316     214     1318     -1
+    1     217     219     313     315     222     318     -1
+    1     217     218     1514     1515     220     1517     -1
+    1     224     225     1994     1995     227     1997     -1
+    1     223     225     1173     1172     228     1175     -1
+    1     223     224     651     649     226     654     -1
+    1     230     231     2156     2157     233     2159     -1
+    1     229     231     2150     2149     234     2152     -1
+    1     229     230     643     644     232     646     -1
+    1     236     237     259     260     239     262     -1
+    1     235     237     548     547     240     550     -1
+    1     235     236     1688     1689     238     1691     -1
+    1     242     243     1351     1352     245     1354     -1
+    1     241     243     1436     1435     246     1438     -1
+    1     241     242     1357     1358     244     1360     -1
+    1     248     249     595     596     251     598     -1
+    1     247     249     608     607     252     610     -1
+    1     247     248     583     584     250     586     -1
+    1     253     255     277     279     258     282     -1
+    1     253     254     381     379     256     384     -1
+    1     260     261     825     823     263     828     -1
+    1     259     261     1676     1675     264     1678     -1
+    1     266     267     1772     1773     269     1775     -1
+    1     265     267     1760     1759     270     1762     -1
+    1     265     266     1827     1825     268     1830     -1
+    1     271     273     565     567     276     570     -1
+    1     271     272     579     577     274     582     -1
+    1     278     279     411     409     281     414     -1
+    1     277     278     1549     1550     280     1552     -1
+    1     284     285     495     493     287     498     -1
+    1     283     285     1593     1592     288     1595     -1
+    1     283     284     1022     1023     286     1025     -1
+    1     290     291     1117     1118     293     1120     -1
+    1     289     291     1131     1130     294     1133     -1
+    1     289     290     1124     1125     292     1127     -1
+    1     296     297     1544     1545     299     1547     -1
+    1     295     297     357     356     300     359     -1
+    1     295     296     363     361     298     366     -1
+    1     302     303     2121     2119     305     2124     -1
+    1     301     303     1495     1497     306     1500     -1
+    1     308     309     619     620     311     622     -1
+    1     307     309     614     613     312     616     -1
+    1     307     308     1653     1651     310     1656     -1
+    1     320     321     333     331     323     336     -1
+    1     319     321     1497     1496     324     1499     -1
+    1     319     320     1533     1531     322     1536     -1
+    1     326     327     1508     1509     329     1511     -1
+    1     325     327     1533     1532     330     1535     -1
+    1     325     326     361     362     328     364     -1
+    1     332     333     350     351     335     353     -1
+    1     331     332     1491     1489     334     1494     -1
+    1     337     339     1460     1459     342     1462     -1
+    1     337     338     1442     1443     340     1445     -1
+    1     344     345     631     632     347     634     -1
+    1     343     345     627     626     348     629     -1
+    1     343     344     1495     1496     346     1498     -1
+    1     355     356     1525     1526     358     1528     -1
+    1     362     363     373     374     365     376     -1
+    1     367     369     1527     1526     372     1529     -1
+    1     374     375     1503     1501     377     1506     -1
+    1     373     375     381     380     378     383     -1
+    1     379     380     1538     1539     382     1541     -1
+    1     385     387     1543     1545     390     1548     -1
+    1     385     386     1537     1538     388     1540     -1
+    1     392     393     1550     1551     395     1553     -1
+    1     391     392     1647     1645     394     1650     -1
+    1     397     399     597     596     402     599     -1
+    1     397     398     1640     1641     400     1643     -1
+    1     404     405     1557     1555     407     1560     -1
+    1     403     404     2132     2133     406     2135     -1
+    1     410     411     415     416     413     418     -1
+    1     409     410     2126     2127     412     2129     -1
+    1     416     417     2131     2132     419     2134     -1
+    1     415     417     1537     1539     420     1542     -1
+    1     422     423     1639     1640     425     1642     -1
+    1     421     423     692     691     426     694     -1
+    1     421     422     1646     1647     424     1649     -1
+    1     428     429     1575     1573     431     1578     -1
+    1     427     429     1581     1580     432     1583     -1
+    1     434     435     1561     1562     437     1564     -1
+    1     433     434     1574     1575     436     1577     -1
+    1     440     441     447     445     443     450     -1
+    1     439     441     1587     1586     444     1589     -1
+    1     439     440     1579     1580     442     1582     -1
+    1     445     446     506     507     448     509     -1
+    1     452     453     1869     1867     455     1872     -1
+    1     451     452     2138     2139     454     2141     -1
+    1     457     459     1137     1136     462     1139     -1
+    1     457     458     1131     1129     460     1134     -1
+    1     463     464     1585     1586     466     1588     -1
+    1     469     471     961     963     474     966     -1
+    1     469     470     1783     1784     472     1786     -1
+    1     476     477     2173     2174     479     2176     -1
+    1     475     477     1059     1058     480     1061     -1
+    1     475     476     980     981     478     983     -1
+    1     482     483     1004     1005     485     1007     -1
+    1     481     483     1041     1040     486     1043     -1
+    1     481     482     1118     1119     484     1121     -1
+    1     488     489     1917     1915     491     1920     -1
+    1     487     489     1081     1083     492     1086     -1
+    1     487     488     1063     1064     490     1066     -1
+    1     494     495     609     607     497     612     -1
+    1     493     494     578     579     496     581     -1
+    1     499     501     1141     1143     504     1146     -1
+    1     499     500     1629     1627     502     1632     -1
+    1     505     507     585     584     510     587     -1
+    1     505     506     577     578     508     580     -1
+    1     512     513     1027     1028     515     1030     -1
+    1     511     513     1597     1599     516     1602     -1
+    1     511     512     536     537     514     539     -1
+    1     518     519     1610     1611     521     1613     -1
+    1     517     519     817     819     522     822     -1
+    1     517     518     1766     1767     520     1769     -1
+    1     524     525     1856     1857     527     1859     -1
+    1     523     525     1850     1849     528     1852     -1
+    1     523     524     775     776     526     778     -1
+    1     530     531     1591     1592     533     1594     -1
+    1     529     530     1598     1599     532     1601     -1
+    1     535     537     1603     1605     540     1608     -1
+    1     535     536     1826     1827     538     1829     -1
+    1     542     543     1615     1616     545     1618     -1
+    1     541     543     1435     1437     546     1440     -1
+    1     541     542     1621     1622     544     1624     -1
+    1     548     549     1681     1682     551     1684     -1
+    1     547     549     1220     1219     552     1222     -1
+    1     554     555     608     609     557     611     -1
+    1     553     555     1670     1669     558     1672     -1
+    1     553     554     1617     1615     556     1620     -1
+    1     560     561     1347     1345     563     1350     -1
+    1     559     561     1663     1665     564     1668     -1
+    1     565     566     1627     1628     568     1630     -1
+    1     571     573     1143     1142     576     1145     -1
+    1     571     572     1135     1136     574     1138     -1
+    1     583     585     590     589     588     592     -1
+    1     589     591     1639     1641     594     1644     -1
+    1     595     597     602     601     600     604     -1
+    1     602     603     2150     2151     605     2153     -1
+    1     601     603     1671     1670     606     1673     -1
+    1     613     615     1503     1502     618     1505     -1
+    1     620     621     1657     1658     623     1660     -1
+    1     625     627     1502     1501     630     1504     -1
+    1     625     626     1531     1532     628     1534     -1
+    1     632     633     1471     1472     635     1474     -1
+    1     631     633     1466     1465     636     1468     -1
+    1     638     639     1659     1657     641     1662     -1
+    1     637     638     1447     1448     640     1450     -1
+    1     644     645     1449     1447     647     1452     -1
+    1     643     645     1665     1664     648     1667     -1
+    1     650     651     1221     1219     653     1224     -1
+    1     656     657     1953     1951     659     1956     -1
+    1     655     657     1695     1694     660     1697     -1
+    1     655     656     1687     1688     658     1690     -1
+    1     662     663     1065     1063     665     1068     -1
+    1     661     663     764     763     666     766     -1
+    1     661     662     1959     1957     664     1962     -1
+    1     668     669     997     998     671     1000     -1
+    1     667     669     800     799     672     802     -1
+    1     674     675     974     975     677     977     -1
+    1     673     675     991     993     678     996     -1
+    1     673     674     1897     1898     676     1900     -1
+    1     679     681     1909     1911     684     1914     -1
+    1     679     680     1903     1904     682     1906     -1
+    1     686     687     991     992     689     994     -1
+    1     685     687     1911     1910     690     1913     -1
+    1     692     693     1563     1561     695     1566     -1
+    1     691     693     1549     1551     696     1554     -1
+    1     698     699     915     913     701     918     -1
+    1     697     699     1706     1705     702     1708     -1
+    1     697     698     1699     1700     700     1702     -1
+    1     704     705     1706     1707     707     1709     -1
+    1     703     705     1899     1898     708     1901     -1
+    1     703     704     1885     1886     706     1888     -1
+    1     709     711     1699     1701     714     1704     -1
+    1     709     710     722     723     712     725     -1
+    1     715     717     919     921     720     924     -1
+    1     715     716     913     914     718     916     -1
+    1     721     722     1711     1712     724     1714     -1
+    1     728     729     746     747     731     749     -1
+    1     740     741     896     897     743     899     -1
+    1     739     741     746     745     744     748     -1
+    1     739     740     1723     1724     742     1726     -1
+    1     745     747     1719     1718     750     1721     -1
+    1     752     753     902     903     755     905     -1
+    1     751     753     2169     2168     756     2171     -1
+    1     758     759     1815     1813     761     1818     -1
+    1     757     759     1807     1809     762     1812     -1
+    1     764     765     956     957     767     959     -1
+    1     763     765     1965     1964     768     1967     -1
+    1     770     771     2163     2161     773     2166     -1
+    1     769     771     871     873     774     876     -1
+    1     769     770     859     860     772     862     -1
+    1     776     777     1843     1844     779     1846     -1
+    1     775     777     1833     1832     780     1835     -1
+    1     782     783     2169     2167     785     2172     -1
+    1     781     783     908     907     786     910     -1
+    1     781     782     1807     1808     784     1810     -1
+    1     788     789     1731     1729     791     1734     -1
+    1     794     795     805     806     797     808     -1
+    1     793     795     1820     1819     798     1822     -1
+    1     793     794     853     854     796     856     -1
+    1     800     801     1899     1897     803     1902     -1
+    1     799     801     915     914     804     917     -1
+    1     806     807     1735     1736     809     1738     -1
+    1     812     813     1947     1945     815     1950     -1
+    1     811     813     973     975     816     978     -1
+    1     811     812     999     997     814     1002     -1
+    1     818     819     1754     1755     821     1757     -1
+    1     817     818     1737     1735     820     1740     -1
+    1     824     825     1609     1610     827     1612     -1
+    1     823     824     1693     1694     826     1696     -1
+    1     829     831     859     861     834     864     -1
+    1     829     830     1755     1753     832     1758     -1
+    1     836     837     1857     1855     839     1860     -1
+    1     835     837     1009     1011     840     1014     -1
+    1     835     836     1388     1389     838     1391     -1
+    1     841     843     867     866     846     869     -1
+    1     847     849     1737     1736     852     1739     -1
+    1     847     848     1729     1730     850     1732     -1
+    1     853     855     1731     1730     858     1733     -1
+    1     860     861     1010     1011     863     1013     -1
+    1     865     867     2168     2167     870     2170     -1
+    1     865     866     2161     2162     868     2164     -1
+    1     872     873     1389     1387     875     1392     -1
+    1     871     872     1778     1779     874     1781     -1
+    1     877     879     1777     1779     882     1782     -1
+    1     877     878     1813     1814     880     1816     -1
+    1     883     885     1797     1796     888     1799     -1
+    1     883     884     961     962     886     964     -1
+    1     890     891     931     932     893     934     -1
+    1     889     891     1856     1855     894     1858     -1
+    1     889     890     1875     1873     892     1878     -1
+    1     895     897     909     908     900     911     -1
+    1     895     896     901     902     898     904     -1
+    1     926     927     1015     1016     929     1018     -1
+    1     925     927     1010     1009     930     1012     -1
+    1     925     926     1831     1832     928     1834     -1
+    1     932     933     1862     1863     935     1865     -1
+    1     931     933     938     937     936     940     -1
+    1     938     939     2145     2143     941     2148     -1
+    1     937     939     1851     1850     942     1853     -1
+    1     944     945     1633     1634     947     1636     -1
+    1     943     944     1839     1837     946     1842     -1
+    1     949     951     1861     1863     954     1866     -1
+    1     949     950     1874     1875     952     1877     -1
+    1     955     957     1034     1033     960     1036     -1
+    1     955     956     1003     1004     958     1006     -1
+    1     962     963     1881     1879     965     1884     -1
+    1     968     969     2181     2179     971     2184     -1
+    1     967     969     1083     1082     972     1085     -1
+    1     973     974     1904     1905     976     1907     -1
+    1     979     981     1052     1051     984     1054     -1
+    1     979     980     2180     2181     982     2183     -1
+    1     986     987     1035     1033     989     1038     -1
+    1     985     987     1058     1057     990     1060     -1
+    1     985     986     1041     1039     988     1044     -1
+    1     992     993     1891     1892     995     1894     -1
+    1     998     999     1821     1819     1001     1824     -1
+    1     1003     1005     1124     1123     1008     1126     -1
+    1     1016     1017     1767     1765     1019     1770     -1
+    1     1015     1017     1754     1753     1020     1756     -1
+    1     1021     1023     1598     1597     1026     1600     -1
+    1     1028     1029     1844     1845     1031     1847     -1
+    1     1027     1029     1838     1837     1032     1840     -1
+    1     1034     1035     1052     1053     1037     1055     -1
+    1     1045     1047     1963     1965     1050     1968     -1
+    1     1045     1046     1125     1123     1048     1128     -1
+    1     1051     1053     1065     1064     1056     1067     -1
+    1     1070     1071     1940     1941     1073     1943     -1
+    1     1069     1070     2174     2175     1072     2177     -1
+    1     1075     1076     1951     1952     1078     1954     -1
+    1     1081     1082     1927     1928     1084     1930     -1
+    1     1088     1089     1905     1903     1091     1908     -1
+    1     1094     1095     1963     1964     1097     1966     -1
+    1     1093     1095     1869     1868     1098     1871     -1
+    1     1093     1094     1861     1862     1096     1864     -1
+    1     1100     1101     1279     1280     1103     1282     -1
+    1     1099     1101     1287     1286     1104     1289     -1
+    1     1099     1100     1369     1370     1102     1372     -1
+    1     1105     1107     1971     1970     1110     1973     -1
+    1     1105     1106     1322     1323     1108     1325     -1
+    1     1129     1130     1141     1142     1132     1144     -1
+    1     1135     1137     1585     1587     1140     1590     -1
+    1     1148     1149     1999     2000     1151     2002     -1
+    1     1147     1149     1994     1993     1152     1996     -1
+    1     1147     1148     2025     2023     1150     2028     -1
+    1     1154     1155     1159     1160     1157     1162     -1
+    1     1153     1154     2000     2001     1156     2003     -1
+    1     1160     1161     2024     2025     1163     2027     -1
+    1     1159     1161     1298     1297     1164     1300     -1
+    1     1165     1167     1183     1185     1170     1188     -1
+    1     1165     1166     1207     1208     1168     1210     -1
+    1     1171     1173     2017     2019     1176     2022     -1
+    1     1171     1172     1177     1178     1174     1180     -1
+    1     1178     1179     2012     2013     1181     2015     -1
+    1     1184     1185     2005     2006     1187     2008     -1
+    1     1183     1184     2018     2019     1186     2021     -1
+    1     1190     1191     1976     1977     1193     1979     -1
+    1     1189     1191     1285     1287     1194     1290     -1
+    1     1189     1190     1303     1304     1192     1306     -1
+    1     1196     1197     1213     1214     1199     1216     -1
+    1     1195     1197     2042     2041     1200     2044     -1
+    1     1195     1196     1375     1376     1198     1378     -1
+    1     1202     1203     2192     2193     1205     2195     -1
+    1     1201     1203     2186     2185     1206     2188     -1
+    1     1201     1202     1214     1215     1204     1217     -1
+    1     1208     1209     1424     1425     1211     1427     -1
+    1     1213     1215     2035     2037     1218     2040     -1
+    1     1220     1221     1353     1351     1223     1356     -1
+    1     1226     1227     2108     2109     1229     2111     -1
+    1     1225     1227     2102     2101     1230     2104     -1
+    1     1225     1226     1431     1429     1228     1434     -1
+    1     1231     1233     1267     1269     1236     1272     -1
+    1     1231     1232     2073     2071     1234     2076     -1
+    1     1238     1239     2186     2187     1241     2189     -1
+    1     1237     1238     1429     1430     1240     1432     -1
+    1     1244     1245     1419     1417     1247     1422     -1
+    1     1243     1244     2096     2097     1246     2099     -1
+    1     1250     1251     2054     2055     1253     2057     -1
+    1     1249     1251     2191     2193     1254     2196     -1
+    1     1256     1257     1412     1413     1259     1415     -1
+    1     1261     1263     2053     2055     1266     2058     -1
+    1     1261     1262     1418     1419     1264     1421     -1
+    1     1268     1269     1280     1281     1271     1283     -1
+    1     1273     1275     1412     1411     1278     1414     -1
+    1     1273     1274     2095     2096     1276     2098     -1
+    1     1279     1281     1291     1293     1284     1296     -1
+    1     1285     1286     2078     2079     1288     2081     -1
+    1     1292     1293     1304     1305     1295     1307     -1
+    1     1291     1292     1376     1377     1294     1379     -1
+    1     1298     1299     2029     2030     1301     2032     -1
+    1     1297     1299     1381     1383     1302     1386     -1
+    1     1303     1305     1382     1381     1308     1384     -1
+    1     1310     1311     2031     2029     1313     2034     -1
+    1     1309     1311     1317     1316     1314     1319     -1
+    1     1315     1317     1329     1328     1320     1331     -1
+    1     1327     1329     1340     1339     1332     1342     -1
+    1     1334     1335     1448     1449     1337     1451     -1
+    1     1333     1335     1442     1441     1338     1444     -1
+    1     1340     1341     2030     2031     1343     2033     -1
+    1     1339     1341     2024     2023     1344     2026     -1
+    1     1346     1347     1352     1353     1349     1355     -1
+    1     1345     1346     1995     1993     1348     1998     -1
+    1     1358     1359     2157     2155     1361     2160     -1
+    1     1357     1359     1664     1663     1362     1666     -1
+    1     1364     1365     1423     1424     1367     1426     -1
+    1     1363     1365     2108     2107     1368     2110     -1
+    1     1370     1371     2083     2084     1373     2086     -1
+    1     1369     1371     2072     2071     1374     2074     -1
+    1     1382     1383     2043     2041     1385     2046     -1
+    1     1387     1388     1789     1790     1390     1792     -1
+    1     1394     1395     1791     1789     1397     1794     -1
+    1     1393     1395     1881     1880     1398     1883     -1
+    1     1393     1394     1873     1874     1396     1876     -1
+    1     1400     1401     1820     1821     1403     1823     -1
+    1     1399     1400     1747     1748     1402     1750     -1
+    1     1406     1407     1749     1747     1409     1752     -1
+    1     1405     1407     1953     1952     1410     1955     -1
+    1     1405     1406     1945     1946     1408     1948     -1
+    1     1411     1413     2091     2090     1416     2093     -1
+    1     1417     1418     2089     2090     1420     2092     -1
+    1     1423     1425     2103     2102     1428     2105     -1
+    1     1436     1437     1683     1681     1439     1686     -1
+    1     1453     1455     1467     1466     1458     1469     -1
+    1     1465     1467     1652     1651     1470     1654     -1
+    1     1472     1473     2120     2121     1475     2123     -1
+    1     1471     1473     2114     2113     1476     2116     -1
+    1     1478     1479     1513     1514     1481     1516     -1
+    1     1477     1479     1485     1484     1482     1487     -1
+    1     1483     1485     2120     2119     1488     2122     -1
+    1     1489     1490     1507     1508     1492     1510     -1
+    1     1507     1509     1519     1521     1512     1524     -1
+    1     1513     1515     2115     2114     1518     2117     -1
+    1     1519     1520     1527     1525     1522     1530     -1
+    1     1556     1557     1569     1567     1559     1572     -1
+    1     1555     1556     2125     2126     1558     2128     -1
+    1     1562     1563     2127     2125     1565     2130     -1
+    1     1567     1568     1573     1574     1570     1576     -1
+    1     1591     1593     1617     1616     1596     1619     -1
+    1     1604     1605     1773     1771     1607     1776     -1
+    1     1603     1604     1676     1677     1606     1679     -1
+    1     1609     1611     1772     1771     1614     1774     -1
+    1     1622     1623     1682     1683     1625     1685     -1
+    1     1628     1629     2139     2137     1631     2142     -1
+    1     1634     1635     2144     2145     1637     2147     -1
+    1     1633     1635     2138     2137     1638     2140     -1
+    1     1652     1653     1658     1659     1655     1661     -1
+    1     1669     1671     2156     2155     1674     2158     -1
+    1     1693     1695     1742     1741     1698     1744     -1
+    1     1711     1713     1717     1719     1716     1722     -1
+    1     1742     1743     1748     1749     1745     1751     -1
+    1     1760     1761     1765     1766     1763     1768     -1
+    1     1759     1761     1831     1833     1764     1836     -1
+    1     1777     1778     1785     1783     1780     1788     -1
+    1     1784     1785     1790     1791     1787     1793     -1
+    1     1795     1797     1929     1928     1800     1931     -1
+    1     1795     1796     1921     1922     1798     1924     -1
+    1     1802     1803     1923     1921     1805     1926     -1
+    1     1801     1802     1879     1880     1804     1882     -1
+    1     1814     1815     2162     2163     1817     2165     -1
+    1     1825     1826     1845     1843     1828     1848     -1
+    1     1838     1839     1851     1849     1841     1854     -1
+    1     1867     1868     2143     2144     1870     2146     -1
+    1     1886     1887     1893     1891     1889     1896     -1
+    1     1916     1917     1922     1923     1919     1925     -1
+    1     1915     1916     1958     1959     1918     1961     -1
+    1     1933     1935     1939     1941     1938     1944     -1
+    1     1933     1934     2179     2180     1936     2182     -1
+    1     1969     1970     1987     1988     1972     1990     -1
+    1     1975     1976     1989     1987     1978     1992     -1
+    1     1981     1983     2085     2084     1986     2087     -1
+    1     1988     1989     2079     2077     1991     2082     -1
+    1     1999     2001     2018     2017     2004     2020     -1
+    1     2036     2037     2042     2043     2039     2045     -1
+    1     2048     2049     2191     2192     2051     2194     -1
+    1     2047     2049     2066     2065     2052     2068     -1
+    1     2053     2054     2059     2060     2056     2062     -1
+    1     2059     2061     2067     2066     2064     2069     -1
+    0     366     0
+    1     1     2     3     4     5     6
+    1     7     8     9     10     11     12
+    1     13     14     15     16     17     18
+    1     19     20     21     22     23     24
+    1     25     26     27     28     29     30
+    1     31     32     33     34     35     36
+    1     37     38     39     40     41     42
+    1     43     44     45     46     47     48
+    1     49     50     51     52     53     54
+    1     55     56     57     58     59     60
+    1     61     62     63     64     65     66
+    1     67     68     69     70     71     72
+    1     73     74     75     76     77     78
+    1     79     80     81     82     83     84
+    1     85     86     87     88     89     90
+    1     91     92     93     94     95     96
+    1     97     98     99     100     101     102
+    1     103     104     105     106     107     108
+    1     109     110     111     112     113     114
+    1     115     116     117     118     119     120
+    1     121     122     123     124     125     126
+    1     127     128     129     130     131     132
+    1     133     134     135     136     137     138
+    1     139     140     141     142     143     144
+    1     145     146     147     148     149     150
+    1     151     152     153     154     155     156
+    1     157     158     159     160     161     162
+    1     163     164     165     166     167     168
+    1     169     170     171     172     173     174
+    1     175     176     177     178     179     180
+    1     181     182     183     184     185     186
+    1     187     188     189     190     191     192
+    1     193     194     195     196     197     198
+    1     199     200     201     202     203     204
+    1     205     206     207     208     209     210
+    1     211     212     213     214     215     216
+    1     217     218     219     220     221     222
+    1     223     224     225     226     227     228
+    1     229     230     231     232     233     234
+    1     235     236     237     238     239     240
+    1     241     242     243     244     245     246
+    1     247     248     249     250     251     252
+    1     253     254     255     256     257     258
+    1     259     260     261     262     263     264
+    1     265     266     267     268     269     270
+    1     271     272     273     274     275     276
+    1     277     278     279     280     281     282
+    1     283     284     285     286     287     288
+    1     289     290     291     292     293     294
+    1     295     296     297     298     299     300
+    1     301     302     303     304     305     306
+    1     307     308     309     310     311     312
+    1     313     314     315     316     317     318
+    1     319     320     321     322     323     324
+    1     325     326     327     328     329     330
+    1     331     332     333     334     335     336
+    1     337     338     339     340     341     342
+    1     343     344     345     346     347     348
+    1     349     350     351     352     353     354
+    1     355     356     357     358     359     360
+    1     361     362     363     364     365     366
+    1     367     368     369     370     371     372
+    1     373     374     375     376     377     378
+    1     379     380     381     382     383     384
+    1     385     386     387     388     389     390
+    1     391     392     393     394     395     396
+    1     397     398     399     400     401     402
+    1     403     404     405     406     407     408
+    1     409     410     411     412     413     414
+    1     415     416     417     418     419     420
+    1     421     422     423     424     425     426
+    1     427     428     429     430     431     432
+    1     433     434     435     436     437     438
+    1     439     440     441     442     443     444
+    1     445     446     447     448     449     450
+    1     451     452     453     454     455     456
+    1     457     458     459     460     461     462
+    1     463     464     465     466     467     468
+    1     469     470     471     472     473     474
+    1     475     476     477     478     479     480
+    1     481     482     483     484     485     486
+    1     487     488     489     490     491     492
+    1     493     494     495     496     497     498
+    1     499     500     501     502     503     504
+    1     505     506     507     508     509     510
+    1     511     512     513     514     515     516
+    1     517     518     519     520     521     522
+    1     523     524     525     526     527     528
+    1     529     530     531     532     533     534
+    1     535     536     537     538     539     540
+    1     541     542     543     544     545     546
+    1     547     548     549     550     551     552
+    1     553     554     555     556     557     558
+    1     559     560     561     562     563     564
+    1     565     566     567     568     569     570
+    1     571     572     573     574     575     576
+    1     577     578     579     580     581     582
+    1     583     584     585     586     587     588
+    1     589     590     591     592     593     594
+    1     595     596     597     598     599     600
+    1     601     602     603     604     605     606
+    1     607     608     609     610     611     612
+    1     613     614     615     616     617     618
+    1     619     620     621     622     623     624
+    1     625     626     627     628     629     630
+    1     631     632     633     634     635     636
+    1     637     638     639     640     641     642
+    1     643     644     645     646     647     648
+    1     649     650     651     652     653     654
+    1     655     656     657     658     659     660
+    1     661     662     663     664     665     666
+    1     667     668     669     670     671     672
+    1     673     674     675     676     677     678
+    1     679     680     681     682     683     684
+    1     685     686     687     688     689     690
+    1     691     692     693     694     695     696
+    1     697     698     699     700     701     702
+    1     703     704     705     706     707     708
+    1     709     710     711     712     713     714
+    1     715     716     717     718     719     720
+    1     721     722     723     724     725     726
+    1     727     728     729     730     731     732
+    1     733     734     735     736     737     738
+    1     739     740     741     742     743     744
+    1     745     746     747     748     749     750
+    1     751     752     753     754     755     756
+    1     757     758     759     760     761     762
+    1     763     764     765     766     767     768
+    1     769     770     771     772     773     774
+    1     775     776     777     778     779     780
+    1     781     782     783     784     785     786
+    1     787     788     789     790     791     792
+    1     793     794     795     796     797     798
+    1     799     800     801     802     803     804
+    1     805     806     807     808     809     810
+    1     811     812     813     814     815     816
+    1     817     818     819     820     821     822
+    1     823     824     825     826     827     828
+    1     829     830     831     832     833     834
+    1     835     836     837     838     839     840
+    1     841     842     843     844     845     846
+    1     847     848     849     850     851     852
+    1     853     854     855     856     857     858
+    1     859     860     861     862     863     864
+    1     865     866     867     868     869     870
+    1     871     872     873     874     875     876
+    1     877     878     879     880     881     882
+    1     883     884     885     886     887     888
+    1     889     890     891     892     893     894
+    1     895     896     897     898     899     900
+    1     901     902     903     904     905     906
+    1     907     908     909     910     911     912
+    1     913     914     915     916     917     918
+    1     919     920     921     922     923     924
+    1     925     926     927     928     929     930
+    1     931     932     933     934     935     936
+    1     937     938     939     940     941     942
+    1     943     944     945     946     947     948
+    1     949     950     951     952     953     954
+    1     955     956     957     958     959     960
+    1     961     962     963     964     965     966
+    1     967     968     969     970     971     972
+    1     973     974     975     976     977     978
+    1     979     980     981     982     983     984
+    1     985     986     987     988     989     990
+    1     991     992     993     994     995     996
+    1     997     998     999     1000     1001     1002
+    1     1003     1004     1005     1006     1007     1008
+    1     1009     1010     1011     1012     1013     1014
+    1     1015     1016     1017     1018     1019     1020
+    1     1021     1022     1023     1024     1025     1026
+    1     1027     1028     1029     1030     1031     1032
+    1     1033     1034     1035     1036     1037     1038
+    1     1039     1040     1041     1042     1043     1044
+    1     1045     1046     1047     1048     1049     1050
+    1     1051     1052     1053     1054     1055     1056
+    1     1057     1058     1059     1060     1061     1062
+    1     1063     1064     1065     1066     1067     1068
+    1     1069     1070     1071     1072     1073     1074
+    1     1075     1076     1077     1078     1079     1080
+    1     1081     1082     1083     1084     1085     1086
+    1     1087     1088     1089     1090     1091     1092
+    1     1093     1094     1095     1096     1097     1098
+    1     1099     1100     1101     1102     1103     1104
+    1     1105     1106     1107     1108     1109     1110
+    1     1111     1112     1113     1114     1115     1116
+    1     1117     1118     1119     1120     1121     1122
+    1     1123     1124     1125     1126     1127     1128
+    1     1129     1130     1131     1132     1133     1134
+    1     1135     1136     1137     1138     1139     1140
+    1     1141     1142     1143     1144     1145     1146
+    1     1147     1148     1149     1150     1151     1152
+    1     1153     1154     1155     1156     1157     1158
+    1     1159     1160     1161     1162     1163     1164
+    1     1165     1166     1167     1168     1169     1170
+    1     1171     1172     1173     1174     1175     1176
+    1     1177     1178     1179     1180     1181     1182
+    1     1183     1184     1185     1186     1187     1188
+    1     1189     1190     1191     1192     1193     1194
+    1     1195     1196     1197     1198     1199     1200
+    1     1201     1202     1203     1204     1205     1206
+    1     1207     1208     1209     1210     1211     1212
+    1     1213     1214     1215     1216     1217     1218
+    1     1219     1220     1221     1222     1223     1224
+    1     1225     1226     1227     1228     1229     1230
+    1     1231     1232     1233     1234     1235     1236
+    1     1237     1238     1239     1240     1241     1242
+    1     1243     1244     1245     1246     1247     1248
+    1     1249     1250     1251     1252     1253     1254
+    1     1255     1256     1257     1258     1259     1260
+    1     1261     1262     1263     1264     1265     1266
+    1     1267     1268     1269     1270     1271     1272
+    1     1273     1274     1275     1276     1277     1278
+    1     1279     1280     1281     1282     1283     1284
+    1     1285     1286     1287     1288     1289     1290
+    1     1291     1292     1293     1294     1295     1296
+    1     1297     1298     1299     1300     1301     1302
+    1     1303     1304     1305     1306     1307     1308
+    1     1309     1310     1311     1312     1313     1314
+    1     1315     1316     1317     1318     1319     1320
+    1     1321     1322     1323     1324     1325     1326
+    1     1327     1328     1329     1330     1331     1332
+    1     1333     1334     1335     1336     1337     1338
+    1     1339     1340     1341     1342     1343     1344
+    1     1345     1346     1347     1348     1349     1350
+    1     1351     1352     1353     1354     1355     1356
+    1     1357     1358     1359     1360     1361     1362
+    1     1363     1364     1365     1366     1367     1368
+    1     1369     1370     1371     1372     1373     1374
+    1     1375     1376     1377     1378     1379     1380
+    1     1381     1382     1383     1384     1385     1386
+    1     1387     1388     1389     1390     1391     1392
+    1     1393     1394     1395     1396     1397     1398
+    1     1399     1400     1401     1402     1403     1404
+    1     1405     1406     1407     1408     1409     1410
+    1     1411     1412     1413     1414     1415     1416
+    1     1417     1418     1419     1420     1421     1422
+    1     1423     1424     1425     1426     1427     1428
+    1     1429     1430     1431     1432     1433     1434
+    1     1435     1436     1437     1438     1439     1440
+    1     1441     1442     1443     1444     1445     1446
+    1     1447     1448     1449     1450     1451     1452
+    1     1453     1454     1455     1456     1457     1458
+    1     1459     1460     1461     1462     1463     1464
+    1     1465     1466     1467     1468     1469     1470
+    1     1471     1472     1473     1474     1475     1476
+    1     1477     1478     1479     1480     1481     1482
+    1     1483     1484     1485     1486     1487     1488
+    1     1489     1490     1491     1492     1493     1494
+    1     1495     1496     1497     1498     1499     1500
+    1     1501     1502     1503     1504     1505     1506
+    1     1507     1508     1509     1510     1511     1512
+    1     1513     1514     1515     1516     1517     1518
+    1     1519     1520     1521     1522     1523     1524
+    1     1525     1526     1527     1528     1529     1530
+    1     1531     1532     1533     1534     1535     1536
+    1     1537     1538     1539     1540     1541     1542
+    1     1543     1544     1545     1546     1547     1548
+    1     1549     1550     1551     1552     1553     1554
+    1     1555     1556     1557     1558     1559     1560
+    1     1561     1562     1563     1564     1565     1566
+    1     1567     1568     1569     1570     1571     1572
+    1     1573     1574     1575     1576     1577     1578
+    1     1579     1580     1581     1582     1583     1584
+    1     1585     1586     1587     1588     1589     1590
+    1     1591     1592     1593     1594     1595     1596
+    1     1597     1598     1599     1600     1601     1602
+    1     1603     1604     1605     1606     1607     1608
+    1     1609     1610     1611     1612     1613     1614
+    1     1615     1616     1617     1618     1619     1620
+    1     1621     1622     1623     1624     1625     1626
+    1     1627     1628     1629     1630     1631     1632
+    1     1633     1634     1635     1636     1637     1638
+    1     1639     1640     1641     1642     1643     1644
+    1     1645     1646     1647     1648     1649     1650
+    1     1651     1652     1653     1654     1655     1656
+    1     1657     1658     1659     1660     1661     1662
+    1     1663     1664     1665     1666     1667     1668
+    1     1669     1670     1671     1672     1673     1674
+    1     1675     1676     1677     1678     1679     1680
+    1     1681     1682     1683     1684     1685     1686
+    1     1687     1688     1689     1690     1691     1692
+    1     1693     1694     1695     1696     1697     1698
+    1     1699     1700     1701     1702     1703     1704
+    1     1705     1706     1707     1708     1709     1710
+    1     1711     1712     1713     1714     1715     1716
+    1     1717     1718     1719     1720     1721     1722
+    1     1723     1724     1725     1726     1727     1728
+    1     1729     1730     1731     1732     1733     1734
+    1     1735     1736     1737     1738     1739     1740
+    1     1741     1742     1743     1744     1745     1746
+    1     1747     1748     1749     1750     1751     1752
+    1     1753     1754     1755     1756     1757     1758
+    1     1759     1760     1761     1762     1763     1764
+    1     1765     1766     1767     1768     1769     1770
+    1     1771     1772     1773     1774     1775     1776
+    1     1777     1778     1779     1780     1781     1782
+    1     1783     1784     1785     1786     1787     1788
+    1     1789     1790     1791     1792     1793     1794
+    1     1795     1796     1797     1798     1799     1800
+    1     1801     1802     1803     1804     1805     1806
+    1     1807     1808     1809     1810     1811     1812
+    1     1813     1814     1815     1816     1817     1818
+    1     1819     1820     1821     1822     1823     1824
+    1     1825     1826     1827     1828     1829     1830
+    1     1831     1832     1833     1834     1835     1836
+    1     1837     1838     1839     1840     1841     1842
+    1     1843     1844     1845     1846     1847     1848
+    1     1849     1850     1851     1852     1853     1854
+    1     1855     1856     1857     1858     1859     1860
+    1     1861     1862     1863     1864     1865     1866
+    1     1867     1868     1869     1870     1871     1872
+    1     1873     1874     1875     1876     1877     1878
+    1     1879     1880     1881     1882     1883     1884
+    1     1885     1886     1887     1888     1889     1890
+    1     1891     1892     1893     1894     1895     1896
+    1     1897     1898     1899     1900     1901     1902
+    1     1903     1904     1905     1906     1907     1908
+    1     1909     1910     1911     1912     1913     1914
+    1     1915     1916     1917     1918     1919     1920
+    1     1921     1922     1923     1924     1925     1926
+    1     1927     1928     1929     1930     1931     1932
+    1     1933     1934     1935     1936     1937     1938
+    1     1939     1940     1941     1942     1943     1944
+    1     1945     1946     1947     1948     1949     1950
+    1     1951     1952     1953     1954     1955     1956
+    1     1957     1958     1959     1960     1961     1962
+    1     1963     1964     1965     1966     1967     1968
+    1     1969     1970     1971     1972     1973     1974
+    1     1975     1976     1977     1978     1979     1980
+    1     1981     1982     1983     1984     1985     1986
+    1     1987     1988     1989     1990     1991     1992
+    1     1993     1994     1995     1996     1997     1998
+    1     1999     2000     2001     2002     2003     2004
+    1     2005     2006     2007     2008     2009     2010
+    1     2011     2012     2013     2014     2015     2016
+    1     2017     2018     2019     2020     2021     2022
+    1     2023     2024     2025     2026     2027     2028
+    1     2029     2030     2031     2032     2033     2034
+    1     2035     2036     2037     2038     2039     2040
+    1     2041     2042     2043     2044     2045     2046
+    1     2047     2048     2049     2050     2051     2052
+    1     2053     2054     2055     2056     2057     2058
+    1     2059     2060     2061     2062     2063     2064
+    1     2065     2066     2067     2068     2069     2070
+    1     2071     2072     2073     2074     2075     2076
+    1     2077     2078     2079     2080     2081     2082
+    1     2083     2084     2085     2086     2087     2088
+    1     2089     2090     2091     2092     2093     2094
+    1     2095     2096     2097     2098     2099     2100
+    1     2101     2102     2103     2104     2105     2106
+    1     2107     2108     2109     2110     2111     2112
+    1     2113     2114     2115     2116     2117     2118
+    1     2119     2120     2121     2122     2123     2124
+    1     2125     2126     2127     2128     2129     2130
+    1     2131     2132     2133     2134     2135     2136
+    1     2137     2138     2139     2140     2141     2142
+    1     2143     2144     2145     2146     2147     2148
+    1     2149     2150     2151     2152     2153     2154
+    1     2155     2156     2157     2158     2159     2160
+    1     2161     2162     2163     2164     2165     2166
+    1     2167     2168     2169     2170     2171     2172
+    1     2173     2174     2175     2176     2177     2178
+    1     2179     2180     2181     2182     2183     2184
+    1     2185     2186     2187     2188     2189     2190
+    1     2191     2192     2193     2194     2195     2196
+    0
+    0
+    0
+    0
+   
\ No newline at end of file
diff --git a/examples/fem/crack2D/driver.C b/examples/fem/crack2D/driver.C
new file mode 100644 (file)
index 0000000..74ea743
--- /dev/null
@@ -0,0 +1,102 @@
+#include <fstream.h>
+#include "crack.h"
+
+extern "C" void
+init(void)
+{
+}
+
+static void
+initNodes(GlobalData *gd)
+{
+  int idx;
+  for(idx=0; idx<gd->nn; idx++) {
+    Node *np = &(gd->nodes[idx]);
+    // Zero out node accumulators, update node positions
+    np->Rco.y = np->Rco.x = np->Rin.y =  np->Rin.x = 0.0;
+    np->disp.x += gd->delta2*np->accel.x + gd->delta*np->vel.x;
+    np->disp.y += gd->delta2*np->accel.y + gd->delta*np->vel.y;
+  }
+}
+
+static void 
+updateNodes(GlobalData *gd, double prop, double slope)
+{
+  for(int idx=0; idx<gd->nn; idx++) {
+    Node *np = &(gd->nodes[idx]);
+    if(!np->isbnd) {
+      double aX, aY;
+      aX = (np->Rco.x-np->Rin.x)*np->xM.x;
+      aY = (np->Rco.y-np->Rin.y)*np->xM.y;
+      np->vel.x += (gd->delta*(np->accel.x+aX)*(double) 0.5);
+      np->vel.y += (gd->delta*(np->accel.y+aY)*(double) 0.5);
+      np->accel.x = aX;
+      np->accel.y = aY;
+    } else {
+      double acc;
+      if (!(np->id1)) {
+        np->vel.x = (np->r.x)*prop;
+        np->accel.x = (np->r.x)*slope;
+      } else {
+        acc = (np->r.x*prop+ np->Rco.x - np->Rin.x)*np->xM.x;
+        np->vel.x += (gd->delta*(np->accel.x+acc)*0.5);
+        np->accel.x = acc;
+      }
+      if (!(np->id2)) {
+        np->vel.y = np->r.y*prop;
+        np->accel.y = np->r.y*slope;
+      } else {
+        acc = (np->r.y*prop+ np->Rco.y - np->Rin.y)*np->xM.y;
+        np->vel.y = (np->vel.y + gd->delta*(np->accel.y+acc)*0.5);
+        np->accel.y = acc;
+      }
+    }
+  }
+}
+
+extern "C" void
+driver(int nn, int *nnums, int ne, int *enums, int npere, int *conn)
+{
+  int myid = FEM_My_Partition();
+  GlobalData *gd = new GlobalData;
+  Node *nodes = new Node[nn];
+  Element *elements = new Element[ne];
+  gd->nn = nn;
+  gd->nnums = nnums;
+  gd->ne = ne;
+  gd->enums = enums;
+  gd->npere = npere;
+  gd->conn = conn;
+  gd->nodes = nodes;
+  gd->elements = elements;
+  readFile(gd);
+  int i;
+  int kk = -1;
+  double prop, slope;
+  for(i=0;i<gd->nTime;i++)
+  {
+    CkPrintf("[%d] Beginning iteration %d\n", myid, i);
+    if (gd->ts_proportion[kk+1] == i)
+    {
+      kk++;
+      prop = gd->proportion[kk];
+      slope = (gd->proportion[kk+1]-prop)/gd->delta;
+      slope /= (double) (gd->ts_proportion[kk+1]- gd->ts_proportion[kk]);
+    }
+    else
+    {
+      prop = (double)(i - gd->ts_proportion[kk])*
+                    slope*gd->delta+gd->proportion[kk];
+    }
+    initNodes(gd);
+    lst_NL(gd);
+    lst_coh2(gd);
+    updateNodes(gd, prop, slope);
+    // FIXME: fem_update here.
+  }
+}
+
+extern "C" void
+finalize()
+{
+}
diff --git a/examples/fem/crack2D/getmesh.C b/examples/fem/crack2D/getmesh.C
new file mode 100644 (file)
index 0000000..5a8651f
--- /dev/null
@@ -0,0 +1,78 @@
+#include <fstream.h>
+
+int 
+main(int, char**)
+{
+  ifstream fi;
+  ofstream fo;
+  fi.open("crck_bar.inp");
+  if (!fi)
+  {
+    fprintf(stderr, "Cannot open crck_bar.inp for reading\n");
+    return 1;
+  }
+  fo.open("crck_bar.mesh");
+  if (!fo)
+  {
+    fprintf(stderr, "Cannot open crck_bar.mesh for writing\n");
+    return 1;
+  }
+  int itmp, i, j, k;
+  double dtmp;
+  fi >> itmp >> itmp;
+  fi >> dtmp;
+  int nump;
+  fi >> nump; // number of props
+  for (i=0;i<nump;i++)
+    fi >> itmp >> dtmp;
+  int numnp;
+  fi >> numnp; // number of nodes
+  for (i=0;i<numnp;i++)
+    fi >> itmp >> dtmp >> dtmp;
+  int numb; // num boundary nodes
+  fi >> numb;
+  for(i=0;i<numb;i++)
+    fi >> itmp >> itmp >> itmp >> dtmp >> dtmp;
+  int numclst; // number of cohesive elements
+  fi >> itmp >> numclst >> itmp >> itmp >> itmp;
+  int *cnodes = new int[numclst*6];
+  k = 0;
+  for (i=0; i< numclst; i++)
+  {
+    fi >> itmp;
+    for(j=0;j<6;j++)
+      fi >> cnodes[k++];
+    fi >> itmp;
+  }
+  int numlst; // number of vol elements
+  fi >> itmp >> numlst >> itmp;
+  int *vnodes = new int[numlst*6];
+  k = 0;
+  for (i=0; i< numlst; i++)
+  {
+    fi >> itmp;
+    for(j=0;j<6;j++)
+      fi >> vnodes[k++];
+  }
+  fo << numlst+numclst << ' ' << numnp << ' ' << 6 << endl;
+  k = 0;
+  for(i=0;i<numclst;i++)
+  {
+    for(j=0;j<6;j++)
+      fo << cnodes[k++]-1 << ' ';
+    fo << endl;
+  }
+  k = 0;
+  for(i=0;i<numlst;i++)
+  {
+    for(j=0;j<6;j++)
+      fo << vnodes[k++]-1 << ' ';
+    fo << endl;
+  }
+  cout << "Total Nodes: " << numnp << endl;
+  cout << "Total Cohesive Elements: " << numclst << endl;
+  cout << "Total Volumetric Elements: " << numlst << endl;
+  fi.close();
+  fo.close();
+  return 0;
+}
diff --git a/examples/fem/crack2D/lst_NL.C b/examples/fem/crack2D/lst_NL.C
new file mode 100644 (file)
index 0000000..2a6047d
--- /dev/null
@@ -0,0 +1,150 @@
+/*Subroutine: 
+  lst_NL:  converted from Fortran on 9/7/99 by Orion Sky Lawlor
+           Computes each node's Rin.
+*/
+
+#include "crack.h"
+
+//BtoR: sum the "B" spatial derivative array into the
+// "R" partial internal force array.
+
+static void BtoR(const Coord *B/*6-array*/,
+    Coord *R/*6-array*/,Node **n/*6-array*/,
+    Element *v,const VolMaterial *vm,
+    int pkCount,double iaa)
+
+{
+  double dudx,dvdy,dudy,dvdx;      //    ! partial derivatives of disp
+  double E11,E22,E12;              //            ! strains
+  int k;                         //Loop index
+
+//C-----Calculate displacement gradient (H)
+  dudx=dvdy=dudy=dvdx=0.0;
+  for (k=0;k<6;k++) {
+    dudx+=B[k].x*n[k]->disp.x *iaa;
+    dvdy+=B[k].y*n[k]->disp.y *iaa;
+    dudy+=B[k].y*n[k]->disp.x *iaa;
+    dvdx+=B[k].x*n[k]->disp.y *iaa;
+  }
+
+//C-----Calculate Lagrangian strain (E)
+  
+  E11 = dudx + 0.5*(dudx*dudx + dvdx*dvdx);
+  E22 = dvdy + 0.5*(dvdy*dvdy + dudy*dudy);
+  E12 = dudy + dvdx + dudx*dudy + dvdy*dvdx;
+  
+//C-----Calculate the pkCount'th Piola-Kirchhoff stress (S)
+  
+  double s11c,s12c,s22c;
+  s11c=v->v.s11l[pkCount] = E11*vm->c[0] + E22*vm->c[1];
+  s22c=v->v.s22l[pkCount] = E11*vm->c[1] + E22*vm->c[2];
+  s12c=v->v.s12l[pkCount] = E12*vm->c[3];
+  
+//Update R
+  for (k=0;k<6;k++) {
+    R[k].x+=s11c*B[k].x*(1.0+dudx)  +  s22c*B[k].y* dudy+
+      s12c*(B[k].y*(1.0+dudx) + B[k].x*dudy);
+    R[k].y+=s11c*B[k].x*dvdx        +  s22c*B[k].y*(1.0+dvdy)+
+      s12c*(B[k].y*     dvdx  + B[k].x*(1.0+dvdy));
+  }
+
+}
+
+
+void
+lst_NL(GlobalData *gd)
+{
+  int idx;
+  // skip over cohesive elements
+  for(idx=0; idx<gd->ne && gd->nnums[idx]<gd->numCLST; idx++);
+  for(; idx<gd->ne; idx++) {
+    Element *v = &(gd->elements[idx]);
+    Node *n[6];
+    int k;
+    for(k=0;k<6;k++)
+      n[k] = &(gd->nodes[gd->conn[idx*6+k]]);
+
+    const double c5v3 = 1.66666666666667;
+    const double c1v3 = 0.33333333333333;
+    const double c2v3 = 0.66666666666667;
+    const double c1v6 = 0.166666666666667;
+  
+    double x21,y21,x31,y31,x32,y32;         //! coor(1,n2)-coor(1,n1)
+    double iaa;                       //! inverse of twice the element area
+      
+    Coord B[6];                           //! spacial derivatives
+    Coord R[6];                           //! partial internal force
+      
+    //Nodes 0, 1, and 2 are the outer corners;
+    //Nodes 4, 5, and 6 are on the midpoints of the sides.
+      
+    VolMaterial *vm = &(gd->volm[v->material]); //Current material
+      
+    for (k=0;k<6;k++) {
+        R[k].x=R[k].y=0.0;
+    }
+      
+    x21 = n[1]->pos.x-n[0]->pos.x;
+    y21 = n[1]->pos.y-n[0]->pos.y;
+    x31 = n[2]->pos.x-n[0]->pos.x;
+    y31 = n[2]->pos.y-n[0]->pos.y;
+    x32 = n[2]->pos.x-n[1]->pos.x;
+    y32 = n[2]->pos.y-n[1]->pos.y;
+    iaa = 1.0/(x21*y31-x31*y21);
+      
+    //Perform the three stages of Piola-Kirchoff stress
+    //First stage:
+    B[0].x = -y32*c5v3;             //                !   B1 = aa*dN1/dx
+    B[0].y = x32*c5v3;              //                !   B2 = aa*dN1/dy
+    B[1].x = -y31*c1v3;             //                !   B3 = aa*dN2/dx
+    B[1].y = x31*c1v3;              //                !   B4 = aa*dN2/dy
+    B[2].x = y21*c1v3;              //                !   B5 = aa*dN3/dx
+    B[2].y = -x21*c1v3;             //                !   B6 = aa*dN3/dy
+    B[3].x = (4.0*y31 - y32)*c2v3;  //                !   B7 = aa*dN4/dx
+    B[3].y = (x32 - 4.0*x31)*c2v3;  //                !   B8 = aa*dN4/dy
+    B[4].x = (y31 - y21)*c2v3;      //                !   B9 = aa*dN5/dx
+    B[4].y = (x21 - x31)*c2v3;      //                !   B10 = aa*dN5/dy
+    B[5].x = -(y32 + 4.0*y21)*c2v3; //                !   B11 = aa*dN6/dx
+    B[5].y = (x32 + 4.0*x21)*c2v3;  //                !   B12 = aa*dN6/dy
+          
+    BtoR(B,R,n,v,vm,0,iaa);
+          
+    //Second stage:
+    B[0].x = y32*c1v3;               //   !   B1 = aa*dN1/dx
+    B[0].y = -x32*c1v3;              //   !   B2 = aa*dN1/dy
+    B[1].x = y31*c5v3;               //   !   B3 = aa*dN2/dx
+    B[1].y = -x31*c5v3;              //   !   B4 = aa*dN2/dy
+    B[2].x = y21*c1v3;               //   !   B5 = aa*dN3/dx
+    B[2].y = -x21*c1v3;              //   !   B6 = aa*dN3/dy
+    B[3].x = (y31 - 4.0*y32)*c2v3;   //   !   B7 = aa*dN4/dx
+    B[3].y = (4.0*x32 - x31)*c2v3;   //   !   B8 = aa*dN4/dy
+    B[4].x = (y31 - 4.0*y21)*c2v3;   //   !   B9 = aa*dN5/dx
+    B[4].y = (4.0*x21 - x31)*c2v3;   //   !   B10 = aa*dN5/dy
+    B[5].x = -(y32 + y21)*c2v3;      //   !   B11 = aa*dN6/dx
+    B[5].y = (x32 + x21)*c2v3;       //   !   B12 = aa*dN6/dy
+    
+    BtoR(B,R,n,v,vm,1,iaa);
+    
+    //Third stage:
+    B[0].x = y32*c1v3;               //  !   B1 = aa*dN1/dx
+    B[0].y = -x32*c1v3;              //  !   B2 = aa*dN1/dy
+    B[1].x = -y31*c1v3;              //  !   B3 = aa*dN2/dx
+    B[1].y = x31*c1v3;               //  !   B4 = aa*dN2/dy
+    B[2].x = -y21*c5v3;              //  !   B5 = aa*dN3/dx
+    B[2].y = x21*c5v3;               //  !   B6 = aa*dN3/dy
+    B[3].x = (y31 - y32)*c2v3;       //  !   B7 = aa*dN4/dx
+    B[3].y = (x32 - x31)*c2v3;       //  !   B8 = aa*dN4/dy
+    B[4].x = (4.0*y31 - y21)*c2v3;   //  !   B9 = aa*dN5/dx
+    B[4].y = (x21 - 4.0*x31)*c2v3;   //  !   B10 = aa*dN5/dy
+    B[5].x = -(4.0*y32 + y21)*c2v3;  //  !   B11 = aa*dN6/dx
+    B[5].y = (4.0*x32 + x21)*c2v3;   //  !   B12 = aa*dN6/dy
+    
+    BtoR(B,R,n,v,vm,2,iaa);
+    
+    //Now we've calculated R.  Update each node`s Rin with the new values.
+    for (k=0;k<6;k++) {
+      n[k]->Rin.x+=c1v6*R[k].x;
+      n[k]->Rin.y+=c1v6*R[k].y;
+    }
+  }
+}
diff --git a/examples/fem/crack2D/lst_coh2.C b/examples/fem/crack2D/lst_coh2.C
new file mode 100644 (file)
index 0000000..73f55ea
--- /dev/null
@@ -0,0 +1,109 @@
+#include "crack.h"
+#include <math.h>
+
+
+//Update cohesive element traction itNo given Dn, Dt
+static void 
+updateTraction(double Dn, double Dt,Element *coh,CohMaterial *m,int itNo)
+{
+  double Tn, Tt;
+  double c=coh->c.sidel[1],s=coh->c.sidel[2];
+  double x = 1.0 - sqrt(Dn*Dn+Dt*Dt);
+  
+  if (x <= 0.0) {
+    coh->c.Sthresh[itNo] = 0.0;
+  } else if (x <= coh->c.Sthresh[itNo])
+    coh->c.Sthresh[itNo] = x;
+  
+  if (Dn > 0.0)
+    Tn = coh->c.Sthresh[itNo]/(1.0-coh->c.Sthresh[itNo])*m->sigmax*Dn;
+  else
+    Tn = m->Sinit/(1.0-m->Sinit)*m->sigmax*Dn;
+  
+  Tt = coh->c.Sthresh[itNo]/(1.0-coh->c.Sthresh[itNo])*m->taumax*Dt;
+  coh->c.T[itNo].x = c*Tt - s*Tn;
+  coh->c.T[itNo].y = s*Tt + c*Tn;
+}       
+
+void
+lst_coh2(GlobalData *gd) 
+{
+  int idx;
+  for(idx=0; idx<gd->ne && gd->enums[idx]<gd->numCLST; idx++) {
+    Element *coh = &(gd->elements[idx]);
+    Node *n[6];
+    int k;
+    for(k=0;k<6;k++)
+      n[k] = &(gd->nodes[gd->conn[idx*6+k]]);
+
+    // local variables:
+    double deltn, deltt;  // the char length of current element
+    double length, c, s;
+    double Dx1, Dx2, Dx3;
+    double Dy1, Dy2, Dy3;
+    double Dn1, Dn2, Dn3;
+    double Dt1, Dt2, Dt3;
+    double Dn, Dt;
+    double x;
+    double Rx1,Rx2,Rx3;  // cohesive forces at nodes 1,2,5
+    double Ry1,Ry2,Ry3;  // cohesive forces at nodes 1,2,5
+  
+    // g1,g3: gauss quadrature points
+    // w1,w2,w3: weights; w1=w2, w3=w4
+    
+    CohMaterial *m = &(gd->cohm[coh->material]);
+      
+    deltn = (double)1.0 / m->deltan;
+    deltt = (double)1.0 / m->deltat;
+    length = coh->c.sidel[0] * (double)0.5;
+    c = coh->c.sidel[1];
+    s = coh->c.sidel[2];
+    
+    Dx1 = n[3]->disp.x - n[0]->disp.x;
+    Dy1 = n[3]->disp.y - n[0]->disp.y;
+    Dt1 =  (c*Dx1 + s*Dy1)*deltt;
+    Dn1 = (-s*Dx1 + c*Dy1)*deltn;
+    Dx2 = n[2]->disp.x - n[1]->disp.x;
+    Dy2 = n[2]->disp.y - n[1]->disp.y;
+    Dt2 =  (c*Dx2 + s*Dy2)*deltt;
+    Dn2 = (-s*Dx2 + c*Dy2)*deltn;
+    Dx3 = n[5]->disp.x - n[4]->disp.x;
+    Dy3 = n[5]->disp.y - n[4]->disp.y;
+    Dt3 =  (c*Dx3 + s*Dy3)*deltt;
+    Dn3 = (-s*Dx3 + c*Dy3)*deltn;
+      
+    // gauss point 1
+    Dt = g1*g1*0.5*(Dt1+Dt2-2.0*Dt3)+g1*0.5*(Dt2-Dt1)+Dt3;
+    Dn = g1*g1*0.5*(Dn1+Dn2-2.0*Dn3)+g1*0.5*(Dn2-Dn1)+Dn3;
+    updateTraction(Dn,Dt,coh,m,0);
+      
+    // gauss point 2
+    updateTraction(Dn3,Dt3,coh,m,1);
+  
+    // gauss point 3
+    Dt = g3*g3*0.5*(Dt1+Dt2-2.0*Dt3)+g3*0.5*(Dt2-Dt1)+Dt3;
+    Dn = g3*g3*0.5*(Dn1+Dn2-2.0*Dn3)+g3*0.5*(Dn2-Dn1)+Dn3;
+    updateTraction(Dn,Dt,coh,m,2);
+  
+    // cohesive forces
+    x = length*w1*g1*0.5;
+    Rx1 = (coh->c.T[0].x*(g1-1.0)+coh->c.T[2].x*(g1+1.0))*x;
+    Ry1 = (coh->c.T[0].y*(g1-1.0)+coh->c.T[2].y*(g1+1.0))*x;
+    Rx2 = (coh->c.T[0].x*(g1+1.0)+coh->c.T[2].x*(g1-1.0))*x;
+    Ry2 = (coh->c.T[0].y*(g1+1.0)+coh->c.T[2].y*(g1-1.0))*x;
+    Rx3 = ((coh->c.T[0].x+coh->c.T[2].x)*w1*(1.0-g1*g1)+coh->c.T[1].x*w2)*length;
+    Ry3 = ((coh->c.T[0].y+coh->c.T[2].y)*w1*(1.0-g1*g1)+coh->c.T[1].y*w2)*length;
+    n[0]->Rco.x += Rx1;
+    n[0]->Rco.y += Ry1;
+    n[3]->Rco.x -= Rx1;
+    n[3]->Rco.y -= Ry1;
+    n[1]->Rco.x += Rx2;
+    n[1]->Rco.y += Ry2;
+    n[2]->Rco.x -= Rx2;
+    n[2]->Rco.y -= Ry2;
+    n[4]->Rco.x += Rx3;
+    n[4]->Rco.y += Ry3;
+    n[5]->Rco.x -= Rx3;
+    n[5]->Rco.y -= Ry3;
+  }
+}
diff --git a/examples/fem/crack2D/readfile.C b/examples/fem/crack2D/readfile.C
new file mode 100644 (file)
index 0000000..77ef202
--- /dev/null
@@ -0,0 +1,167 @@
+#include "crack.h"
+#include <fstream.h>
+
+#define MAXLINE 1024
+#define cl(fd, buffer, n) do {\
+                               for(int i=0; i<n; i++) \
+                                 fd.getline(buffer, MAXLINE);\
+                             } while(0)
+void
+readFile(GlobalData *gd)
+{
+  ifstream fg; // global parameter file
+  ifstream fm; // mesh file
+  char buf[MAXLINE];
+  int i, itmp;
+  double dtmp;
+
+  fg.open("cohesive.inp");
+  fm.open("crck_bar.inp");
+  if (!fg || !fm)
+  {
+    CkAbort("Cannot open input files for reading.\n");
+    return;
+  }
+  cl(fg,buf,3); // ignore first three lines
+  fg >> gd->nTime >> gd->steps >> gd->tsintFull 
+     >> gd->tsintEnergy >> gd->restart;
+  cl(fg, buf, 2);
+  fg >> gd->nplane >> gd->ncoh >> gd->lin;
+  cl(fg, buf, 2);
+  //read volumetric material properties
+  fg >> gd->numMatVol;
+  gd->volm = new VolMaterial[gd->numMatVol];
+  VolMaterial *v;
+  cl(fg, buf, 2);
+  for (i=0; i<gd->numMatVol; i++)
+  {
+    v = &(gd->volm[i]);
+    fg >> v->e1 >> v->e2 >> v->g12 >> v->g23;
+    fg >> v->xnu12 >> v->xnu23 >> v->rho;
+    fg >> v->alpha1 >> v->alpha2 >> v->theta;
+    cl(fg, buf, 2);
+  }
+  //read cohesive material properties
+  fg >> gd->numMatCoh;
+  gd->cohm = new CohMaterial[gd->numMatCoh];
+  CohMaterial *c;
+  cl(fg, buf, 2);
+  for (i=0; i<gd->numMatCoh; i++)
+  {
+    c = &(gd->cohm[i]);
+    fg >> c->deltan >> c->deltat >> c->sigmax
+       >> c->taumax >> c->mu;
+    if (gd->ncoh)
+      fg >> c->Sinit;
+    cl(fg, buf, 2);
+  }
+  //read impact data
+  fg >> gd->imp >> gd->voImp >> gd->massImp >> gd->radiusImp;
+  cl(fg, buf, 2);
+  fg >> gd->eImp >> gd->xnuImp;
+  cl(fg, buf, 2);
+  gd->voImp = 0; gd->massImp = 1.0; gd->radiusImp = 0; 
+  gd->eImp = 1.0; gd->xnuImp = 0.3;
+  //read thermal load
+  fg >> dtmp;
+  fg.close();
+  fm >> itmp >> itmp >> gd->delta >> gd->numProp;
+  gd->delta /= (double) gd->steps;
+  gd->delta2 = gd->delta*gd->delta*0.5;
+  gd->ts_proportion = new int[gd->numProp];
+  gd->proportion = new double[gd->numProp];
+  for (i=0; i< gd->numProp; i++) {
+    fm >> gd->ts_proportion[i] >> gd->proportion[i];
+  }
+  //read nodal co-ordinates
+  fm >> gd->numNP;
+  cl(fm,buf,1);
+  int curline = 0;
+  for(i=0; i<gd->nn; i++)
+  {
+    cl(fg, buf, gd->nnums[i]-curline);
+    curline = gd->nnums[i];
+    Node *np = &(gd->nodes[i]);
+    fm >> itmp >> np->pos.x >> np->pos.y;
+    if(itmp != gd->nnums[i]+1)
+    {
+      CkError("Expected to read node %d got %d\n", gd->nnums[i]+1, itmp);
+      CkAbort("");
+    }
+    np->xM.x = np->xM.y = 0;
+    np->disp.x = np->disp.y = 0.0;
+    np->vel.x = np->vel.y = 0.0;
+    np->accel.x = np->accel.y = 0.0;
+  }
+  cl(fg, buf, gd->numNP-curline);
+  //read nodal boundary conditions
+  fm >> gd->numBound;
+  for (i=0; i<gd->numBound; i++)
+  {
+    int j, id0;
+    fm >> id0;
+    id0--;
+    for (j=0; j<gd->nn && id0!=gd->nnums[j]; j++);
+    if(j==gd->nn)
+    {
+      cl(fm, buf, 1);
+      continue;
+    }
+    Node *np = &(gd->nodes[j]);
+    np->isbnd = 1;
+    fm >> np->id1 >> np->id2 >> np->r.x >> np->r.y;
+  }
+  //read cohesive element data, determine the length and angle
+  fm >> itmp >> gd->numCLST >> itmp >> itmp >> itmp;
+  cl(fm,buf,1);
+  curline = 0;
+  for(i=0;i<gd->ne && gd->enums[i]<gd->numCLST;i++);
+  gd->scoh = 0;
+  gd->ecoh = i;
+  gd->svol = i;
+  gd->evol = gd->ne;
+  for(i=0; i<gd->ne && gd->enums[i]<gd->numCLST; i++)
+  {
+    cl(fg, buf, gd->enums[i]-curline);
+    curline = gd->enums[i];
+    Element *ep = &(gd->elements[i]);
+    fm >> ep->material; ep->material--;
+    int k;
+    for(k=0;k<6;k++) {
+      fm >> itmp;
+    }
+    fm >> itmp;
+    Node *np1 = &(gd->nodes[gd->conn[6*gd->enums[i]+1]]);
+    Node *np2 = &(gd->nodes[gd->conn[6*gd->enums[i]]]);
+    Coh *coh = &(ep->c);
+    coh->Sthresh[2] = coh->Sthresh[1] =
+      coh->Sthresh[0] = gd->cohm[ep->material].Sinit;
+    double x = np1->pos.x - np2->pos.x;
+    double y = np1->pos.y - np2->pos.y;
+    coh->sidel[0] = sqrt(x*x+y*y);
+    coh->sidel[1] = x/coh->sidel[0];
+    coh->sidel[2] = y/coh->sidel[0];
+  }
+  cl(fm, buf, gd->numCLST-curline-1);
+  fm >> itmp >> gd->numLST >> itmp;
+  cl(fm,buf,1);
+  curline = gd->numCLST;
+  for(gd->svol;i<gd->ne;i++)
+  {
+    cl(fm, buf, gd->enums[i]-curline);
+    curline = gd->enums[i];
+    Element *ep = &(gd->elements[i]);
+    fm >> ep->material; ep->material--;
+    int k;
+    for(k=0;k<6;k++) fm >> itmp;
+    for(k=0;k<3;k++)
+    {
+      ep->v.s11l[k] = 0.0;
+      ep->v.s22l[k] = 0.0;
+      ep->v.s12l[k] = 0.0;
+    }
+  }
+  cl(fm, buf, gd->numCLST-curline);
+  fm.close();
+  vol_elem(gd);
+}
diff --git a/examples/fem/crack2D/volume.C b/examples/fem/crack2D/volume.C
new file mode 100644 (file)
index 0000000..a8ff307
--- /dev/null
@@ -0,0 +1,83 @@
+#include "crack.h"
+#include <math.h>
+
+void 
+vol_elem(GlobalData *gd)
+{
+  double x, x1, x2, x3;
+
+  //Compute the elastic stiffness constants for each material type
+  int matNo;
+  for (matNo=0;matNo<gd->numMatVol;matNo++)
+  {
+    VolMaterial *vm=&(gd->volm[matNo]);
+    switch (gd->nplane)
+    {
+      case 1://Orthotropic plane strain
+        double sT,cT,xnu21;
+        cT = cos(vm->theta*1.74532925199e-2);
+        sT = sin(vm->theta*1.74532925199e-2);
+        xnu21 = vm->xnu12*vm->e2/vm->e1;
+        x = 1.0 - vm->xnu23*vm->xnu23 - 
+          2.0*vm->xnu12*xnu21*(1.0 + vm->xnu23);
+        x1 = vm->e1*(1.0-vm->xnu23*vm->xnu23) / x;
+        x2 = xnu21*vm->e1*(1.0+vm->xnu23) / x;
+        x3 = vm->e2*(vm->xnu23+vm->xnu12*xnu21) / x;
+        vm->c[2] = vm->e2*(1.0-vm->xnu12*xnu21) / x;
+        vm->c[0] = x1*cT*cT*cT*cT + 2.0*(x2+2.0*vm->g12)*cT*cT*sT*sT +
+          vm->c[2]*sT*sT*sT*sT;
+        vm->c[1] = x2*cT*cT + x3*sT*sT;
+        vm->c[3] = vm->g12*cT*cT + vm->g23*sT*sT;
+        break;
+      case 0: //Plane stress (isotropic)
+        vm->c[0] = vm->e1 / (1.0 - vm->xnu12*vm->xnu12);
+        vm->c[1] = vm->e1*vm->xnu12 / (1.0 - vm->xnu12*vm->xnu12);
+        vm->c[2] = vm->c[0];
+        vm->c[3] = vm->e1/ (2.0 * (1.0 + vm->xnu12));
+        break;
+      case 2: //Axisymmetric (isotropic)
+        vm->c[0] = vm->e1 * (1.0 - vm->xnu12) / ((1.0 + vm->xnu12)*
+                                                 (1.0 - 2.0*vm->xnu12));
+        vm->c[1] = vm->e1 * vm->xnu12 / ((1.0 + vm->xnu12)*
+                                         (1.0 - 2.0*vm->xnu12));
+        vm->c[2] = vm->e1 / (2.0*(1.0 + vm->xnu12));
+        break;
+      default:
+        CkAbort("Unknown planar analysis type passed to vol_elem!\n");
+    }
+  }
+
+  //Update the node-by-node mass of each element
+  int volNo;
+  for (volNo=gd->svol;volNo<gd->ne;volNo++)
+  {
+    Node *n[6];              //Pointers to each of the triangle's nodes
+    int k;                  //Loop index
+    for (k=0;k<6;k++)
+      n[k]=&(gd->nodes[gd->conn[volNo*6+k]]);
+    //Compute the mass of this element
+    double area=((n[1]->pos.x-n[0]->pos.x)*(n[2]->pos.y-n[0]->pos.y)-
+                 (n[2]->pos.x-n[0]->pos.x)*(n[1]->pos.y-n[0]->pos.y));
+    double mass=gd->volm[gd->elements[volNo].material].rho*area/114.0;
+    //Divide the mass among the element's nodes
+    for (k=0;k<3;k++) {
+      n[k]->xM.x+=mass*3.0;
+      n[k]->xM.y+=mass*3.0;
+    }
+    for (k=3;k<6;k++) {
+      n[k]->xM.x+=mass*16.0;
+      n[k]->xM.y+=mass*16.0;
+    }
+  }  
+  // FIXME: Do an FEM_Update here
+}
+
+// FIXME: incorporate this
+void updateElementMass(Element *v, Node **n)
+{
+  for(int k=0;k<3;k++) {
+    v->v.s11l[k] = 0.0;
+    v->v.s22l[k] = 0.0;
+    v->v.s12l[k] = 0.0;
+  }
+}