diff --git a/fpp-ptc-sandbox/code/CMakeLists.txt b/fpp-ptc-sandbox/code/CMakeLists.txt index 8887750..53cc249 100644 --- a/fpp-ptc-sandbox/code/CMakeLists.txt +++ b/fpp-ptc-sandbox/code/CMakeLists.txt @@ -33,6 +33,7 @@ set (EXE_SPECS cmake_files/cmake.z_canonize_matt_fake_maps cmake_files/cmake.z_resonance cmake_files/cmake.z_spin1 + cmake_files/cmake.z_spin_res ) diff --git a/fpp-ptc-sandbox/code/EQUATION.TXT b/fpp-ptc-sandbox/code/EQUATION.TXT new file mode 100644 index 0000000..96ce0b9 --- /dev/null +++ b/fpp-ptc-sandbox/code/EQUATION.TXT @@ -0,0 +1,20 @@ + + etall 1, NO = 2, NV = 6, INA = 56 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 6 + 0 -0.1604272270583351E-13 0 0 0 0 0 0 + 1 3.532841697202252 1 0 0 0 0 0 + 1 0.4306530717865838 0 1 0 0 0 0 + -3 0.000000000000000 0 0 0 0 0 0 + + etall 1, NO = 2, NV = 6, INA = 57 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 6 + 0 0.1842970220877760E-13 0 0 0 0 0 0 + 1 -2.342063456623922 1 0 0 0 0 0 + 1 -3.255712043070225 0 1 0 0 0 0 + -3 0.000000000000000 0 0 0 0 0 0 diff --git a/fpp-ptc-sandbox/code/Makefile b/fpp-ptc-sandbox/code/Makefile index b9f3122..e9e94b1 100644 --- a/fpp-ptc-sandbox/code/Makefile +++ b/fpp-ptc-sandbox/code/Makefile @@ -428,6 +428,19 @@ z_spin1-exe/fast: $(MAKE) $(MAKESILENT) -f CMakeFiles/z_spin1-exe.dir/build.make CMakeFiles/z_spin1-exe.dir/build .PHONY : z_spin1-exe/fast +#============================================================================= +# Target rules for targets named z_spin_res-exe + +# Build rule for target. +z_spin_res-exe: cmake_check_build_system + $(MAKE) $(MAKESILENT) -f CMakeFiles/Makefile2 z_spin_res-exe +.PHONY : z_spin_res-exe + +# fast build rule for target. +z_spin_res-exe/fast: + $(MAKE) $(MAKESILENT) -f CMakeFiles/z_spin_res-exe.dir/build.make CMakeFiles/z_spin_res-exe.dir/build +.PHONY : z_spin_res-exe/fast + normal_simple.o: normal_simple.f90.o .PHONY : normal_simple.o @@ -860,6 +873,30 @@ z_spin1.f90.s: $(MAKE) $(MAKESILENT) -f CMakeFiles/z_spin1-exe.dir/build.make CMakeFiles/z_spin1-exe.dir/z_spin1.f90.s .PHONY : z_spin1.f90.s +z_spin_res.o: z_spin_res.f90.o +.PHONY : z_spin_res.o + +# target to build an object file +z_spin_res.f90.o: + $(MAKE) $(MAKESILENT) -f CMakeFiles/z_spin_res-exe.dir/build.make CMakeFiles/z_spin_res-exe.dir/z_spin_res.f90.o +.PHONY : z_spin_res.f90.o + +z_spin_res.i: z_spin_res.f90.i +.PHONY : z_spin_res.i + +# target to preprocess a source file +z_spin_res.f90.i: + $(MAKE) $(MAKESILENT) -f CMakeFiles/z_spin_res-exe.dir/build.make CMakeFiles/z_spin_res-exe.dir/z_spin_res.f90.i +.PHONY : z_spin_res.f90.i + +z_spin_res.s: z_spin_res.f90.s +.PHONY : z_spin_res.s + +# target to generate assembly for a file +z_spin_res.f90.s: + $(MAKE) $(MAKESILENT) -f CMakeFiles/z_spin_res-exe.dir/build.make CMakeFiles/z_spin_res-exe.dir/z_spin_res.f90.s +.PHONY : z_spin_res.f90.s + z_sub_i.o: z_sub_i.f90.o .PHONY : z_sub_i.o @@ -1030,6 +1067,7 @@ help: @echo "... z_shift-exe" @echo "... z_simple_operations-exe" @echo "... z_spin1-exe" + @echo "... z_spin_res-exe" @echo "... z_sub_i-exe" @echo "... z_sub_j-exe" @echo "... z_track_da_tpsa-exe" @@ -1090,6 +1128,9 @@ help: @echo "... z_spin1.o" @echo "... z_spin1.i" @echo "... z_spin1.s" + @echo "... z_spin_res.o" + @echo "... z_spin_res.i" + @echo "... z_spin_res.s" @echo "... z_sub_i.o" @echo "... z_sub_i.i" @echo "... z_sub_i.s" diff --git a/fpp-ptc-sandbox/code/cmake_files/cmake.z_spin_res b/fpp-ptc-sandbox/code/cmake_files/cmake.z_spin_res new file mode 100644 index 0000000..e65fb9e --- /dev/null +++ b/fpp-ptc-sandbox/code/cmake_files/cmake.z_spin_res @@ -0,0 +1,14 @@ +set (EXENAME z_spin_res) + +FILE (GLOB SRC_FILES z_spin_res.f90) + +set (INC_DIRS +) + +set (LINK_LIBS + forest + bsim + bmad + sim_utils + ${ACC_BMAD_LINK_LIBS} +) diff --git a/fpp-ptc-sandbox/code/fit_tune.txt b/fpp-ptc-sandbox/code/fit_tune.txt index 24646cf..7ce82ff 100644 --- a/fpp-ptc-sandbox/code/fit_tune.txt +++ b/fpp-ptc-sandbox/code/fit_tune.txt @@ -1,11 +1,13 @@ select layout 1 - L MAX - 1.0000000000000000E-002 - CUTTING ALGORITHM - 2 - LIMIT FOR CUTTING - 10000 10000 - MISALIGN EVERYTHING - 0.001 0.001 0.001 0.001 0.001 0.001 5 + set families + 2 + 1 qf1 + 2 1 + 1 qd1 + 2 2 + FITTUNE + 0.0000000001 + 0.367269834767751 0.404042919 !2293877 !6057 0.336 0.218992743820551 + deallocate families return diff --git a/fpp-ptc-sandbox/code/plot.dat b/fpp-ptc-sandbox/code/plot.dat new file mode 100644 index 0000000..e69de29 diff --git a/fpp-ptc-sandbox/code/spinmap.txt b/fpp-ptc-sandbox/code/spinmap.txt index 154d263..d6322e9 100644 --- a/fpp-ptc-sandbox/code/spinmap.txt +++ b/fpp-ptc-sandbox/code/spinmap.txt @@ -5,41 +5,41 @@ I COEFFICIENT ORDER EXPONENTS NO = 3 NV = 4 - 0 0.1583591029206294E-03 0.000000000000000 0 0 0 0 - 1 -0.7133082123844400 0.000000000000000 1 0 0 0 - 1 7.760252162396446 0.000000000000000 0 1 0 0 - 1 -0.2061942199606380 0.000000000000000 0 0 1 0 - 1 -0.6710032459762710 0.000000000000000 0 0 0 1 - 2 11.08435544404162 0.000000000000000 2 0 0 0 - 2 147.5115941047682 0.000000000000000 1 1 0 0 - 2 2788.949729070081 0.000000000000000 0 2 0 0 - 2 -21.14422257782249 0.000000000000000 1 0 1 0 - 2 -399.7035173568318 0.000000000000000 0 1 1 0 - 2 -214.6804470606113 0.000000000000000 1 0 0 1 - 2 -2378.403918890698 0.000000000000000 0 1 0 1 - 2 29.00610402969301 0.000000000000000 0 0 2 0 - 2 190.2349421247553 0.000000000000000 0 0 1 1 - 2 272.1124669076372 0.000000000000000 0 0 0 2 - 3 1287.165015612577 0.000000000000000 3 0 0 0 - 3 15920.05359617360 0.000000000000000 2 1 0 0 - 3 -67574.32808238995 0.000000000000000 1 2 0 0 - 3 557763.4052516703 0.000000000000000 0 3 0 0 - 3 -1789.229112963190 0.000000000000000 2 0 1 0 - 3 22544.70047301857 0.000000000000000 1 1 1 0 - 3 63510.17314151523 0.000000000000000 0 2 1 0 - 3 -3454.628075444478 0.000000000000000 2 0 0 1 - 3 19466.55113407743 0.000000000000000 1 1 0 1 - 3 869586.4712777771 0.000000000000000 0 2 0 1 - 3 7529.324986373725 0.000000000000000 1 0 2 0 - 3 83073.47401893517 0.000000000000000 0 1 2 0 - 3 8187.948348912192 0.000000000000000 1 0 1 1 - 3 63730.40584072565 0.000000000000000 0 1 1 1 - 3 105249.1316853605 0.000000000000000 1 0 0 2 - 3 478263.1809695730 0.000000000000000 0 1 0 2 - 3 -477.5206034254909 0.000000000000000 0 0 3 0 - 3 -3158.707728920645 0.000000000000000 0 0 2 1 - 3 -33913.71680906291 0.000000000000000 0 0 1 2 - 3 21146.59285884549 0.000000000000000 0 0 0 3 + 0 0.4867812994196517E-04 0.000000000000000 0 0 0 0 + 1 -0.7009790553862907 0.000000000000000 1 0 0 0 + 1 7.904434314383582 0.000000000000000 0 1 0 0 + 1 0.1066153604074146 0.000000000000000 0 0 1 0 + 1 0.2299414479930226E-01 0.000000000000000 0 0 0 1 + 2 11.79497559214498 0.000000000000000 2 0 0 0 + 2 158.3575925664820 0.000000000000000 1 1 0 0 + 2 2673.674442760778 0.000000000000000 0 2 0 0 + 2 10.48669129302925 0.000000000000000 1 0 1 0 + 2 -15.65130750712133 0.000000000000000 0 1 1 0 + 2 -99.02923431614474 0.000000000000000 1 0 0 1 + 2 -1180.000693870767 0.000000000000000 0 1 0 1 + 2 38.41036578566355 0.000000000000000 0 0 2 0 + 2 27.97031603167960 0.000000000000000 0 0 1 1 + 2 49.99458464053005 0.000000000000000 0 0 0 2 + 3 1427.339372465137 0.000000000000000 3 0 0 0 + 3 16816.85430015757 0.000000000000000 2 1 0 0 + 3 144.8299485483585 0.000000000000000 1 2 0 0 + 3 561954.5739163718 0.000000000000000 0 3 0 0 + 3 -2037.929988219854 0.000000000000000 2 0 1 0 + 3 -13901.70576711066 0.000000000000000 1 1 1 0 + 3 35475.75559255313 0.000000000000000 0 2 1 0 + 3 -7234.501404451950 0.000000000000000 2 0 0 1 + 3 132146.4025601656 0.000000000000000 1 1 0 1 + 3 553799.3263696884 0.000000000000000 0 2 0 1 + 3 6495.962490463056 0.000000000000000 1 0 2 0 + 3 27483.89047271314 0.000000000000000 0 1 2 0 + 3 7053.617749403057 0.000000000000000 1 0 1 1 + 3 214200.6547345958 0.000000000000000 0 1 1 1 + 3 125616.6591170086 0.000000000000000 1 0 0 2 + 3 1390133.288887875 0.000000000000000 0 1 0 2 + 3 47.62893353950619 0.000000000000000 0 0 3 0 + 3 6261.804641697765 0.000000000000000 0 0 2 1 + 3 -68996.22011253830 0.000000000000000 0 0 1 2 + 3 -79792.34464947744 0.000000000000000 0 0 0 3 -35 0.000000000000000 0.000000000000000 0 0 0 0 1, NO = 3, NV = 4, INA = 205 @@ -47,41 +47,41 @@ I COEFFICIENT ORDER EXPONENTS NO = 3 NV = 4 - 0 -0.1888090075450990E-04 0.000000000000000 0 0 0 0 - 1 -0.7203636799033739E-01 0.000000000000000 1 0 0 0 - 1 -0.6361260999793767 0.000000000000000 0 1 0 0 - 1 -0.3093068624609216E-02 0.000000000000000 0 0 1 0 - 1 0.5188745160021657E-01 0.000000000000000 0 0 0 1 - 2 -0.8730028407742259 0.000000000000000 2 0 0 0 - 2 4.849108186508680 0.000000000000000 1 1 0 0 - 2 203.5862568544221 0.000000000000000 0 2 0 0 - 2 1.218891753462678 0.000000000000000 1 0 1 0 - 2 -25.06647677039867 0.000000000000000 0 1 1 0 - 2 27.50503236210943 0.000000000000000 1 0 0 1 - 2 -315.2853072006382 0.000000000000000 0 1 0 1 - 2 2.602136919108817 0.000000000000000 0 0 2 0 - 2 0.8197056136191350 0.000000000000000 0 0 1 1 - 2 -27.19275023767941 0.000000000000000 0 0 0 2 - 3 -86.47398763814783 0.000000000000000 3 0 0 0 - 3 736.2694927450048 0.000000000000000 2 1 0 0 - 3 -15080.66355770527 0.000000000000000 1 2 0 0 - 3 179101.2474619584 0.000000000000000 0 3 0 0 - 3 70.83190561331241 0.000000000000000 2 0 1 0 - 3 4895.090439509268 0.000000000000000 1 1 1 0 - 3 -31405.00678947513 0.000000000000000 0 2 1 0 - 3 812.7123304378706 0.000000000000000 2 0 0 1 - 3 -16646.44214646829 0.000000000000000 1 1 0 1 - 3 -89141.18929828267 0.000000000000000 0 2 0 1 - 3 -684.8463600691666 0.000000000000000 1 0 2 0 - 3 9165.502047414531 0.000000000000000 0 1 2 0 - 3 604.8412998372014 0.000000000000000 1 0 1 1 - 3 19901.42394506867 0.000000000000000 0 1 1 1 - 3 -9228.993731061710 0.000000000000000 1 0 0 2 - 3 219246.0293633156 0.000000000000000 0 1 0 2 - 3 -85.12190858984732 0.000000000000000 0 0 3 0 - 3 -1014.891958238282 0.000000000000000 0 0 2 1 - 3 -3524.614363378083 0.000000000000000 0 0 1 2 - 3 1799.255703016682 0.000000000000000 0 0 0 3 + 0 0.4079952402053416E-05 0.000000000000000 0 0 0 0 + 1 -0.6950431993446177E-01 0.000000000000000 1 0 0 0 + 1 -0.6462030576832420 0.000000000000000 0 1 0 0 + 1 0.8046445421895778E-02 0.000000000000000 0 0 1 0 + 1 -0.2046721326935349E-01 0.000000000000000 0 0 0 1 + 2 -0.2785683768333875 0.000000000000000 2 0 0 0 + 2 2.099235149605022 0.000000000000000 1 1 0 0 + 2 179.8213181320455 0.000000000000000 0 2 0 0 + 2 -2.619366061310086 0.000000000000000 1 0 1 0 + 2 14.47203346047138 0.000000000000000 0 1 1 0 + 2 6.214871289366293 0.000000000000000 1 0 0 1 + 2 -94.91918593386094 0.000000000000000 0 1 0 1 + 2 1.986744420528817 0.000000000000000 0 0 2 0 + 2 -7.568461716908912 0.000000000000000 0 0 1 1 + 2 8.506437651304697 0.000000000000000 0 0 0 2 + 3 -85.07192985189931 0.000000000000000 3 0 0 0 + 3 777.0739203572730 0.000000000000000 2 1 0 0 + 3 -18800.90327233213 0.000000000000000 1 2 0 0 + 3 151245.0278252386 0.000000000000000 0 3 0 0 + 3 -17.23120377380008 0.000000000000000 2 0 1 0 + 3 4729.248840907421 0.000000000000000 1 1 1 0 + 3 -6236.397535203525 0.000000000000000 0 2 1 0 + 3 766.6717124172396 0.000000000000000 2 0 0 1 + 3 -1731.648292483020 0.000000000000000 1 1 0 1 + 3 -64564.33287938734 0.000000000000000 0 2 0 1 + 3 -1064.169138978218 0.000000000000000 1 0 2 0 + 3 7629.698042606490 0.000000000000000 0 1 2 0 + 3 136.1495662169655 0.000000000000000 1 0 1 1 + 3 6661.426994793539 0.000000000000000 0 1 1 1 + 3 -5793.228572617023 0.000000000000000 1 0 0 2 + 3 121071.5804719628 0.000000000000000 0 1 0 2 + 3 124.4945153667090 0.000000000000000 0 0 3 0 + 3 -1240.439126355940 0.000000000000000 0 0 2 1 + 3 -4608.512684806366 0.000000000000000 0 0 1 2 + 3 6940.331747395168 0.000000000000000 0 0 0 3 -35 0.000000000000000 0.000000000000000 0 0 0 0 1, NO = 3, NV = 4, INA = 204 @@ -89,41 +89,41 @@ I COEFFICIENT ORDER EXPONENTS NO = 3 NV = 4 - 0 -0.4126698713745131E-02 0.000000000000000 0 0 0 0 - 1 0.6945630430047658E-02 0.000000000000000 1 0 0 0 - 1 -0.4825859270554137 0.000000000000000 0 1 0 0 - 1 0.5621490866329165 0.000000000000000 0 0 1 0 - 1 3.237416273678870 0.000000000000000 0 0 0 1 - 2 5.053111573697564 0.000000000000000 2 0 0 0 - 2 -56.01404633628826 0.000000000000000 1 1 0 0 - 2 304.8213108038958 0.000000000000000 0 2 0 0 - 2 -4.416282678529211 0.000000000000000 1 0 1 0 - 2 199.7745096981538 0.000000000000000 0 1 1 0 - 2 -4.500624373991428 0.000000000000000 1 0 0 1 - 2 324.6999425254253 0.000000000000000 0 1 0 1 - 2 -8.160587998580747 0.000000000000000 0 0 2 0 - 2 -122.8527580157144 0.000000000000000 0 0 1 1 - 2 366.0036403511324 0.000000000000000 0 0 0 2 - 3 -94.51904188010491 0.000000000000000 3 0 0 0 - 3 -3954.555250550836 0.000000000000000 2 1 0 0 - 3 55038.11541286786 0.000000000000000 1 2 0 0 - 3 248677.4430619288 0.000000000000000 0 3 0 0 - 3 3674.287221807684 0.000000000000000 2 0 1 0 - 3 -624.0580368642341 0.000000000000000 1 1 1 0 - 3 338522.3058582309 0.000000000000000 0 2 1 0 - 3 -9327.433339453610 0.000000000000000 2 0 0 1 - 3 70405.01747205549 0.000000000000000 1 1 0 1 - 3 -1196057.979318736 0.000000000000000 0 2 0 1 - 3 -54.43857690995839 0.000000000000000 1 0 2 0 - 3 -13750.96098015015 0.000000000000000 0 1 2 0 - 3 -4651.988927235364 0.000000000000000 1 0 1 1 - 3 17887.72225623039 0.000000000000000 0 1 1 1 - 3 15566.21140084440 0.000000000000000 1 0 0 2 - 3 -9064.795633102127 0.000000000000000 0 1 0 2 - 3 3100.162059063046 0.000000000000000 0 0 3 0 - 3 -2969.401825660620 0.000000000000000 0 0 2 1 - 3 53180.69159010795 0.000000000000000 0 0 1 2 - 3 -246554.7157853577 0.000000000000000 0 0 0 3 + 0 -0.1881736176817633E-02 0.000000000000000 0 0 0 0 + 1 0.1643878192699564E-01 0.000000000000000 1 0 0 0 + 1 0.1493855966153149 0.000000000000000 0 1 0 0 + 1 -0.8084495364632538 0.000000000000000 0 0 1 0 + 1 2.026592039121844 0.000000000000000 0 0 0 1 + 2 -1.951964557166611 0.000000000000000 2 0 0 0 + 2 -37.66634032326170 0.000000000000000 1 1 0 0 + 2 -533.8580278189322 0.000000000000000 0 2 0 0 + 2 11.67739061655207 0.000000000000000 1 0 1 0 + 2 128.4727431544466 0.000000000000000 0 1 1 0 + 2 8.384769157953315 0.000000000000000 1 0 0 1 + 2 62.97219966619150 0.000000000000000 0 1 0 1 + 2 -2.032881058943547 0.000000000000000 0 0 2 0 + 2 20.12523675149727 0.000000000000000 0 0 1 1 + 2 -146.5222148685612 0.000000000000000 0 0 0 2 + 3 -316.3773522568532 0.000000000000000 3 0 0 0 + 3 -123.4292062587774 0.000000000000000 2 1 0 0 + 3 81739.39826328863 0.000000000000000 1 2 0 0 + 3 93251.98243610471 0.000000000000000 0 3 0 0 + 3 2552.217051829479 0.000000000000000 2 0 1 0 + 3 -11220.32318963180 0.000000000000000 1 1 1 0 + 3 211094.2433088830 0.000000000000000 0 2 1 0 + 3 10424.11056195601 0.000000000000000 2 0 0 1 + 3 56149.20469275284 0.000000000000000 1 1 0 1 + 3 1615081.217774376 0.000000000000000 0 2 0 1 + 3 1259.920024982692 0.000000000000000 1 0 2 0 + 3 1456.852971214741 0.000000000000000 0 1 2 0 + 3 352.9416413601475 0.000000000000000 1 0 1 1 + 3 -152411.5703108387 0.000000000000000 0 1 1 1 + 3 -28971.12394811244 0.000000000000000 1 0 0 2 + 3 -105101.4208378543 0.000000000000000 0 1 0 2 + 3 1934.339423639165 0.000000000000000 0 0 3 0 + 3 19259.04091146962 0.000000000000000 0 0 2 1 + 3 26930.09344611266 0.000000000000000 0 0 1 2 + 3 131107.1894905354 0.000000000000000 0 0 0 3 -35 0.000000000000000 0.000000000000000 0 0 0 0 1, NO = 3, NV = 4, INA = 203 @@ -131,41 +131,41 @@ I COEFFICIENT ORDER EXPONENTS NO = 3 NV = 4 - 0 -0.3601273821355091E-03 0.000000000000000 0 0 0 0 - 1 -0.2494691594887007E-01 0.000000000000000 1 0 0 0 - 1 -0.1058681631046830 0.000000000000000 0 1 0 0 - 1 -0.1982130454057530 0.000000000000000 0 0 1 0 - 1 0.6601028475669869 0.000000000000000 0 0 0 1 - 2 -4.376269461405133 0.000000000000000 2 0 0 0 - 2 11.01664518299118 0.000000000000000 1 1 0 0 - 2 -457.9230187108817 0.000000000000000 0 2 0 0 - 2 2.060325443216622 0.000000000000000 1 0 1 0 - 2 68.53168350076965 0.000000000000000 0 1 1 0 - 2 30.40484154667094 0.000000000000000 1 0 0 1 - 2 9.854935514329327 0.000000000000000 0 1 0 1 - 2 -6.727843738258257 0.000000000000000 0 0 2 0 - 2 -27.32116687289197 0.000000000000000 0 0 1 1 - 2 -156.3340257896515 0.000000000000000 0 0 0 2 - 3 -71.81469861741031 0.000000000000000 3 0 0 0 - 3 1160.533026190067 0.000000000000000 2 1 0 0 - 3 25460.86606623653 0.000000000000000 1 2 0 0 - 3 1897.817977684681 0.000000000000000 0 3 0 0 - 3 518.2278078062291 0.000000000000000 2 0 1 0 - 3 1736.023000627927 0.000000000000000 1 1 1 0 - 3 71987.97423937204 0.000000000000000 0 2 1 0 - 3 3147.931219162976 0.000000000000000 2 0 0 1 - 3 -26203.75767485784 0.000000000000000 1 1 0 1 - 3 350900.1819036597 0.000000000000000 0 2 0 1 - 3 148.0121694885408 0.000000000000000 1 0 2 0 - 3 -1232.760490002174 0.000000000000000 0 1 2 0 - 3 -162.3009674600773 0.000000000000000 1 0 1 1 - 3 -28614.95172149240 0.000000000000000 0 1 1 1 - 3 -9626.069464687000 0.000000000000000 1 0 0 2 - 3 86773.78587599265 0.000000000000000 0 1 0 2 - 3 654.5148985310880 0.000000000000000 0 0 3 0 - 3 5340.101751440514 0.000000000000000 0 0 2 1 - 3 7911.504523563914 0.000000000000000 0 0 1 2 - 3 56059.49444207293 0.000000000000000 0 0 0 3 + 0 0.2169309298279741E-03 0.000000000000000 0 0 0 0 + 1 0.1089831076000199E-02 0.000000000000000 1 0 0 0 + 1 -0.1340935961490033 0.000000000000000 0 1 0 0 + 1 -0.1612596266934681 0.000000000000000 0 0 1 0 + 1 -0.8356237888601501 0.000000000000000 0 0 0 1 + 2 -2.109642153719256 0.000000000000000 2 0 0 0 + 2 6.702969933420714 0.000000000000000 1 1 0 0 + 2 -177.9506367930843 0.000000000000000 0 2 0 0 + 2 -1.213411706449580 0.000000000000000 1 0 1 0 + 2 -71.23633107352856 0.000000000000000 0 1 1 0 + 2 -5.741433603522974 0.000000000000000 1 0 0 1 + 2 57.55400546403282 0.000000000000000 0 1 0 1 + 2 -1.765419941567711 0.000000000000000 0 0 2 0 + 2 13.73339712575856 0.000000000000000 0 0 1 1 + 2 -56.46272202393999 0.000000000000000 0 0 0 2 + 3 -26.12920442163664 0.000000000000000 3 0 0 0 + 3 2056.933264341425 0.000000000000000 2 1 0 0 + 3 1602.776374712306 0.000000000000000 1 2 0 0 + 3 -49011.23128974235 0.000000000000000 0 3 0 0 + 3 -1037.826133533944 0.000000000000000 2 0 1 0 + 3 6998.503672639559 0.000000000000000 1 1 1 0 + 3 -53169.98408589326 0.000000000000000 0 2 1 0 + 3 1957.465253320000 0.000000000000000 2 0 0 1 + 3 1277.912528791556 0.000000000000000 1 1 0 1 + 3 260159.9398063142 0.000000000000000 0 2 0 1 + 3 151.7248934854840 0.000000000000000 1 0 2 0 + 3 -1066.836334195409 0.000000000000000 0 1 2 0 + 3 -2251.648465699132 0.000000000000000 1 0 1 1 + 3 -30775.92467424635 0.000000000000000 0 1 1 1 + 3 -6506.364083884773 0.000000000000000 1 0 0 2 + 3 65774.41023640917 0.000000000000000 0 1 0 2 + 3 -481.4346121063184 0.000000000000000 0 0 3 0 + 3 386.0481962047043 0.000000000000000 0 0 2 1 + 3 -15340.11089393405 0.000000000000000 0 0 1 2 + 3 33257.14667513099 0.000000000000000 0 0 0 3 -35 0.000000000000000 0.000000000000000 0 0 0 0 No Spin Matrix Quaternion @@ -176,41 +176,41 @@ I COEFFICIENT ORDER EXPONENTS NO = 3 NV = 4 - 0 -0.2969128299686966 0.000000000000000 0 0 0 0 - 1 0.1683509718624019 0.000000000000000 1 0 0 0 - 1 3.686791826057608 0.000000000000000 0 1 0 0 - 1 0.2163044561922609E-01 0.000000000000000 0 0 1 0 - 1 0.6760006361666161E-01 0.000000000000000 0 0 0 1 - 2 14.72273476865321 0.000000000000000 2 0 0 0 - 2 -28.69129992621258 0.000000000000000 1 1 0 0 - 2 937.6946113845058 0.000000000000000 0 2 0 0 - 2 -2.145618232199157 0.000000000000000 1 0 1 0 - 2 42.19140043402259 0.000000000000000 0 1 1 0 - 2 -66.26003012375647 0.000000000000000 1 0 0 1 - 2 647.3761542316439 0.000000000000000 0 1 0 1 - 2 -12.69777140507599 0.000000000000000 0 0 2 0 - 2 -6.093424204488350 0.000000000000000 0 0 1 1 - 2 -12.60820781523897 0.000000000000000 0 0 0 2 - 3 219.0139495868201 0.000000000000000 3 0 0 0 - 3 888.7851266193368 0.000000000000000 2 1 0 0 - 3 -28102.60508271792 0.000000000000000 1 2 0 0 - 3 -432960.1804460597 0.000000000000000 0 3 0 0 - 3 -128.3263004465716 0.000000000000000 2 0 1 0 - 3 -6455.347203699675 0.000000000000000 1 1 1 0 - 3 38047.71909784441 0.000000000000000 0 2 1 0 - 3 -1970.413560641252 0.000000000000000 2 0 0 1 - 3 52511.62192976443 0.000000000000000 1 1 0 1 - 3 220936.5476906177 0.000000000000000 0 2 0 1 - 3 971.2670858686553 0.000000000000000 1 0 2 0 - 3 -18251.67779067983 0.000000000000000 0 1 2 0 - 3 -1453.478589491812 0.000000000000000 1 0 1 1 - 3 -26902.94646973653 0.000000000000000 0 1 1 1 - 3 21990.83431476812 0.000000000000000 1 0 0 2 - 3 -476238.1708054012 0.000000000000000 0 1 0 2 - 3 -25.12515866064648 0.000000000000000 0 0 3 0 - 3 1388.858803540460 0.000000000000000 0 0 2 1 - 3 2918.810554919819 0.000000000000000 0 0 1 2 - 3 -15611.90528092934 0.000000000000000 0 0 0 3 + 0 -0.2967996746715293 0.000000000000000 0 0 0 0 + 1 0.1561012093764574 0.000000000000000 1 0 0 0 + 1 3.722883276874034 0.000000000000000 0 1 0 0 + 1 0.7641205223954218E-03 0.000000000000000 0 0 1 0 + 1 0.4864449934371465E-01 0.000000000000000 0 0 0 1 + 2 13.21880030461130 0.000000000000000 2 0 0 0 + 2 -14.94621012856408 0.000000000000000 1 1 0 0 + 2 1007.362437331714 0.000000000000000 0 2 0 0 + 2 3.828320578099249 0.000000000000000 1 0 1 0 + 2 -51.71278156687806 0.000000000000000 0 1 1 0 + 2 -16.65630677712585 0.000000000000000 1 0 0 1 + 2 212.6376876499664 0.000000000000000 0 1 0 1 + 2 1.338473254453042 0.000000000000000 0 0 2 0 + 2 -10.32682115008832 0.000000000000000 0 0 1 1 + 2 -85.94626785873777 0.000000000000000 0 0 0 2 + 3 199.6521334388494 0.000000000000000 3 0 0 0 + 3 -515.7965012541285 0.000000000000000 2 1 0 0 + 3 -16162.99820050433 0.000000000000000 1 2 0 0 + 3 -335079.2856413586 0.000000000000000 0 3 0 0 + 3 109.8694266611851 0.000000000000000 2 0 1 0 + 3 -8124.952059859799 0.000000000000000 1 1 1 0 + 3 -5281.774484877198 0.000000000000000 0 2 1 0 + 3 -1588.876188464621 0.000000000000000 2 0 0 1 + 3 7835.633751390963 0.000000000000000 1 1 0 1 + 3 117859.4620509410 0.000000000000000 0 2 0 1 + 3 1568.816943522881 0.000000000000000 1 0 2 0 + 3 -21238.99041101434 0.000000000000000 0 1 2 0 + 3 -948.1986771149096 0.000000000000000 1 0 1 1 + 3 9601.635012672821 0.000000000000000 0 1 1 1 + 3 14406.69961249195 0.000000000000000 1 0 0 2 + 3 -231052.7544266776 0.000000000000000 0 1 0 2 + 3 -294.9104355931647 0.000000000000000 0 0 3 0 + 3 2497.396182369676 0.000000000000000 0 0 2 1 + 3 11638.41008922243 0.000000000000000 0 0 1 2 + 3 -22477.70617443228 0.000000000000000 0 0 0 3 -35 0.000000000000000 0.000000000000000 0 0 0 0 1, NO = 3, NV = 4, INA = 185 @@ -218,41 +218,41 @@ I COEFFICIENT ORDER EXPONENTS NO = 3 NV = 4 - 0 -0.1639748677798469E-01 0.000000000000000 0 0 0 0 - 1 -0.3008061215458598E-01 0.000000000000000 1 0 0 0 - 1 -1.295553605858698 0.000000000000000 0 1 0 0 - 1 4.626504860806483 0.000000000000000 0 0 1 0 - 1 8.719670307722497 0.000000000000000 0 0 0 1 - 2 39.35264616845731 0.000000000000000 2 0 0 0 - 2 -80.86751438913197 0.000000000000000 1 1 0 0 - 2 5547.935754774854 0.000000000000000 0 2 0 0 - 2 -49.34011516073829 0.000000000000000 1 0 1 0 - 2 452.1427490796136 0.000000000000000 0 1 1 0 - 2 39.02815553835716 0.000000000000000 1 0 0 1 - 2 298.9463122929935 0.000000000000000 0 1 0 1 - 2 25.20180266458752 0.000000000000000 0 0 2 0 - 2 21.71139328065176 0.000000000000000 0 0 1 1 - 2 2236.103272947546 0.000000000000000 0 0 0 2 - 3 640.9862784863635 0.000000000000000 3 0 0 0 - 3 -16584.61291160005 0.000000000000000 2 1 0 0 - 3 -355487.2726889580 0.000000000000000 1 2 0 0 - 3 1705894.705090054 0.000000000000000 0 3 0 0 - 3 3246.128535696794 0.000000000000000 2 0 1 0 - 3 -30561.22635820419 0.000000000000000 1 1 1 0 - 3 -215811.2089756964 0.000000000000000 0 2 1 0 - 3 -35779.30432043518 0.000000000000000 2 0 0 1 - 3 177666.8433217286 0.000000000000000 1 1 0 1 - 3 -5951521.149866415 0.000000000000000 0 2 0 1 - 3 -2537.074376518879 0.000000000000000 1 0 2 0 - 3 49700.66719624455 0.000000000000000 0 1 2 0 - 3 15548.56420941157 0.000000000000000 1 0 1 1 - 3 324164.4846352939 0.000000000000000 0 1 1 1 - 3 25342.22686092175 0.000000000000000 1 0 0 2 - 3 590190.0848609661 0.000000000000000 0 1 0 2 - 3 37.27461925625222 0.000000000000000 0 0 3 0 - 3 -41747.67462195644 0.000000000000000 0 0 2 1 - 3 -47485.77005419776 0.000000000000000 0 0 1 2 - 3 -951850.9994747553 0.000000000000000 0 0 0 3 + 0 -0.9378109818871268E-02 0.000000000000000 0 0 0 0 + 1 -0.1258080225944594 0.000000000000000 1 0 0 0 + 1 1.454695886838768 0.000000000000000 0 1 0 0 + 1 2.165374026944721 0.000000000000000 0 0 1 0 + 1 23.21846595999425 0.000000000000000 0 0 0 1 + 2 22.03083501149388 0.000000000000000 2 0 0 0 + 2 -138.2128837365694 0.000000000000000 1 1 0 0 + 2 2008.454246825898 0.000000000000000 0 2 0 0 + 2 -45.92789927553395 0.000000000000000 1 0 1 0 + 2 1368.199884666971 0.000000000000000 0 1 1 0 + 2 200.6658755188116 0.000000000000000 1 0 0 1 + 2 1059.706922082436 0.000000000000000 0 1 0 1 + 2 24.11612789441240 0.000000000000000 0 0 2 0 + 2 -81.70273938238094 0.000000000000000 0 0 1 1 + 2 489.6805613774040 0.000000000000000 0 0 0 2 + 3 -259.5390573040842 0.000000000000000 3 0 0 0 + 3 -23587.01078240259 0.000000000000000 2 1 0 0 + 3 -58189.18272880526 0.000000000000000 1 2 0 0 + 3 748180.4065702861 0.000000000000000 0 3 0 0 + 3 14522.12177141163 0.000000000000000 2 0 1 0 + 3 -90083.95892733682 0.000000000000000 1 1 1 0 + 3 1165562.509352523 0.000000000000000 0 2 1 0 + 3 -10110.14911353981 0.000000000000000 2 0 0 1 + 3 51786.86696377716 0.000000000000000 1 1 0 1 + 3 -1710650.753515610 0.000000000000000 0 2 0 1 + 3 -4832.914175220283 0.000000000000000 1 0 2 0 + 3 41999.25789822237 0.000000000000000 0 1 2 0 + 3 34757.31247922957 0.000000000000000 1 0 1 1 + 3 141403.7331520373 0.000000000000000 0 1 1 1 + 3 17790.35390684681 0.000000000000000 1 0 0 2 + 3 -825601.0919344509 0.000000000000000 0 1 0 2 + 3 11156.88470526157 0.000000000000000 0 0 3 0 + 3 24760.27019748771 0.000000000000000 0 0 2 1 + 3 186100.3386721462 0.000000000000000 0 0 1 2 + 3 -240750.2675011277 0.000000000000000 0 0 0 3 -35 0.000000000000000 0.000000000000000 0 0 0 0 1, NO = 3, NV = 4, INA = 184 @@ -260,41 +260,41 @@ I COEFFICIENT ORDER EXPONENTS NO = 3 NV = 4 - 0 0.9546649636982508 0.000000000000000 0 0 0 0 - 1 0.5001225667162666E-01 0.000000000000000 1 0 0 0 - 1 1.118510277439809 0.000000000000000 0 1 0 0 - 1 0.6051051633064486E-01 0.000000000000000 0 0 1 0 - 1 0.3674925494957817 0.000000000000000 0 0 0 1 - 2 5.288917504058189 0.000000000000000 2 0 0 0 - 2 -13.70211930915342 0.000000000000000 1 1 0 0 - 2 357.0365375968584 0.000000000000000 0 2 0 0 - 2 -1.318252230439930 0.000000000000000 1 0 1 0 - 2 40.22979595461808 0.000000000000000 0 1 1 0 - 2 -15.90379471531404 0.000000000000000 1 0 0 1 - 2 230.8850476490815 0.000000000000000 0 1 0 1 - 2 -17.01256218881513 0.000000000000000 0 0 2 0 - 2 -23.11879927918381 0.000000000000000 0 0 1 1 - 2 -100.0193522176275 0.000000000000000 0 0 0 2 - 3 75.61892731202389 0.000000000000000 3 0 0 0 - 3 -499.3947503263357 0.000000000000000 2 1 0 0 - 3 -10579.50287513764 0.000000000000000 1 2 0 0 - 3 -89105.88873359577 0.000000000000000 0 3 0 0 - 3 2.295342477294064 0.000000000000000 2 0 1 0 - 3 -2665.178232846630 0.000000000000000 1 1 1 0 - 3 -5099.048978746289 0.000000000000000 0 2 1 0 - 3 -1843.616358771547 0.000000000000000 2 0 0 1 - 3 27022.51012616824 0.000000000000000 1 1 0 1 - 3 -75442.34171377755 0.000000000000000 0 2 0 1 - 3 609.0177808834723 0.000000000000000 1 0 2 0 - 3 -5772.902989115783 0.000000000000000 0 1 2 0 - 3 -271.8340142302644 0.000000000000000 1 0 1 1 - 3 -27092.00057815559 0.000000000000000 0 1 1 1 - 3 4997.673230002756 0.000000000000000 1 0 0 2 - 3 -152936.6630823335 0.000000000000000 0 1 0 2 - 3 -60.38827595500342 0.000000000000000 0 0 3 0 - 3 -442.8709133961884 0.000000000000000 0 0 2 1 - 3 -2948.706628588041 0.000000000000000 0 0 1 2 - 3 -48222.96219295739 0.000000000000000 0 0 0 3 + 0 0.9548199704114351 0.000000000000000 0 0 0 0 + 1 0.4791277590665655E-01 0.000000000000000 1 0 0 0 + 1 1.192876545001144 0.000000000000000 0 1 0 0 + 1 -0.5979764743139625E-01 0.000000000000000 0 0 1 0 + 1 0.3222169398816387 0.000000000000000 0 0 0 1 + 2 4.271212982714482 0.000000000000000 2 0 0 0 + 2 -8.171019169182360 0.000000000000000 1 1 0 0 + 2 305.3082523842008 0.000000000000000 0 2 0 0 + 2 2.219590679479423 0.000000000000000 1 0 1 0 + 2 16.78458143258052 0.000000000000000 0 1 1 0 + 2 1.114140659593143 0.000000000000000 1 0 0 1 + 2 19.36559943479044 0.000000000000000 0 1 0 1 + 2 -24.15383806588803 0.000000000000000 0 0 2 0 + 2 -12.50225827341920 0.000000000000000 0 0 1 1 + 2 -331.5542469702108 0.000000000000000 0 0 0 2 + 3 47.99157589903309 0.000000000000000 3 0 0 0 + 3 -568.1116104102456 0.000000000000000 2 1 0 0 + 3 -1059.339287743586 0.000000000000000 1 2 0 0 + 3 -96031.59597441691 0.000000000000000 0 3 0 0 + 3 249.8709511362710 0.000000000000000 2 0 1 0 + 3 -4344.988696793234 0.000000000000000 1 1 1 0 + 3 4918.572818855383 0.000000000000000 0 2 1 0 + 3 -668.8585155362409 0.000000000000000 2 0 0 1 + 3 9356.922274684650 0.000000000000000 1 1 0 1 + 3 44950.09109969577 0.000000000000000 0 2 0 1 + 3 1075.572438533631 0.000000000000000 1 0 2 0 + 3 -3318.733755449175 0.000000000000000 0 1 2 0 + 3 1187.825534275374 0.000000000000000 1 0 1 1 + 3 -48722.62663894740 0.000000000000000 0 1 1 1 + 3 -1902.800290632470 0.000000000000000 1 0 0 2 + 3 -103539.1437517130 0.000000000000000 0 1 0 2 + 3 112.5499786383581 0.000000000000000 0 0 3 0 + 3 1716.372965986014 0.000000000000000 0 0 2 1 + 3 4001.423802867926 0.000000000000000 0 0 1 2 + 3 -12469.50320389639 0.000000000000000 0 0 0 3 -35 0.000000000000000 0.000000000000000 0 0 0 0 1, NO = 3, NV = 4, INA = 183 @@ -302,40 +302,40 @@ I COEFFICIENT ORDER EXPONENTS NO = 3 NV = 4 - 0 -0.1373684513865455E-01 0.000000000000000 0 0 0 0 - 1 -0.1272029950140894 0.000000000000000 1 0 0 0 - 1 -0.4083468942955703 0.000000000000000 0 1 0 0 - 1 -1.784844990389734 0.000000000000000 0 0 1 0 - 1 13.66982409386089 0.000000000000000 0 0 0 1 - 2 4.109768766835080 0.000000000000000 2 0 0 0 - 2 -179.7039967947391 0.000000000000000 1 1 0 0 - 2 -1469.852627613200 0.000000000000000 0 2 0 0 - 2 20.54063730362542 0.000000000000000 1 0 1 0 - 2 971.6338173048675 0.000000000000000 0 1 1 0 - 2 136.8085503868049 0.000000000000000 1 0 0 1 - 2 515.6243942970827 0.000000000000000 0 1 0 1 - 2 -42.74911083982582 0.000000000000000 0 0 2 0 - 2 -338.5584312754533 0.000000000000000 0 0 1 1 - 2 226.4183888686971 0.000000000000000 0 0 0 2 - 3 -168.2558299474317 0.000000000000000 3 0 0 0 - 3 -32131.88773240398 0.000000000000000 2 1 0 0 - 3 314925.9090240216 0.000000000000000 1 2 0 0 - 3 930473.6171353194 0.000000000000000 0 3 0 0 - 3 11711.36926503535 0.000000000000000 2 0 1 0 - 3 -19168.87380198314 0.000000000000000 1 1 1 0 - 3 1086488.893119412 0.000000000000000 0 2 1 0 - 3 -15765.64103980689 0.000000000000000 2 0 0 1 - 3 276744.2374849044 0.000000000000000 1 1 0 1 - 3 -691992.4089474851 0.000000000000000 0 2 0 1 - 3 5187.422356754444 0.000000000000000 1 0 2 0 - 3 -45648.16087356248 0.000000000000000 0 1 2 0 - 3 -18837.41570575458 0.000000000000000 1 0 1 1 - 3 -392860.4453980401 0.000000000000000 0 1 1 1 - 3 -5594.999229623592 0.000000000000000 1 0 0 2 - 3 -556435.7079250492 0.000000000000000 0 1 0 2 - 3 10249.14386602788 0.000000000000000 0 0 3 0 - 3 13165.56139303168 0.000000000000000 0 0 2 1 - 3 188126.0311651550 0.000000000000000 0 0 1 2 - 3 -235706.1670886399 0.000000000000000 0 0 0 3 + 0 -0.1186710894051485E-01 0.000000000000000 0 0 0 0 + 1 0.5031794767687248E-01 0.000000000000000 1 0 0 0 + 1 1.718152610595122 0.000000000000000 0 1 0 0 + 1 -6.541601200960002 0.000000000000000 0 0 1 0 + 1 6.360114683690653 0.000000000000000 0 0 0 1 + 2 -2.459940195208099 0.000000000000000 2 0 0 0 + 2 -128.7511789586227 0.000000000000000 1 1 0 0 + 2 -1359.255258515615 0.000000000000000 0 2 0 0 + 2 68.20995150526528 0.000000000000000 1 0 1 0 + 2 875.1471484527590 0.000000000000000 0 1 1 0 + 2 130.4036202523218 0.000000000000000 1 0 0 1 + 2 -782.7741343302763 0.000000000000000 0 1 0 1 + 2 4.760601453542805 0.000000000000000 0 0 2 0 + 2 46.00169979835894 0.000000000000000 0 0 1 1 + 2 -491.3467018616159 0.000000000000000 0 0 0 2 + 3 -979.7420749397230 0.000000000000000 3 0 0 0 + 3 -6559.359203140742 0.000000000000000 2 1 0 0 + 3 311326.1226868047 0.000000000000000 1 2 0 0 + 3 458644.9266322721 0.000000000000000 0 3 0 0 + 3 12071.04317801135 0.000000000000000 2 0 1 0 + 3 -35141.43208086607 0.000000000000000 1 1 1 0 + 3 1000912.122522561 0.000000000000000 0 2 1 0 + 3 34078.90718138096 0.000000000000000 2 0 0 1 + 3 202936.6062126463 0.000000000000000 1 1 0 1 + 3 5319594.090665693 0.000000000000000 0 2 0 1 + 3 4815.941688151748 0.000000000000000 1 0 2 0 + 3 -231.9273177052596 0.000000000000000 0 1 2 0 + 3 4192.027018295088 0.000000000000000 1 0 1 1 + 3 -508942.7259145931 0.000000000000000 0 1 1 1 + 3 -74754.69691117053 0.000000000000000 1 0 0 2 + 3 -315730.3314054444 0.000000000000000 0 1 0 2 + 3 9512.649240258916 0.000000000000000 0 0 3 0 + 3 64952.13203115657 0.000000000000000 0 0 2 1 + 3 110087.7128129228 0.000000000000000 0 0 1 2 + 3 434530.8143226430 0.000000000000000 0 0 0 3 -35 0.000000000000000 0.000000000000000 0 0 0 0 No Stochastic Radiation diff --git a/fpp-ptc-sandbox/code/z_resonance.f90 b/fpp-ptc-sandbox/code/z_resonance.f90 index e198fc3..cef3406 100644 --- a/fpp-ptc-sandbox/code/z_resonance.f90 +++ b/fpp-ptc-sandbox/code/z_resonance.f90 @@ -178,9 +178,37 @@ program Resonance !kall c_linear_a(T,A) !call print(A) ! stop - + !call print(T) + !stopS + res=.true. + + normal_form%positive=.false. +if(res) normal_form%nres=(c_%no+1)/3 + write(6,*) " # of resonance terms ",normal_form%nres + do i=1, (c_%no+1)/3 + if(res) normal_form%m(1,i)=3*i ! 1/3 order + enddo + +mr2=0 +if(res) then + do i=1,c_%nd +mr2=mr2+normal_form%m(i,1)**2 +enddo + +s=0 +do i=1,c_%nd +mu(i)=phat(i) +s=s+normal_form%m(i,1)*mu(i) +mu(i)=mu(i)*twopi +enddo +p_res=nint(s) + endif + call c_normal(T,normal_form,phase=phat,canonize=.true.) - write(*,*) "this should not have been reached" + call print(normal_form%atot) + write(*,*) "+=========================================" + stop + !write(*,*) "this should not have been reached" T = ci_phasor()*normal_form%atot**(-1)*T*normal_form%atot*c_phasor() call print(T) stop @@ -192,29 +220,7 @@ program Resonance - res=.true. - - normal_form%positive=.false. - if(res) normal_form%nres=(c_%no+1)/3 - write(6,*) " # of resonance terms ",normal_form%nres - do i=1, (c_%no+1)/3 - if(res) normal_form%m(1,i)=3*i ! 1/3 order - enddo - - mr2=0 - if(res) then - do i=1,c_%nd - mr2=mr2+normal_form%m(i,1)**2 - enddo - - s=0 - do i=1,c_%nd - mu(i)=phat(i) - s=s+normal_form%m(i,1)*mu(i) - mu(i)=mu(i)*twopi - enddo - p_res=nint(s) - endif + call c_normal(T,normal_form,canonize=.true.) diff --git a/fpp-ptc-sandbox/code/z_spin1.f90 b/fpp-ptc-sandbox/code/z_spin1.f90 index 1d69793..79f6232 100644 --- a/fpp-ptc-sandbox/code/z_spin1.f90 +++ b/fpp-ptc-sandbox/code/z_spin1.f90 @@ -52,11 +52,6 @@ program Guignard_normal_form_average_x ALS=>m_u%start call build_lattice(ALS,mis,exact=.false.,thin=.false.,onecell=.false.) - - - - !call read_ptc_command77("fit_tune.txt") - p=>als%start misa=0 call move_to(als,p,"BEND") if(state%nocavity) then @@ -69,6 +64,41 @@ program Guignard_normal_form_average_x misa=0.1d0*misa call MISALIGN_FIBRE(p,misa) + call kanalnummer(mf,file="fit_tune.txt") + write(mf,*) "select layout" + write(mf,*) " 1 " + write(mf,*) "set families " + write(mf,*) " 2 " + write(mf,*) " 1 qf1 " + write(mf,*) "2 1 " + write(mf,*) " 1 qd1 " + write(mf,*) "2 2 " + write(mf,*) "FITTUNE " + write(mf,*) "0.0000000001 " + write(mf,*) "0.367269834767751 0.404042919 !2293877 !6057 0.336 0.218992743820551 " + write(mf,*) "deallocate families " + write(mf,*) "return" + close(mf) + + + + call read_ptc_command77("fit_tune.txt") + + + !call read_ptc_command77("fit_tune.txt") + ! p=>als%start + ! misa=0 + ! call move_to(als,p,"BEND") + ! if(state%nocavity) then + ! misa(4)=0.03d0 + ! misa(6)=0.03d0 + ! else + ! misa(4)=0.003d0 + ! misa(6)=0.003d0 + ! endif + ! misa=0.1d0*misa + ! call MISALIGN_FIBRE(p,misa) +! courant_snyder_teng_edwards=.true. time_lie_choice=.true. @@ -85,13 +115,14 @@ program Guignard_normal_form_average_x MY_DELTA=>mondelta my_fix=>closed_orbit - !call phase_advance_n(6) + !call phase_advance_n() + !write(*,*) "hi" !write(6,format3) als%end%tm%lf%k(1,2,2) write(6,format3) als%end%t2%lf%phase write(6,format3) als%end%t2%lf%damping write(6,format3) als%end%t2%lf%spin write(6,format1) als%end%t2%lf%phase(1)*5 - + !stop !goto 1001 call kanalnummer(mf,"C:\document\my_tex_papers\fpp_handbook\julia\lat.txt") @@ -101,7 +132,6 @@ program Guignard_normal_form_average_x map_order=3 skew=.true. - call init_all(state,map_order,0) call alloc(one_turn_map, id_s,U_c,A,U,a_cs,N_c ) @@ -131,16 +161,24 @@ program Guignard_normal_form_average_x call print(one_turn_map,i) close(I) !stop + normal_form%nres=1 + normal_form%m(2,1)=1 + normal_form%ms(1)=-1 + write(*,*) "entering NF" call c_normal(one_turn_map,normal_form,dospin=state%spin,phase=phase) ! (6) - +write(*,*) "exiting NF" id_s=one_turn_map one_turn_map=ci_phasor()*normal_form%atot**(-1)*id_s*normal_form%atot*c_phasor() + one_turn_map%e_ij = 0 call clean(one_turn_map, one_turn_map, prec=1.e-7_dp) - one_turn_map%e_ij = 0 - !call c_q0_to_qr(one_turn_map%q, one_turn_map%q) ! only necessary when leaving spin resonance in so we make i, k parts 0 call print(one_turn_map) + stop + + + !call c_q0_to_qr(one_turn_map%q, one_turn_map%q) ! only necessary when leaving spin resonance in so we make i, k parts 0 + call print(one_turn_map%q) write(*,*) "hi" stop diff --git a/fpp-ptc-sandbox/code/z_spin_res.f90 b/fpp-ptc-sandbox/code/z_spin_res.f90 new file mode 100644 index 0000000..7a5aa9f --- /dev/null +++ b/fpp-ptc-sandbox/code/z_spin_res.f90 @@ -0,0 +1,358 @@ +program spin_res + use madx_ptc_module +use pointer_lattice +implicit none +! lesson 1 +type(probe) xs0,xs1,xst +type(probe_8) xs +type(layout), pointer :: als +integer mf,mf1,i,k,pos,no,kp,nturn,mfa,ks,kn,j +type(fibre),pointer:: p , fib +type(integration_node),pointer:: tt +type(internal_state),target :: state +real(dp) prec,cut, closed(6),x(6), ph(4),stu(2),damp(3),thi,nu_mod,circ,epsr,f3r,f3i,emi,xturn +real(dp) spin1,spin2,periodicity +complex(dp) epsb +logical first,breaksym +logical :: mis=.false.,thin=.false. +type(c_damap) id,Q,U_0,U_1,U_2,fp,N,Nc +type(c_taylor) phase(4),nu_spin +type(c_normal_form) c_n +type(c_spinor) isf +type(c_vector_field) fq +type(c_quaternion) qf +type(quaternion) q0,qr +TYPE(c_spinor) N_spin +TYPE(work) werk +real(dp) nus_in,nus_out,dnudn,dnus,pol_rat,dnuda,a_in,a_out,dnda,del_a,nus_shift +character(255) filename + + + integer jj(lnv) +first=.true.;lmax = 10.d0;use_info = .true.;prec=1.d-7;thin=.false. + +!ALWAYS_EXACT_PATCHING=.false. +use_quaternion=.true. +! check_excessive_cutting=.false. +call ptc_ini_no_append + call append_empty_layout(m_u) + als=>m_u%start + n_cai=-i_ +call build_lu(als) + + +breaksym=.true. + + + + +! call read_lattice_append(M_U,'C:\document\etienne_programs_new\programs_for_learning\fu\Fu_flat_new.txt') +periodicity=1.d0 + +als=>m_u%start +call MAKE_node_LAYOUT(als) +call survey(als) + +!goto 1001 + +state=only_4d0+spin0 !+modulation0 + +werk=als%start + + +call get_length(als,circ) +write(6,*) "circumference ",circ + +p=>als%start + +if(breaksym) then +call move_to(als,p,"QF2SPIN") +call add(p,2,0,0.01d0) +call add(p,2,0,40.0d0) +!call add(p,2,0,-20.d0) +!call add(p,3,0,-100.d0) +endif + +write(6,*) p%ag/werk%gamma0I/periodicity + + p=>als%start + +do i=1,als%n + + + p%ag=periodicity*(2.17557d0)*werk%gamma0I + if(breaksym) p%ag=periodicity*(2.d0+ 0.139870246560962d0)*werk%gamma0I + + + + +p=>p%next + +enddo + + +closed=0 +xs0=closed + + +p=>als%start !%next%next + +werk=p + closed=0 +call FIND_ORBIT_x(closed,state, 1.0e-7_dp, fibre1=p) + +write(6,*) check_stable +write(6,format6) closed + + +!stop +!goto 1000 + +write(6,format6) closed +call init(state,3,0) + + +call alloc(id,Q,U_0,U_1,U_2,fp,N,Nc) +call alloc(c_n) +call alloc(xs) +call alloc(isf) +call alloc(phase) +call alloc(nu_spin) +call alloc(fq) +call alloc(qf) +call alloc(N_spin) + + +xs0=closed +id=1 +xs=xs0+id +! call print(xs%ac(1),6 ) ! (F2) ! differs from the first edition + +call propagate(xs,state,fibre1=p) +id=xs + + kp=-1 +c_n%M=0 +c_n%NRES=1 +c_n%M(2,1)=1 +c_n%ms(1)=kp +c_n%positive=.false. + + +call c_normal(id,c_n, phase=phase, nu_spin=nu_spin,dospin=state%spin) + +write(6,*) c_n%tune(1:2) +write(6,*) c_n%spin_tune,c_n%quaternion_angle/pi + !call print(phase(2)) +! call print(nu_spin) + +!call print(c_n%h_l) + +!call kanalnummer(i,"spinmap.txt") +!call print(id,i) +!close(i) + +!stop + +N=c_n%atot**(-1)*id*c_n%atot +call clean(N,N,prec=1.d-7) + +U_0=N +U_0%q=1.d0 + fq%q%x(2)=-c_n%tune(2)*pi + + U_0=u_0*exp(fq) + +Nc=Ci_phasor()*N*U_0**(-1)*c_phasor() +fq=c_logf_spin(Nc) + +call c_q0_to_qr(fq%q,qf) + +call clean(qf,qf,prec=1.d-7) + call print(qf) + + + emi=1.d-5 + u_2=c_n%atot*c_phasor() + +u_1=0 +u_1%v(3)=emi +u_1%v(4)=emi + +u_2=u_2.o.u_1 + +x=0 +do i=1,4 +x(i)=u_2%v(i) +enddo + + +do i=0,3 +qf%x(i)=qf%x(i).o.u_1 +enddo +nus_shift=2*qf%x(2) +write(6,*) "nus_shift ", nus_shift + + +epsb=qf%x(1)*2 + +epsr=abs(epsb) + + +epsr=epsr**2 +write(6,*) "square ", epsr +!call print(qf) +!STOP +pause 888 + +use_quaternion=.true. +!use_quaternion=.false. + + +write(6,*) " pol_rat " +read(5,*) pol_rat +dnus=.02d0 + + +nus_in= (A_ELECTRON/werk%gamma0I)/periodicity- dnus/2.d0 +nus_out= (A_ELECTRON/werk%gamma0I)/periodicity +dnus/2.d0 + +a_in= nus_in*werk%gamma0I *periodicity +a_out= nus_out*werk%gamma0I*periodicity +write(6,*) a_in,a_out +write(6,*) nus_in,nus_out + + dnuda=(nus_out-nus_in)/(a_out-a_in) + dnudn= -epsr/log((pol_rat+1.d0)/2.d0)/4.d0 +write(6,*) dnuda,dnudn + +dnda=dnuda/dnudn +xturn=dnda*(a_out-a_in) +cut=dnda*(a_out-a_in) + + +del_a=1.d0/dnda +nus_shift=nus_shift/del_a +nturn=xturn +write(6,*) xturn,nturn + +pause 777 +write(6,format6) x + + +xs0=closed+x +xs0%q=1.d0 + +xs=xs0 + +fib=>als%start + + do j=1,als%n + fib%ag=fib%ag-(a_out-a_in)/2.d0 + fib=>fib%next + enddo + +call kanalnummer(mf,"plot.dat") + +p=>als%start +do i=1,nturn +fib=>als%start + + do j=1,als%n + fib%ag=fib%ag+del_a + fib=>fib%next + enddo + + + +call propagate(xs0,state,fibre1=p) +if(.not.check_stable) stop 777 +if(mod(i,nturn/100)==0) then + + + if(use_quaternion) then + q0=2 +qr=xs0%q + q0=qr*q0*qr**(-1) + + write(mf,format6)float(i), float(i)/nturn ,fib%ag/werk%gamma0I/periodicity, q0%x(2) + write(6,format2)float(i)/nturn, q0%x(2) + + + endif +endif +enddo + +my_estate=>state +close(mf) + + + + +1001 call ptc_end(graphics_maybe=1,flat_file=.false.) + + +contains + +subroutine build_lu(ALS) +use madx_ptc_module +use pointer_lattice +implicit none + + +type(layout), target :: ALS + +real(dp) :: alpha,lbend, cut, ksd, ksf,sig(6) + +type(fibre) L27h +type(fibre) QF2SPIN,QZ1,b2,beg,ende + + + +call make_states(.true.) + + exact_model=.false. + +call update_states +madlength = .false. + +!old_integrator=-100 +call set_mad(p0c = 0.94999986256844970d0, method = 2, step = 1) + +madkind2 = matrix_kick_matrix +!madkind2 = drift_kick_drift + + +ksd=0.d0 +QF2SPIN=QUADRUPOLE("QF2SPIN", 1.d-002, K1= ksd) +QZ1=QUADRUPOLE("QZ1", 1.d-002, K1= -125.66366121076938d0) +Lbend=1.d0 +alpha=pi/2 +B2=sbend("B2", LBEND, ANGLE=ALPHA) + + + + + + beg=marker("start"); + ende=marker("end"); + + + + + + ALS =beg+ QF2SPIN+4*(QZ1 + B2 +QZ1) +ende; + + +ALS = .ring.ALS + +call survey(ALS) + + +! sig=1.d-5; cut=4.d0; +! call MESS_UP_ALIGNMENT(ALS,SIG,cut); + +end subroutine build_lu +end program spin_res + + \ No newline at end of file diff --git a/src/normal.jl b/src/normal.jl index b91ddfe..5ed61e4 100644 --- a/src/normal.jl +++ b/src/normal.jl @@ -100,7 +100,7 @@ function normal(m::DAMap; res=nothing, spin_res=nothing) while idx > 0 # Tune shifts should be left in the map (kernel) because we cannot remove them # if there is damping, we technically could remove them - if !is_tune_shift(j, ords, nhv) && !is_resonance(j, ords, nhv, res) # then remove it + if !is_tune_shift(j, ords, nhv) && !is_orbital_resonance(j, ords, nhv, res, spin_res) # then remove it je = convert(Vector{Int}, ords) je[j] -= 1 lam = 1 @@ -109,11 +109,13 @@ function normal(m::DAMap; res=nothing, spin_res=nothing) end F.x[j] += mono(ords,use=getdesc(m1))*v[]/(1-lam) else # cannot remove it - add to kernel (if close to res and res is specified it will nearly blow up) + println("keeping monomial ", Vector{Int}(ords), ", v = ", v[]) Fker.x[j] += mono(ords,use=getdesc(m1))*v[] end idx = GTPSA.cycle!(nonl.x[j], idx, nn, ords, v) end + println("================================") end kert = exp(Fker,one(m)) #I + Fker @@ -155,7 +157,7 @@ function normal(m::DAMap; res=nothing, spin_res=nothing) Qr_inv = DAMap(Q=inv(Quaternion(scalar.(m1.Q)))) # Now store analogous to eg -> egspin egspin = SVector(cos(nu0)+im*sin(nu0), 1, cos(nu0)-im*sin(nu0)) - + for i =1:mo # get rid of ℛ: linandnonl = m1 ∘ Qr_inv @@ -183,33 +185,40 @@ function normal(m::DAMap; res=nothing, spin_res=nothing) for j=1:3 v = Ref{ComplexF64}() ords = Vector{UInt8}(undef, nn) - idx = GTPSA.cycle!(nr_s[j], 0, nn, ords, v) + idx = GTPSA.cycle!(nr_s[j], 0, nn, ords, v) while idx > 0 # We remove every term in x and z, tune shifts will only be left in y component # because of how we defined everything - if j != 2 || (!is_tune_shift(j, ords, nhv, true) && !is_spin_resonance(j, ords, nhv, res, spin_res)) # then remove it, note spin components are like hamiltonian + # NOTE SPIN RESONANCES WILL HAPPEN WHEN j != 2 SO WE CHECK THAT FIRST!! + if !is_spin_resonance(j, ords, nhv, res, spin_res) && (j != 2 || !is_tune_shift(j, ords, nhv, true)) # then remove it, note spin components are like hamiltonian lam = egspin[j] for k = 1:nhv # ignore coasting plane lam *= eg[k]^ords[k] end na[j] += mono(ords,use=getdesc(m1))*v[]/(1-lam) + else + println("keeping monomial ", Vector{Int}(ords), ", v = ", v[]) end idx = GTPSA.cycle!(nr_s[j], idx, nn, ords, v) end + println("==================================") end - # Now exit the basis na = [(na[1]+na[3])/2, na[2], im*(na[1]-na[3])/2] + # Exponentiate this part now Qnr = DAMap(Q=exp(Quaternion(0,na...))) as = as*Qnr # put in normalizing map m1 = inv(Qnr)*m1*Qnr # kill the terms in m1 end + return a*c*as*c^-1 end - #return a - return a*c*as*c^-1 + + #return as + return a + @@ -252,7 +261,7 @@ function is_tune_shift(varidx, ords, nhv, hamiltonian=false) end """ - is_resonance(varidx, ords, nhv, res) + is_orbital_resonance(varidx, ords, nhv, res) Checks if the monomial corresponds to a particular resonance (and resonance in the same family) @@ -273,14 +282,18 @@ m = 3 6 9 ...]; """ -function is_resonance(varidx, ords, nhv, res) +function is_orbital_resonance(varidx, ords, nhv, res, spin_res) if isnothing(res) return false end + je = convert(Vector{Int}, ords) je[varidx] -= 1 for curresidx=1:size(res, 2) # for each res in the family + if !isnothing(spin_res) && spin_res[curresidx] != 0 + return false # spin res not orbital res + end t1 = 0 t2 = 0 @@ -297,12 +310,18 @@ end """ +nu_s + j dot Q = n + +0*Qx + 1*Qy + 2*Qs = n + +IMPORTANT:::: +For spin resonances, there is only one resonance +in each resonance family: -0*Qx + 1*Qy + 1*Qs = n -m = [0 0 0 ...; - 1 2 3 ...] +m = [0 ; + 1 ] -ms = [1, 2, 3, ...] +ms = [1]; In the code, check - of m + ms """ @@ -310,6 +329,7 @@ function is_spin_resonance(spinidx, ords, nhv, res, spin_res) if isnothing(res) && isnothing(spin_res) return false end + je = convert(Vector{Int}, ords) @assert size(res, 2) == length(spin_res) "Number of resonances in spin_res != number of resonances in res" @@ -318,8 +338,8 @@ function is_spin_resonance(spinidx, ords, nhv, res, spin_res) t2 = 0 for k = 1:2:nhv # ignore coasting plane - t1 += abs(ords[k]-ords[k+1]+res[Int((k+1)/2),curresidx]) - t2 += abs(ords[k]-ords[k+1]-res[Int((k+1)/2),curresidx]) + t1 += abs(je[k]-je[k+1]+res[Int((k+1)/2),curresidx]) + t2 += abs(je[k]-je[k+1]-res[Int((k+1)/2),curresidx]) end m = spin_res[curresidx] diff --git a/test/3nux_leavein/a.map b/test/3nux_leavein/a.map new file mode 100644 index 0000000..386d2e5 --- /dev/null +++ b/test/3nux_leavein/a.map @@ -0,0 +1,448 @@ + 4 Dimensional TPSA/DA map + + 1, NO = 6, NV = 4, INA = 364 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 6 NV = 4 + 1 3.089277205302563 0.000000000000000 1 0 0 0 + 2 -745.8872808832762 0.000000000000000 2 0 0 0 + 2 -2901.782463392580 0.000000000000000 1 1 0 0 + 2 -2237.661842649829 0.000000000000000 0 2 0 0 + 2 1566.142068618037 0.000000000000000 0 0 2 0 + 2 -33.19221198338952 0.000000000000000 0 0 1 1 + 2 -20.10814225498314 0.000000000000000 0 0 0 2 + 3 -152011341.3539437 0.1026769916193888E-06 3 0 0 0 + 3 -526176180.9536611 -0.4603381760869268E-07 2 1 0 0 + 3 -920461186.9808520 0.1438556800271646E-07 1 2 0 0 + 3 -531781604.9508072 -0.1150845440217317E-07 0 3 0 0 + 3 631860548.1592329 -0.7329446897384038E-06 1 0 2 0 + 3 479418861.4823626 0.000000000000000 0 1 2 0 + 3 -83865.52873414631 0.000000000000000 1 0 1 1 + 3 -126178049984.0884 -0.1529905157088896E-05 0 1 1 1 + 3 -47671.95697284423 0.7329446897384038E-06 1 0 0 2 + 3 496892772.9645286 0.000000000000000 0 1 0 2 + 4 832431097020.4404 -0.1885545169252052E-03 4 0 0 0 + 4 617200980073.0580 -0.4713862923130131E-03 3 1 0 0 + 4 -1183941266642.556 0.000000000000000 2 2 0 0 + 4 -1010205290178.882 -0.2828317753878079E-03 1 3 0 0 + 4 -470651479957.5957 -0.2474778034643319E-03 0 4 0 0 + 4 55136775357759.99 0.1960966976022134E-01 2 0 2 0 + 4 292835070489.3641 -0.1930798253314102 1 1 2 0 + 4 -182277500239554.9 -0.1960966976022134E-01 0 2 2 0 + 4 58949247172012.60 -0.9653991266570507 2 0 1 1 + 4 4633949939451547. 0.1327423799153445 1 1 1 1 + 4 -57997103659044.98 0.1930798253314102 0 2 1 1 + 4 -120655484728222.0 -0.2262654203102463E-01 2 0 0 2 + 4 -330489041533.4475 0.1930798253314102 1 1 0 2 + 4 118987352924502.6 0.2262654203102463E-01 0 2 0 2 + 4 11842671496076.10 -0.5656635507756158E-03 0 0 4 0 + 4 156968466496.4189 -0.2262654203102463E-02 0 0 3 1 + 4 22765080954745.46 -0.3393981304653694E-02 0 0 2 2 + 4 -276639234882.5569 -0.1885545169252052E-02 0 0 1 3 + 4 96605488379.43138 0.1696990652326847E-02 0 0 0 4 + 5 0.1643515266434393E+18 6.178554410605125 5 0 0 0 + 5 0.7054211454615286E+17 92.67831615907687 4 1 0 0 + 5 0.5604115985549324E+18 197.7137411393640 3 2 0 0 + 5 -0.1150557863980276E+18 135.9281970333128 2 3 0 0 + 5 0.1603137525544804E+18 284.2135028878358 1 4 0 0 + 5 0.3350740538133411E+18 142.1067514439179 0 5 0 0 + 5 -0.2731769074862176E+20 30384.58595275335 3 0 2 0 + 5 0.6123694404746239E+18 -7696600.515073162 2 1 2 0 + 5 -0.2399501832332957E+20 44874.84068422503 1 2 2 0 + 5 0.1410486954678125E+17 -7691855.385285817 0 3 2 0 + 5 -0.2360664492479624E+20 -502983.7574585420 3 0 1 1 + 5 0.1024006845830624E+23 -30114.27419728938 2 1 1 1 + 5 -0.2361072980378984E+20 -926882.0184613385 1 2 1 1 + 5 0.1034374256203457E+23 -25817.08960471351 0 3 1 1 + 5 0.1308384767192962E+18 10264.12351461776 3 0 0 2 + 5 -0.4040764593012182E+20 10565822.32648761 2 1 0 2 + 5 -0.2897858997247002E+18 846.4619542529022 1 2 0 2 + 5 -0.4078655927231005E+20 10515207.60875593 0 3 0 2 + 5 0.5587241146902801E+19 -484670.5221855085 1 0 4 0 + 5 0.2129438928635072E+17 316341.9858229824 0 1 4 0 + 5 0.1953437297809184E+20 -1062909.072365221 1 0 3 1 + 5 -0.2052239621374630E+22 739350.5349906516 0 1 3 1 + 5 -0.2567887017181306E+22 2923444.804921921 1 0 2 2 + 5 0.9459265118363787E+19 -2328277.015657151 0 1 2 2 + 5 0.2024713052625054E+20 860450.2014385122 1 0 1 3 + 5 -0.2412567882078905E+22 121890.5214124179 0 1 1 3 + 5 -0.3973204633598150E+17 -524608.6978956599 1 0 0 4 + 5 0.9543518921260462E+19 620030.2922130455 0 1 0 4 + 6 -0.8761364208999748E+21 -376446.9631293491 6 0 0 0 + 6 -0.2382129738218097E+22 404917.7418534175 5 1 0 0 + 6 -0.4774555982387246E+21 227766.2297925473 4 2 0 0 + 6 0.4850443874435822E+20 -1012294.354633544 3 3 0 0 + 6 0.4004649202527737E+22 1796822.479474540 2 4 0 0 + 6 0.1874745865923159E+22 1214753.225560253 1 5 0 0 + 6 -0.4374542559199139E+21 376446.9631293491 0 6 0 0 + 6 0.4683282090070992E+24 5908559689.125068 4 0 2 0 + 6 -0.7475525922010925E+24 255415632877.2661 3 1 2 0 + 6 0.1242167115748260E+26 31512318342.00036 2 2 2 0 + 6 -0.8709202891664976E+24 580490074721.0593 1 3 2 0 + 6 0.5857371080695238E+25 -16896407532.05940 0 4 2 0 + 6 -0.3939348340607505E+25 -41463576765.78995 4 0 1 1 + 6 -0.2265570668709163E+27 29024503736.05296 3 1 1 1 + 6 -0.2903069080487569E+26 -26536689130.10557 2 2 1 1 + 6 -0.6919245943906311E+27 102000398843.8433 1 3 1 1 + 6 0.1363139273887431E+26 24878146059.47397 0 4 1 1 + 6 0.3174337384299433E+25 -1762202012.546073 4 0 0 2 + 6 0.8121934565016717E+24 -182439737769.4758 3 1 0 2 + 6 0.1603450609231868E+26 -24048874524.15817 2 2 0 2 + 6 0.8819979374484302E+24 -792783587761.9038 1 3 0 2 + 6 -0.8534759566239648E+25 10262235249.53301 0 4 0 2 + 6 -0.2151948976673188E+26 2455421186.599123 2 0 4 0 + 6 0.2946324590357735E+25 -15341523403.34228 1 1 4 0 + 6 0.1753455310326708E+26 239711303.1772231 0 2 4 0 + 6 0.2382107315507839E+24 15756159171.00018 2 0 3 1 + 6 0.2392210128494570E+26 15341523403.34228 1 1 3 1 + 6 -0.2740219134932185E+24 -22390331453.52657 0 2 3 1 + 6 0.1055243789660128E+26 -24670828175.64502 2 0 2 2 + 6 0.7932206629743500E+25 139317617933.0542 1 1 2 2 + 6 -0.1795813120431623E+26 23841556640.32922 0 2 2 2 + 6 0.4159156073613521E+24 -19902516847.57918 2 0 1 3 + 6 0.7847274546204752E+26 -27106813310.63518 1 1 1 3 + 6 -0.3335773590385476E+24 13268344565.05278 0 2 1 3 + 6 -0.1843595602058684E+26 6912755688.921543 2 0 0 4 + 6 0.2456075995175294E+25 -41048940998.13205 1 1 0 4 + 6 0.1841990600168386E+26 -5461530502.118895 0 2 0 4 + 6 0.2382473590781546E+24 -61927107.14470704 0 0 6 0 + 6 0.4627510079014749E+22 -135242525.7790414 0 0 5 1 + 6 0.8042898795249592E+24 -853566599.8270041 0 0 4 2 + 6 0.2720819823745253E+21 116616309.6537842 0 0 3 3 + 6 0.6248406997172752E+24 -994882891.7338468 0 0 2 4 + 6 -0.7609681382109178E+22 46160622.57128959 0 0 1 5 + 6 0.8789182813793232E+21 16576320.05712428 0 0 0 6 + -108 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 6, NV = 4, INA = 365 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 6 NV = 4 + 1 0.3578838508394021E-01 0.000000000000000 1 0 0 0 + 1 0.3237003135502243 0.000000000000000 0 1 0 0 + 2 447.4405004096221 0.000000000000000 2 0 0 0 + 2 122.6946499037922 0.000000000000000 1 1 0 0 + 2 126.1044629618516 0.000000000000000 0 2 0 0 + 2 -263.7946810376746 0.000000000000000 0 0 2 0 + 2 -4.391605013839287 0.000000000000000 0 0 1 1 + 2 -2.429117597520883 0.000000000000000 0 0 0 2 + 3 -94042149.45299715 -0.1523741361935533E-08 3 0 0 0 + 3 40907188.04324088 -0.5332884950032861E-09 2 1 0 0 + 3 43883162.53099637 -0.1039224982989203E-08 1 2 0 0 + 3 26008404.92319223 -0.1460385377588358E-07 0 3 0 0 + 3 161723211.6391698 -0.8490952756380446E-08 1 0 2 0 + 3 -60643493.82233886 -0.9164670046350750E-07 0 1 2 0 + 3 13218924241.62147 0.1254112743184840E-06 1 0 1 1 + 3 -1461715721.679328 -0.1772350982612484E-07 0 1 1 1 + 3 -52074141.52112622 0.8490952756380446E-08 1 0 0 2 + 3 5761359.537310023 0.9164670046350750E-07 0 1 0 2 + 4 39120171532.80966 0.4226912355981839E-04 4 0 0 0 + 4 -191294730936.8736 0.5381042345830215E-04 3 1 0 0 + 4 -141960028494.6699 0.000000000000000 2 2 0 0 + 4 -79797920572.22917 -0.3276524513300190E-05 1 3 0 0 + 4 26386223940.17202 -0.4732043218448952E-04 0 4 0 0 + 4 12162562844090.88 0.2271723662554798E-03 2 0 2 0 + 4 -25023305257340.35 -0.3501228417551827E-02 1 1 2 0 + 4 9227893203707.139 -0.2271723662554798E-03 0 2 2 0 + 4 256149502653084.1 0.1157630795776883E-01 2 0 1 1 + 4 -8240151752085.648 0.8246286055913164E-01 1 1 1 1 + 4 -256305108275100.6 0.9823500166579645E-02 0 2 1 1 + 4 -1552895000007.370 -0.2621219610640152E-03 2 0 0 2 + 4 25303622473695.79 0.9723197179406992E-03 1 1 0 2 + 4 1355698194918.629 0.2621219610640152E-03 0 2 0 2 + 4 -1988135549998.032 -0.6553049026600380E-05 0 0 4 0 + 4 50032091416.30510 0.5269865819335326E-03 0 0 3 1 + 4 -3696556945066.317 -0.4344602784738409E-03 0 0 2 2 + 4 -1166015892.568950 0.3732984875589041E-03 0 0 1 3 + 4 11684597192.98300 -0.5936924978304660E-04 0 0 0 4 + 5 0.2061206177848168E+17 11.56293790120084 5 0 0 0 + 5 -0.3138953919378784E+17 -13.16916224369166 4 1 0 0 + 5 0.8174207959260483E+17 -5.478350879833211 3 2 0 0 + 5 -0.1246845495373780E+17 -3.604516073110220 2 3 0 0 + 5 3937338320690977. -19.36649052079321 1 4 0 0 + 5 683834867764928.0 -11.30174682814772 0 5 0 0 + 5 -0.8634636343308325E+19 808691.2404540965 3 0 2 0 + 5 0.2789097437043286E+19 -83859.39035384108 2 1 2 0 + 5 -0.9222846826262718E+19 715167.2871203558 1 2 2 0 + 5 0.2706552297800586E+19 -82395.67562978159 0 3 2 0 + 5 -0.1077624115400185E+22 -5329.718024213664 3 0 1 1 + 5 0.1360440296053725E+21 10589.61581769093 2 1 1 1 + 5 -0.1079581572890231E+22 -2782.401847694435 1 2 1 1 + 5 0.1173227341001585E+21 52653.10855765340 0 3 1 1 + 5 0.4202704537650100E+19 -1094547.507032132 3 0 0 2 + 5 -0.5236159527211512E+18 119750.2501600855 2 1 0 2 + 5 0.4304203866687380E+19 -1022406.697939650 1 2 0 2 + 5 -0.4575830449180864E+18 118086.6186153751 0 3 0 2 + 5 0.1742047876717777E+19 -28154.66839217843 1 0 4 0 + 5 -0.5476186342189297E+18 -47174.34581234855 0 1 4 0 + 5 0.2098560313347979E+21 -91203.14574147308 1 0 3 1 + 5 -0.3007657023855131E+20 -102808.4620559751 0 1 3 1 + 5 -0.2780536312442937E+20 336167.1029934262 1 0 2 2 + 5 -0.2707175791265065E+21 279967.9886599452 0 1 2 2 + 5 0.2530023386519223E+21 -2793.493840744346 1 0 1 3 + 5 -0.2581730374071935E+20 123392.7030776301 0 1 1 3 + 5 -0.1000326338561290E+19 -92259.41188102611 1 0 0 4 + 5 0.1063459027991018E+18 -47053.76289608005 0 1 0 4 + 6 -0.3943848165845980E+20 -35187.65771280359 6 0 0 0 + 6 0.2437095477435031E+21 100153.9620794460 5 1 0 0 + 6 0.7575663382356460E+21 -7968.405818945007 4 2 0 0 + 6 0.1271758512262509E+21 115557.0044686595 3 3 0 0 + 6 0.6884730313816796E+20 -85254.44875099520 2 4 0 0 + 6 -0.8347015584671007E+18 -166246.6362358671 1 5 0 0 + 6 -0.1740975226000987E+20 -33757.91947088580 0 6 0 0 + 6 -0.8679265583333774E+24 25180422287.75233 4 0 2 0 + 6 0.3412958415637280E+24 -4948313943.433553 3 1 2 0 + 6 -0.9341840907948567E+24 -24660019571.41544 2 2 2 0 + 6 0.6545795647127197E+24 3683567584.025851 1 3 2 0 + 6 -0.5828213076297702E+24 -18182516750.15716 0 4 2 0 + 6 -0.3120522080130922E+26 3234316841.039136 4 0 1 1 + 6 0.4553826775052128E+25 24666180058.48817 3 1 1 1 + 6 0.3667763321338986E+26 -5694763646.086596 2 2 1 1 + 6 -0.4141854362025551E+25 21340737966.59567 1 3 1 1 + 6 0.1900030582923017E+26 -2340296254.489390 0 4 1 1 + 6 0.1054014119201071E+24 -33821652057.92452 4 0 0 2 + 6 -0.2353862956652947E+25 5315809105.736942 3 1 0 2 + 6 0.1664760850650690E+24 47686138711.29317 2 2 0 2 + 6 -0.8944550450090421E+23 -5491231870.656709 1 3 0 2 + 6 -0.1445907748422901E+24 17758091213.30028 0 4 0 2 + 6 0.2721420942665621E+24 -1524760616.802399 2 0 4 0 + 6 0.4133936417383665E+25 -335220034.4972505 1 1 4 0 + 6 0.3977603947535611E+24 1208412383.977520 0 2 4 0 + 6 0.1358179948645335E+25 845086947.6440912 2 0 3 1 + 6 -0.2566893383885499E+24 -5470294560.701879 1 1 3 1 + 6 -0.1371585882917893E+25 -617817674.9369888 0 2 3 1 + 6 0.1713715690738765E+24 8490352342.576988 2 0 2 2 + 6 -0.3954187515379120E+25 5958586470.637735 1 1 2 2 + 6 0.1052843997132208E+25 -7457247518.744129 0 2 2 2 + 6 0.4366881060829424E+25 -480381258.9344587 2 0 1 3 + 6 -0.1112220913910485E+24 3943714811.325031 1 1 1 3 + 6 -0.4368927843861186E+25 2358610715.870313 0 2 1 3 + 6 -0.8677070091358546E+23 -2255157454.130554 2 0 0 4 + 6 0.3894014815820837E+25 -1567128943.760467 1 1 0 4 + 6 0.8481798143616642E+23 1924398914.101408 0 2 0 4 + 6 -0.4002314973548620E+23 14526857.20962212 0 0 6 0 + 6 0.1120862774233748E+22 29999721.69820810 0 0 5 1 + 6 -0.1336183006815960E+24 124725264.9160905 0 0 4 2 + 6 0.1404010142226543E+22 -8152916.339074723 0 0 3 3 + 6 -0.1022303912820946E+24 131510123.8761820 0 0 2 4 + 6 -0.2722785917262275E+20 -3198910.685885308 0 0 1 5 + 6 0.1057956523906774E+21 3276683.525875952 0 0 0 6 + -109 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 6, NV = 4, INA = 366 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 6 NV = 4 + 1 2.098259839140808 0.000000000000000 0 0 1 0 + 1 0.3639900601656879E-17 0.000000000000000 0 0 0 1 + 2 25.97433672618529 -0.2329536385060402E-15 1 0 1 0 + 2 -22.54439493401693 -0.2981806572877315E-13 0 1 1 0 + 2 28.47162215561574 0.2981806572877315E-13 1 0 0 1 + 2 -27.31519674630692 0.2329536385060402E-15 0 1 0 1 + 3 -42843301821.10007 -0.1656313877361721E-05 2 0 1 0 + 3 -121393.4625260416 -0.2442695944501096E-09 1 1 1 0 + 3 -42850546037.68111 -0.1657935980137366E-05 0 2 1 0 + 3 337501180.1590168 -0.5303007401153823E-04 2 0 0 1 + 3 -104724.3121258802 -0.2211403172256149E-07 1 1 0 1 + 3 337459780.8098767 -0.4302063883975608E-04 0 2 0 1 + 3 8292859686.724035 -0.1193676807254247E-07 0 0 3 0 + 3 -78845770.27078886 0.000000000000000 0 0 2 1 + 3 14925742639.91135 0.4362788541813423E-07 0 0 1 2 + 3 -78857795.97280347 -0.3751980970753684E-06 0 0 0 3 + 4 -538907080151990.6 0.9220868433724254E-01 3 0 1 0 + 4 120623530921980.7 0.2622824798926010 2 1 1 0 + 4 1614539193126740. -0.2786751348858886 1 2 1 0 + 4 -39621402023483.74 -0.6557061997315025E-01 0 3 1 0 + 4 763240910521.3949 -0.6229208897449274 3 0 0 1 + 4 -163981794295798.1 -0.8606143871475970E-01 2 1 0 1 + 4 22047019807.89220 1.639265499328756 1 2 0 1 + 4 53894661743773.37 0.3073622811241418E-01 0 3 0 1 + 4 127350778706.0363 -0.2562853135089483E-02 1 0 3 0 + 4 -72104110365.41389 0.5122704685402363E-03 0 1 3 0 + 4 -1249652493468.322 0.9733138902264491E-02 1 0 2 1 + 4 -716547668690.5276 -0.4442345469372362E-03 0 1 2 1 + 4 -926292495054.4225 -0.2621384038233241E-02 1 0 1 2 + 4 -1320114102881.019 -0.1536811405620709E-02 0 1 1 2 + 4 -64093422984.11141 -0.5122704685402363E-03 1 0 0 3 + 4 70194218958.59415 -0.5137712609285378E-03 0 1 0 3 + 5 0.2186737084700588E+22 -68383.07500503861 4 0 1 0 + 5 5082458796733196. 34377.36463952321 3 1 1 0 + 5 0.4373657192182458E+22 -138270.0777497814 2 2 1 0 + 5 -0.1706478999961719E+17 103134.1921784088 1 3 1 0 + 5 0.2186895309152161E+22 -71024.25957755711 0 4 1 0 + 5 -0.1375322868046026E+20 992747.1532238215 4 0 0 1 + 5 -0.1276876274714857E+17 -2185.337622465152 3 1 0 1 + 5 -0.2754167635907618E+20 2191580.732334035 2 2 0 1 + 5 0.2012574797577315E+17 -1614.610946218852 1 3 0 1 + 5 -0.1376365684755492E+20 1164468.738459105 0 4 0 1 + 5 -0.1196201828209930E+22 177110.4410671571 2 0 3 0 + 5 3437420190974634. 381.8832907236271 1 1 3 0 + 5 -0.1196241401683758E+22 186811.7454334246 0 2 3 0 + 5 0.7312530858580010E+19 -1179910.258824372 2 0 2 1 + 5 -0.5460456298997391E+18 -1898.925154422431 1 1 2 1 + 5 0.7352747719189513E+19 -1157786.207080472 0 2 2 1 + 5 -0.2153102096209258E+22 406750.8172072044 2 0 1 2 + 5 0.4327772185837292E+17 369.2937316887822 1 1 1 2 + 5 -0.2153104314558345E+22 409724.0513992669 0 2 1 2 + 5 0.1053212282767368E+20 -2700733.186753308 2 0 0 3 + 5 -1868767836260464. -1898.925154422431 1 1 0 3 + 5 0.1053466961382467E+20 -2662217.529146039 0 2 0 3 + 5 0.1487625961959540E+21 -56918.40556935658 0 0 5 0 + 5 -0.8717888471185571E+18 15040.32652696131 0 0 4 1 + 5 0.4862420430668817E+21 -530505.5350098569 0 0 3 2 + 5 -0.2491192868139096E+19 12891.70845168112 0 0 2 3 + 5 0.3584521755243155E+21 -566346.3580220697 0 0 1 4 + 5 -0.1619587976018748E+19 38272.25946592834 0 0 0 5 + 6 0.5076306684212753E+26 -11918410246.99610 5 0 1 0 + 6 -0.1145477166038418E+26 -11068580231.38898 4 1 1 0 + 6 -0.1014151566448009E+27 24558584277.84033 3 2 1 0 + 6 -0.7667710966178774E+25 17631181768.96797 2 3 1 0 + 6 -0.1521856073828644E+27 37708908995.53448 1 4 1 0 + 6 0.3787773362980227E+25 -12577460166.46295 0 5 1 0 + 6 -0.1410744100791403E+24 19452934873.69193 5 0 0 1 + 6 0.7651056933188323E+25 18364952030.62075 4 1 0 1 + 6 0.2096440371960648E+24 -36820094453.56947 3 2 0 1 + 6 0.5123971468734493E+25 13399108907.39469 2 3 0 1 + 6 0.3507364367559913E+24 -64950558608.70019 1 4 0 1 + 6 -0.2527935130056373E+25 -5001020298.354550 0 5 0 1 + 6 -0.1141989559030685E+26 -5984502952.716400 3 0 3 0 + 6 0.2626494139378692E+25 -34358087540.30122 2 1 3 0 + 6 0.3419618277956311E+26 13130703537.43069 1 2 3 0 + 6 -0.8597454859495736E+24 14644430754.88249 0 3 3 0 + 6 -0.1063373197778210E+25 -35484582213.75372 3 0 2 1 + 6 -0.4733021902846878E+26 -75615954955.49899 2 1 2 1 + 6 0.3704302707082552E+25 110396477998.3449 1 2 2 1 + 6 0.1588455359990528E+26 25662956779.58974 0 3 2 1 + 6 -0.2357839660271941E+26 24606868023.22802 3 0 1 2 + 6 0.5388882681466103E+25 -92372563223.10490 2 1 1 2 + 6 0.7089532014236735E+26 -73926212945.32024 1 2 1 2 + 6 -0.1670899476478152E+25 34921334877.02747 0 3 1 2 + 6 -0.1169175950930870E+25 -23093140805.77623 3 0 0 3 + 6 -0.5124571763983474E+26 27211886955.58693 2 1 0 3 + 6 0.3534979892407798E+25 74911895784.59117 1 2 0 3 + 6 0.1707248713920268E+26 -8149484903.257922 0 3 0 3 + 6 0.2092191509283146E+22 8813631.344799329 1 0 5 0 + 6 -0.1416231080606560E+22 38434480.13061199 0 1 5 0 + 6 -0.2390592141785819E+23 209842635.7041642 1 0 4 1 + 6 -0.1650196958715096E+23 -115578463.5054718 0 1 4 1 + 6 -0.1892935893268614E+23 194578852.8973738 1 0 3 2 + 6 0.8729331383448138E+22 -207917473.9087132 0 1 3 2 + 6 -0.5592612787525970E+23 36853097.22720578 1 0 2 3 + 6 0.2858892986514514E+22 -196779037.8064607 0 1 2 3 + 6 -0.1390075376869367E+23 44622500.18741893 1 0 1 4 + 6 -0.2338998111080779E+23 -54179553.38626520 0 1 1 4 + 6 -0.6312466404317366E+21 20833000.85791670 1 0 0 5 + 6 0.7515575120036807E+21 61867308.85961772 0 1 0 5 + -100 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 6, NV = 4, INA = 367 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 6 NV = 4 + 1 0.1376883188117553E-01 0.000000000000000 0 0 1 0 + 1 0.4765853977405764 0.000000000000000 0 0 0 1 + 2 -830.0266307116445 0.2708921473929470E-13 1 0 1 0 + 2 -483.3680502366790 0.2689507651725094E-13 0 1 1 0 + 2 -5.712814207872496 0.1956668694614269E-15 1 0 0 1 + 2 4.941347530577305 0.2709227203413004E-13 0 1 0 1 + 3 -54290533.42833087 -0.1159370205762819E-04 2 0 1 0 + 3 -194394743.8205418 0.1420174876573803E-07 1 1 1 0 + 3 -354772280.6976443 -0.1024439429337628E-04 0 2 1 0 + 3 9733358648.508766 0.6391256841998741E-07 2 0 0 1 + 3 31009.19487164752 0.1405823885700024E-07 1 1 0 1 + 3 9735021319.769005 0.7278131775592868E-07 0 2 0 1 + 3 24564775.75113460 0.1135484839886784E-06 0 0 3 0 + 3 -5651264739.457032 -0.7101675834366328E-08 0 0 2 1 + 3 115860293.7962667 0.2862872407163143E-09 0 0 1 2 + 3 -1130563556.160437 -0.3349768479044413E-08 0 0 0 3 + 4 13346548864539.53 -0.1557745080138484 3 0 1 0 + 4 -26537036561623.46 -0.2061883653394258E-01 2 1 1 0 + 4 27305500812673.80 0.4151835500387857 1 2 1 0 + 4 21918606359401.09 0.7481786270890803E-02 0 3 1 0 + 4 122037554755245.6 -0.1760544827318826E-02 3 0 0 1 + 4 -28347880192628.72 -0.1197110866802704 2 1 0 1 + 4 -367088948290105.8 0.9453167685375408E-01 1 2 0 1 + 4 9485398715408.410 0.4488157291143844E-01 0 3 0 1 + 4 -3073643569331.809 0.1844844205583177E-02 1 0 3 0 + 4 -1608015085289.590 0.2360692449615059E-03 0 1 3 0 + 4 -19660767373028.77 0.9946999481608755E-03 1 0 2 1 + 4 33854686125061.96 0.2324162059549456E-02 0 1 2 1 + 4 -5628764374038.941 -0.3740525005393164E-02 1 0 1 2 + 4 -3340565271617.783 -0.9409154486250337E-03 0 1 1 2 + 4 43967058034.96068 0.6947616100005572E-03 1 0 0 3 + 4 5531129402.212612 -0.9342022344190524E-03 0 1 0 3 + 5 0.8341046824785501E+19 240316.4925194468 4 0 1 0 + 5 0.4142053007329714E+19 134.0807029670189 3 1 1 0 + 5 0.2167143390915808E+20 559450.3410149341 2 2 1 0 + 5 0.3641522299834058E+19 768.2734631977004 1 3 1 0 + 5 0.1436938665669512E+20 311184.5667206588 0 4 1 0 + 5 -0.2980431309482988E+21 30160.69195661116 4 0 0 1 + 5 0.3611045474259032E+17 7656.778323628073 3 1 0 1 + 5 -0.5960557992305767E+21 68185.79090572287 2 2 0 1 + 5 -0.5539627068338006E+17 23307.77522451836 1 3 0 1 + 5 -0.2980693099673572E+21 39583.93176163451 0 4 0 1 + 5 -0.3810749247408925E+19 -738742.3460133835 2 0 3 0 + 5 -0.9867112024333270E+18 1222.564545618250 1 1 3 0 + 5 -0.5117172977636930E+19 -742460.8675356501 0 2 3 0 + 5 0.5760003043902553E+21 -290796.1991201608 2 0 2 1 + 5 -0.2576199666114370E+20 475.5626544338864 1 1 2 1 + 5 0.5927159656849820E+21 -293579.1612405238 0 2 2 1 + 5 -0.1428930283274378E+20 293775.0950197216 2 0 1 2 + 5 -0.1197058654552479E+19 2.423314411086892 1 1 1 2 + 5 -0.1656346077496378E+20 273785.6441157568 0 2 1 2 + 5 0.1169262304785450E+21 6312.878199635295 2 0 0 3 + 5 7083864041890351. -500.4842401388141 1 1 0 3 + 5 0.1169331478099869E+21 3881.489917571227 0 2 0 3 + 5 0.4386516036279700E+18 4628.740378859817 0 0 5 0 + 5 -0.1019475978703260E+21 160414.3974204903 0 0 4 1 + 5 0.2737169939775908E+19 -18121.89358738047 0 0 3 2 + 5 -0.1372514448100977E+21 246810.9497717834 0 0 2 3 + 5 0.2889967459079623E+19 21172.81729419631 0 0 1 4 + 5 -0.1147126570733387E+20 -18453.88019700950 0 0 0 5 + 6 -0.6339776822363041E+24 3225046129.553318 5 0 1 0 + 6 -0.2918865279736006E+24 2362081549.280466 4 1 1 0 + 6 -0.2117709291549715E+25 -8302375040.915841 3 2 1 0 + 6 -0.1064128564385389E+25 1738838802.317907 2 3 1 0 + 6 -0.1776262355429691E+25 -10616814083.77897 1 4 1 0 + 6 -0.5503867227986099E+24 -628245321.1628896 0 5 1 0 + 6 -0.6521909694433115E+25 1880724552.916775 5 0 0 1 + 6 0.1508290465625368E+25 -5214670382.989807 4 1 0 1 + 6 0.1309584473303808E+26 -4663278568.200044 3 2 0 1 + 6 0.1029072252479868E+25 -6640520361.768922 2 3 0 1 + 6 0.1959596219453924E+26 -7186635882.163349 1 4 0 1 + 6 -0.5103389663276995E+24 3403367988.590848 0 5 0 1 + 6 0.7363895191797020E+24 -9218421485.381306 3 0 3 0 + 6 0.1276575000790429E+26 -609255858.2375510 2 1 3 0 + 6 -0.3086867018794760E+24 24137458684.92014 1 2 3 0 + 6 -0.3909761116898010E+25 2842646220.356518 0 3 3 0 + 6 0.5225914945324512E+25 62993030.06209031 3 0 2 1 + 6 -0.4307600442951916E+25 -22500569720.97846 2 1 2 1 + 6 -0.8196231392946305E+25 212694688.0019633 1 2 2 1 + 6 -0.2842377096818493E+25 9507467499.182934 0 3 2 1 + 6 0.9403613873021098E+24 13850239650.41151 3 0 1 2 + 6 0.1146822333621040E+26 12730803638.45024 2 1 1 2 + 6 0.4421910671091353E+24 -34515128938.13533 1 2 1 2 + 6 -0.3167952673974035E+25 -8926009058.382936 0 3 1 2 + 6 0.6927846803534563E+24 -635282206.9495819 3 0 0 3 + 6 -0.4950502934730613E+24 -8520839401.215440 2 1 0 3 + 6 -0.2100161624747597E+25 6312498719.120423 1 2 0 3 + 6 0.1634368739078485E+24 4264242009.279778 0 3 0 3 + 6 -0.5726512974477809E+23 179962798.8699778 1 0 5 0 + 6 -0.2978371820676109E+23 -6744096.150854078 0 1 5 0 + 6 -0.3811132696977060E+24 473127785.9176015 1 0 4 1 + 6 0.6691339332915373E+24 15233123.64401768 0 1 4 1 + 6 0.5064601576165629E+24 -982203636.5766430 1 0 3 2 + 6 0.2968975116029379E+24 286483582.2729832 0 1 3 2 + 6 -0.5162798010614871E+24 -207648349.2522964 1 0 2 3 + 6 0.8812709384286934E+24 -161206792.0185730 0 1 2 3 + 6 -0.1356680656880316E+24 480039383.5472942 1 0 1 4 + 6 -0.8341704897438544E+23 -40334408.34316671 0 1 1 4 + 6 0.7618140803233211E+21 -38842702.12548141 1 0 0 5 + 6 0.2438386986667018E+21 -7589801.381268516 0 1 0 5 + -100 0.000000000000000 0.000000000000000 0 0 0 0 + No Spin Matrix + c_quaternion is identity \ No newline at end of file diff --git a/test/3nux_leavein/test.map b/test/3nux_leavein/test.map new file mode 100644 index 0000000..de985da --- /dev/null +++ b/test/3nux_leavein/test.map @@ -0,0 +1,450 @@ + 4 Dimensional TPSA/DA map + + 1, NO = 6, NV = 4, INA = 209 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 6 NV = 4 + 1 -0.4504032720436589 0.000000000000000 1 0 0 0 + 1 -8.012693649976276 0.000000000000000 0 1 0 0 + 2 13.51648819944545 0.000000000000000 2 0 0 0 + 2 -211.3243061376784 0.000000000000000 1 1 0 0 + 2 1533.605381176273 0.000000000000000 0 2 0 0 + 2 36.46768293637611 0.000000000000000 0 0 2 0 + 2 -105.0842516254893 0.000000000000000 0 0 1 1 + 2 -212.5484183849933 0.000000000000000 0 0 0 2 + 3 200.6500469002602 0.000000000000000 3 0 0 0 + 3 457.8340062500000 0.000000000000000 2 1 0 0 + 3 -232432.4109763060 0.000000000000000 1 2 0 0 + 3 740567.6228758973 0.000000000000000 0 3 0 0 + 3 -3833.398768698483 0.000000000000000 1 0 2 0 + 3 29710.42783986440 0.000000000000000 0 1 2 0 + 3 9028.565526505678 0.000000000000000 1 0 1 1 + 3 -5226.992427299499 0.000000000000000 0 1 1 1 + 3 -51172.35172244036 0.000000000000000 1 0 0 2 + 3 200914.1940148843 0.000000000000000 0 1 0 2 + 4 -602.3970021897733 0.000000000000000 4 0 0 0 + 4 -1167806.206589617 0.000000000000000 3 1 0 0 + 4 4199036.924561564 0.000000000000000 2 2 0 0 + 4 36967582.95864053 0.000000000000000 1 3 0 0 + 4 296987447.2780804 0.000000000000000 0 4 0 0 + 4 87412.37732071156 0.000000000000000 2 0 2 0 + 4 1024372.034353754 0.000000000000000 1 1 2 0 + 4 -3502057.048161834 0.000000000000000 0 2 2 0 + 4 602339.5868991513 0.000000000000000 2 0 1 1 + 4 -28692586.36683895 0.000000000000000 1 1 1 1 + 4 -34372579.28765996 0.000000000000000 0 2 1 1 + 4 -3572931.157409074 0.000000000000000 2 0 0 2 + 4 -31701559.07935280 0.000000000000000 1 1 0 2 + 4 372470143.4355974 0.000000000000000 0 2 0 2 + 4 -20277.25356687677 0.000000000000000 0 0 4 0 + 4 -90954.40051067325 0.000000000000000 0 0 3 1 + 4 4936550.143950218 0.000000000000000 0 0 2 2 + 4 -12065425.16042644 0.000000000000000 0 0 1 3 + 4 15844169.87036688 0.000000000000000 0 0 0 4 + 5 -1489551.086723091 0.000000000000000 5 0 0 0 + 5 4639373.250666798 0.000000000000000 4 1 0 0 + 5 238307480.9429767 0.000000000000000 3 2 0 0 + 5 3204784906.780087 0.000000000000000 2 3 0 0 + 5 -3657072551.559564 0.000000000000000 1 4 0 0 + 5 -100798941270.2409 0.000000000000000 0 5 0 0 + 5 5441984.195880674 0.000000000000000 3 0 2 0 + 5 -204136896.1057594 0.000000000000000 2 1 2 0 + 5 1147042664.066144 0.000000000000000 1 2 2 0 + 5 -10841223147.04718 0.000000000000000 0 3 2 0 + 5 -86012270.22508575 0.000000000000000 3 0 1 1 + 5 -853090533.9657284 0.000000000000000 2 1 1 1 + 5 2905414444.545380 0.000000000000000 1 2 1 1 + 5 80769423508.14064 0.000000000000000 0 3 1 1 + 5 -280938040.9906927 0.000000000000000 3 0 0 2 + 5 2635390687.961698 0.000000000000000 2 1 0 2 + 5 20674762226.86228 0.000000000000000 1 2 0 2 + 5 117635951333.7445 0.000000000000000 0 3 0 2 + 5 8468356.626042787 0.000000000000000 1 0 4 0 + 5 -94456011.35247424 0.000000000000000 0 1 4 0 + 5 31201767.50763385 0.000000000000000 1 0 3 1 + 5 485072837.6581241 0.000000000000000 0 1 3 1 + 5 -206430122.1613539 0.000000000000000 1 0 2 2 + 5 -7177483107.352360 0.000000000000000 0 1 2 2 + 5 -2920752406.910351 0.000000000000000 1 0 1 3 + 5 6520911295.903579 0.000000000000000 0 1 1 3 + 5 489719082.9885437 0.000000000000000 1 0 0 4 + 5 41889862436.02633 0.000000000000000 0 1 0 4 + 6 -1704807.758678638 0.000000000000000 6 0 0 0 + 6 608407397.3195839 0.000000000000000 5 1 0 0 + 6 12407878119.11213 0.000000000000000 4 2 0 0 + 6 -11402727408.13223 0.000000000000000 3 3 0 0 + 6 -1282789313348.400 0.000000000000000 2 4 0 0 + 6 934999445722.6575 0.000000000000000 1 5 0 0 + 6 6574701441265.635 0.000000000000000 0 6 0 0 + 6 -524658435.6997423 0.000000000000000 4 0 2 0 + 6 4606534392.618453 0.000000000000000 3 1 2 0 + 6 -52621737464.55686 0.000000000000000 2 2 2 0 + 6 328532361149.8947 0.000000000000000 1 3 2 0 + 6 -147928702105.3727 0.000000000000000 0 4 2 0 + 6 -3209499639.683682 0.000000000000000 4 0 1 1 + 6 29631404975.83156 0.000000000000000 3 1 1 1 + 6 605504298864.0657 0.000000000000000 2 2 1 1 + 6 1862569620645.511 0.000000000000000 1 3 1 1 + 6 -4311882736765.268 0.000000000000000 0 4 1 1 + 6 -262534210.6771260 0.000000000000000 4 0 0 2 + 6 180586417323.1026 0.000000000000000 3 1 0 2 + 6 1744670425333.262 0.000000000000000 2 2 0 2 + 6 -5427062051765.764 0.000000000000000 1 3 0 2 + 6 -11661468519234.95 0.000000000000000 0 4 0 2 + 6 -217903649.1700615 0.000000000000000 2 0 4 0 + 6 2316762499.631013 0.000000000000000 1 1 4 0 + 6 -3846274582.731316 0.000000000000000 0 2 4 0 + 6 -4805592489.583565 0.000000000000000 2 0 3 1 + 6 7078357709.337065 0.000000000000000 1 1 3 1 + 6 -523551570299.3011 0.000000000000000 0 2 3 1 + 6 -57578165781.81641 0.000000000000000 2 0 2 2 + 6 -33398770302.67816 0.000000000000000 1 1 2 2 + 6 2402354367561.075 0.000000000000000 0 2 2 2 + 6 -94838504495.50714 0.000000000000000 2 0 1 3 + 6 1667979890345.480 0.000000000000000 1 1 1 3 + 6 20078961292000.69 0.000000000000000 0 2 1 3 + 6 311765345517.8109 0.000000000000000 2 0 0 4 + 6 3046958401287.689 0.000000000000000 1 1 0 4 + 6 21066469314549.52 0.000000000000000 0 2 0 4 + 6 -37669240.08588751 0.000000000000000 0 0 6 0 + 6 -141774292.0521045 0.000000000000000 0 0 5 1 + 6 1976016545.441733 0.000000000000000 0 0 4 2 + 6 -88736706671.56216 0.000000000000000 0 0 3 3 + 6 -217427642711.3979 0.000000000000000 0 0 2 4 + 6 408665842866.8107 0.000000000000000 0 0 1 5 + 6 1660200142523.079 0.000000000000000 0 0 0 6 + -109 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 6, NV = 4, INA = 208 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 6 NV = 4 + 1 0.8904867536452157E-01 0.000000000000000 1 0 0 0 + 1 -0.6360527601587327 0.000000000000000 0 1 0 0 + 2 25.85974841314479 0.000000000000000 2 0 0 0 + 2 923.6971840642350 0.000000000000000 1 1 0 0 + 2 8062.313669542660 0.000000000000000 0 2 0 0 + 2 -129.7068722904625 0.000000000000000 0 0 2 0 + 2 -0.5893694639479996E-01 0.000000000000000 0 0 1 1 + 2 -9.532844578208875 0.000000000000000 0 0 0 2 + 3 912128.6310134458 0.000000000000000 3 0 0 0 + 3 48759188.62400510 0.000000000000000 2 1 0 0 + 3 867781170.6836473 0.000000000000000 1 2 0 0 + 3 5141299221.573460 0.000000000000000 0 3 0 0 + 3 -13518138.71159396 0.000000000000000 1 0 2 0 + 3 -240549265.4404155 0.000000000000000 0 1 2 0 + 3 945032.8089683487 0.000000000000000 1 0 1 1 + 3 16552532.22537735 0.000000000000000 0 1 1 1 + 3 4495.253818063657 0.000000000000000 1 0 0 2 + 3 116905.7097859786 0.000000000000000 0 1 0 2 + 4 -82261410.42997341 0.000000000000000 4 0 0 0 + 4 -1641935703.988032 0.000000000000000 3 1 0 0 + 4 10429558009.34576 0.000000000000000 2 2 0 0 + 4 75264578498.13565 0.000000000000000 1 3 0 0 + 4 -2955112619989.771 0.000000000000000 0 4 0 0 + 4 539369402.2465667 0.000000000000000 2 0 2 0 + 4 -7750761488.311388 0.000000000000000 1 1 2 0 + 4 -21366238158.59440 0.000000000000000 0 2 2 0 + 4 1863870149.050254 0.000000000000000 2 0 1 1 + 4 42202243792.19951 0.000000000000000 1 1 1 1 + 4 141582632751.9200 0.000000000000000 0 2 1 1 + 4 1255531315.021000 0.000000000000000 2 0 0 2 + 4 45458822893.45134 0.000000000000000 1 1 0 2 + 4 410925323007.8293 0.000000000000000 0 2 0 2 + 4 1094544537.460690 0.000000000000000 0 0 4 0 + 4 -3228587634.841901 0.000000000000000 0 0 3 1 + 4 -6163501640.969437 0.000000000000000 0 0 2 2 + 4 451939640.5618356 0.000000000000000 0 0 1 3 + 4 -2556446.470707916 0.000000000000000 0 0 0 4 + 5 1248221490.486264 0.000000000000000 5 0 0 0 + 5 -79380822993.64948 0.000000000000000 4 1 0 0 + 5 720514644813.5026 0.000000000000000 3 2 0 0 + 5 56886164680052.62 0.000000000000000 2 3 0 0 + 5 163054047661450.2 0.000000000000000 1 4 0 0 + 5 -861370893917143.1 0.000000000000000 0 5 0 0 + 5 8469402253.546633 0.000000000000000 3 0 2 0 + 5 669688799152.4155 0.000000000000000 2 1 2 0 + 5 -4673860622679.846 0.000000000000000 1 2 2 0 + 5 526213922537.3816 0.000000000000000 0 3 2 0 + 5 15298130633.38715 0.000000000000000 3 0 1 1 + 5 1412104692056.720 0.000000000000000 2 1 1 1 + 5 -3155262918701.357 0.000000000000000 1 2 1 1 + 5 117050563576973.6 0.000000000000000 0 3 1 1 + 5 200146857851.8719 0.000000000000000 3 0 0 2 + 5 9206039537940.430 0.000000000000000 2 1 0 2 + 5 70141802633452.73 0.000000000000000 1 2 0 2 + 5 -554185394181437.6 0.000000000000000 0 3 0 2 + 5 -134578218756.7285 0.000000000000000 1 0 4 0 + 5 1043499729610.371 0.000000000000000 0 1 4 0 + 5 233062067931.4333 0.000000000000000 1 0 3 1 + 5 -400799390797.5639 0.000000000000000 0 1 3 1 + 5 -1136685487128.029 0.000000000000000 1 0 2 2 + 5 4538675707256.246 0.000000000000000 0 1 2 2 + 5 3093559755257.843 0.000000000000000 1 0 1 3 + 5 40926771689587.35 0.000000000000000 0 1 1 3 + 5 527305094101.7097 0.000000000000000 1 0 0 4 + 5 9785971106118.014 0.000000000000000 0 1 0 4 + 6 52440303443.25721 0.000000000000000 6 0 0 0 + 6 8719570195578.993 0.000000000000000 5 1 0 0 + 6 105383241300391.5 0.000000000000000 4 2 0 0 + 6 1678776070204416. 0.000000000000000 3 3 0 0 + 6 -6004933609452629. 0.000000000000000 2 4 0 0 + 6 -0.3365469546554607E+18 0.000000000000000 1 5 0 0 + 6 -0.6167778245254751E+17 0.000000000000000 0 6 0 0 + 6 -1079136448124.284 0.000000000000000 4 0 2 0 + 6 -31900531585909.16 0.000000000000000 3 1 2 0 + 6 562493048202285.4 0.000000000000000 2 2 2 0 + 6 -7574946597871215. 0.000000000000000 1 3 2 0 + 6 0.2932558395648182E+17 0.000000000000000 0 4 2 0 + 6 -1861054615830.121 0.000000000000000 4 0 1 1 + 6 150332726759590.4 0.000000000000000 3 1 1 1 + 6 3870360264967518. 0.000000000000000 2 2 1 1 + 6 0.4787896972248809E+17 0.000000000000000 1 3 1 1 + 6 7596001548756252. 0.000000000000000 0 4 1 1 + 6 -3607443625736.787 0.000000000000000 4 0 0 2 + 6 795958181003882.6 0.000000000000000 3 1 0 2 + 6 0.1659266524150417E+17 0.000000000000000 2 2 0 2 + 6 -0.6063773241398414E+17 0.000000000000000 1 3 0 2 + 6 -0.6075807924565674E+18 0.000000000000000 0 4 0 2 + 6 4905471704208.505 0.000000000000000 2 0 4 0 + 6 -20034010135951.78 0.000000000000000 1 1 4 0 + 6 311285148838893.3 0.000000000000000 0 2 4 0 + 6 5684342870749.781 0.000000000000000 2 0 3 1 + 6 -801946567693064.5 0.000000000000000 1 1 3 1 + 6 826689413961356.2 0.000000000000000 0 2 3 1 + 6 -41827790579255.98 0.000000000000000 2 0 2 2 + 6 -282041231483890.3 0.000000000000000 1 1 2 2 + 6 0.1618994000849038E+17 0.000000000000000 0 2 2 2 + 6 194328508242905.6 0.000000000000000 2 0 1 3 + 6 2928439027889932. 0.000000000000000 1 1 1 3 + 6 0.2522442939349130E+17 0.000000000000000 0 2 1 3 + 6 78457686139508.14 0.000000000000000 2 0 0 4 + 6 -268414241438016.2 0.000000000000000 1 1 0 4 + 6 -0.4997530151795554E+17 0.000000000000000 0 2 0 4 + 6 -388362536738.9982 0.000000000000000 0 0 6 0 + 6 -6701501013303.357 0.000000000000000 0 0 5 1 + 6 157273056989763.8 0.000000000000000 0 0 4 2 + 6 -515856954314578.8 0.000000000000000 0 0 3 3 + 6 953963944772813.2 0.000000000000000 0 0 2 4 + 6 947141782588311.8 0.000000000000000 0 0 1 5 + 6 66207944589668.50 0.000000000000000 0 0 0 6 + -109 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 6, NV = 4, INA = 207 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 6 NV = 4 + 1 1.000193609724560 0.000000000000000 0 0 1 0 + 1 -0.3406661955050205E-01 0.000000000000000 0 0 0 1 + 2 -13.14509092409438 0.000000000000000 1 0 1 0 + 2 -6.011614804723042 0.000000000000000 0 1 1 0 + 2 -46.77023710100325 0.000000000000000 1 0 0 1 + 2 119.9001570308618 0.000000000000000 0 1 0 1 + 3 783.8299280980463 0.000000000000000 2 0 1 0 + 3 -8749.689880753658 0.000000000000000 1 1 1 0 + 3 -19096.71666945654 0.000000000000000 0 2 1 0 + 3 -6012.823149106189 0.000000000000000 2 0 0 1 + 3 -3588.839508586664 0.000000000000000 1 1 0 1 + 3 -363588.8007774823 0.000000000000000 0 2 0 1 + 3 322.2939549643885 0.000000000000000 0 0 3 0 + 3 -2811.474793640616 0.000000000000000 0 0 2 1 + 3 -642.4509821365777 0.000000000000000 0 0 1 2 + 3 -66728.50862305475 0.000000000000000 0 0 0 3 + 4 -1750.816131563152 0.000000000000000 3 0 1 0 + 4 -959726.8798164810 0.000000000000000 2 1 1 0 + 4 -3999400.642509552 0.000000000000000 1 2 1 0 + 4 34524537.14517007 0.000000000000000 0 3 1 0 + 4 -169817.4029644272 0.000000000000000 3 0 0 1 + 4 -671934.9954707411 0.000000000000000 2 1 0 1 + 4 57583590.31716622 0.000000000000000 1 2 0 1 + 4 -7339014.498590179 0.000000000000000 0 3 0 1 + 4 -43288.27016541988 0.000000000000000 1 0 3 0 + 4 126594.9489489179 0.000000000000000 0 1 3 0 + 4 699809.9130055548 0.000000000000000 1 0 2 1 + 4 -7614373.744621934 0.000000000000000 0 1 2 1 + 4 -4126029.227410547 0.000000000000000 1 0 1 2 + 4 -22935749.27084871 0.000000000000000 0 1 1 2 + 4 3091338.020641199 0.000000000000000 1 0 0 3 + 4 -31490675.46997997 0.000000000000000 0 1 0 3 + 5 -2196971.521133668 0.000000000000000 4 0 1 0 + 5 -15403983.93482250 0.000000000000000 3 1 1 0 + 5 91266133.01032060 0.000000000000000 2 2 1 0 + 5 7619970028.405921 0.000000000000000 1 3 1 0 + 5 -2530378317.655665 0.000000000000000 0 4 1 0 + 5 -4899568.876740783 0.000000000000000 4 0 0 1 + 5 336196336.3108041 0.000000000000000 3 1 0 1 + 5 338287312.2019672 0.000000000000000 2 2 0 1 + 5 -11205775714.83471 0.000000000000000 1 3 0 1 + 5 -131804683130.4673 0.000000000000000 0 4 0 1 + 5 411175.4158843628 0.000000000000000 2 0 3 0 + 5 45956658.71560873 0.000000000000000 1 1 3 0 + 5 -419205751.0552785 0.000000000000000 0 2 3 0 + 5 -37008516.23440759 0.000000000000000 2 0 2 1 + 5 -1250814098.049803 0.000000000000000 1 1 2 1 + 5 975588431.7816904 0.000000000000000 0 2 2 1 + 5 -324015906.0133275 0.000000000000000 2 0 1 2 + 5 5063714963.191343 0.000000000000000 1 1 1 2 + 5 22018494934.20858 0.000000000000000 0 2 1 2 + 5 547127446.2336122 0.000000000000000 2 0 0 3 + 5 8361829872.178164 0.000000000000000 1 1 0 3 + 5 -38022998347.49440 0.000000000000000 0 2 0 3 + 5 -157566.3896060446 0.000000000000000 0 0 5 0 + 5 -14381338.55084014 0.000000000000000 0 0 4 1 + 5 22850051.24461664 0.000000000000000 0 0 3 2 + 5 -1537045138.638453 0.000000000000000 0 0 2 3 + 5 631393386.6922882 0.000000000000000 0 0 1 4 + 5 -2668278455.428359 0.000000000000000 0 0 0 5 + 6 -1826169.746091389 0.000000000000000 5 0 1 0 + 6 -337489388.9044386 0.000000000000000 4 1 1 0 + 6 45483563684.28324 0.000000000000000 3 2 1 0 + 6 -35737806340.70488 0.000000000000000 2 3 1 0 + 6 4946503059.937771 0.000000000000000 1 4 1 0 + 6 -17779000571119.54 0.000000000000000 0 5 1 0 + 6 344898707.7331856 0.000000000000000 5 0 0 1 + 6 2224148846.655696 0.000000000000000 4 1 0 1 + 6 6796887080.286690 0.000000000000000 3 2 0 1 + 6 -1885872489229.926 0.000000000000000 2 3 0 1 + 6 3018738185163.818 0.000000000000000 1 4 0 1 + 6 26982194088314.40 0.000000000000000 0 5 0 1 + 6 39178634.06589998 0.000000000000000 3 0 3 0 + 6 -7200902179.262613 0.000000000000000 2 1 3 0 + 6 28311636930.47040 0.000000000000000 1 2 3 0 + 6 -319878410493.8215 0.000000000000000 0 3 3 0 + 6 -3484833465.522489 0.000000000000000 3 0 2 1 + 6 2094435213.892307 0.000000000000000 2 1 2 1 + 6 796301061842.4338 0.000000000000000 1 2 2 1 + 6 4065607273712.190 0.000000000000000 0 3 2 1 + 6 6801280232.621964 0.000000000000000 3 0 1 2 + 6 372837652752.2585 0.000000000000000 2 1 1 2 + 6 1053971977949.150 0.000000000000000 1 2 1 2 + 6 -23392800720935.33 0.000000000000000 0 3 1 2 + 6 61566487123.31325 0.000000000000000 3 0 0 3 + 6 -249932053233.1758 0.000000000000000 2 1 0 3 + 6 -3744557211070.133 0.000000000000000 1 2 0 3 + 6 -32582144824872.22 0.000000000000000 0 3 0 3 + 6 68322527.82914048 0.000000000000000 1 0 5 0 + 6 -1355423661.648436 0.000000000000000 0 1 5 0 + 6 603294700.2358246 0.000000000000000 1 0 4 1 + 6 5071084629.074301 0.000000000000000 0 1 4 1 + 6 -26104987238.43517 0.000000000000000 1 0 3 2 + 6 -313726469542.2896 0.000000000000000 0 1 3 2 + 6 -31277582711.92221 0.000000000000000 1 0 2 3 + 6 1360318674269.239 0.000000000000000 0 1 2 3 + 6 486363295227.9961 0.000000000000000 1 0 1 4 + 6 -3506039777.749512 0.000000000000000 0 1 1 4 + 6 92232937817.15804 0.000000000000000 1 0 0 5 + 6 -4669879832358.473 0.000000000000000 0 1 0 5 + -100 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 6, NV = 4, INA = 206 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 6 NV = 4 + 1 0.1758953291135623E-02 0.000000000000000 0 0 1 0 + 1 0.9997465177594933 0.000000000000000 0 0 0 1 + 2 110.3015246054108 0.000000000000000 1 0 1 0 + 2 2124.571692806800 0.000000000000000 0 1 1 0 + 2 9.300091551782643 0.000000000000000 1 0 0 1 + 2 -66.14317986083039 0.000000000000000 0 1 0 1 + 3 6082386.583224361 0.000000000000000 2 0 1 0 + 3 216601295.2679267 0.000000000000000 1 1 1 0 + 3 1926107656.225572 0.000000000000000 0 2 1 0 + 3 -212995.8970383699 0.000000000000000 2 0 0 1 + 3 -7455644.376630494 0.000000000000000 1 1 0 1 + 3 -65330536.91814293 0.000000000000000 0 2 0 1 + 3 -10014894.34484767 0.000000000000000 0 0 3 0 + 3 1049610.587404526 0.000000000000000 0 0 2 1 + 3 22881.78002614812 0.000000000000000 0 0 1 2 + 3 -822.4352749921588 0.000000000000000 0 0 0 3 + 4 -445264815.4249681 0.000000000000000 3 0 1 0 + 4 -3670618707.371643 0.000000000000000 2 1 1 0 + 4 33589026328.15363 0.000000000000000 1 2 1 0 + 4 -749247030707.7306 0.000000000000000 0 3 1 0 + 4 -272722105.9314457 0.000000000000000 3 0 0 1 + 4 -9385716673.792971 0.000000000000000 2 1 0 1 + 4 -66193058322.51290 0.000000000000000 1 2 0 1 + 4 255214507664.4940 0.000000000000000 0 3 0 1 + 4 -590046108.1978782 0.000000000000000 1 0 3 0 + 4 -17362285875.78894 0.000000000000000 0 1 3 0 + 4 4248104371.295781 0.000000000000000 1 0 2 1 + 4 47512843083.16444 0.000000000000000 0 1 2 1 + 4 5564054697.475570 0.000000000000000 1 0 1 2 + 4 100690562428.9566 0.000000000000000 0 1 1 2 + 4 -203024752.7600704 0.000000000000000 1 0 0 3 + 4 -3613309810.766264 0.000000000000000 0 1 0 3 + 5 9628195120.488598 0.000000000000000 4 0 1 0 + 5 -150981282066.3891 0.000000000000000 3 1 1 0 + 5 7354207519274.849 0.000000000000000 2 2 1 0 + 5 60634689055229.01 0.000000000000000 1 3 1 0 + 5 -317889280062520.4 0.000000000000000 0 4 1 0 + 5 -19505928741.90552 0.000000000000000 4 0 0 1 + 5 -1321971630858.799 0.000000000000000 3 1 0 1 + 5 -17768171117195.07 0.000000000000000 2 2 0 1 + 5 -46296962913365.02 0.000000000000000 1 3 0 1 + 5 -778890153902803.9 0.000000000000000 0 4 0 1 + 5 119352348472.8810 0.000000000000000 2 0 3 0 + 5 1141689432093.531 0.000000000000000 1 1 3 0 + 5 -9630886549203.328 0.000000000000000 0 2 3 0 + 5 -196825147394.1696 0.000000000000000 2 0 2 1 + 5 -3299512249739.670 0.000000000000000 1 1 2 1 + 5 -3693399918834.738 0.000000000000000 0 2 2 1 + 5 933334941649.7539 0.000000000000000 2 0 1 2 + 5 18754478232431.60 0.000000000000000 1 1 1 2 + 5 -113054072634939.6 0.000000000000000 0 2 1 2 + 5 -713783375065.7948 0.000000000000000 2 0 0 3 + 5 -19293376668851.97 0.000000000000000 1 1 0 3 + 5 -112366779945212.1 0.000000000000000 0 2 0 3 + 5 30233909358.97941 0.000000000000000 0 0 5 0 + 5 -146241503651.0735 0.000000000000000 0 0 4 1 + 5 -113750427771.5203 0.000000000000000 0 0 3 2 + 5 3350056846235.519 0.000000000000000 0 0 2 3 + 5 1168330859821.891 0.000000000000000 0 0 1 4 + 5 -47683652175.39589 0.000000000000000 0 0 0 5 + 6 -118583081145.1229 0.000000000000000 5 0 1 0 + 6 29715597476703.57 0.000000000000000 4 1 1 0 + 6 -14345031767634.37 0.000000000000000 3 2 1 0 + 6 -4023940338744820. 0.000000000000000 2 3 1 0 + 6 -0.4817259529701473E+17 0.000000000000000 1 4 1 0 + 6 7652331795798241. 0.000000000000000 0 5 1 0 + 6 1154376813420.622 0.000000000000000 5 0 0 1 + 6 -22707520084976.76 0.000000000000000 4 1 0 1 + 6 -795074211385608.9 0.000000000000000 3 2 0 1 + 6 0.1327378969412048E+17 0.000000000000000 2 3 0 1 + 6 0.1139661441117794E+18 0.000000000000000 1 4 0 1 + 6 0.2220479472496324E+18 0.000000000000000 0 5 0 1 + 6 -7243692363723.579 0.000000000000000 3 0 3 0 + 6 3892525997716.647 0.000000000000000 2 1 3 0 + 6 -1258664491487304. 0.000000000000000 1 2 3 0 + 6 5404398133531487. 0.000000000000000 0 3 3 0 + 6 4269323675093.062 0.000000000000000 3 0 2 1 + 6 606726880591276.4 0.000000000000000 2 1 2 1 + 6 0.1527135667569204E+17 0.000000000000000 1 2 2 1 + 6 3006102012504674. 0.000000000000000 0 3 2 1 + 6 -1315910676250.283 0.000000000000000 3 0 1 2 + 6 1879288180080792. 0.000000000000000 2 1 1 2 + 6 -0.1433413291602494E+17 0.000000000000000 1 2 1 2 + 6 -0.2320051045174506E+18 0.000000000000000 0 3 1 2 + 6 -48223282690285.11 0.000000000000000 3 0 0 3 + 6 -1104338260484292. 0.000000000000000 2 1 0 3 + 6 751803838024012.6 0.000000000000000 1 2 0 3 + 6 -0.5689709767631211E+17 0.000000000000000 0 3 0 3 + 6 -7131474606658.672 0.000000000000000 1 0 5 0 + 6 65218182935846.08 0.000000000000000 0 1 5 0 + 6 29067716117372.28 0.000000000000000 1 0 4 1 + 6 139507222842585.6 0.000000000000000 0 1 4 1 + 6 -130983968396715.0 0.000000000000000 1 0 3 2 + 6 -1704623984223982. 0.000000000000000 0 1 3 2 + 6 423750436132733.7 0.000000000000000 1 0 2 3 + 6 6378526856201166. 0.000000000000000 0 1 2 3 + 6 -246448111787994.3 0.000000000000000 1 0 1 4 + 6 -0.1321716685988446E+17 0.000000000000000 0 1 1 4 + 6 -448914898839678.6 0.000000000000000 1 0 0 5 + 6 -6330036205718577. 0.000000000000000 0 1 0 5 + -100 0.000000000000000 0.000000000000000 0 0 0 0 + No Spin Matrix + c_quaternion is identity + No Stochastic Radiation diff --git a/test/runtests.jl b/test/runtests.jl index 1a7720f..796176e 100644 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -75,6 +75,11 @@ end @test norm(Σ - Σ_fpp) < tol + m = read_fpp_map("spin_res/test.map") + R_fpp = read_fpp_map("spin_res/R.map") + a = normal(m,res=[0; 1], spin_res=[-1]) # Q_y - Q_s resonance + + end diff --git a/test/spin1/a.map b/test/spin1/a.map deleted file mode 100644 index 3ad74dc..0000000 --- a/test/spin1/a.map +++ /dev/null @@ -1,336 +0,0 @@ - 4 Dimensional TPSA/DA map - - 1, NO = 3, NV = 4, INA = 347 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 1 -3.165461932401119 0.000000000000000 1 0 0 0 - 1 -0.6252284612438701 0.000000000000000 0 1 0 0 - 1 0.1328279888522940 0.000000000000000 0 0 1 0 - 1 0.2123268890064651 0.000000000000000 0 0 0 1 - 2 42.68562022110221 0.3030231144566080E-31 2 0 0 0 - 2 -200.6042486864551 -0.3701873875581051E-30 1 1 0 0 - 2 203.0876493467034 -0.2405361097240284E-30 0 2 0 0 - 2 -186.2877061563023 0.4774177842943975E-30 1 0 1 0 - 2 716.7530951733407 0.7105427357601002E-14 0 1 1 0 - 2 30.65142118654810 -0.7105427357601001E-14 1 0 0 1 - 2 -27.91265683058521 -0.4268976172681256E-30 0 1 0 1 - 2 70.45805802755858 0.1065814103640150E-13 0 0 2 0 - 2 89.74189823159145 -0.8881784197001256E-15 0 0 1 1 - 2 58.70942104212858 -0.3552713678800500E-14 0 0 0 2 - 3 -8546.211803976572 0.4547473508864647E-12 3 0 0 0 - 3 8085.723521504641 0.3994411562050484E-27 2 1 0 0 - 3 23355.43235312018 0.6110274232647549E-28 1 2 0 0 - 3 -15998.18126711272 0.3122444467003882E-27 0 3 0 0 - 3 54688.96193234755 -0.1091393642127514E-10 2 0 1 0 - 3 89245.90208894022 -0.1455191522836685E-10 1 1 1 0 - 3 -25740.15732727661 -0.3637978807091712E-11 0 2 1 0 - 3 21497.51081891500 -0.7275957614183429E-11 2 0 0 1 - 3 10112.89744737610 -0.7275957614183429E-11 1 1 0 1 - 3 29.21473707343725 0.1091393642127514E-10 0 2 0 1 - 3 97270.34949381967 0.3637978807091713E-11 1 0 2 0 - 3 -27715.78379896367 0.7275957614183426E-11 0 1 2 0 - 3 -125010.9473231767 0.2182787284255028E-10 1 0 1 1 - 3 -61295.88721892839 -0.2182787284255027E-10 0 1 1 1 - 3 -18182.58734266024 -0.3637978807091711E-11 1 0 0 2 - 3 28253.28853397982 -0.1091393642127515E-10 0 1 0 2 - 3 13167.48639131377 -0.4547473508864635E-12 0 0 3 0 - 3 -8135.451800469542 0.2728484105318788E-11 0 0 2 1 - 3 -6267.625892942094 0.5456968210637572E-11 0 0 1 2 - 3 11778.77042910094 0.1364242052659393E-11 0 0 0 3 - -34 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 348 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 1 0.4275905604281688E-01 0.000000000000000 1 0 0 0 - 1 -0.3081487940500551 0.000000000000000 0 1 0 0 - 1 -0.2945662483720033E-02 0.000000000000000 0 0 1 0 - 1 -0.2102643629566357E-01 0.000000000000000 0 0 0 1 - 2 -10.97204508551239 -0.2388601202501189E-32 2 0 0 0 - 2 -4.365942352798449 0.7627066728731646E-32 1 1 0 0 - 2 4.287577267795555 -0.7386443224127784E-31 0 2 0 0 - 2 -64.79543558893457 -0.4440892098500627E-15 1 0 1 0 - 2 -5.313871027224780 0.1776356839400251E-14 0 1 1 0 - 2 -1.022575875896327 0.4654677436365622E-31 1 0 0 1 - 2 -4.137241217989332 -0.6661338147750943E-15 0 1 0 1 - 2 1.440354095949309 -0.1665334536937735E-15 0 0 2 0 - 2 1.842337677119144 -0.8881784197001251E-15 0 0 1 1 - 2 -8.556188294465398 0.5551115123125777E-16 0 0 0 2 - 3 -1247.348148546697 0.2036381039580450E-28 3 0 0 0 - 3 386.6505504654470 -0.5325528101547588E-28 2 1 0 0 - 3 263.1638091815568 -0.1284313329208480E-27 1 2 0 0 - 3 -821.0591489034543 -0.2273736754432321E-12 0 3 0 0 - 3 579.8053391014437 -0.6821210263296963E-12 2 0 1 0 - 3 -9681.973973853046 0.9094947017729286E-12 1 1 1 0 - 3 -56.45270832571323 0.6821210263296966E-12 0 2 1 0 - 3 -1379.832885501070 -0.2273736754432320E-12 2 0 0 1 - 3 -626.1587334136013 -0.2799165464037778E-27 1 1 0 1 - 3 -2330.849716187910 0.7958078640513123E-12 0 2 0 1 - 3 -6239.950951165020 0.5790923296444820E-12 1 0 2 0 - 3 4190.021061320819 -0.1364242052659393E-11 0 1 2 0 - 3 -2804.127846144846 -0.4547473508864641E-11 1 0 1 1 - 3 11889.23054163142 -0.7034373084024993E-12 0 1 1 1 - 3 3631.378773285719 -0.1033839680530946E-11 1 0 0 2 - 3 1497.976699084921 0.1818989403545857E-11 0 1 0 2 - 3 308.8942386633769 0.5684341886080804E-12 0 0 3 0 - 3 -3226.922095826239 -0.2842170943040403E-12 0 0 2 1 - 3 -801.1148490851450 -0.2273736754432322E-12 0 0 1 2 - 3 -1059.282317292492 -0.1136868377216160E-12 0 0 0 3 - -34 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 349 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 1 0.6943408875947243E-01 0.000000000000000 1 0 0 0 - 1 0.3696011290279373E-01 0.000000000000000 0 1 0 0 - 1 -1.308057866358130 0.000000000000000 0 0 1 0 - 1 -1.536966698285183 0.000000000000000 0 0 0 1 - 2 -162.7602439364891 0.2233576285734730E-31 2 0 0 0 - 2 44.45112862418956 0.9941875077202168E-31 1 1 0 0 - 2 -180.4514139347281 -0.1288092018379228E-30 0 2 0 0 - 2 41.76283679950394 -0.2295727204663649E-30 1 0 1 0 - 2 -11.39654077210529 0.5329070518200753E-14 0 1 1 0 - 2 72.30754143393978 -0.5329070518200751E-14 1 0 0 1 - 2 8.414719135902047 0.5982440249417986E-31 0 1 0 1 - 2 -191.8530549632474 -0.7105427357601002E-14 0 0 2 0 - 2 91.77538343289859 0.4123651178881686E-30 0 0 1 1 - 2 -67.06861700041779 -0.7105427357601001E-14 0 0 0 2 - 3 12376.15968638786 0.1818989403545856E-11 3 0 0 0 - 3 -36186.01983002596 -0.4547473508864642E-11 2 1 0 0 - 3 485.9393218717656 -0.1455191522836685E-10 1 2 0 0 - 3 10947.44381290157 -0.9094947017729280E-12 0 3 0 0 - 3 -5117.094072788933 0.2182787284255028E-10 2 0 1 0 - 3 -76904.77971426502 0.2644382160027355E-26 1 1 1 0 - 3 7375.235929785207 -0.2182787284255028E-10 0 2 1 0 - 3 -38738.43704727563 0.9094947017729339E-12 2 0 0 1 - 3 40654.96202837688 -0.4001776687800885E-10 1 1 0 1 - 3 10319.91298595206 0.2728484105318785E-11 0 2 0 1 - 3 14335.10169939692 0.4547473508864640E-12 1 0 2 0 - 3 -13126.89440902276 -0.1818989403545856E-11 0 1 2 0 - 3 1334.860300846432 -0.3637978807091712E-11 1 0 1 1 - 3 -4992.863499130713 0.3637978807091712E-11 0 1 1 1 - 3 2630.505730128729 0.3183231456205248E-11 1 0 0 2 - 3 10275.54715117433 0.1818989403545856E-11 0 1 0 2 - 3 5825.814150002814 0.5684341886080800E-12 0 0 3 0 - 3 -3952.048119102576 -0.4547473508864642E-12 0 0 2 1 - 3 -521.6505457911165 0.1023181539494544E-11 0 0 1 2 - 3 -4798.697405839369 -0.2273736754432323E-12 0 0 0 3 - -34 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 350 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 1 -0.2175801382902183E-01 0.000000000000000 1 0 0 0 - 1 -0.4279791951483362E-01 0.000000000000000 0 1 0 0 - 1 0.3574030933540097 0.000000000000000 0 0 1 0 - 1 -0.3462008946181011 0.000000000000000 0 0 0 1 - 2 -38.45102617425694 -0.5551115123125781E-16 2 0 0 0 - 2 -1.610657819485401 -0.2220446049250312E-15 1 1 0 0 - 2 -35.13284270826650 0.5551115123125783E-16 0 2 0 0 - 2 4.168006404569821 -0.1248778997469990E-30 1 0 1 0 - 2 -12.61878139530729 -0.8881784197001252E-15 0 1 1 0 - 2 -12.43518496266496 -0.2664535259100375E-14 1 0 0 1 - 2 -25.25182559050909 -0.1776356839400250E-14 0 1 0 1 - 2 -43.49919301538070 -0.1776356839400250E-14 0 0 2 0 - 2 -28.08788714854483 -0.4875272756633948E-30 0 0 1 1 - 2 -11.92877696015589 -0.1776356839400250E-14 0 0 0 2 - 3 -260.4025828774706 -0.3579290961003561E-28 3 0 0 0 - 3 -11688.98238440336 0.1818989403545856E-11 2 1 0 0 - 3 -4201.376672459543 -0.2990076439485527E-27 1 2 0 0 - 3 -69.13322642267048 -0.7958078640513126E-12 0 3 0 0 - 3 -11099.68999782649 0.1136868377216162E-12 2 0 1 0 - 3 7114.159094835281 -0.7275957614183426E-11 1 1 1 0 - 3 5327.546231476639 0.5684341886080816E-12 0 2 1 0 - 3 -1131.596541937585 -0.3637978807091713E-11 2 0 0 1 - 3 16444.91180593867 0.2273736754432319E-12 1 1 0 1 - 3 -404.7714832423368 0.3637978807091714E-11 0 2 0 1 - 3 -5676.148455085884 -0.2273736754432320E-12 1 0 2 0 - 3 -4464.936167118038 -0.4547473508864638E-12 0 1 2 0 - 3 -11686.34096684730 -0.1364242052659392E-11 1 0 1 1 - 3 8507.151651365233 -0.9094947017729277E-12 0 1 1 1 - 3 -1957.942406578657 -0.4547473508864639E-12 1 0 0 2 - 3 -3316.219635749429 0.4547473508864641E-12 0 1 0 2 - 3 -2682.985710313675 0.1136868377216160E-12 0 0 3 0 - 3 209.2739804527577 -0.3193254993482035E-28 0 0 2 1 - 3 -19.04197541296079 -0.1136868377216160E-12 0 0 1 2 - 3 -618.9295968078433 0.2500449022466576E-28 0 0 0 3 - -34 0.000000000000000 0.000000000000000 0 0 0 0 - No Spin Matrix - Quaternion - c_quaternion - - 1, NO = 3, NV = 4, INA = 360 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 0.5269975716308092 0.000000000000000 0 0 0 0 - 1 0.6618296125517910E-01 0.000000000000000 1 0 0 0 - 1 -0.5378957350051947 0.000000000000000 0 1 0 0 - 1 -0.1032943032044960 0.000000000000000 0 0 1 0 - 1 -0.2337677165057460E-01 0.000000000000000 0 0 0 1 - 2 -15.62899676848601 -0.8881784197001252E-15 2 0 0 0 - 2 -8.435195078094683 0.1776356839400250E-14 1 1 0 0 - 2 15.04855283234293 0.8881784197001249E-15 0 2 0 0 - 2 -117.7509855061457 -0.2664535259100374E-14 1 0 1 0 - 2 -8.390950760287767 0.2502168183747897E-29 0 1 1 0 - 2 -8.447739281253293 -0.3552713678800501E-14 1 0 0 1 - 2 -11.55246120615549 -0.2664535259100376E-14 0 1 0 1 - 2 7.009831431700046 -0.3885780586188050E-15 0 0 2 0 - 2 1.029476086559878 -0.2664535259100375E-14 0 0 1 1 - 2 -8.253414313076640 0.4996003610813201E-15 0 0 0 2 - 3 -8582.359369915690 0.3410605131648485E-12 3 0 0 0 - 3 -4510.370709467620 -0.7275957614183428E-11 2 1 0 0 - 3 16221.54255947597 -0.4547473508864640E-12 1 2 0 0 - 3 -356.6457305494384 0.1364242052659393E-11 0 3 0 0 - 3 -13724.13945814193 -0.4547473508864641E-11 2 0 1 0 - 3 -27888.91313889303 -0.1136868377216150E-12 1 1 1 0 - 3 -13950.30528439172 -0.9094947017729278E-12 0 2 1 0 - 3 -8952.095251187297 -0.1023181539494546E-11 2 0 0 1 - 3 -1236.935415391401 0.1818989403545855E-11 1 1 0 1 - 3 3348.144655738395 0.2842170943040401E-11 0 2 0 1 - 3 -10517.22112572823 0.2629008122312370E-11 1 0 2 0 - 3 7255.031962390278 0.3637978807091708E-11 0 1 2 0 - 3 -3818.955651436335 -0.5456968210637570E-11 1 0 1 1 - 3 25776.14878082439 0.1335820343228987E-11 0 1 1 1 - 3 4616.874073323601 0.9947598300641374E-13 1 0 0 2 - 3 2627.072078790962 0.2728484105318784E-11 0 1 0 2 - 3 -8398.233413557815 0.8526512829121207E-12 0 0 3 0 - 3 -5006.075938650246 -0.4547473508864643E-12 0 0 2 1 - 3 -9843.238317753156 -0.4547473508864645E-12 0 0 1 2 - 3 -1906.670261476162 -0.1421085471520197E-12 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 361 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 0.1539220957562536E-02 0.000000000000000 0 0 0 0 - 1 -0.4774720224267424E-01 0.000000000000000 1 0 0 0 - 1 0.2056523162158802 0.000000000000000 0 1 0 0 - 1 1.094311025597547 0.000000000000000 0 0 1 0 - 1 1.070601058997268 0.000000000000000 0 0 0 1 - 2 92.31730625964889 -0.9769962616701373E-14 2 0 0 0 - 2 -23.64549761836914 0.5329070518200751E-14 1 1 0 0 - 2 120.8352587548097 -0.2664535259100375E-14 0 2 0 0 - 2 1.253972413732481 -0.5773159728050812E-14 1 0 1 0 - 2 40.36481996537252 0.7105427357600999E-14 0 1 1 0 - 2 -55.92742034574643 0.7105427357600999E-14 1 0 0 1 - 2 167.0950241186235 -0.8437694987151191E-14 0 1 0 1 - 2 128.9626044810322 0.2486899575160350E-13 0 0 2 0 - 2 -81.10931204136922 0.7105427357601005E-14 0 0 1 1 - 2 31.09173237418942 -0.1065814103640151E-13 0 0 0 2 - 3 -10938.34805122015 0.4547473508864649E-12 3 0 0 0 - 3 31659.30484001551 -0.1421527351208263E-26 2 1 0 0 - 3 9916.510840689503 0.1091393642127514E-10 1 2 0 0 - 3 9398.123981895818 0.3183231456205248E-11 0 3 0 0 - 3 -14813.82737034286 -0.7275957614183427E-10 2 0 1 0 - 3 54885.44211883304 -0.1728039933368563E-10 1 1 1 0 - 3 -17249.94595176070 -0.5093170329928395E-10 0 2 1 0 - 3 15705.01130630442 -0.9094947017729272E-12 2 0 0 1 - 3 -25579.43118830498 0.2182787284255029E-10 1 1 0 1 - 3 -9168.084714709228 0.2091837814077735E-10 0 2 0 1 - 3 -8837.603442355921 -0.3183231456205243E-11 1 0 2 0 - 3 16196.30960089367 -0.8077935669463161E-27 0 1 2 0 - 3 -15347.97049277173 -0.2728484105318779E-11 1 0 1 1 - 3 -20188.37024968313 -0.9094947017729274E-12 0 1 1 1 - 3 -5715.990499986716 0.3183231456205257E-11 1 0 0 2 - 3 29172.30219051577 0.1364242052659390E-11 0 1 0 2 - 3 -14675.78534966607 -0.3120703695458358E-10 0 0 3 0 - 3 -1402.910603710945 -0.1000444171950220E-10 0 0 2 1 - 3 -10599.48710207822 -0.3115019353572279E-10 0 0 1 2 - 3 1174.719348914304 -0.1000444171950222E-10 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 362 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 -0.8498653953974559 0.000000000000000 0 0 0 0 - 1 0.4095326925763497E-01 0.000000000000000 1 0 0 0 - 1 -0.3331741747771645 0.000000000000000 0 1 0 0 - 1 -0.6207037111177373E-01 0.000000000000000 0 0 1 0 - 1 -0.1255682413122932E-01 0.000000000000000 0 0 0 1 - 2 -9.502616341057701 -0.4440892098500625E-15 2 0 0 0 - 2 -5.313455422623803 0.1776356839400250E-14 1 1 0 0 - 2 9.823781992346504 0.4440892098500624E-15 0 2 0 0 - 2 -73.18248459190309 -0.1332267629550188E-14 1 0 1 0 - 2 -4.859570211080615 -0.7105427357600998E-14 0 1 1 0 - 2 -5.585951728947380 -0.3552713678800503E-14 1 0 0 1 - 2 -6.743941239798468 -0.1332267629550188E-14 0 1 0 1 - 2 5.429274165673431 -0.3330669073875470E-15 0 0 2 0 - 2 2.396937634230391 -0.8881784197001247E-15 0 0 1 1 - 2 -3.883206305100393 0.3330669073875470E-15 0 0 0 2 - 3 -5368.558725933373 0.2273736754432321E-12 3 0 0 0 - 3 -2714.705308603199 -0.1818989403545856E-11 2 1 0 0 - 3 10056.96401679240 -0.6821210263296961E-12 1 2 0 0 - 3 -207.4543416992766 0.6821210263296959E-12 0 3 0 0 - 3 -8369.770702993063 -0.1364242052659393E-11 2 0 1 0 - 3 -17151.13378831840 -0.1421085471520193E-12 1 1 1 0 - 3 -8460.492768804030 -0.1364242052659392E-11 0 2 1 0 - 3 -5273.908228608519 -0.1847411112976260E-12 2 0 0 1 - 3 -853.6065462194812 -0.4078410879992631E-27 1 1 0 1 - 3 2399.965834307547 0.2003730514843482E-11 0 2 0 1 - 3 -6547.790508736306 0.1236344360222574E-11 1 0 2 0 - 3 4623.339033652176 0.4547473508864634E-12 0 1 2 0 - 3 -2517.799256615574 -0.6366462912410498E-11 1 0 1 1 - 3 16204.11539896878 0.4831690603168681E-12 0 1 1 1 - 3 2659.067097074627 0.3552713678800504E-12 1 0 0 2 - 3 1779.887095859442 0.2273736754432321E-11 0 1 0 2 - 3 -5011.891670696306 0.2842170943040402E-12 0 0 3 0 - 3 -2950.248080752798 -0.2273736754432318E-12 0 0 2 1 - 3 -6175.191725652178 -0.2273736754432323E-12 0 0 1 2 - 3 -1080.372716265485 0.8526512829121209E-13 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 363 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 1 -0.1687201381062306 -0.3851859888774472E-33 1 0 0 0 - 1 -0.1485566660458297 0.3851859888774472E-33 0 1 0 0 - 1 0.4805551321764081 -0.5551115123125781E-16 0 0 1 0 - 1 0.9252288422627845 -0.2775557561562891E-16 0 0 0 1 - 2 100.7690816572896 0.9769962616701382E-14 2 0 0 0 - 2 -30.51421247046459 0.8881784197001254E-14 1 1 0 0 - 2 109.7567011315633 -0.2664535259100377E-14 0 2 0 0 - 2 -8.347045729998380 0.7105427357601002E-14 1 0 1 0 - 2 61.90574126935259 -0.1065814103640150E-13 0 1 1 0 - 2 -100.8190889330763 0.3552713678800501E-14 1 0 0 1 - 2 -107.9038282490171 -0.7105427357601002E-14 0 1 0 1 - 2 101.8415252953052 0.1421085471520200E-13 0 0 2 0 - 2 -20.83578357891845 0.2842170943040400E-13 0 0 1 1 - 2 55.44846474977768 -0.5329070518200756E-14 0 0 0 2 - 3 -15217.27735362319 0.2728484105318784E-11 3 0 0 0 - 3 26972.27615613576 -0.1818989403545858E-11 2 1 0 0 - 3 -1746.421127839207 -0.5048709793414476E-28 1 2 0 0 - 3 -21739.06049795085 0.6821210263296950E-12 0 3 0 0 - 3 -3247.328413053327 -0.3637978807091712E-10 2 0 1 0 - 3 26332.58551617702 0.1182343112304807E-10 1 1 1 0 - 3 5307.281930918076 -0.9094947017729266E-11 0 2 1 0 - 3 5655.569645500581 -0.5138645065017044E-10 2 0 0 1 - 3 -28502.48322547848 0.3092281986027955E-10 1 1 0 1 - 3 -12493.01216624389 -0.5593392415903505E-10 0 2 0 1 - 3 -45046.38477962303 0.9094947017729250E-12 1 0 2 0 - 3 -22411.84566451363 -0.1615587133892632E-26 0 1 2 0 - 3 -25220.26592321765 0.3637978807091706E-11 1 0 1 1 - 3 56978.77914565820 0.3637978807091715E-11 0 1 1 1 - 3 9179.938226329125 0.3637978807091715E-11 1 0 0 2 - 3 -7483.827955893589 0.1818989403545866E-11 0 1 0 2 - 3 -2149.782625976277 0.1261923898709937E-10 0 0 3 0 - 3 -7909.850557967234 -0.2523847797419875E-10 0 0 2 1 - 3 2659.135756612739 0.1273292582482100E-10 0 0 1 2 - 3 -3115.973367070202 -0.2410160959698260E-10 0 0 0 3 - -34 0.000000000000000 0.000000000000000 0 0 0 0 - No Stochastic Radiation \ No newline at end of file diff --git a/test/spin1/test.map b/test/spin1/test.map deleted file mode 100644 index c10b74d..0000000 --- a/test/spin1/test.map +++ /dev/null @@ -1,341 +0,0 @@ -4 Dimensional TPSA/DA map - - 1, NO = 3, NV = 4, INA = 206 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 0.1583591029206141E-03 0.000000000000000 0 0 0 0 - 1 -0.7133082123844348 0.000000000000000 1 0 0 0 - 1 7.760252162396466 0.000000000000000 0 1 0 0 - 1 -0.2061942199606399 0.000000000000000 0 0 1 0 - 1 -0.6710032459762874 0.000000000000000 0 0 0 1 - 2 11.08435544404222 0.000000000000000 2 0 0 0 - 2 147.5115941047698 0.000000000000000 1 1 0 0 - 2 2788.949729070067 0.000000000000000 0 2 0 0 - 2 -21.14422257782299 0.000000000000000 1 0 1 0 - 2 -399.7035173568319 0.000000000000000 0 1 1 0 - 2 -214.6804470606147 0.000000000000000 1 0 0 1 - 2 -2378.403918890681 0.000000000000000 0 1 0 1 - 2 29.00610402969371 0.000000000000000 0 0 2 0 - 2 190.2349421247565 0.000000000000000 0 0 1 1 - 2 272.1124669076335 0.000000000000000 0 0 0 2 - 3 1287.165015612604 0.000000000000000 3 0 0 0 - 3 15920.05359617343 0.000000000000000 2 1 0 0 - 3 -67574.32808239023 0.000000000000000 1 2 0 0 - 3 557763.4052516725 0.000000000000000 0 3 0 0 - 3 -1789.229112963170 0.000000000000000 2 0 1 0 - 3 22544.70047301858 0.000000000000000 1 1 1 0 - 3 63510.17314152337 0.000000000000000 0 2 1 0 - 3 -3454.628075444616 0.000000000000000 2 0 0 1 - 3 19466.55113408004 0.000000000000000 1 1 0 1 - 3 869586.4712777968 0.000000000000000 0 2 0 1 - 3 7529.324986373755 0.000000000000000 1 0 2 0 - 3 83073.47401893395 0.000000000000000 0 1 2 0 - 3 8187.948348911816 0.000000000000000 1 0 1 1 - 3 63730.40584072094 0.000000000000000 0 1 1 1 - 3 105249.1316853604 0.000000000000000 1 0 0 2 - 3 478263.1809695527 0.000000000000000 0 1 0 2 - 3 -477.5206034254476 0.000000000000000 0 0 3 0 - 3 -3158.707728920650 0.000000000000000 0 0 2 1 - 3 -33913.71680906254 0.000000000000000 0 0 1 2 - 3 21146.59285884735 0.000000000000000 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 205 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 -0.1888090075451473E-04 0.000000000000000 0 0 0 0 - 1 -0.7203636799033827E-01 0.000000000000000 1 0 0 0 - 1 -0.6361260999793732 0.000000000000000 0 1 0 0 - 1 -0.3093068624609299E-02 0.000000000000000 0 0 1 0 - 1 0.5188745160021785E-01 0.000000000000000 0 0 0 1 - 2 -0.8730028407742543 0.000000000000000 2 0 0 0 - 2 4.849108186508623 0.000000000000000 1 1 0 0 - 2 203.5862568544241 0.000000000000000 0 2 0 0 - 2 1.218891753462692 0.000000000000000 1 0 1 0 - 2 -25.06647677039923 0.000000000000000 0 1 1 0 - 2 27.50503236210960 0.000000000000000 1 0 0 1 - 2 -315.2853072006434 0.000000000000000 0 1 0 1 - 2 2.602136919108823 0.000000000000000 0 0 2 0 - 2 0.8197056136191634 0.000000000000000 0 0 1 1 - 2 -27.19275023767935 0.000000000000000 0 0 0 2 - 3 -86.47398763814840 0.000000000000000 3 0 0 0 - 3 736.2694927450048 0.000000000000000 2 1 0 0 - 3 -15080.66355770427 0.000000000000000 1 2 0 0 - 3 179101.2474619581 0.000000000000000 0 3 0 0 - 3 70.83190561331128 0.000000000000000 2 0 1 0 - 3 4895.090439509293 0.000000000000000 1 1 1 0 - 3 -31405.00678947525 0.000000000000000 0 2 1 0 - 3 812.7123304378517 0.000000000000000 2 0 0 1 - 3 -16646.44214646858 0.000000000000000 1 1 0 1 - 3 -89141.18929828558 0.000000000000000 0 2 0 1 - 3 -684.8463600691530 0.000000000000000 1 0 2 0 - 3 9165.502047414513 0.000000000000000 0 1 2 0 - 3 604.8412998371864 0.000000000000000 1 0 1 1 - 3 19901.42394506920 0.000000000000000 0 1 1 1 - 3 -9228.993731061499 0.000000000000000 1 0 0 2 - 3 219246.0293633172 0.000000000000000 0 1 0 2 - 3 -85.12190858985120 0.000000000000000 0 0 3 0 - 3 -1014.891958238344 0.000000000000000 0 0 2 1 - 3 -3524.614363378187 0.000000000000000 0 0 1 2 - 3 1799.255703016103 0.000000000000000 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 204 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 -0.4126698713745159E-02 0.000000000000000 0 0 0 0 - 1 0.6945630430047534E-02 0.000000000000000 1 0 0 0 - 1 -0.4825859270554154 0.000000000000000 0 1 0 0 - 1 0.5621490866329152 0.000000000000000 0 0 1 0 - 1 3.237416273678845 0.000000000000000 0 0 0 1 - 2 5.053111573697766 0.000000000000000 2 0 0 0 - 2 -56.01404633628905 0.000000000000000 1 1 0 0 - 2 304.8213108038999 0.000000000000000 0 2 0 0 - 2 -4.416282678528900 0.000000000000000 1 0 1 0 - 2 199.7745096981564 0.000000000000000 0 1 1 0 - 2 -4.500624373992011 0.000000000000000 1 0 0 1 - 2 324.6999425254193 0.000000000000000 0 1 0 1 - 2 -8.160587998580692 0.000000000000000 0 0 2 0 - 2 -122.8527580157155 0.000000000000000 0 0 1 1 - 2 366.0036403511446 0.000000000000000 0 0 0 2 - 3 -94.51904188010047 0.000000000000000 3 0 0 0 - 3 -3954.555250550721 0.000000000000000 2 1 0 0 - 3 55038.11541286981 0.000000000000000 1 2 0 0 - 3 248677.4430619434 0.000000000000000 0 3 0 0 - 3 3674.287221807707 0.000000000000000 2 0 1 0 - 3 -624.0580368639139 0.000000000000000 1 1 1 0 - 3 338522.3058582250 0.000000000000000 0 2 1 0 - 3 -9327.433339453630 0.000000000000000 2 0 0 1 - 3 70405.01747205325 0.000000000000000 1 1 0 1 - 3 -1196057.979318757 0.000000000000000 0 2 0 1 - 3 -54.43857690995458 0.000000000000000 1 0 2 0 - 3 -13750.96098014984 0.000000000000000 0 1 2 0 - 3 -4651.988927235378 0.000000000000000 1 0 1 1 - 3 17887.72225623483 0.000000000000000 0 1 1 1 - 3 15566.21140084499 0.000000000000000 1 0 0 2 - 3 -9064.795633094589 0.000000000000000 0 1 0 2 - 3 3100.162059062996 0.000000000000000 0 0 3 0 - 3 -2969.401825660421 0.000000000000000 0 0 2 1 - 3 53180.69159010693 0.000000000000000 0 0 1 2 - 3 -246554.7157853621 0.000000000000000 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 203 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 -0.3601273821355099E-03 0.000000000000000 0 0 0 0 - 1 -0.2494691594887103E-01 0.000000000000000 1 0 0 0 - 1 -0.1058681631046877 0.000000000000000 0 1 0 0 - 1 -0.1982130454057517 0.000000000000000 0 0 1 0 - 1 0.6601028475670008 0.000000000000000 0 0 0 1 - 2 -4.376269461405176 0.000000000000000 2 0 0 0 - 2 11.01664518299111 0.000000000000000 1 1 0 0 - 2 -457.9230187108863 0.000000000000000 0 2 0 0 - 2 2.060325443216720 0.000000000000000 1 0 1 0 - 2 68.53168350076911 0.000000000000000 0 1 1 0 - 2 30.40484154667142 0.000000000000000 1 0 0 1 - 2 9.854935514327153 0.000000000000000 0 1 0 1 - 2 -6.727843738258255 0.000000000000000 0 0 2 0 - 2 -27.32116687289210 0.000000000000000 0 0 1 1 - 2 -156.3340257896522 0.000000000000000 0 0 0 2 - 3 -71.81469861740725 0.000000000000000 3 0 0 0 - 3 1160.533026190064 0.000000000000000 2 1 0 0 - 3 25460.86606623748 0.000000000000000 1 2 0 0 - 3 1897.817977686078 0.000000000000000 0 3 0 0 - 3 518.2278078062360 0.000000000000000 2 0 1 0 - 3 1736.023000627993 0.000000000000000 1 1 1 0 - 3 71987.97423937160 0.000000000000000 0 2 1 0 - 3 3147.931219162960 0.000000000000000 2 0 0 1 - 3 -26203.75767485805 0.000000000000000 1 1 0 1 - 3 350900.1819036526 0.000000000000000 0 2 0 1 - 3 148.0121694885487 0.000000000000000 1 0 2 0 - 3 -1232.760490002236 0.000000000000000 0 1 2 0 - 3 -162.3009674600962 0.000000000000000 1 0 1 1 - 3 -28614.95172149233 0.000000000000000 0 1 1 1 - 3 -9626.069464686956 0.000000000000000 1 0 0 2 - 3 86773.78587599380 0.000000000000000 0 1 0 2 - 3 654.5148985310881 0.000000000000000 0 0 3 0 - 3 5340.101751440428 0.000000000000000 0 0 2 1 - 3 7911.504523563732 0.000000000000000 0 0 1 2 - 3 56059.49444207287 0.000000000000000 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - No Spin Matrix - Quaternion - c_quaternion - - 1, NO = 3, NV = 4, INA = 186 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 -0.8314146548813478 0.000000000000000 0 0 0 0 - 1 -0.8548546920666077E-01 0.000000000000000 1 0 0 0 - 1 -1.858780119873707 0.000000000000000 0 1 0 0 - 1 0.9503147677543166E-02 0.000000000000000 0 0 1 0 - 1 0.1550731793655365 0.000000000000000 0 0 0 1 - 2 -7.469631516119694 0.000000000000000 2 0 0 0 - 2 14.98268513244695 0.000000000000000 1 1 0 0 - 2 -486.1939105883108 0.000000000000000 0 2 0 0 - 2 0.6423143581375751 0.000000000000000 1 0 1 0 - 2 -11.39635129555268 0.000000000000000 0 1 1 0 - 2 35.28524549600866 0.000000000000000 1 0 0 1 - 2 -327.3367606341089 0.000000000000000 0 1 0 1 - 2 -0.4143168856506030 0.000000000000000 0 0 2 0 - 2 -5.438079808215720 0.000000000000000 0 0 1 1 - 2 -86.50816673132395 0.000000000000000 0 0 0 2 - 3 -112.5453765879508 0.000000000000000 3 0 0 0 - 3 -566.4409170590384 0.000000000000000 2 1 0 0 - 3 15316.52974328784 0.000000000000000 1 2 0 0 - 3 223421.6775648115 0.000000000000000 0 3 0 0 - 3 27.39529045497964 0.000000000000000 2 0 1 0 - 3 3303.897011291280 0.000000000000000 1 1 1 0 - 3 -27213.83411547573 0.000000000000000 0 2 1 0 - 3 1142.605642821450 0.000000000000000 2 0 0 1 - 3 -25996.41018279795 0.000000000000000 1 1 0 1 - 3 -67860.88807468773 0.000000000000000 0 2 0 1 - 3 -362.5328079761680 0.000000000000000 1 0 2 0 - 3 9552.563296390306 0.000000000000000 0 1 2 0 - 3 1110.415677491335 0.000000000000000 1 0 1 1 - 3 2335.129636726653 0.000000000000000 0 1 1 1 - 3 -12551.31975790701 0.000000000000000 1 0 0 2 - 3 237306.3603909586 0.000000000000000 0 1 0 2 - 3 -30.87460160152715 0.000000000000000 0 0 3 0 - 3 28.60070020265874 0.000000000000000 0 0 2 1 - 3 -3555.652153837332 0.000000000000000 0 0 1 2 - 3 16624.27387403020 0.000000000000000 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 185 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 0.1453732038693338E-02 0.000000000000000 0 0 0 0 - 1 -0.8773365936382678E-01 0.000000000000000 1 0 0 0 - 1 -0.4028364807869175 0.000000000000000 0 1 0 0 - 1 -0.6940778729981379 0.000000000000000 0 0 1 0 - 1 -1.532887447219051 0.000000000000000 0 0 0 1 - 2 -4.275904109976432 0.000000000000000 2 0 0 0 - 2 23.14000404808969 0.000000000000000 1 1 0 0 - 2 234.0890216445948 0.000000000000000 0 2 0 0 - 2 -5.806186055780102 0.000000000000000 1 0 1 0 - 2 162.3709104370310 0.000000000000000 0 1 1 0 - 2 105.6797522604061 0.000000000000000 1 0 0 1 - 2 198.3571939390740 0.000000000000000 0 1 0 1 - 2 -1.215863077294948 0.000000000000000 0 0 2 0 - 2 56.11329075560710 0.000000000000000 0 0 1 1 - 2 -76.36813338253387 0.000000000000000 0 0 0 2 - 3 56.58630546500015 0.000000000000000 3 0 0 0 - 3 2193.606932765343 0.000000000000000 2 1 0 0 - 3 -127155.9026493125 0.000000000000000 1 2 0 0 - 3 424860.5065652152 0.000000000000000 0 3 0 0 - 3 -1766.396431186215 0.000000000000000 2 0 1 0 - 3 -978.3141867898171 0.000000000000000 1 1 1 0 - 3 -180703.5732627107 0.000000000000000 0 2 1 0 - 3 7960.231479458607 0.000000000000000 2 0 0 1 - 3 -58802.52223400104 0.000000000000000 1 1 0 1 - 3 106066.0446020791 0.000000000000000 0 2 0 1 - 3 -563.3543861407430 0.000000000000000 1 0 2 0 - 3 31606.51278796261 0.000000000000000 0 1 2 0 - 3 11604.99127689512 0.000000000000000 1 0 1 1 - 3 -41864.36532169023 0.000000000000000 0 1 1 1 - 3 -38205.85480987165 0.000000000000000 1 0 0 2 - 3 154209.7946058578 0.000000000000000 0 1 0 2 - 3 -1529.906452723396 0.000000000000000 0 0 3 0 - 3 5275.662254130144 0.000000000000000 0 0 2 1 - 3 -39497.42007533793 0.000000000000000 0 0 1 2 - 3 86177.00458696159 0.000000000000000 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 184 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 -0.5556498409026547 0.000000000000000 0 0 0 0 - 1 0.1277041819260989 0.000000000000000 1 0 0 0 - 1 2.781305250579574 0.000000000000000 0 1 0 0 - 1 -0.1556473184270137E-01 0.000000000000000 0 0 1 0 - 1 -0.2391299850942822 0.000000000000000 0 0 0 1 - 2 11.18739276352362 0.000000000000000 2 0 0 0 - 2 -21.24209034225953 0.000000000000000 1 1 0 0 - 2 737.6294948521365 0.000000000000000 0 2 0 0 - 2 -0.8543934233995063 0.000000000000000 1 0 1 0 - 2 17.97758051414892 0.000000000000000 0 1 1 0 - 2 -52.49330792491032 0.000000000000000 1 0 0 1 - 2 486.8628675218655 0.000000000000000 0 1 0 1 - 2 1.130496630307132 0.000000000000000 0 0 2 0 - 2 9.309927303735270 0.000000000000000 0 0 1 1 - 2 134.1153585557756 0.000000000000000 0 0 0 2 - 3 173.2169896106523 0.000000000000000 3 0 0 0 - 3 935.5591754634897 0.000000000000000 2 1 0 0 - 3 -23214.17450967109 0.000000000000000 1 2 0 0 - 3 -330013.1380374780 0.000000000000000 0 3 0 0 - 3 -47.17398601156798 0.000000000000000 2 0 1 0 - 3 -4971.382644178220 0.000000000000000 1 1 1 0 - 3 39061.01096856797 0.000000000000000 0 2 1 0 - 3 -1695.441120019737 0.000000000000000 2 0 0 1 - 3 37969.87762993237 0.000000000000000 1 1 0 1 - 3 108098.0101626962 0.000000000000000 0 2 0 1 - 3 549.7723738972493 0.000000000000000 1 0 2 0 - 3 -14460.98039658978 0.000000000000000 0 1 2 0 - 3 -1781.561573940901 0.000000000000000 1 0 1 1 - 3 -3784.536865453733 0.000000000000000 0 1 1 1 - 3 18705.56419054104 0.000000000000000 1 0 0 2 - 3 -354605.5490990755 0.000000000000000 0 1 0 2 - 3 40.77276943509479 0.000000000000000 0 0 3 0 - 3 -62.17111831599973 0.000000000000000 0 0 2 1 - 3 4730.253944103382 0.000000000000000 0 0 1 2 - 3 -23018.13828583136 0.000000000000000 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - - 1, NO = 3, NV = 4, INA = 183 - ********************************************* - - I COEFFICIENT ORDER EXPONENTS - NO = 3 NV = 4 - 0 -0.9014524633456889E-03 0.000000000000000 0 0 0 0 - 1 -0.1384179866484397E-01 0.000000000000000 1 0 0 0 - 1 -0.6660411325475244 0.000000000000000 0 1 0 0 - 1 -0.2901081604197308 0.000000000000000 0 0 1 0 - 1 1.901394410800662 0.000000000000000 0 0 0 1 - 2 4.020775579060298 0.000000000000000 2 0 0 0 - 2 -68.08564894216218 0.000000000000000 1 1 0 0 - 2 669.6291578676124 0.000000000000000 0 2 0 0 - 2 -6.231156375135741 0.000000000000000 1 0 1 0 - 2 148.3903477897647 0.000000000000000 0 1 1 0 - 2 54.61268127545821 0.000000000000000 1 0 0 1 - 2 347.6820854972664 0.000000000000000 0 1 0 1 - 2 -2.594675672029847 0.000000000000000 0 0 2 0 - 2 -58.41438141732203 0.000000000000000 0 0 1 1 - 2 349.5457495533942 0.000000000000000 0 0 0 2 - 3 -229.7938308753919 0.000000000000000 3 0 0 0 - 3 -7481.149604844648 0.000000000000000 2 1 0 0 - 3 38565.48391334945 0.000000000000000 1 2 0 0 - 3 719199.5801649451 0.000000000000000 0 3 0 0 - 3 3167.523459743854 0.000000000000000 2 0 1 0 - 3 8953.521511852123 0.000000000000000 1 1 1 0 - 3 214254.8614415547 0.000000000000000 0 2 1 0 - 3 -6342.058788455323 0.000000000000000 2 0 0 1 - 3 55731.83544880304 0.000000000000000 1 1 0 1 - 3 -1304801.386837105 0.000000000000000 0 2 0 1 - 3 1436.151522363004 0.000000000000000 1 0 2 0 - 3 -12152.99673179859 0.000000000000000 0 1 2 0 - 3 -10612.34182622812 0.000000000000000 1 0 1 1 - 3 -76745.90710982753 0.000000000000000 0 1 1 1 - 3 -30706.06634584456 0.000000000000000 1 0 0 2 - 3 535042.3309686909 0.000000000000000 0 1 0 2 - 3 2623.817126474333 0.000000000000000 0 0 3 0 - 3 -7949.022430507368 0.000000000000000 0 0 2 1 - 3 21048.09727979810 0.000000000000000 0 0 1 2 - 3 -188761.7996327568 0.000000000000000 0 0 0 3 - -35 0.000000000000000 0.000000000000000 0 0 0 0 - No Stochastic Radiation \ No newline at end of file diff --git a/test/spin_res/R.map b/test/spin_res/R.map new file mode 100644 index 0000000..9808620 --- /dev/null +++ b/test/spin_res/R.map @@ -0,0 +1,90 @@ + 4 Dimensional TPSA/DA map + + 1, NO = 3, NV = 4, INA = 206 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 -0.6719421772137812 -0.7406036122523244 1 0 0 0 + 3 27915.25617558173 -25327.22998076627 2 1 0 0 + 3 61805.95414685494 -56075.91795010498 1 0 1 1 + -3 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 205 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 -0.6719421772137812 0.7406036122523244 0 1 0 0 + 3 27915.25617558174 25327.22998076627 1 2 0 0 + 3 61805.95414685493 56075.91795010498 0 1 1 1 + -3 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 204 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 -0.8236855419826870 -0.5670468480898976 0 0 1 0 + 3 47322.03693360092 -68739.43091417720 1 1 1 0 + 3 16703.14333386590 -24262.78836768822 0 0 2 1 + -3 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 203 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 -0.8236855419826870 0.5670468480898976 0 0 0 1 + 3 47322.03693360092 68739.43091417720 1 1 0 1 + 3 16703.14333386590 24262.78836768822 0 0 1 2 + -3 0.000000000000000 0.000000000000000 0 0 0 0 + No Spin Matrix + Quaternion + c_quaternion + + 1, NO = 3, NV = 4, INA = 186 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 -0.2967996746715294 0.000000000000000 0 0 0 0 + 2 265.5897515029528 0.000000000000000 1 1 0 0 + 2 -270.2377358853063 0.000000000000000 0 0 1 1 + -3 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 185 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 0.1222875968392615 0.1065321708673125 0 0 1 0 + 1 0.1222875968392615 -0.1065321708673125 0 0 0 1 + 3 -19998.40360102746 -28291.04654197930 1 1 1 0 + 3 -19998.40360102732 28291.04654197941 1 1 0 1 + 3 15276.30544215209 -34693.60823067735 0 0 2 1 + 3 15276.30544215199 34693.60823067739 0 0 1 2 + -6 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 184 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 0.9549397641290652 0.000000000000000 0 0 0 0 + 2 82.54651738590167 0.000000000000000 1 1 0 0 + 2 -84.04622137341740 0.000000000000000 0 0 1 1 + -3 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 183 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 -0.1065321708673100 0.1222875968392678 0 0 1 0 + 1 -0.1065321708673100 -0.1222875968392678 0 0 0 1 + 3 28291.04654197940 -19998.40360102728 1 1 1 0 + 3 28291.04654197929 19998.40360102741 1 1 0 1 + 3 34693.60823067732 15276.30544215194 0 0 2 1 + 3 34693.60823067728 -15276.30544215204 0 0 1 2 + -6 0.000000000000000 0.000000000000000 0 0 0 0 \ No newline at end of file diff --git a/test/spin_res/m1.map b/test/spin_res/m1.map new file mode 100644 index 0000000..a1b8e15 --- /dev/null +++ b/test/spin_res/m1.map @@ -0,0 +1,165 @@ + 4 Dimensional TPSA/DA map + + 1, NO = 3, NV = 4, INA = 1469 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 1.000000000000000 0.000000000000000 1 0 0 0 + -1 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 1468 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 1.000000000000000 0.000000000000000 0 1 0 0 + -1 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 1467 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 1.000000000000000 0.000000000000000 0 0 1 0 + -1 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 1364 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 1.000000000000000 0.000000000000000 0 0 0 1 + -1 0.000000000000000 0.000000000000000 0 0 0 0 + Spin Matrix is identity + Quaternion + c_quaternion + + 1, NO = 3, NV = 4, INA = 1337 + ********************************************* + + ALL COMPONENTS 0.0_dp + NO = 3 NV = 4 + -1 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 1336 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 1.804291595517040 -0.7895290673527218 1 0 0 0 + 1 1.804291595517040 0.7895290673527218 0 1 0 0 + 1 -0.1758396404783295 -7.992402111212650 0 0 1 0 + 1 -0.1758396404783295 7.992402111212650 0 0 0 1 + 2 -34.75231970264403 60.84313848300265 2 0 0 0 + 2 59.21359339099943 0.000000000000000 1 1 0 0 + 2 -34.75231970264403 -60.84313848300265 0 2 0 0 + 2 61.13320801443875 -116.7409753541305 1 0 1 0 + 2 0.6064822360506228 9.050286557304629 0 1 1 0 + 2 0.6064822360506232 -9.050286557304624 1 0 0 1 + 2 61.13320801443871 116.7409753541305 0 1 0 1 + 2 45.87541738610921 11.49315627422870 0 0 2 0 + 2 21.91227965433630 0.000000000000000 0 0 1 1 + 2 45.87541738610921 -11.49315627422870 0 0 0 2 + 3 -6401.802964764841 9495.646355768655 3 0 0 0 + 3 -273165.4571038964 245527.6245380901 2 1 0 0 + 3 -273165.4571038962 -245527.6245380902 1 2 0 0 + 3 -6401.802964764842 -9495.646355768657 0 3 0 0 + 3 -7129.176736115559 21032.49149962682 2 0 1 0 + 3 332867.7591318095 118716.7144834473 1 1 1 0 + 3 45946.36085290540 119438.7285734195 0 2 1 0 + 3 45946.36085290539 -119438.7285734196 2 0 0 1 + 3 332867.7591318096 -118716.7144834472 1 1 0 1 + 3 -7129.176736115548 -21032.49149962682 0 2 0 1 + 3 32792.98133753835 -46052.25536964472 1 0 2 0 + 3 307783.5201201316 -4469.041559306614 0 1 2 0 + 3 5120.035559836484 138988.6842617372 1 0 1 1 + 3 5120.035559836193 -138988.6842617373 0 1 1 1 + 3 307783.5201201315 4469.041559306585 1 0 0 2 + 3 32792.98133753839 46052.25536964470 0 1 0 2 + 3 -5511.799733281073 -14023.65889657745 0 0 3 0 + 3 124315.5782203694 5001.477123418837 0 0 2 1 + 3 124315.5782203695 -5001.477123418837 0 0 1 2 + 3 -5511.799733281076 14023.65889657745 0 0 0 3 + -34 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 1335 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 0.1925253553910997E-01 -0.5085601734430411 1 0 0 0 + 1 0.1925253553910997E-01 0.5085601734430411 0 1 0 0 + 1 -0.6036548025990877E-01 -0.1625487285276035E-01 0 0 1 0 + 1 -0.6036548025990877E-01 0.1625487285276035E-01 0 0 0 1 + 2 -7.649203632630429 4.507734058446177 2 0 0 0 + 2 -7.649203632630429 -4.507734058446177 0 2 0 0 + 2 -25.66656080888161 2.422887075856805 1 0 1 0 + 2 -179.0138995422793 -52.08553089322935 0 1 1 0 + 2 -179.0138995422793 52.08553089322935 1 0 0 1 + 2 -25.66656080888161 -2.422887075856806 0 1 0 1 + 2 4.366402114609000 32.49460280436217 0 0 2 0 + 2 4.366402114609002 -32.49460280436217 0 0 0 2 + 3 2488.606802970196 759.5840312421527 3 0 0 0 + 3 9296.781587247087 4280.393527109613 2 1 0 0 + 3 9296.781587247087 -4280.393527109613 1 2 0 0 + 3 2488.606802970196 -759.5840312421520 0 3 0 0 + 3 13603.80822903459 -10566.02024204493 2 0 1 0 + 3 -9213.943060506343 -21675.90904738754 1 1 1 0 + 3 -7535.912657945539 -6080.319127977014 0 2 1 0 + 3 -7535.912657945532 6080.319127977010 2 0 0 1 + 3 -9213.943060506339 21675.90904738753 1 1 0 1 + 3 13603.80822903460 10566.02024204493 0 2 0 1 + 3 -2061.448159198135 -3736.265370486371 1 0 2 0 + 3 4945.489701386839 14163.30578130331 0 1 2 0 + 3 19776.61076055888 2867.768636873776 1 0 1 1 + 3 19776.61076055888 -2867.768636873776 0 1 1 1 + 3 4945.489701386834 -14163.30578130331 1 0 0 2 + 3 -2061.448159198135 3736.265370486373 0 1 0 2 + 3 -745.1444119828938 1239.070971228067 0 0 3 0 + 3 12759.15691456897 2847.425387468326 0 0 2 1 + 3 12759.15691456897 -2847.425387468335 0 0 1 2 + 3 -745.1444119828925 -1239.070971228068 0 0 0 3 + -32 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 1334 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 1 1.324047224444345 3.177374492276109 1 0 0 0 + 1 1.324047224444345 -3.177374492276109 0 1 0 0 + 1 -7.992402111212650 0.1758396404783295 0 0 1 0 + 1 -7.992402111212650 -0.1758396404783295 0 0 0 1 + 2 -5.336185713529815 -86.41635729788548 2 0 0 0 + 2 68.52718287435025 0.7105427357601002E-14 1 1 0 0 + 2 -5.336185713529815 86.41635729788548 0 2 0 0 + 2 15.90337298612287 42.20555585042437 1 0 1 0 + 2 3.050249160614152 -2.378236625363715 0 1 1 0 + 2 3.050249160614156 2.378236625363714 1 0 0 1 + 2 15.90337298612285 -42.20555585042435 0 1 0 1 + 2 56.61986808131009 20.31901020924470 0 0 2 0 + 2 57.74249061105725 0.1776356839400250E-14 0 0 1 1 + 2 56.61986808131009 -20.31901020924470 0 0 0 2 + 3 -3897.855136607241 -2568.794973616555 3 0 0 0 + 3 -75677.54709539931 -55124.16743197650 2 1 0 0 + 3 -75677.54709539925 55124.16743197636 1 2 0 0 + 3 -3897.855136607241 2568.794973616554 0 3 0 0 + 3 -14443.65505745811 8297.739774543916 2 0 1 0 + 3 118716.7144834473 -332867.7591318095 1 1 1 0 + 3 -176427.6297036265 78637.05922161807 0 2 1 0 + 3 -176427.6297036264 -78637.05922161805 2 0 0 1 + 3 118716.7144834472 332867.7591318096 1 1 0 1 + 3 -14443.65505745811 -8297.739774543930 0 2 0 1 + 3 14630.12285265951 2592.698601018272 1 0 2 0 + 3 147204.3906348685 -168831.7150092767 0 1 2 0 + 3 -447943.1259363448 -533098.6438489691 1 0 1 1 + 3 -447943.1259363450 533098.6438489690 0 1 1 1 + 3 147204.3906348683 168831.7150092766 1 0 0 2 + 3 14630.12285265948 -2592.698601018305 0 1 0 2 + 3 5517.260162036746 -3462.071535057766 0 0 3 0 + 3 5001.477123418837 -124315.5782203694 0 0 2 1 + 3 5001.477123418837 124315.5782203695 0 0 1 2 + 3 5517.260162036740 3462.071535057763 0 0 0 3 + -34 0.000000000000000 0.000000000000000 0 0 0 0 + No Stochastic Radiation \ No newline at end of file diff --git a/test/spin_res/test.map b/test/spin_res/test.map new file mode 100644 index 0000000..d6322e9 --- /dev/null +++ b/test/spin_res/test.map @@ -0,0 +1,341 @@ + 4 Dimensional TPSA/DA map + + 1, NO = 3, NV = 4, INA = 206 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 0.4867812994196517E-04 0.000000000000000 0 0 0 0 + 1 -0.7009790553862907 0.000000000000000 1 0 0 0 + 1 7.904434314383582 0.000000000000000 0 1 0 0 + 1 0.1066153604074146 0.000000000000000 0 0 1 0 + 1 0.2299414479930226E-01 0.000000000000000 0 0 0 1 + 2 11.79497559214498 0.000000000000000 2 0 0 0 + 2 158.3575925664820 0.000000000000000 1 1 0 0 + 2 2673.674442760778 0.000000000000000 0 2 0 0 + 2 10.48669129302925 0.000000000000000 1 0 1 0 + 2 -15.65130750712133 0.000000000000000 0 1 1 0 + 2 -99.02923431614474 0.000000000000000 1 0 0 1 + 2 -1180.000693870767 0.000000000000000 0 1 0 1 + 2 38.41036578566355 0.000000000000000 0 0 2 0 + 2 27.97031603167960 0.000000000000000 0 0 1 1 + 2 49.99458464053005 0.000000000000000 0 0 0 2 + 3 1427.339372465137 0.000000000000000 3 0 0 0 + 3 16816.85430015757 0.000000000000000 2 1 0 0 + 3 144.8299485483585 0.000000000000000 1 2 0 0 + 3 561954.5739163718 0.000000000000000 0 3 0 0 + 3 -2037.929988219854 0.000000000000000 2 0 1 0 + 3 -13901.70576711066 0.000000000000000 1 1 1 0 + 3 35475.75559255313 0.000000000000000 0 2 1 0 + 3 -7234.501404451950 0.000000000000000 2 0 0 1 + 3 132146.4025601656 0.000000000000000 1 1 0 1 + 3 553799.3263696884 0.000000000000000 0 2 0 1 + 3 6495.962490463056 0.000000000000000 1 0 2 0 + 3 27483.89047271314 0.000000000000000 0 1 2 0 + 3 7053.617749403057 0.000000000000000 1 0 1 1 + 3 214200.6547345958 0.000000000000000 0 1 1 1 + 3 125616.6591170086 0.000000000000000 1 0 0 2 + 3 1390133.288887875 0.000000000000000 0 1 0 2 + 3 47.62893353950619 0.000000000000000 0 0 3 0 + 3 6261.804641697765 0.000000000000000 0 0 2 1 + 3 -68996.22011253830 0.000000000000000 0 0 1 2 + 3 -79792.34464947744 0.000000000000000 0 0 0 3 + -35 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 205 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 0.4079952402053416E-05 0.000000000000000 0 0 0 0 + 1 -0.6950431993446177E-01 0.000000000000000 1 0 0 0 + 1 -0.6462030576832420 0.000000000000000 0 1 0 0 + 1 0.8046445421895778E-02 0.000000000000000 0 0 1 0 + 1 -0.2046721326935349E-01 0.000000000000000 0 0 0 1 + 2 -0.2785683768333875 0.000000000000000 2 0 0 0 + 2 2.099235149605022 0.000000000000000 1 1 0 0 + 2 179.8213181320455 0.000000000000000 0 2 0 0 + 2 -2.619366061310086 0.000000000000000 1 0 1 0 + 2 14.47203346047138 0.000000000000000 0 1 1 0 + 2 6.214871289366293 0.000000000000000 1 0 0 1 + 2 -94.91918593386094 0.000000000000000 0 1 0 1 + 2 1.986744420528817 0.000000000000000 0 0 2 0 + 2 -7.568461716908912 0.000000000000000 0 0 1 1 + 2 8.506437651304697 0.000000000000000 0 0 0 2 + 3 -85.07192985189931 0.000000000000000 3 0 0 0 + 3 777.0739203572730 0.000000000000000 2 1 0 0 + 3 -18800.90327233213 0.000000000000000 1 2 0 0 + 3 151245.0278252386 0.000000000000000 0 3 0 0 + 3 -17.23120377380008 0.000000000000000 2 0 1 0 + 3 4729.248840907421 0.000000000000000 1 1 1 0 + 3 -6236.397535203525 0.000000000000000 0 2 1 0 + 3 766.6717124172396 0.000000000000000 2 0 0 1 + 3 -1731.648292483020 0.000000000000000 1 1 0 1 + 3 -64564.33287938734 0.000000000000000 0 2 0 1 + 3 -1064.169138978218 0.000000000000000 1 0 2 0 + 3 7629.698042606490 0.000000000000000 0 1 2 0 + 3 136.1495662169655 0.000000000000000 1 0 1 1 + 3 6661.426994793539 0.000000000000000 0 1 1 1 + 3 -5793.228572617023 0.000000000000000 1 0 0 2 + 3 121071.5804719628 0.000000000000000 0 1 0 2 + 3 124.4945153667090 0.000000000000000 0 0 3 0 + 3 -1240.439126355940 0.000000000000000 0 0 2 1 + 3 -4608.512684806366 0.000000000000000 0 0 1 2 + 3 6940.331747395168 0.000000000000000 0 0 0 3 + -35 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 204 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 -0.1881736176817633E-02 0.000000000000000 0 0 0 0 + 1 0.1643878192699564E-01 0.000000000000000 1 0 0 0 + 1 0.1493855966153149 0.000000000000000 0 1 0 0 + 1 -0.8084495364632538 0.000000000000000 0 0 1 0 + 1 2.026592039121844 0.000000000000000 0 0 0 1 + 2 -1.951964557166611 0.000000000000000 2 0 0 0 + 2 -37.66634032326170 0.000000000000000 1 1 0 0 + 2 -533.8580278189322 0.000000000000000 0 2 0 0 + 2 11.67739061655207 0.000000000000000 1 0 1 0 + 2 128.4727431544466 0.000000000000000 0 1 1 0 + 2 8.384769157953315 0.000000000000000 1 0 0 1 + 2 62.97219966619150 0.000000000000000 0 1 0 1 + 2 -2.032881058943547 0.000000000000000 0 0 2 0 + 2 20.12523675149727 0.000000000000000 0 0 1 1 + 2 -146.5222148685612 0.000000000000000 0 0 0 2 + 3 -316.3773522568532 0.000000000000000 3 0 0 0 + 3 -123.4292062587774 0.000000000000000 2 1 0 0 + 3 81739.39826328863 0.000000000000000 1 2 0 0 + 3 93251.98243610471 0.000000000000000 0 3 0 0 + 3 2552.217051829479 0.000000000000000 2 0 1 0 + 3 -11220.32318963180 0.000000000000000 1 1 1 0 + 3 211094.2433088830 0.000000000000000 0 2 1 0 + 3 10424.11056195601 0.000000000000000 2 0 0 1 + 3 56149.20469275284 0.000000000000000 1 1 0 1 + 3 1615081.217774376 0.000000000000000 0 2 0 1 + 3 1259.920024982692 0.000000000000000 1 0 2 0 + 3 1456.852971214741 0.000000000000000 0 1 2 0 + 3 352.9416413601475 0.000000000000000 1 0 1 1 + 3 -152411.5703108387 0.000000000000000 0 1 1 1 + 3 -28971.12394811244 0.000000000000000 1 0 0 2 + 3 -105101.4208378543 0.000000000000000 0 1 0 2 + 3 1934.339423639165 0.000000000000000 0 0 3 0 + 3 19259.04091146962 0.000000000000000 0 0 2 1 + 3 26930.09344611266 0.000000000000000 0 0 1 2 + 3 131107.1894905354 0.000000000000000 0 0 0 3 + -35 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 203 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 0.2169309298279741E-03 0.000000000000000 0 0 0 0 + 1 0.1089831076000199E-02 0.000000000000000 1 0 0 0 + 1 -0.1340935961490033 0.000000000000000 0 1 0 0 + 1 -0.1612596266934681 0.000000000000000 0 0 1 0 + 1 -0.8356237888601501 0.000000000000000 0 0 0 1 + 2 -2.109642153719256 0.000000000000000 2 0 0 0 + 2 6.702969933420714 0.000000000000000 1 1 0 0 + 2 -177.9506367930843 0.000000000000000 0 2 0 0 + 2 -1.213411706449580 0.000000000000000 1 0 1 0 + 2 -71.23633107352856 0.000000000000000 0 1 1 0 + 2 -5.741433603522974 0.000000000000000 1 0 0 1 + 2 57.55400546403282 0.000000000000000 0 1 0 1 + 2 -1.765419941567711 0.000000000000000 0 0 2 0 + 2 13.73339712575856 0.000000000000000 0 0 1 1 + 2 -56.46272202393999 0.000000000000000 0 0 0 2 + 3 -26.12920442163664 0.000000000000000 3 0 0 0 + 3 2056.933264341425 0.000000000000000 2 1 0 0 + 3 1602.776374712306 0.000000000000000 1 2 0 0 + 3 -49011.23128974235 0.000000000000000 0 3 0 0 + 3 -1037.826133533944 0.000000000000000 2 0 1 0 + 3 6998.503672639559 0.000000000000000 1 1 1 0 + 3 -53169.98408589326 0.000000000000000 0 2 1 0 + 3 1957.465253320000 0.000000000000000 2 0 0 1 + 3 1277.912528791556 0.000000000000000 1 1 0 1 + 3 260159.9398063142 0.000000000000000 0 2 0 1 + 3 151.7248934854840 0.000000000000000 1 0 2 0 + 3 -1066.836334195409 0.000000000000000 0 1 2 0 + 3 -2251.648465699132 0.000000000000000 1 0 1 1 + 3 -30775.92467424635 0.000000000000000 0 1 1 1 + 3 -6506.364083884773 0.000000000000000 1 0 0 2 + 3 65774.41023640917 0.000000000000000 0 1 0 2 + 3 -481.4346121063184 0.000000000000000 0 0 3 0 + 3 386.0481962047043 0.000000000000000 0 0 2 1 + 3 -15340.11089393405 0.000000000000000 0 0 1 2 + 3 33257.14667513099 0.000000000000000 0 0 0 3 + -35 0.000000000000000 0.000000000000000 0 0 0 0 + No Spin Matrix + Quaternion + c_quaternion + + 1, NO = 3, NV = 4, INA = 186 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 -0.2967996746715293 0.000000000000000 0 0 0 0 + 1 0.1561012093764574 0.000000000000000 1 0 0 0 + 1 3.722883276874034 0.000000000000000 0 1 0 0 + 1 0.7641205223954218E-03 0.000000000000000 0 0 1 0 + 1 0.4864449934371465E-01 0.000000000000000 0 0 0 1 + 2 13.21880030461130 0.000000000000000 2 0 0 0 + 2 -14.94621012856408 0.000000000000000 1 1 0 0 + 2 1007.362437331714 0.000000000000000 0 2 0 0 + 2 3.828320578099249 0.000000000000000 1 0 1 0 + 2 -51.71278156687806 0.000000000000000 0 1 1 0 + 2 -16.65630677712585 0.000000000000000 1 0 0 1 + 2 212.6376876499664 0.000000000000000 0 1 0 1 + 2 1.338473254453042 0.000000000000000 0 0 2 0 + 2 -10.32682115008832 0.000000000000000 0 0 1 1 + 2 -85.94626785873777 0.000000000000000 0 0 0 2 + 3 199.6521334388494 0.000000000000000 3 0 0 0 + 3 -515.7965012541285 0.000000000000000 2 1 0 0 + 3 -16162.99820050433 0.000000000000000 1 2 0 0 + 3 -335079.2856413586 0.000000000000000 0 3 0 0 + 3 109.8694266611851 0.000000000000000 2 0 1 0 + 3 -8124.952059859799 0.000000000000000 1 1 1 0 + 3 -5281.774484877198 0.000000000000000 0 2 1 0 + 3 -1588.876188464621 0.000000000000000 2 0 0 1 + 3 7835.633751390963 0.000000000000000 1 1 0 1 + 3 117859.4620509410 0.000000000000000 0 2 0 1 + 3 1568.816943522881 0.000000000000000 1 0 2 0 + 3 -21238.99041101434 0.000000000000000 0 1 2 0 + 3 -948.1986771149096 0.000000000000000 1 0 1 1 + 3 9601.635012672821 0.000000000000000 0 1 1 1 + 3 14406.69961249195 0.000000000000000 1 0 0 2 + 3 -231052.7544266776 0.000000000000000 0 1 0 2 + 3 -294.9104355931647 0.000000000000000 0 0 3 0 + 3 2497.396182369676 0.000000000000000 0 0 2 1 + 3 11638.41008922243 0.000000000000000 0 0 1 2 + 3 -22477.70617443228 0.000000000000000 0 0 0 3 + -35 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 185 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 -0.9378109818871268E-02 0.000000000000000 0 0 0 0 + 1 -0.1258080225944594 0.000000000000000 1 0 0 0 + 1 1.454695886838768 0.000000000000000 0 1 0 0 + 1 2.165374026944721 0.000000000000000 0 0 1 0 + 1 23.21846595999425 0.000000000000000 0 0 0 1 + 2 22.03083501149388 0.000000000000000 2 0 0 0 + 2 -138.2128837365694 0.000000000000000 1 1 0 0 + 2 2008.454246825898 0.000000000000000 0 2 0 0 + 2 -45.92789927553395 0.000000000000000 1 0 1 0 + 2 1368.199884666971 0.000000000000000 0 1 1 0 + 2 200.6658755188116 0.000000000000000 1 0 0 1 + 2 1059.706922082436 0.000000000000000 0 1 0 1 + 2 24.11612789441240 0.000000000000000 0 0 2 0 + 2 -81.70273938238094 0.000000000000000 0 0 1 1 + 2 489.6805613774040 0.000000000000000 0 0 0 2 + 3 -259.5390573040842 0.000000000000000 3 0 0 0 + 3 -23587.01078240259 0.000000000000000 2 1 0 0 + 3 -58189.18272880526 0.000000000000000 1 2 0 0 + 3 748180.4065702861 0.000000000000000 0 3 0 0 + 3 14522.12177141163 0.000000000000000 2 0 1 0 + 3 -90083.95892733682 0.000000000000000 1 1 1 0 + 3 1165562.509352523 0.000000000000000 0 2 1 0 + 3 -10110.14911353981 0.000000000000000 2 0 0 1 + 3 51786.86696377716 0.000000000000000 1 1 0 1 + 3 -1710650.753515610 0.000000000000000 0 2 0 1 + 3 -4832.914175220283 0.000000000000000 1 0 2 0 + 3 41999.25789822237 0.000000000000000 0 1 2 0 + 3 34757.31247922957 0.000000000000000 1 0 1 1 + 3 141403.7331520373 0.000000000000000 0 1 1 1 + 3 17790.35390684681 0.000000000000000 1 0 0 2 + 3 -825601.0919344509 0.000000000000000 0 1 0 2 + 3 11156.88470526157 0.000000000000000 0 0 3 0 + 3 24760.27019748771 0.000000000000000 0 0 2 1 + 3 186100.3386721462 0.000000000000000 0 0 1 2 + 3 -240750.2675011277 0.000000000000000 0 0 0 3 + -35 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 184 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 0.9548199704114351 0.000000000000000 0 0 0 0 + 1 0.4791277590665655E-01 0.000000000000000 1 0 0 0 + 1 1.192876545001144 0.000000000000000 0 1 0 0 + 1 -0.5979764743139625E-01 0.000000000000000 0 0 1 0 + 1 0.3222169398816387 0.000000000000000 0 0 0 1 + 2 4.271212982714482 0.000000000000000 2 0 0 0 + 2 -8.171019169182360 0.000000000000000 1 1 0 0 + 2 305.3082523842008 0.000000000000000 0 2 0 0 + 2 2.219590679479423 0.000000000000000 1 0 1 0 + 2 16.78458143258052 0.000000000000000 0 1 1 0 + 2 1.114140659593143 0.000000000000000 1 0 0 1 + 2 19.36559943479044 0.000000000000000 0 1 0 1 + 2 -24.15383806588803 0.000000000000000 0 0 2 0 + 2 -12.50225827341920 0.000000000000000 0 0 1 1 + 2 -331.5542469702108 0.000000000000000 0 0 0 2 + 3 47.99157589903309 0.000000000000000 3 0 0 0 + 3 -568.1116104102456 0.000000000000000 2 1 0 0 + 3 -1059.339287743586 0.000000000000000 1 2 0 0 + 3 -96031.59597441691 0.000000000000000 0 3 0 0 + 3 249.8709511362710 0.000000000000000 2 0 1 0 + 3 -4344.988696793234 0.000000000000000 1 1 1 0 + 3 4918.572818855383 0.000000000000000 0 2 1 0 + 3 -668.8585155362409 0.000000000000000 2 0 0 1 + 3 9356.922274684650 0.000000000000000 1 1 0 1 + 3 44950.09109969577 0.000000000000000 0 2 0 1 + 3 1075.572438533631 0.000000000000000 1 0 2 0 + 3 -3318.733755449175 0.000000000000000 0 1 2 0 + 3 1187.825534275374 0.000000000000000 1 0 1 1 + 3 -48722.62663894740 0.000000000000000 0 1 1 1 + 3 -1902.800290632470 0.000000000000000 1 0 0 2 + 3 -103539.1437517130 0.000000000000000 0 1 0 2 + 3 112.5499786383581 0.000000000000000 0 0 3 0 + 3 1716.372965986014 0.000000000000000 0 0 2 1 + 3 4001.423802867926 0.000000000000000 0 0 1 2 + 3 -12469.50320389639 0.000000000000000 0 0 0 3 + -35 0.000000000000000 0.000000000000000 0 0 0 0 + + 1, NO = 3, NV = 4, INA = 183 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 3 NV = 4 + 0 -0.1186710894051485E-01 0.000000000000000 0 0 0 0 + 1 0.5031794767687248E-01 0.000000000000000 1 0 0 0 + 1 1.718152610595122 0.000000000000000 0 1 0 0 + 1 -6.541601200960002 0.000000000000000 0 0 1 0 + 1 6.360114683690653 0.000000000000000 0 0 0 1 + 2 -2.459940195208099 0.000000000000000 2 0 0 0 + 2 -128.7511789586227 0.000000000000000 1 1 0 0 + 2 -1359.255258515615 0.000000000000000 0 2 0 0 + 2 68.20995150526528 0.000000000000000 1 0 1 0 + 2 875.1471484527590 0.000000000000000 0 1 1 0 + 2 130.4036202523218 0.000000000000000 1 0 0 1 + 2 -782.7741343302763 0.000000000000000 0 1 0 1 + 2 4.760601453542805 0.000000000000000 0 0 2 0 + 2 46.00169979835894 0.000000000000000 0 0 1 1 + 2 -491.3467018616159 0.000000000000000 0 0 0 2 + 3 -979.7420749397230 0.000000000000000 3 0 0 0 + 3 -6559.359203140742 0.000000000000000 2 1 0 0 + 3 311326.1226868047 0.000000000000000 1 2 0 0 + 3 458644.9266322721 0.000000000000000 0 3 0 0 + 3 12071.04317801135 0.000000000000000 2 0 1 0 + 3 -35141.43208086607 0.000000000000000 1 1 1 0 + 3 1000912.122522561 0.000000000000000 0 2 1 0 + 3 34078.90718138096 0.000000000000000 2 0 0 1 + 3 202936.6062126463 0.000000000000000 1 1 0 1 + 3 5319594.090665693 0.000000000000000 0 2 0 1 + 3 4815.941688151748 0.000000000000000 1 0 2 0 + 3 -231.9273177052596 0.000000000000000 0 1 2 0 + 3 4192.027018295088 0.000000000000000 1 0 1 1 + 3 -508942.7259145931 0.000000000000000 0 1 1 1 + 3 -74754.69691117053 0.000000000000000 1 0 0 2 + 3 -315730.3314054444 0.000000000000000 0 1 0 2 + 3 9512.649240258916 0.000000000000000 0 0 3 0 + 3 64952.13203115657 0.000000000000000 0 0 2 1 + 3 110087.7128129228 0.000000000000000 0 0 1 2 + 3 434530.8143226430 0.000000000000000 0 0 0 3 + -35 0.000000000000000 0.000000000000000 0 0 0 0 + No Stochastic Radiation