aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorrhaas <rhaas@8e189c6b-2ab8-4400-aa02-70a9cfce18b9>2012-06-22 17:18:12 +0000
committerrhaas <rhaas@8e189c6b-2ab8-4400-aa02-70a9cfce18b9>2012-06-22 17:18:12 +0000
commit323497aef47512b4fdf138ac6d0e36b1bef7adfb (patch)
tree5a97128fc0ebd803b9e10ea2e17c303166dcaed4
parentb7f6e4f941dabfa1f05b79e8ed04eef8a12638e6 (diff)
EOS_Omni: add table reader test
This adds a test (plus sample hdf5 table) using a small routine EOS_OMNI_dumptable to output the read in data as ASCII into a user selected file. Code and table kindly provided by Evan O'Connor. Small changes afterwards by Roland Haas. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/trunk@63 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
Diffstat
-rw-r--r--param.ccl15
-rw-r--r--schedule.ccl10
-rw-r--r--src/nuc_eos/dumpASCIItable.F90102
-rw-r--r--src/nuc_eos/make.code.defn2
-rw-r--r--test/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.h5bin0 -> 21384 bytes
-rw-r--r--test/eos_omni.par11
-rw-r--r--test/eos_omni/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.asc1268
-rw-r--r--test/test.ccl1
8 files changed, 1407 insertions, 2 deletions
diff --git a/param.ccl b/param.ccl
index 7a2db1c..cc98523 100644
--- a/param.ccl
+++ b/param.ccl
@@ -66,7 +66,20 @@ BOOLEAN do_energy_shift "shift energies around?" STEERABLE=ALWAYS
{
} "yes"
-STRING nuceos_table_name "nuceos table name (hdf5)" STEERABLE=ALWAYS
+BOOLEAN dump_nuceos_table "Dump table in ASCII at beginning" STEERABLE=RECOVER
+{
+} "No"
+
+STRING dump_nuceos_table_name "nuceos dump table name (ASCII)" STEERABLE=RECOVER
+{
+ .* :: "Can be anything"
+} "blah.asc"
+
+STRING nuceos_table_name "nuceos table name (hdf5)"
{
.* :: "Can be anything"
} "blah.h5"
+
+shares: IO
+
+USES STRING out_dir \ No newline at end of file
diff --git a/schedule.ccl b/schedule.ccl
index 74b5d33..9920825 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -14,5 +14,15 @@ if (nuceos_read_table)
LANG: C
OPTIONS: global
} "Read EOS HDF5 table"
+
+ if (dump_nuceos_table)
+ {
+ SCHEDULE EOS_OMNI_dumptable after EOS_OMNI_ReadTable AT CCTK_BASEGRID
+ {
+ LANG: FORTRAN
+ OPTIONS: global
+ } "Dump EOS HDF5 table in ASCII"
+ }
+
}
diff --git a/src/nuc_eos/dumpASCIItable.F90 b/src/nuc_eos/dumpASCIItable.F90
new file mode 100644
index 0000000..ee89acf
--- /dev/null
+++ b/src/nuc_eos/dumpASCIItable.F90
@@ -0,0 +1,102 @@
+#include "cctk.h"
+#include "cctk_Parameters.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+
+subroutine EOS_OMNI_dumptable(CCTK_ARGUMENTS)
+
+ use eosmodule
+ implicit none
+
+ DECLARE_CCTK_PARAMETERS
+ DECLARE_CCTK_ARGUMENTS
+ DECLARE_CCTK_FUNCTIONS
+
+ integer :: irho,itemp,iye,n
+ integer :: strlength1,strlength2
+ integer, parameter :: maxstrlength=200
+ character(len=200) :: fname, dirname, fullpath
+
+ if (CCTK_MyProc(cctkGH) .ne. 0) then
+ return
+ end if
+
+ call CCTK_FortranString(strlength1,out_dir,dirname)
+ if (strlength1 .gt. maxstrlength) then
+ call CCTK_WARN(0,"The output directory string is too long")
+ end if
+
+ call CCTK_FortranString(strlength2,dump_nuceos_table_name,fname)
+ if (strlength2 .gt. maxstrlength) then
+ call CCTK_WARN(0,"The output filename string is too long")
+ end if
+
+ if (strlength1+strlength2+1 .gt. maxstrlength) then
+ call CCTK_WARN(0,"The full output path string is too long")
+ end if
+
+ fullpath = trim(dirname)//'/'//trim(fname)
+
+ call CCTK_Info(CCTK_THORNSTRING,"*******************************");
+ call CCTK_Info(CCTK_THORNSTRING,"Dumping nuc_eos table file in ASCII:");
+ call CCTK_Info(CCTK_THORNSTRING,trim(fullpath));
+ call CCTK_Info(CCTK_THORNSTRING,"*******************************");
+
+ open(unit=473,file=trim(fullpath),status='unknown',form='formatted',position='rewind')
+
+ write(473,"('# ',a20,/,i4)") "nrho:",nrho
+ write(473,"('# ',a20,/,i4)") "ntemp:",ntemp
+ write(473,"('# ',a20,/,i4)") "nye:",nye
+
+ write(473,"('# ',a20,/,1P1E18.9)") "energy shift:",energy_shift
+
+ write(473,"('# ',a20,/,1P2E18.9)") "rho min and max:",eos_rhomin,eos_rhomax
+ write(473,"('# ',a20,/,1P2E18.9)") "ye min and max:",eos_yemin,eos_yemax
+ write(473,"('# ',a20,/,1P2E18.9)") "temp min and max:",eos_tempmin,eos_tempmax
+
+ write(473,"('# ',a20,/)") "log rho points:"
+ do irho=1,nrho
+ write(473,"(E18.9)") logrho(irho)
+ enddo
+ write(473,"('# ',a20)") "log temp points:"
+ do itemp=1,ntemp
+ write(473,"(E18.9)") logtemp(itemp)
+ enddo
+ write(473,"('#',a20)") "ye points:"
+ do iye=1,nye
+ write(473,"(E18.9)") ye(iye)
+ enddo
+
+ write(473,"('# ',a20,/,i4)") "nvars:",nvars
+ write(473,"('# ',a20)") "table mappings:"
+ write(473,"('# ',a20)") " 1 -> logpress"
+ write(473,"('# ',a20)") " 2 -> logenergy"
+ write(473,"('# ',a20)") " 3 -> entropy"
+ write(473,"('# ',a20)") " 4 -> munu"
+ write(473,"('# ',a20)") " 5 -> cs2"
+ write(473,"('# ',a20)") " 6 -> dedT"
+ write(473,"('# ',a20)") " 7 -> dpdrhoe"
+ write(473,"('# ',a20)") " 8 -> dpderho"
+ write(473,"('# ',a20)") " 9 -> muhat"
+ write(473,"('# ',a20)") "10 -> mu_e"
+ write(473,"('# ',a20)") "11 -> mu_p"
+ write(473,"('# ',a20)") "12 -> mu_n"
+ write(473,"('# ',a20)") "13 -> xa"
+ write(473,"('# ',a20)") "14 -> xh"
+ write(473,"('# ',a20)") "15 -> xn"
+ write(473,"('# ',a20)") "16 -> xp"
+ write(473,"('# ',a20)") "17 -> abar"
+ write(473,"('# ',a20)") "18 -> zbar"
+ write(473,"('# ',a20)") "19 -> gamma"
+
+ do irho=1,nrho
+ do itemp=1,ntemp
+ do iye=1,nye
+ do n=1,nvars
+ write(473,"(i4,i4,i4,i4,E18.9)") irho,itemp,iye,n,alltables(irho,itemp,iye,n)
+ enddo
+ enddo
+ enddo
+ enddo
+
+end subroutine EOS_OMNI_dumptable
diff --git a/src/nuc_eos/make.code.defn b/src/nuc_eos/make.code.defn
index 734ab63..ebc57ad 100644
--- a/src/nuc_eos/make.code.defn
+++ b/src/nuc_eos/make.code.defn
@@ -1,6 +1,6 @@
SRCS = eosmodule.F90 nuc_eos.F90 bisection.F90 \
findtemp.F90 findrho.F90 linterp_many.F90 readtable.c \
- linterp.F
+ linterp.F dumpASCIItable.F90
SUBDIRS =
diff --git a/test/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.h5 b/test/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.h5
new file mode 100644
index 0000000..a013049
--- /dev/null
+++ b/test/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.h5
Binary files differ
diff --git a/test/eos_omni.par b/test/eos_omni.par
new file mode 100644
index 0000000..8d97d25
--- /dev/null
+++ b/test/eos_omni.par
@@ -0,0 +1,11 @@
+ActiveThorns = "eos_omni IOUtil PUGH"
+
+EOS_Omni::nuceos_read_table = yes
+EOS_Omni::nuceos_table_name = "../../../arrangements/EinsteinEOS/EOS_Omni/test/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.h5" #$ENV{'EOS_TABLE_PATH'}/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.h5
+EOS_Omni::dump_nuceos_table = yes
+EOS_Omni::dump_nuceos_table_name = "TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.asc"
+
+IO::out_dir = $parfile
+
+#IOASCII::out1D_every = 1
+#IOASCII::out1D_vars = ""
diff --git a/test/eos_omni/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.asc b/test/eos_omni/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.asc
new file mode 100644
index 0000000..2ba12f5
--- /dev/null
+++ b/test/eos_omni/TEST_TABLEEOS_rho4_temp4_ye4_version1.0_20120514.asc
@@ -0,0 +1,1268 @@
+# nrho:
+ 4
+# ntemp:
+ 4
+# nye:
+ 4
+# energy shift:
+ 0.000000000E+00
+# rho min and max:
+ 1.000000000E+13 1.584893192E+13
+# ye min and max:
+ 2.000000000E-01 2.500000000E-01
+# temp min and max:
+ 1.000000000E+00 1.258925412E+00
+# log rho points:
+
+ 0.130000000E+02
+ 0.130666667E+02
+ 0.131333333E+02
+ 0.132000000E+02
+# log temp points:
+ 0.000000000E+00
+ 0.333333333E-01
+ 0.666666667E-01
+ 0.100000000E+00
+# ye points:
+ 0.200000000E+00
+ 0.216666667E+00
+ 0.233333333E+00
+ 0.250000000E+00
+# nvars:
+ 19
+# table mappings:
+# 1 -> logpress
+# 2 -> logenergy
+# 3 -> entropy
+# 4 -> munu
+# 5 -> cs2
+# 6 -> dedT
+# 7 -> dpdrhoe
+# 8 -> dpderho
+# 9 -> muhat
+# 10 -> mu_e
+# 11 -> mu_p
+# 12 -> mu_n
+# 13 -> xa
+# 14 -> xh
+# 15 -> xn
+# 16 -> xp
+# 17 -> abar
+# 18 -> zbar
+# 19 -> gamma
+ 1 1 1 1 0.315194459E+02
+ 1 1 1 2 0.190401865E+02
+ 1 1 1 3 0.626265020E+00
+ 1 1 1 4 -0.337628160E+02
+ 1 1 1 5 0.434856444E+19
+ 1 1 1 6 0.474795459E+18
+ 1 1 1 7 0.206663980E+19
+ 1 1 1 8 0.578756290E+13
+ 1 1 1 9 0.337628160E+02
+ 1 1 1 10 0.000000000E+00
+ 1 1 1 11 -0.309085200E+02
+ 1 1 1 12 0.285429600E+01
+ 1 1 1 13 0.243180600E-13
+ 1 1 1 14 0.752350100E+00
+ 1 1 1 15 0.247649900E+00
+ 1 1 1 16 0.205969600E-14
+ 1 1 1 17 0.175091400E+03
+ 1 1 1 18 0.465452000E+02
+ 1 1 1 19 0.131492219E+01
+ 1 1 2 1 0.315564264E+02
+ 1 1 2 2 0.190596440E+02
+ 1 1 2 3 0.599554076E+00
+ 1 1 2 4 -0.331028000E+02
+ 1 1 2 5 0.468308289E+19
+ 1 1 2 6 0.478635864E+18
+ 1 1 2 7 0.226237031E+19
+ 1 1 2 8 0.566244862E+13
+ 1 1 2 9 0.331028000E+02
+ 1 1 2 10 0.000000000E+00
+ 1 1 2 11 -0.304587144E+02
+ 1 1 2 12 0.264408553E+01
+ 1 1 2 13 0.390064804E-13
+ 1 1 2 14 0.799733996E+00
+ 1 1 2 15 0.200266004E+00
+ 1 1 2 16 0.320859170E-14
+ 1 1 2 17 0.176534081E+03
+ 1 1 2 18 0.478275544E+02
+ 1 1 2 19 0.130048527E+01
+ 1 1 3 1 0.315920997E+02
+ 1 1 3 2 0.190796823E+02
+ 1 1 3 3 0.564935397E+00
+ 1 1 3 4 -0.323601052E+02
+ 1 1 3 5 0.495008356E+19
+ 1 1 3 6 0.425729974E+18
+ 1 1 3 7 0.264081594E+19
+ 1 1 3 8 0.507088192E+13
+ 1 1 3 9 0.323601052E+02
+ 1 1 3 10 0.000000000E+00
+ 1 1 3 11 -0.299534000E+02
+ 1 1 3 12 0.240670520E+01
+ 1 1 3 13 0.668891830E-13
+ 1 1 3 14 0.837877252E+00
+ 1 1 3 15 0.162122748E+00
+ 1 1 3 16 0.532925967E-14
+ 1 1 3 17 0.177770770E+03
+ 1 1 3 18 0.495058748E+02
+ 1 1 3 19 0.126623077E+01
+ 1 1 4 1 0.316260359E+02
+ 1 1 4 2 0.191000254E+02
+ 1 1 4 3 0.516822724E+00
+ 1 1 4 4 -0.314919366E+02
+ 1 1 4 5 0.522834497E+19
+ 1 1 4 6 0.374311992E+18
+ 1 1 4 7 0.297805957E+19
+ 1 1 4 8 0.470964643E+13
+ 1 1 4 9 0.314919366E+02
+ 1 1 4 10 0.000000000E+00
+ 1 1 4 11 -0.293611800E+02
+ 1 1 4 12 0.213075660E+01
+ 1 1 4 13 0.125220900E-12
+ 1 1 4 14 0.875112600E+00
+ 1 1 4 15 0.124887400E+00
+ 1 1 4 16 0.959418900E-14
+ 1 1 4 17 0.179103700E+03
+ 1 1 4 18 0.511659100E+02
+ 1 1 4 19 0.123688231E+01
+ 1 2 1 1 0.315222156E+02
+ 1 2 1 2 0.190416293E+02
+ 1 2 1 3 0.660711914E+00
+ 1 2 1 4 -0.336092014E+02
+ 1 2 1 5 0.434846621E+19
+ 1 2 1 6 0.468913161E+18
+ 1 2 1 7 0.209524430E+19
+ 1 2 1 8 0.571465738E+13
+ 1 2 1 9 0.336092014E+02
+ 1 2 1 10 0.000000000E+00
+ 1 2 1 11 -0.308083719E+02
+ 1 2 1 12 0.280082948E+01
+ 1 2 1 13 0.351987318E-12
+ 1 2 1 14 0.750281132E+00
+ 1 2 1 15 0.249718868E+00
+ 1 2 1 16 0.292622427E-13
+ 1 2 1 17 0.174424797E+03
+ 1 2 1 18 0.464958533E+02
+ 1 2 1 19 0.130653344E+01
+ 1 2 2 1 0.315589352E+02
+ 1 2 2 2 0.190610347E+02
+ 1 2 2 3 0.633692510E+00
+ 1 2 2 4 -0.329480451E+02
+ 1 2 2 5 0.468101344E+19
+ 1 2 2 6 0.467722998E+18
+ 1 2 2 7 0.230683494E+19
+ 1 2 2 8 0.555549172E+13
+ 1 2 2 9 0.329480451E+02
+ 1 2 2 10 0.000000000E+00
+ 1 2 2 11 -0.303588848E+02
+ 1 2 2 12 0.258916031E+01
+ 1 2 2 13 0.540236557E-12
+ 1 2 2 14 0.793895117E+00
+ 1 2 2 15 0.206104883E+00
+ 1 2 2 16 0.439033773E-13
+ 1 2 2 17 0.175968136E+03
+ 1 2 2 18 0.480245169E+02
+ 1 2 2 19 0.129242308E+01
+ 1 2 3 1 0.315939404E+02
+ 1 2 3 2 0.190808635E+02
+ 1 2 3 3 0.594557026E+00
+ 1 2 3 4 -0.321767168E+02
+ 1 2 3 5 0.495965366E+19
+ 1 2 3 6 0.422328531E+18
+ 1 2 3 7 0.266223158E+19
+ 1 2 3 8 0.505373306E+13
+ 1 2 3 9 0.321767168E+02
+ 1 2 3 10 0.000000000E+00
+ 1 2 3 11 -0.298341150E+02
+ 1 2 3 12 0.234260179E+01
+ 1 2 3 13 0.902817830E-12
+ 1 2 3 14 0.835082063E+00
+ 1 2 3 15 0.164917937E+00
+ 1 2 3 16 0.711626882E-13
+ 1 2 3 17 0.176810347E+03
+ 1 2 3 18 0.494032258E+02
+ 1 2 3 19 0.126331288E+01
+ 1 2 4 1 0.316274260E+02
+ 1 2 4 2 0.191010164E+02
+ 1 2 4 3 0.541914008E+00
+ 1 2 4 4 -0.313019673E+02
+ 1 2 4 5 0.521419270E+19
+ 1 2 4 6 0.399578678E+18
+ 1 2 4 7 0.299734779E+19
+ 1 2 4 8 0.468927295E+13
+ 1 2 4 9 0.313019673E+02
+ 1 2 4 10 0.000000000E+00
+ 1 2 4 11 -0.292387772E+02
+ 1 2 4 12 0.206319014E+01
+ 1 2 4 13 0.159787662E-11
+ 1 2 4 14 0.873955930E+00
+ 1 2 4 15 0.126044070E+00
+ 1 2 4 16 0.122006787E-12
+ 1 2 4 17 0.177048951E+03
+ 1 2 4 18 0.506458464E+02
+ 1 2 4 19 0.122959209E+01
+ 1 3 1 1 0.315252416E+02
+ 1 3 1 2 0.190432535E+02
+ 1 3 1 3 0.696275848E+00
+ 1 3 1 4 -0.334275753E+02
+ 1 3 1 5 0.439131729E+19
+ 1 3 1 6 0.501571155E+18
+ 1 3 1 7 0.213753832E+19
+ 1 3 1 8 0.561803121E+13
+ 1 3 1 9 0.334275753E+02
+ 1 3 1 10 0.000000000E+00
+ 1 3 1 11 -0.306898872E+02
+ 1 3 1 12 0.273768813E+01
+ 1 3 1 13 0.384741861E-11
+ 1 3 1 14 0.748328001E+00
+ 1 3 1 15 0.251671999E+00
+ 1 3 1 16 0.315111706E-12
+ 1 3 1 17 0.173750622E+03
+ 1 3 1 18 0.464370368E+02
+ 1 3 1 19 0.131024714E+01
+ 1 3 2 1 0.315615908E+02
+ 1 3 2 2 0.190625695E+02
+ 1 3 2 3 0.668217453E+00
+ 1 3 2 4 -0.327737010E+02
+ 1 3 2 5 0.470824884E+19
+ 1 3 2 6 0.477858826E+18
+ 1 3 2 7 0.237631620E+19
+ 1 3 2 8 0.537880382E+13
+ 1 3 2 9 0.327737010E+02
+ 1 3 2 10 0.000000000E+00
+ 1 3 2 11 -0.302479143E+02
+ 1 3 2 12 0.252578666E+01
+ 1 3 2 13 0.563665790E-11
+ 1 3 2 14 0.789471721E+00
+ 1 3 2 15 0.210528279E+00
+ 1 3 2 16 0.455271302E-12
+ 1 3 2 17 0.175218143E+03
+ 1 3 2 18 0.480877655E+02
+ 1 3 2 19 0.129201817E+01
+ 1 3 3 1 0.315960132E+02
+ 1 3 3 2 0.190822059E+02
+ 1 3 3 3 0.625596908E+00
+ 1 3 3 4 -0.319633736E+02
+ 1 3 3 5 0.498670822E+19
+ 1 3 3 6 0.452532618E+18
+ 1 3 3 7 0.268667794E+19
+ 1 3 3 8 0.501567081E+13
+ 1 3 3 9 0.319633736E+02
+ 1 3 3 10 0.000000000E+00
+ 1 3 3 11 -0.296950340E+02
+ 1 3 3 12 0.226833967E+01
+ 1 3 3 13 0.928242674E-11
+ 1 3 3 14 0.831988340E+00
+ 1 3 3 15 0.168011660E+00
+ 1 3 3 16 0.725559079E-12
+ 1 3 3 17 0.176314517E+03
+ 1 3 3 18 0.494478733E+02
+ 1 3 3 19 0.126415638E+01
+ 1 3 4 1 0.316292166E+02
+ 1 3 4 2 0.191023048E+02
+ 1 3 4 3 0.574208050E+00
+ 1 3 4 4 -0.310891248E+02
+ 1 3 4 5 0.523257005E+19
+ 1 3 4 6 0.446288391E+18
+ 1 3 4 7 0.301882615E+19
+ 1 3 4 8 0.466750863E+13
+ 1 3 4 9 0.310891248E+02
+ 1 3 4 10 0.000000000E+00
+ 1 3 4 11 -0.291026963E+02
+ 1 3 4 12 0.198642845E+01
+ 1 3 4 13 0.155991450E-10
+ 1 3 4 14 0.868837055E+00
+ 1 3 4 15 0.131162945E+00
+ 1 3 4 16 0.119188258E-11
+ 1 3 4 17 0.176730296E+03
+ 1 3 4 18 0.508525433E+02
+ 1 3 4 19 0.122884890E+01
+ 1 4 1 1 0.315287506E+02
+ 1 4 1 2 0.190451583E+02
+ 1 4 1 3 0.734301740E+00
+ 1 4 1 4 -0.332363629E+02
+ 1 4 1 5 0.453946149E+19
+ 1 4 1 6 0.503665796E+18
+ 1 4 1 7 0.177378507E+19
+ 1 4 1 8 0.691215987E+13
+ 1 4 1 9 0.332363629E+02
+ 1 4 1 10 0.000000000E+00
+ 1 4 1 11 -0.305679368E+02
+ 1 4 1 12 0.266842616E+01
+ 1 4 1 13 0.322264710E-10
+ 1 4 1 14 0.744361828E+00
+ 1 4 1 15 0.255638172E+00
+ 1 4 1 16 0.262655539E-11
+ 1 4 1 17 0.172444348E+03
+ 1 4 1 18 0.463334742E+02
+ 1 4 1 19 0.134354956E+01
+ 1 4 2 1 0.315644176E+02
+ 1 4 2 2 0.190642507E+02
+ 1 4 2 3 0.702828246E+00
+ 1 4 2 4 -0.325584892E+02
+ 1 4 2 5 0.484401203E+19
+ 1 4 2 6 0.464497743E+18
+ 1 4 2 7 0.200315166E+19
+ 1 4 2 8 0.657441480E+13
+ 1 4 2 9 0.325584892E+02
+ 1 4 2 10 0.000000000E+00
+ 1 4 2 11 -0.301091763E+02
+ 1 4 2 12 0.244931295E+01
+ 1 4 2 13 0.470536636E-10
+ 1 4 2 14 0.786671061E+00
+ 1 4 2 15 0.213328939E+00
+ 1 4 2 16 0.376649741E-11
+ 1 4 2 17 0.174391667E+03
+ 1 4 2 18 0.480313349E+02
+ 1 4 2 19 0.132064971E+01
+ 1 4 3 1 0.315984128E+02
+ 1 4 3 2 0.190837743E+02
+ 1 4 3 3 0.659494730E+00
+ 1 4 3 4 -0.317355493E+02
+ 1 4 3 5 0.510293069E+19
+ 1 4 3 6 0.453270285E+18
+ 1 4 3 7 0.227243833E+19
+ 1 4 3 8 0.620891423E+13
+ 1 4 3 9 0.317355493E+02
+ 1 4 3 10 0.000000000E+00
+ 1 4 3 11 -0.295489506E+02
+ 1 4 3 12 0.218659878E+01
+ 1 4 3 13 0.740981732E-10
+ 1 4 3 14 0.827535649E+00
+ 1 4 3 15 0.172464351E+00
+ 1 4 3 16 0.578803586E-11
+ 1 4 3 17 0.175389298E+03
+ 1 4 3 18 0.494530594E+02
+ 1 4 3 19 0.128649126E+01
+ 1 4 4 1 0.316312284E+02
+ 1 4 4 2 0.191037573E+02
+ 1 4 4 3 0.607816102E+00
+ 1 4 4 4 -0.308402339E+02
+ 1 4 4 5 0.533700058E+19
+ 1 4 4 6 0.439520318E+18
+ 1 4 4 7 0.259039830E+19
+ 1 4 4 8 0.580466950E+13
+ 1 4 4 9 0.308402339E+02
+ 1 4 4 10 0.000000000E+00
+ 1 4 4 11 -0.289432036E+02
+ 1 4 4 12 0.189703033E+01
+ 1 4 4 13 0.121692782E-09
+ 1 4 4 14 0.863663057E+00
+ 1 4 4 15 0.136336943E+00
+ 1 4 4 16 0.930546051E-11
+ 1 4 4 17 0.175946308E+03
+ 1 4 4 18 0.509302517E+02
+ 1 4 4 19 0.124758119E+01
+ 2 1 1 1 0.315961745E+02
+ 2 1 1 2 0.190584537E+02
+ 2 1 1 3 0.506953221E+00
+ 2 1 1 4 0.338646779E+02
+ 2 1 1 5 0.438432816E+19
+ 2 1 1 6 0.460796338E+18
+ 2 1 1 7 0.250773832E+19
+ 2 1 1 8 0.601452614E+13
+ 2 1 1 9 0.344502627E+02
+ 2 1 1 10 0.683149406E+02
+ 2 1 1 11 -0.314615622E+02
+ 2 1 1 12 0.298870048E+01
+ 2 1 1 13 0.865631956E-14
+ 2 1 1 14 0.762747493E+00
+ 2 1 1 15 0.237252507E+00
+ 2 1 1 16 0.985680548E-15
+ 2 1 1 17 0.183434274E+03
+ 2 1 1 18 0.480979396E+02
+ 2 1 1 19 0.129537426E+01
+ 2 1 2 1 0.316330616E+02
+ 2 1 2 2 0.190789082E+02
+ 2 1 2 3 0.480122511E+00
+ 2 1 2 4 0.364038559E+02
+ 2 1 2 5 0.475461081E+19
+ 2 1 2 6 0.441756753E+18
+ 2 1 2 7 0.272369466E+19
+ 2 1 2 8 0.599360640E+13
+ 2 1 2 9 0.337607887E+02
+ 2 1 2 10 0.701646446E+02
+ 2 1 2 11 -0.309899645E+02
+ 2 1 2 12 0.277082423E+01
+ 2 1 2 13 0.142583490E-13
+ 2 1 2 14 0.808157885E+00
+ 2 1 2 15 0.191842115E+00
+ 2 1 2 16 0.157783992E-14
+ 2 1 2 17 0.184478964E+03
+ 2 1 2 18 0.494587048E+02
+ 2 1 2 19 0.129038721E+01
+ 2 1 3 1 0.316684893E+02
+ 2 1 3 2 0.190998028E+02
+ 2 1 3 3 0.452094633E+00
+ 2 1 3 4 0.389194042E+02
+ 2 1 3 5 0.512012780E+19
+ 2 1 3 6 0.396083175E+18
+ 2 1 3 7 0.298964758E+19
+ 2 1 3 8 0.584555975E+13
+ 2 1 3 9 0.330023023E+02
+ 2 1 3 10 0.719217065E+02
+ 2 1 3 11 -0.304731997E+02
+ 2 1 3 12 0.252910261E+01
+ 2 1 3 13 0.247493666E-13
+ 2 1 3 14 0.846722905E+00
+ 2 1 3 15 0.153277095E+00
+ 2 1 3 16 0.264646295E-14
+ 2 1 3 17 0.186039982E+03
+ 2 1 3 18 0.512674559E+02
+ 2 1 3 19 0.128073170E+01
+ 2 1 4 1 0.317023250E+02
+ 2 1 4 2 0.191209278E+02
+ 2 1 4 3 0.418761608E+00
+ 2 1 4 4 0.415122817E+02
+ 2 1 4 5 0.548187562E+19
+ 2 1 4 6 0.399000298E+18
+ 2 1 4 7 0.339931787E+19
+ 2 1 4 8 0.534782251E+13
+ 2 1 4 9 0.320846385E+02
+ 2 1 4 10 0.735969202E+02
+ 2 1 4 11 -0.298461081E+02
+ 2 1 4 12 0.223853035E+01
+ 2 1 4 13 0.482410265E-13
+ 2 1 4 14 0.883766644E+00
+ 2 1 4 15 0.116233356E+00
+ 2 1 4 16 0.492899585E-14
+ 2 1 4 17 0.187155507E+03
+ 2 1 4 18 0.529427448E+02
+ 2 1 4 19 0.126844259E+01
+ 2 2 1 1 0.315985155E+02
+ 2 2 1 2 0.190597963E+02
+ 2 2 1 3 0.542291550E+00
+ 2 2 1 4 0.339965865E+02
+ 2 2 1 5 0.439986228E+19
+ 2 2 1 6 0.452978742E+18
+ 2 2 1 7 0.250321183E+19
+ 2 2 1 8 0.607824629E+13
+ 2 2 1 9 0.343103695E+02
+ 2 2 1 10 0.683069560E+02
+ 2 2 1 11 -0.313704320E+02
+ 2 2 1 12 0.293993749E+01
+ 2 2 1 13 0.133189427E-12
+ 2 2 1 14 0.759445400E+00
+ 2 2 1 15 0.240554600E+00
+ 2 2 1 16 0.146138950E-13
+ 2 2 1 17 0.182550244E+03
+ 2 2 1 18 0.480741763E+02
+ 2 2 1 19 0.129297528E+01
+ 2 2 2 1 0.316351160E+02
+ 2 2 2 2 0.190801361E+02
+ 2 2 2 3 0.513893558E+00
+ 2 2 2 4 0.365389043E+02
+ 2 2 2 5 0.476113622E+19
+ 2 2 2 6 0.440274015E+18
+ 2 2 2 7 0.275527548E+19
+ 2 2 2 8 0.592248836E+13
+ 2 2 2 9 0.336179660E+02
+ 2 2 2 10 0.701568703E+02
+ 2 2 2 11 -0.308978491E+02
+ 2 2 2 12 0.272011686E+01
+ 2 2 2 13 0.210190568E-12
+ 2 2 2 14 0.804173706E+00
+ 2 2 2 15 0.195826294E+00
+ 2 2 2 16 0.223779937E-13
+ 2 2 2 17 0.184405340E+03
+ 2 2 2 18 0.496839053E+02
+ 2 2 2 19 0.128606011E+01
+ 2 2 3 1 0.316701451E+02
+ 2 2 3 2 0.191008520E+02
+ 2 2 3 3 0.482458803E+00
+ 2 2 3 4 0.390702921E+02
+ 2 2 3 5 0.512597685E+19
+ 2 2 3 6 0.405421073E+18
+ 2 2 3 7 0.306900028E+19
+ 2 2 3 8 0.564568386E+13
+ 2 2 3 9 0.328438298E+02
+ 2 2 3 10 0.719141218E+02
+ 2 2 3 11 -0.303713939E+02
+ 2 2 3 12 0.247243584E+01
+ 2 2 3 13 0.351031186E-12
+ 2 2 3 14 0.844673648E+00
+ 2 2 3 15 0.155326352E+00
+ 2 2 3 16 0.363833824E-13
+ 2 2 3 17 0.185543797E+03
+ 2 2 3 18 0.512547690E+02
+ 2 2 3 19 0.127731561E+01
+ 2 2 4 1 0.317037366E+02
+ 2 2 4 2 0.191219346E+02
+ 2 2 4 3 0.449290317E+00
+ 2 2 4 4 0.416738181E+02
+ 2 2 4 5 0.548741548E+19
+ 2 2 4 6 0.392549942E+18
+ 2 2 4 7 0.343316957E+19
+ 2 2 4 8 0.528089410E+13
+ 2 2 4 9 0.319156899E+02
+ 2 2 4 10 0.735895080E+02
+ 2 2 4 11 -0.297376192E+02
+ 2 2 4 12 0.217807064E+01
+ 2 2 4 13 0.649860601E-12
+ 2 2 4 14 0.879343062E+00
+ 2 2 4 15 0.120656938E+00
+ 2 2 4 16 0.647523277E-13
+ 2 2 4 17 0.186443917E+03
+ 2 2 4 18 0.530065923E+02
+ 2 2 4 19 0.126560420E+01
+ 2 3 1 1 0.316011700E+02
+ 2 3 1 2 0.190612951E+02
+ 2 3 1 3 0.578931490E+00
+ 2 3 1 4 0.341459158E+02
+ 2 3 1 5 0.444110359E+19
+ 2 3 1 6 0.474262651E+18
+ 2 3 1 7 0.248695703E+19
+ 2 3 1 8 0.620338910E+13
+ 2 3 1 9 0.341517309E+02
+ 2 3 1 10 0.682976466E+02
+ 2 3 1 11 -0.312680079E+02
+ 2 3 1 12 0.288372297E+01
+ 2 3 1 13 0.152551770E-11
+ 2 3 1 14 0.757142311E+00
+ 2 3 1 15 0.242857689E+00
+ 2 3 1 16 0.162263336E-12
+ 2 3 1 17 0.181772437E+03
+ 2 3 1 18 0.480153953E+02
+ 2 3 1 19 0.129714214E+01
+ 2 3 2 1 0.316374102E+02
+ 2 3 2 2 0.190815435E+02
+ 2 3 2 3 0.549927823E+00
+ 2 3 2 4 0.366870165E+02
+ 2 3 2 5 0.478962883E+19
+ 2 3 2 6 0.459334690E+18
+ 2 3 2 7 0.276854554E+19
+ 2 3 2 8 0.592266312E+13
+ 2 3 2 9 0.334607895E+02
+ 2 3 2 10 0.701478060E+02
+ 2 3 2 11 -0.307983384E+02
+ 2 3 2 12 0.266245116E+01
+ 2 3 2 13 0.230710215E-11
+ 2 3 2 14 0.799416188E+00
+ 2 3 2 15 0.200583812E+00
+ 2 3 2 16 0.238619006E-12
+ 2 3 2 17 0.183797776E+03
+ 2 3 2 18 0.498149174E+02
+ 2 3 2 19 0.128693999E+01
+ 2 3 3 1 0.316720107E+02
+ 2 3 3 2 0.191021165E+02
+ 2 3 3 3 0.516363739E+00
+ 2 3 3 4 0.392525611E+02
+ 2 3 3 5 0.514585437E+19
+ 2 3 3 6 0.437101405E+18
+ 2 3 3 7 0.312817016E+19
+ 2 3 3 8 0.550314469E+13
+ 2 3 3 9 0.326527177E+02
+ 2 3 3 10 0.719052788E+02
+ 2 3 3 11 -0.302478324E+02
+ 2 3 3 12 0.240488531E+01
+ 2 3 3 13 0.378330720E-11
+ 2 3 3 14 0.840373469E+00
+ 2 3 3 15 0.159626531E+00
+ 2 3 3 16 0.382206541E-12
+ 2 3 3 17 0.184735327E+03
+ 2 3 3 18 0.512925637E+02
+ 2 3 3 19 0.127677218E+01
+ 2 3 4 1 0.317052770E+02
+ 2 3 4 2 0.191230619E+02
+ 2 3 4 3 0.481082237E+00
+ 2 3 4 4 0.418587613E+02
+ 2 3 4 5 0.551125058E+19
+ 2 3 4 6 0.407543083E+18
+ 2 3 4 7 0.344914972E+19
+ 2 3 4 8 0.527401338E+13
+ 2 3 4 9 0.317221047E+02
+ 2 3 4 10 0.735808660E+02
+ 2 3 4 11 -0.296144647E+02
+ 2 3 4 12 0.210764000E+01
+ 2 3 4 13 0.663521290E-11
+ 2 3 4 14 0.875579576E+00
+ 2 3 4 15 0.124420424E+00
+ 2 3 4 16 0.647728215E-12
+ 2 3 4 17 0.185947312E+03
+ 2 3 4 18 0.530926364E+02
+ 2 3 4 19 0.126660108E+01
+ 2 4 1 1 0.316042442E+02
+ 2 4 1 2 0.190629974E+02
+ 2 4 1 3 0.617584301E+00
+ 2 4 1 4 0.343088169E+02
+ 2 4 1 5 0.452888444E+19
+ 2 4 1 6 0.468985565E+18
+ 2 4 1 7 0.208919974E+19
+ 2 4 1 8 0.776010155E+13
+ 2 4 1 9 0.339779758E+02
+ 2 4 1 10 0.682867927E+02
+ 2 4 1 11 -0.311577433E+02
+ 2 4 1 12 0.282023246E+01
+ 2 4 1 13 0.134256291E-10
+ 2 4 1 14 0.754773328E+00
+ 2 4 1 15 0.245226672E+00
+ 2 4 1 16 0.139040158E-11
+ 2 4 1 17 0.180882573E+03
+ 2 4 1 18 0.479303034E+02
+ 2 4 1 19 0.131345051E+01
+ 2 4 2 1 0.316399892E+02
+ 2 4 2 2 0.190830952E+02
+ 2 4 2 3 0.586897258E+00
+ 2 4 2 4 0.368609157E+02
+ 2 4 2 5 0.487382766E+19
+ 2 4 2 6 0.447735446E+18
+ 2 4 2 7 0.234487621E+19
+ 2 4 2 8 0.742341535E+13
+ 2 4 2 9 0.332763222E+02
+ 2 4 2 10 0.701372379E+02
+ 2 4 2 11 -0.306815293E+02
+ 2 4 2 12 0.259479293E+01
+ 2 4 2 13 0.198111773E-10
+ 2 4 2 14 0.796207800E+00
+ 2 4 2 15 0.203792200E+00
+ 2 4 2 16 0.201230765E-11
+ 2 4 2 17 0.182906043E+03
+ 2 4 2 18 0.497729896E+02
+ 2 4 2 19 0.130181016E+01
+ 2 4 3 1 0.316741239E+02
+ 2 4 3 2 0.191035376E+02
+ 2 4 3 3 0.551848075E+00
+ 2 4 3 4 0.394563678E+02
+ 2 4 3 5 0.521542734E+19
+ 2 4 3 6 0.429836760E+18
+ 2 4 3 7 0.269969028E+19
+ 2 4 3 8 0.687457982E+13
+ 2 4 3 9 0.324386007E+02
+ 2 4 3 10 0.718949685E+02
+ 2 4 3 11 -0.301105181E+02
+ 2 4 3 12 0.232808251E+01
+ 2 4 3 13 0.317498529E-10
+ 2 4 3 14 0.836421243E+00
+ 2 4 3 15 0.163578757E+00
+ 2 4 3 16 0.313481318E-11
+ 2 4 3 17 0.183755024E+03
+ 2 4 3 18 0.512614579E+02
+ 2 4 3 19 0.128775317E+01
+ 2 4 4 1 0.317070298E+02
+ 2 4 4 2 0.191243207E+02
+ 2 4 4 3 0.514065938E+00
+ 2 4 4 4 0.420707598E+02
+ 2 4 4 5 0.557572415E+19
+ 2 4 4 6 0.399393855E+18
+ 2 4 4 7 0.298559936E+19
+ 2 4 4 8 0.663401019E+13
+ 2 4 4 9 0.315000303E+02
+ 2 4 4 10 0.735707901E+02
+ 2 4 4 11 -0.294734468E+02
+ 2 4 4 12 0.202658341E+01
+ 2 4 4 13 0.534406439E-10
+ 2 4 4 14 0.872417040E+00
+ 2 4 4 15 0.127582960E+00
+ 2 4 4 16 0.514339223E-11
+ 2 4 4 17 0.185124683E+03
+ 2 4 4 18 0.530493660E+02
+ 2 4 4 19 0.127625713E+01
+ 3 1 1 1 0.316822976E+02
+ 3 1 1 2 0.190769247E+02
+ 3 1 1 3 0.471456404E+00
+ 3 1 1 4 0.367559307E+02
+ 3 1 1 5 0.470106935E+19
+ 3 1 1 6 0.437471724E+18
+ 3 1 1 7 0.280072188E+19
+ 3 1 1 8 0.726680628E+13
+ 3 1 1 9 0.351502953E+02
+ 3 1 1 10 0.719062260E+02
+ 3 1 1 11 -0.320282000E+02
+ 3 1 1 12 0.312209530E+01
+ 3 1 1 13 0.307201557E-14
+ 3 1 1 14 0.772451330E+00
+ 3 1 1 15 0.227548670E+00
+ 3 1 1 16 0.455654279E-15
+ 3 1 1 17 0.192669111E+03
+ 3 1 1 18 0.498851135E+02
+ 3 1 1 19 0.132810489E+01
+ 3 1 2 1 0.317196878E+02
+ 3 1 2 2 0.190983142E+02
+ 3 1 2 3 0.448245362E+00
+ 3 1 2 4 0.394180026E+02
+ 3 1 2 5 0.512382990E+19
+ 3 1 2 6 0.411468494E+18
+ 3 1 2 7 0.316117368E+19
+ 3 1 2 8 0.688875217E+13
+ 3 1 2 9 0.344349128E+02
+ 3 1 2 10 0.738529155E+02
+ 3 1 2 11 -0.315375846E+02
+ 3 1 2 12 0.289732822E+01
+ 3 1 2 13 0.516769591E-14
+ 3 1 2 14 0.817085039E+00
+ 3 1 2 15 0.182914961E+00
+ 3 1 2 16 0.742748315E-15
+ 3 1 2 17 0.193814836E+03
+ 3 1 2 18 0.513939339E+02
+ 3 1 2 19 0.132812905E+01
+ 3 1 3 1 0.317555037E+02
+ 3 1 3 2 0.191200434E+02
+ 3 1 3 3 0.423163846E+00
+ 3 1 3 4 0.420486388E+02
+ 3 1 3 5 0.556377009E+19
+ 3 1 3 6 0.376403451E+18
+ 3 1 3 7 0.348076149E+19
+ 3 1 3 8 0.673908671E+13
+ 3 1 3 9 0.336534724E+02
+ 3 1 3 10 0.757021112E+02
+ 3 1 3 11 -0.310043548E+02
+ 3 1 3 12 0.264911755E+01
+ 3 1 3 13 0.910479221E-14
+ 3 1 3 14 0.854931144E+00
+ 3 1 3 15 0.145068856E+00
+ 3 1 3 16 0.126510306E-14
+ 3 1 3 17 0.196087874E+03
+ 3 1 3 18 0.535175700E+02
+ 3 1 3 19 0.132800215E+01
+ 3 1 4 1 0.317896561E+02
+ 3 1 4 2 0.191419233E+02
+ 3 1 4 3 0.392780938E+00
+ 3 1 4 4 0.447668925E+02
+ 3 1 4 5 0.601406817E+19
+ 3 1 4 6 0.404612606E+18
+ 3 1 4 7 0.406474283E+19
+ 3 1 4 8 0.583021850E+13
+ 3 1 4 9 0.326982794E+02
+ 3 1 4 10 0.774651719E+02
+ 3 1 4 11 -0.303508278E+02
+ 3 1 4 12 0.234745162E+01
+ 3 1 4 13 0.182269183E-13
+ 3 1 4 14 0.891876489E+00
+ 3 1 4 15 0.108123511E+00
+ 3 1 4 16 0.244181359E-14
+ 3 1 4 17 0.196733081E+03
+ 3 1 4 18 0.551459033E+02
+ 3 1 4 19 0.132692257E+01
+ 3 2 1 1 0.316844999E+02
+ 3 2 1 2 0.190781463E+02
+ 3 2 1 3 0.505002726E+00
+ 3 2 1 4 0.368698445E+02
+ 3 2 1 5 0.472259245E+19
+ 3 2 1 6 0.430754424E+18
+ 3 2 1 7 0.280921588E+19
+ 3 2 1 8 0.731171699E+13
+ 3 2 1 9 0.350287952E+02
+ 3 2 1 10 0.718986397E+02
+ 3 2 1 11 -0.319500848E+02
+ 3 2 1 12 0.307871037E+01
+ 3 2 1 13 0.497301741E-13
+ 3 2 1 14 0.769501332E+00
+ 3 2 1 15 0.230498668E+00
+ 3 2 1 16 0.702412109E-14
+ 3 2 1 17 0.191769944E+03
+ 3 2 1 18 0.498426541E+02
+ 3 2 1 19 0.132743692E+01
+ 3 2 2 1 0.317214894E+02
+ 3 2 2 2 0.190994079E+02
+ 3 2 2 3 0.479673401E+00
+ 3 2 2 4 0.395402529E+02
+ 3 2 2 5 0.513667845E+19
+ 3 2 2 6 0.415328808E+18
+ 3 2 2 7 0.319359963E+19
+ 3 2 2 8 0.681690553E+13
+ 3 2 2 9 0.343052759E+02
+ 3 2 2 10 0.738455289E+02
+ 3 2 2 11 -0.314546326E+02
+ 3 2 2 12 0.285064337E+01
+ 3 2 2 13 0.802209634E-13
+ 3 2 2 14 0.814307583E+00
+ 3 2 2 15 0.185692417E+00
+ 3 2 2 16 0.110255214E-13
+ 3 2 2 17 0.193508307E+03
+ 3 2 2 18 0.514876697E+02
+ 3 2 2 19 0.132594745E+01
+ 3 2 3 1 0.317569884E+02
+ 3 2 3 2 0.191209951E+02
+ 3 2 3 3 0.452031537E+00
+ 3 2 3 4 0.421805815E+02
+ 3 2 3 5 0.557374776E+19
+ 3 2 3 6 0.384882089E+18
+ 3 2 3 7 0.355729135E+19
+ 3 2 3 8 0.652007342E+13
+ 3 2 3 9 0.335143233E+02
+ 3 2 3 10 0.756949048E+02
+ 3 2 3 11 -0.309160752E+02
+ 3 2 3 12 0.259824809E+01
+ 3 2 3 13 0.135499488E-12
+ 3 2 3 14 0.853091249E+00
+ 3 2 3 15 0.146908751E+00
+ 3 2 3 16 0.180065867E-13
+ 3 2 3 17 0.195657394E+03
+ 3 2 3 18 0.535152507E+02
+ 3 2 3 19 0.132584344E+01
+ 3 2 4 1 0.317909324E+02
+ 3 2 4 2 0.191428961E+02
+ 3 2 4 3 0.423602426E+00
+ 3 2 4 4 0.449120736E+02
+ 3 2 4 5 0.602554197E+19
+ 3 2 4 6 0.381631911E+18
+ 3 2 4 7 0.402465571E+19
+ 3 2 4 8 0.598303922E+13
+ 3 2 4 9 0.325460558E+02
+ 3 2 4 10 0.774581294E+02
+ 3 2 4 11 -0.302534919E+02
+ 3 2 4 12 0.229256395E+01
+ 3 2 4 13 0.258042854E-12
+ 3 2 4 14 0.885636930E+00
+ 3 2 4 15 0.114363070E+00
+ 3 2 4 16 0.332054646E-13
+ 3 2 4 17 0.196751684E+03
+ 3 2 4 18 0.555396003E+02
+ 3 2 4 19 0.132555291E+01
+ 3 3 1 1 0.316869857E+02
+ 3 3 1 2 0.190795149E+02
+ 3 3 1 3 0.539895521E+00
+ 3 3 1 4 0.369973145E+02
+ 3 3 1 5 0.475188480E+19
+ 3 3 1 6 0.451041349E+18
+ 3 3 1 7 0.286436667E+19
+ 3 3 1 8 0.717137068E+13
+ 3 3 1 9 0.348924801E+02
+ 3 3 1 10 0.718897947E+02
+ 3 3 1 11 -0.318634119E+02
+ 3 3 1 12 0.302906819E+01
+ 3 3 1 13 0.595059421E-12
+ 3 3 1 14 0.767388423E+00
+ 3 3 1 15 0.232611577E+00
+ 3 3 1 16 0.804027991E-13
+ 3 3 1 17 0.191033003E+03
+ 3 3 1 18 0.497877084E+02
+ 3 3 1 19 0.132804741E+01
+ 3 3 2 1 0.317235579E+02
+ 3 3 2 2 0.191006923E+02
+ 3 3 2 3 0.514017995E+00
+ 3 3 2 4 0.396749699E+02
+ 3 3 2 5 0.516333130E+19
+ 3 3 2 6 0.438446151E+18
+ 3 3 2 7 0.320683832E+19
+ 3 3 2 8 0.683303230E+13
+ 3 3 2 9 0.341619469E+02
+ 3 3 2 10 0.738369168E+02
+ 3 3 2 11 -0.313641234E+02
+ 3 3 2 12 0.279782344E+01
+ 3 3 2 13 0.923618919E-12
+ 3 3 2 14 0.809557707E+00
+ 3 3 2 15 0.190442293E+00
+ 3 3 2 16 0.122102386E-12
+ 3 3 2 17 0.193131558E+03
+ 3 3 2 18 0.516889290E+02
+ 3 3 2 19 0.132649427E+01
+ 3 3 3 1 0.317586648E+02
+ 3 3 3 2 0.191221403E+02
+ 3 3 3 3 0.484225576E+00
+ 3 3 3 4 0.423385426E+02
+ 3 3 3 5 0.559649015E+19
+ 3 3 3 6 0.415407742E+18
+ 3 3 3 7 0.361986152E+19
+ 3 3 3 8 0.636706791E+13
+ 3 3 3 9 0.333479602E+02
+ 3 3 3 10 0.756865028E+02
+ 3 3 3 11 -0.308100265E+02
+ 3 3 3 12 0.253793375E+01
+ 3 3 3 13 0.151986729E-11
+ 3 3 3 14 0.849312396E+00
+ 3 3 3 15 0.150687604E+00
+ 3 3 3 16 0.194578057E-12
+ 3 3 3 17 0.194712318E+03
+ 3 3 3 18 0.534937136E+02
+ 3 3 3 19 0.132612444E+01
+ 3 3 4 1 0.317923161E+02
+ 3 3 4 2 0.191439002E+02
+ 3 3 4 3 0.453310615E+00
+ 3 3 4 4 0.450766447E+02
+ 3 3 4 5 0.604916538E+19
+ 3 3 4 6 0.378996321E+18
+ 3 3 4 7 0.396541830E+19
+ 3 3 4 8 0.621121371E+13
+ 3 3 4 9 0.323732737E+02
+ 3 3 4 10 0.774499184E+02
+ 3 3 4 11 -0.301445698E+02
+ 3 3 4 12 0.222870392E+01
+ 3 3 4 13 0.273766425E-11
+ 3 3 4 14 0.882787585E+00
+ 3 3 4 15 0.117212415E+00
+ 3 3 4 16 0.341458729E-12
+ 3 3 4 17 0.196028206E+03
+ 3 3 4 18 0.555139790E+02
+ 3 3 4 19 0.132651671E+01
+ 3 4 1 1 0.316896689E+02
+ 3 4 1 2 0.190810656E+02
+ 3 4 1 3 0.576624264E+00
+ 3 4 1 4 0.371440933E+02
+ 3 4 1 5 0.480439848E+19
+ 3 4 1 6 0.445352944E+18
+ 3 4 1 7 0.251464409E+19
+ 3 4 1 8 0.869719575E+13
+ 3 4 1 9 0.347353889E+02
+ 3 4 1 10 0.718794822E+02
+ 3 4 1 11 -0.317641843E+02
+ 3 4 1 12 0.297120460E+01
+ 3 4 1 13 0.546928533E-11
+ 3 4 1 14 0.765327167E+00
+ 3 4 1 15 0.234672833E+00
+ 3 4 1 16 0.714368954E-12
+ 3 4 1 17 0.190206474E+03
+ 3 4 1 18 0.497059250E+02
+ 3 4 1 19 0.133445344E+01
+ 3 4 2 1 0.317258904E+02
+ 3 4 2 2 0.191021106E+02
+ 3 4 2 3 0.549327738E+00
+ 3 4 2 4 0.398255288E+02
+ 3 4 2 5 0.521877730E+19
+ 3 4 2 6 0.427574406E+18
+ 3 4 2 7 0.276066543E+19
+ 3 4 2 8 0.858227330E+13
+ 3 4 2 9 0.340013470E+02
+ 3 4 2 10 0.738268758E+02
+ 3 4 2 11 -0.312639217E+02
+ 3 4 2 12 0.273742523E+01
+ 3 4 2 13 0.819914010E-11
+ 3 4 2 14 0.806403404E+00
+ 3 4 2 15 0.193596596E+00
+ 3 4 2 16 0.104751912E-11
+ 3 4 2 17 0.192600751E+03
+ 3 4 2 18 0.517484953E+02
+ 3 4 2 19 0.133355748E+01
+ 3 4 3 1 0.317605720E+02
+ 3 4 3 2 0.191234324E+02
+ 3 4 3 3 0.518013479E+00
+ 3 4 3 4 0.425245073E+02
+ 3 4 3 5 0.564342357E+19
+ 3 4 3 6 0.409117279E+18
+ 3 4 3 7 0.317038726E+19
+ 3 4 3 8 0.796495407E+13
+ 3 4 3 9 0.331521996E+02
+ 3 4 3 10 0.756767069E+02
+ 3 4 3 11 -0.306850211E+02
+ 3 4 3 12 0.246717844E+01
+ 3 4 3 13 0.132780284E-10
+ 3 4 3 14 0.845749215E+00
+ 3 4 3 15 0.154250785E+00
+ 3 4 3 16 0.165078940E-11
+ 3 4 3 17 0.193646535E+03
+ 3 4 3 18 0.534250469E+02
+ 3 4 3 19 0.133138586E+01
+ 3 4 4 1 0.317938841E+02
+ 3 4 4 2 0.191450111E+02
+ 3 4 4 3 0.483797066E+00
+ 3 4 4 4 0.452623763E+02
+ 3 4 4 5 0.609080360E+19
+ 3 4 4 6 0.369633250E+18
+ 3 4 4 7 0.346007468E+19
+ 3 4 4 8 0.784068719E+13
+ 3 4 4 9 0.321779688E+02
+ 3 4 4 10 0.774403451E+02
+ 3 4 4 11 -0.300221835E+02
+ 3 4 4 12 0.215578530E+01
+ 3 4 4 13 0.227356937E-10
+ 3 4 4 14 0.880868175E+00
+ 3 4 4 15 0.119131825E+00
+ 3 4 4 16 0.276184313E-11
+ 3 4 4 17 0.195224215E+03
+ 3 4 4 18 0.554067623E+02
+ 3 4 4 19 0.133083391E+01
+ 4 1 1 1 0.317686152E+02
+ 4 1 1 2 0.190955844E+02
+ 4 1 1 3 0.438604960E+00
+ 4 1 1 4 0.398206229E+02
+ 4 1 1 5 0.490779700E+19
+ 4 1 1 6 0.406400979E+18
+ 4 1 1 7 0.295943179E+19
+ 4 1 1 8 0.834134835E+13
+ 4 1 1 9 0.358651960E+02
+ 4 1 1 10 0.756858189E+02
+ 4 1 1 11 -0.326102500E+02
+ 4 1 1 12 0.325494600E+01
+ 4 1 1 13 0.106531400E-14
+ 4 1 1 14 0.781882900E+00
+ 4 1 1 15 0.218117100E+00
+ 4 1 1 16 0.222025800E-15
+ 4 1 1 17 0.202845100E+03
+ 4 1 1 18 0.518863200E+02
+ 4 1 1 19 0.132516940E+01
+ 4 1 2 1 0.318064535E+02
+ 4 1 2 2 0.191178446E+02
+ 4 1 2 3 0.417872028E+00
+ 4 1 2 4 0.426099736E+02
+ 4 1 2 5 0.535514602E+19
+ 4 1 2 6 0.393036953E+18
+ 4 1 2 7 0.339458251E+19
+ 4 1 2 8 0.768896413E+13
+ 4 1 2 9 0.351246181E+02
+ 4 1 2 10 0.777345917E+02
+ 4 1 2 11 -0.321013111E+02
+ 4 1 2 12 0.302330700E+01
+ 4 1 2 13 0.183786437E-14
+ 4 1 2 14 0.826879622E+00
+ 4 1 2 15 0.173120378E+00
+ 4 1 2 16 0.365949541E-15
+ 4 1 2 17 0.204751825E+03
+ 4 1 2 18 0.536509719E+02
+ 4 1 2 19 0.132531113E+01
+ 4 1 3 1 0.318425955E+02
+ 4 1 3 2 0.191403981E+02
+ 4 1 3 3 0.396346724E+00
+ 4 1 3 4 0.453657349E+02
+ 4 1 3 5 0.582005063E+19
+ 4 1 3 6 0.365577379E+18
+ 4 1 3 7 0.373080143E+19
+ 4 1 3 8 0.753633620E+13
+ 4 1 3 9 0.343150290E+02
+ 4 1 3 10 0.796807639E+02
+ 4 1 3 11 -0.315483026E+02
+ 4 1 3 12 0.276672643E+01
+ 4 1 3 13 0.332166481E-14
+ 4 1 3 14 0.862454652E+00
+ 4 1 3 15 0.137545348E+00
+ 4 1 3 16 0.635733189E-15
+ 4 1 3 17 0.208022347E+03
+ 4 1 3 18 0.562795359E+02
+ 4 1 3 19 0.132535196E+01
+ 4 1 4 1 0.318771259E+02
+ 4 1 4 2 0.191630171E+02
+ 4 1 4 3 0.369330084E+00
+ 4 1 4 4 0.482008515E+02
+ 4 1 4 5 0.630320018E+19
+ 4 1 4 6 0.379829463E+18
+ 4 1 4 7 0.432114649E+19
+ 4 1 4 8 0.659816173E+13
+ 4 1 4 9 0.333354373E+02
+ 4 1 4 10 0.815362888E+02
+ 4 1 4 11 -0.308774400E+02
+ 4 1 4 12 0.245799730E+01
+ 4 1 4 13 0.682153300E-14
+ 4 1 4 14 0.898932300E+00
+ 4 1 4 15 0.101067700E+00
+ 4 1 4 16 0.123731200E-14
+ 4 1 4 17 0.208270900E+03
+ 4 1 4 18 0.579217400E+02
+ 4 1 4 19 0.132566938E+01
+ 4 2 1 1 0.317705403E+02
+ 4 2 1 2 0.190966714E+02
+ 4 2 1 3 0.469746941E+00
+ 4 2 1 4 0.399143565E+02
+ 4 2 1 5 0.492312047E+19
+ 4 2 1 6 0.403249548E+18
+ 4 2 1 7 0.296805636E+19
+ 4 2 1 8 0.836816673E+13
+ 4 2 1 9 0.357642545E+02
+ 4 2 1 10 0.756786110E+02
+ 4 2 1 11 -0.325471876E+02
+ 4 2 1 12 0.321706686E+01
+ 4 2 1 13 0.184174232E-13
+ 4 2 1 14 0.780392837E+00
+ 4 2 1 15 0.219607163E+00
+ 4 2 1 16 0.350812921E-14
+ 4 2 1 17 0.202261684E+03
+ 4 2 1 18 0.518358723E+02
+ 4 2 1 19 0.132342753E+01
+ 4 2 2 1 0.318080265E+02
+ 4 2 2 2 0.191188433E+02
+ 4 2 2 3 0.447941097E+00
+ 4 2 2 4 0.427158768E+02
+ 4 2 2 5 0.536658602E+19
+ 4 2 2 6 0.394282480E+18
+ 4 2 2 7 0.338525439E+19
+ 4 2 2 8 0.776873463E+13
+ 4 2 2 9 0.350116968E+02
+ 4 2 2 10 0.777275736E+02
+ 4 2 2 11 -0.320305726E+02
+ 4 2 2 12 0.298112422E+01
+ 4 2 2 13 0.305048790E-13
+ 4 2 2 14 0.824080983E+00
+ 4 2 2 15 0.175919017E+00
+ 4 2 2 16 0.558063206E-14
+ 4 2 2 17 0.203672237E+03
+ 4 2 2 18 0.535493303E+02
+ 4 2 2 19 0.132334058E+01
+ 4 2 3 1 0.318439151E+02
+ 4 2 3 2 0.191412801E+02
+ 4 2 3 3 0.424360270E+00
+ 4 2 3 4 0.454844614E+02
+ 4 2 3 5 0.583063082E+19
+ 4 2 3 6 0.360962205E+18
+ 4 2 3 7 0.372376808E+19
+ 4 2 3 8 0.759797864E+13
+ 4 2 3 9 0.341894556E+02
+ 4 2 3 10 0.796739171E+02
+ 4 2 3 11 -0.314694290E+02
+ 4 2 3 12 0.272002667E+01
+ 4 2 3 13 0.529814349E-13
+ 4 2 3 14 0.860374106E+00
+ 4 2 3 15 0.139625894E+00
+ 4 2 3 16 0.934614619E-14
+ 4 2 3 17 0.207043698E+03
+ 4 2 3 18 0.561502209E+02
+ 4 2 3 19 0.132373321E+01
+ 4 2 4 1 0.318782344E+02
+ 4 2 4 2 0.191638870E+02
+ 4 2 4 3 0.398152680E+00
+ 4 2 4 4 0.483347396E+02
+ 4 2 4 5 0.631456877E+19
+ 4 2 4 6 0.361516596E+18
+ 4 2 4 7 0.426905350E+19
+ 4 2 4 8 0.681066530E+13
+ 4 2 4 9 0.331948580E+02
+ 4 2 4 10 0.815295976E+02
+ 4 2 4 11 -0.307880784E+02
+ 4 2 4 12 0.240677960E+01
+ 4 2 4 13 0.103455321E-12
+ 4 2 4 14 0.893481146E+00
+ 4 2 4 15 0.106518854E+00
+ 4 2 4 16 0.173764851E-13
+ 4 2 4 17 0.208024369E+03
+ 4 2 4 18 0.582061440E+02
+ 4 2 4 19 0.132467473E+01
+ 4 3 1 1 0.317727353E+02
+ 4 3 1 2 0.190979079E+02
+ 4 3 1 3 0.502636336E+00
+ 4 3 1 4 0.400240561E+02
+ 4 3 1 5 0.494717330E+19
+ 4 3 1 6 0.430892704E+18
+ 4 3 1 7 0.300325293E+19
+ 4 3 1 8 0.827972401E+13
+ 4 3 1 9 0.356461510E+02
+ 4 3 1 10 0.756702072E+02
+ 4 3 1 11 -0.324733848E+02
+ 4 3 1 12 0.317276625E+01
+ 4 3 1 13 0.234539991E-12
+ 4 3 1 14 0.778968513E+00
+ 4 3 1 15 0.221031487E+00
+ 4 3 1 16 0.413328033E-13
+ 4 3 1 17 0.201693007E+03
+ 4 3 1 18 0.517846827E+02
+ 4 3 1 19 0.132318890E+01
+ 4 3 2 1 0.318098841E+02
+ 4 3 2 2 0.191200033E+02
+ 4 3 2 3 0.480329435E+00
+ 4 3 2 4 0.428368699E+02
+ 4 3 2 5 0.539364534E+19
+ 4 3 2 6 0.414490713E+18
+ 4 3 2 7 0.335875576E+19
+ 4 3 2 8 0.794807231E+13
+ 4 3 2 9 0.348825211E+02
+ 4 3 2 10 0.777193911E+02
+ 4 3 2 11 -0.319495224E+02
+ 4 3 2 12 0.293299873E+01
+ 4 3 2 13 0.375492953E-12
+ 4 3 2 14 0.820043829E+00
+ 4 3 2 15 0.179956171E+00
+ 4 3 2 16 0.638159472E-13
+ 4 3 2 17 0.203445979E+03
+ 4 3 2 18 0.537531783E+02
+ 4 3 2 19 0.132433631E+01
+ 4 3 3 1 0.318454236E+02
+ 4 3 3 2 0.191422749E+02
+ 4 3 3 3 0.453684286E+00
+ 4 3 3 4 0.456189465E+02
+ 4 3 3 5 0.585563779E+19
+ 4 3 3 6 0.386754791E+18
+ 4 3 3 7 0.375731168E+19
+ 4 3 3 8 0.754383152E+13
+ 4 3 3 9 0.340469876E+02
+ 4 3 3 10 0.796659342E+02
+ 4 3 3 11 -0.313802742E+02
+ 4 3 3 12 0.266671342E+01
+ 4 3 3 13 0.625487105E-12
+ 4 3 3 14 0.858665901E+00
+ 4 3 3 15 0.141334099E+00
+ 4 3 3 16 0.103018131E-12
+ 4 3 3 17 0.206281704E+03
+ 4 3 3 18 0.560548608E+02
+ 4 3 3 19 0.132480091E+01
+ 4 3 4 1 0.318794719E+02
+ 4 3 4 2 0.191648015E+02
+ 4 3 4 3 0.426523746E+00
+ 4 3 4 4 0.484778916E+02
+ 4 3 4 5 0.633466989E+19
+ 4 3 4 6 0.364052665E+18
+ 4 3 4 7 0.420815222E+19
+ 4 3 4 8 0.706148171E+13
+ 4 3 4 9 0.330439046E+02
+ 4 3 4 10 0.815217962E+02
+ 4 3 4 11 -0.306942640E+02
+ 4 3 4 12 0.234964062E+01
+ 4 3 4 13 0.115469822E-11
+ 4 3 4 14 0.890724381E+00
+ 4 3 4 15 0.109275619E+00
+ 4 3 4 16 0.182551471E-12
+ 4 3 4 17 0.207301995E+03
+ 4 3 4 18 0.581835414E+02
+ 4 3 4 19 0.132511060E+01
+ 4 4 1 1 0.317752201E+02
+ 4 4 1 2 0.190993447E+02
+ 4 4 1 3 0.538110719E+00
+ 4 4 1 4 0.401533823E+02
+ 4 4 1 5 0.500158943E+19
+ 4 4 1 6 0.430381399E+18
+ 4 4 1 7 0.261456349E+19
+ 4 4 1 8 0.101665648E+14
+ 4 4 1 9 0.355070268E+02
+ 4 4 1 10 0.756604091E+02
+ 4 4 1 11 -0.323863387E+02
+ 4 4 1 12 0.312068814E+01
+ 4 4 1 13 0.228103651E-11
+ 4 4 1 14 0.776178278E+00
+ 4 4 1 15 0.223821722E+00
+ 4 4 1 16 0.376743011E-12
+ 4 4 1 17 0.200635817E+03
+ 4 4 1 18 0.516983030E+02
+ 4 4 1 19 0.133011124E+01
+ 4 4 2 1 0.318119935E+02
+ 4 4 2 2 0.191212873E+02
+ 4 4 2 3 0.513708469E+00
+ 4 4 2 4 0.429721604E+02
+ 4 4 2 5 0.544632648E+19
+ 4 4 2 6 0.404527599E+18
+ 4 4 2 7 0.288168339E+19
+ 4 4 2 8 0.100181898E+14
+ 4 4 2 9 0.347376905E+02
+ 4 4 2 10 0.777098509E+02
+ 4 4 2 11 -0.318596499E+02
+ 4 4 2 12 0.287804059E+01
+ 4 4 2 13 0.350824862E-11
+ 4 4 2 14 0.817249066E+00
+ 4 4 2 15 0.182750934E+00
+ 4 4 2 16 0.562708328E-12
+ 4 4 2 17 0.203398908E+03
+ 4 4 2 18 0.539245198E+02
+ 4 4 2 19 0.133079192E+01
+ 4 4 3 1 0.318471524E+02
+ 4 4 3 2 0.191434465E+02
+ 4 4 3 3 0.485737746E+00
+ 4 4 3 4 0.457812812E+02
+ 4 4 3 5 0.589828289E+19
+ 4 4 3 6 0.388442614E+18
+ 4 4 3 7 0.331412313E+19
+ 4 4 3 8 0.929236726E+13
+ 4 4 3 9 0.338753456E+02
+ 4 4 3 10 0.796566268E+02
+ 4 4 3 11 -0.312719740E+02
+ 4 4 3 12 0.260337163E+01
+ 4 4 3 13 0.573412555E-11
+ 4 4 3 14 0.855598255E+00
+ 4 4 3 15 0.144401745E+00
+ 4 4 3 16 0.891048470E-12
+ 4 4 3 17 0.205233491E+03
+ 4 4 3 18 0.559699769E+02
+ 4 4 3 19 0.132914770E+01
+ 4 4 4 1 0.318808718E+02
+ 4 4 4 2 0.191658231E+02
+ 4 4 4 3 0.455960115E+00
+ 4 4 4 4 0.486390021E+02
+ 4 4 4 5 0.636940094E+19
+ 4 4 4 6 0.356651602E+18
+ 4 4 4 7 0.369723920E+19
+ 4 4 4 8 0.887367983E+13
+ 4 4 4 9 0.328736983E+02
+ 4 4 4 10 0.815127004E+02
+ 4 4 4 11 -0.305893962E+02
+ 4 4 4 12 0.228430212E+01
+ 4 4 4 13 0.100055378E-10
+ 4 4 4 14 0.888767175E+00
+ 4 4 4 15 0.111232825E+00
+ 4 4 4 16 0.150679787E-11
+ 4 4 4 17 0.206444236E+03
+ 4 4 4 18 0.580703929E+02
+ 4 4 4 19 0.132808789E+01
diff --git a/test/test.ccl b/test/test.ccl
new file mode 100644
index 0000000..f6277fe
--- /dev/null
+++ b/test/test.ccl
@@ -0,0 +1 @@
+RELTOL 1e-14
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-03 10:39:38 +0000