From afc4b2395a043b017278122f0d159b08d3a8cee4 Mon Sep 17 00:00:00 2001 From: mattsignorelli Date: Fri, 19 Jul 2024 16:55:11 -0400 Subject: [PATCH] working on factorize --- fpp-ptc-sandbox/code/CMakeLists.txt | 1 + fpp-ptc-sandbox/code/Makefile | 41 + .../code/cmake_files/cmake.z_benggston | 14 + fpp-ptc-sandbox/code/map.txt | 239 +++++ fpp-ptc-sandbox/code/plot.dat | 100 +++ fpp-ptc-sandbox/code/z_benggston.f90 | 578 ++++++++++++ fpp-ptc-sandbox/code/z_spin_res.f90 | 2 +- src/NonlinearNormalForm.jl | 3 +- src/normal.jl | 80 +- test/factorize/a.map | 203 +++++ test/factorize/a4.map | 820 ++++++++++++++++++ 11 files changed, 2076 insertions(+), 5 deletions(-) create mode 100644 fpp-ptc-sandbox/code/cmake_files/cmake.z_benggston create mode 100644 fpp-ptc-sandbox/code/map.txt create mode 100644 fpp-ptc-sandbox/code/z_benggston.f90 create mode 100644 test/factorize/a.map create mode 100644 test/factorize/a4.map diff --git a/fpp-ptc-sandbox/code/CMakeLists.txt b/fpp-ptc-sandbox/code/CMakeLists.txt index 53cc249..e3db20c 100644 --- a/fpp-ptc-sandbox/code/CMakeLists.txt +++ b/fpp-ptc-sandbox/code/CMakeLists.txt @@ -34,6 +34,7 @@ set (EXE_SPECS cmake_files/cmake.z_resonance cmake_files/cmake.z_spin1 cmake_files/cmake.z_spin_res + cmake_files/cmake.z_benggston ) diff --git a/fpp-ptc-sandbox/code/Makefile b/fpp-ptc-sandbox/code/Makefile index e9e94b1..e1d56af 100644 --- a/fpp-ptc-sandbox/code/Makefile +++ b/fpp-ptc-sandbox/code/Makefile @@ -441,6 +441,19 @@ 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 +#============================================================================= +# Target rules for targets named z_benggston-exe + +# Build rule for target. +z_benggston-exe: cmake_check_build_system + $(MAKE) $(MAKESILENT) -f CMakeFiles/Makefile2 z_benggston-exe +.PHONY : z_benggston-exe + +# fast build rule for target. +z_benggston-exe/fast: + $(MAKE) $(MAKESILENT) -f CMakeFiles/z_benggston-exe.dir/build.make CMakeFiles/z_benggston-exe.dir/build +.PHONY : z_benggston-exe/fast + normal_simple.o: normal_simple.f90.o .PHONY : normal_simple.o @@ -465,6 +478,30 @@ normal_simple.f90.s: $(MAKE) $(MAKESILENT) -f CMakeFiles/normal_simple-exe.dir/build.make CMakeFiles/normal_simple-exe.dir/normal_simple.f90.s .PHONY : normal_simple.f90.s +z_benggston.o: z_benggston.f90.o +.PHONY : z_benggston.o + +# target to build an object file +z_benggston.f90.o: + $(MAKE) $(MAKESILENT) -f CMakeFiles/z_benggston-exe.dir/build.make CMakeFiles/z_benggston-exe.dir/z_benggston.f90.o +.PHONY : z_benggston.f90.o + +z_benggston.i: z_benggston.f90.i +.PHONY : z_benggston.i + +# target to preprocess a source file +z_benggston.f90.i: + $(MAKE) $(MAKESILENT) -f CMakeFiles/z_benggston-exe.dir/build.make CMakeFiles/z_benggston-exe.dir/z_benggston.f90.i +.PHONY : z_benggston.f90.i + +z_benggston.s: z_benggston.f90.s +.PHONY : z_benggston.s + +# target to generate assembly for a file +z_benggston.f90.s: + $(MAKE) $(MAKESILENT) -f CMakeFiles/z_benggston-exe.dir/build.make CMakeFiles/z_benggston-exe.dir/z_benggston.f90.s +.PHONY : z_benggston.f90.s + z_canonize_matt_fake_maps.o: z_canonize_matt_fake_maps.f90.o .PHONY : z_canonize_matt_fake_maps.o @@ -1050,6 +1087,7 @@ help: @echo "... edit_cache" @echo "... rebuild_cache" @echo "... normal_simple-exe" + @echo "... z_benggston-exe" @echo "... z_canonize_matt_fake_maps-exe" @echo "... z_coast_matt_fake_maps-exe" @echo "... z_cut-exe" @@ -1077,6 +1115,9 @@ help: @echo "... normal_simple.o" @echo "... normal_simple.i" @echo "... normal_simple.s" + @echo "... z_benggston.o" + @echo "... z_benggston.i" + @echo "... z_benggston.s" @echo "... z_canonize_matt_fake_maps.o" @echo "... z_canonize_matt_fake_maps.i" @echo "... z_canonize_matt_fake_maps.s" diff --git a/fpp-ptc-sandbox/code/cmake_files/cmake.z_benggston b/fpp-ptc-sandbox/code/cmake_files/cmake.z_benggston new file mode 100644 index 0000000..43d2bd4 --- /dev/null +++ b/fpp-ptc-sandbox/code/cmake_files/cmake.z_benggston @@ -0,0 +1,14 @@ +set (EXENAME z_benggston) + +FILE (GLOB SRC_FILES z_benggston.f90) + +set (INC_DIRS +) + +set (LINK_LIBS + forest + bsim + bmad + sim_utils + ${ACC_BMAD_LINK_LIBS} +) diff --git a/fpp-ptc-sandbox/code/map.txt b/fpp-ptc-sandbox/code/map.txt new file mode 100644 index 0000000..11db5f8 --- /dev/null +++ b/fpp-ptc-sandbox/code/map.txt @@ -0,0 +1,239 @@ + 6 Dimensional TPSA/DA map + + 1, NO = 2, NV = 8, INA = 526 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 3.248739248164878 0.000000000000000 1 0 0 0 0 0 0 0 + 1 -0.1666181892190699E-04 0.000000000000000 0 0 0 0 0 1 0 0 + 2 9.481060110485970 0.8881784197001250E-15 2 0 0 0 0 0 0 0 + 2 -1.415180496021328 -0.1232595164407831E-31 1 1 0 0 0 0 0 0 + 2 210.3141529350647 -0.8881784197001252E-15 0 2 0 0 0 0 0 0 + 2 102.1100983383984 -0.1998401444325282E-14 0 0 2 0 0 0 0 0 + 2 -3.415504135280438 0.000000000000000 0 0 1 1 0 0 0 0 + 2 -8.850217497580280 -0.1998401444325282E-14 0 0 0 2 0 0 0 0 + 2 4.752917252847052 0.1213973082665156E-35 1 0 0 0 0 1 0 0 + 2 0.2113546110192892 0.6321614590947053E-37 0 1 0 0 0 1 0 0 + 2 0.6261305504023522 0.000000000000000 1 0 0 0 0 0 1 0 + 2 0.5772675543170591 0.000000000000000 0 1 0 0 0 0 1 0 + 2 -0.6261492747210906 0.000000000000000 1 0 0 0 0 0 0 1 + 2 -0.5772697454811313 0.000000000000000 0 1 0 0 0 0 0 1 + 2 1.066755907013854 0.000000000000000 0 0 0 0 0 2 0 0 + 2 -0.1057805123175021 0.000000000000000 0 0 0 0 0 1 1 0 + 2 0.1057781125343333 0.000000000000000 0 0 0 0 0 1 0 1 + -17 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 527 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 0.1262958087429899E-01 0.000000000000000 1 0 0 0 0 0 0 0 + 1 0.3078117151337622 0.000000000000000 0 1 0 0 0 0 0 0 + 1 0.2176579745711674E-09 0.000000000000000 0 0 0 0 0 1 0 0 + 2 -1.177767725300218 -0.5081713830180124E-33 2 0 0 0 0 0 0 0 + 2 -1.802125513034857 0.000000000000000 1 1 0 0 0 0 0 0 + 2 0.8846459970520613 0.000000000000000 0 2 0 0 0 0 0 0 + 2 0.5771097021168961 -0.4440892098500626E-15 0 0 2 0 0 0 0 0 + 2 -7.959512334345833 -0.3146056920670404E-33 0 0 1 1 0 0 0 0 + 2 -0.4013131322035950 -0.4440892098500626E-15 0 0 0 2 0 0 0 0 + 2 0.3850256010038856E-01 0.000000000000000 1 0 0 0 0 1 0 0 + 2 -0.4495080028843736 0.000000000000000 0 1 0 0 0 1 0 0 + 2 0.5712907938500170E-01 0.000000000000000 1 0 0 0 0 0 1 0 + 2 -0.5708050329320462E-01 0.000000000000000 0 1 0 0 0 0 1 0 + 2 -0.5712935978493515E-01 0.000000000000000 1 0 0 0 0 0 0 1 + 2 0.5708226886770709E-01 0.000000000000000 0 1 0 0 0 0 0 1 + 2 0.3349377514389761E-02 0.000000000000000 0 0 0 0 0 2 0 0 + 2 0.3128260423381617E-01 0.000000000000000 0 0 0 0 0 1 1 0 + 2 -0.3128261177046555E-01 0.000000000000000 0 0 0 0 0 1 0 1 + -18 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 528 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 1.998246462619670 0.000000000000000 0 0 1 0 0 0 0 0 + 2 51.58521946741487 -0.1110223024625157E-15 1 0 1 0 0 0 0 0 + 2 -2.100820821566178 0.3552713678800501E-14 0 1 1 0 0 0 0 0 + 2 4.763768354838248 0.5329070518200751E-14 1 0 0 1 0 0 0 0 + 2 -10.88724853368484 0.3219646771412954E-14 0 1 0 1 0 0 0 0 + 2 -6.943564711368939 0.3295974604355933E-15 0 0 1 0 0 1 0 0 + 2 0.2366158864055927E-01 -0.4440892098500626E-15 0 0 0 1 0 1 0 0 + 2 0.4448179007408046E-01 -0.5551115123125783E-16 0 0 1 0 0 0 1 0 + 2 0.5063300944771596 -0.2775557561562891E-16 0 0 0 1 0 0 1 0 + 2 -0.9502844475768888 0.2220446049250313E-15 0 0 1 0 0 0 0 1 + 2 -1.142662996648280 0.000000000000000 0 0 0 1 0 0 0 1 + -11 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 529 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 -0.1593841720719695E-02 0.000000000000000 0 0 1 0 0 0 0 0 + 1 0.5004387690440424 0.000000000000000 0 0 0 1 0 0 0 0 + 2 0.5446393799924837 0.1776356839400250E-14 1 0 1 0 0 0 0 0 + 2 -31.45659046770351 0.8053453737222327E-15 0 1 1 0 0 0 0 0 + 2 -12.92274847319755 0.2818925648462312E-16 1 0 0 1 0 0 0 0 + 2 0.5348112742120492 0.000000000000000 0 1 0 1 0 0 0 0 + 2 0.1146411108386130E-01 -0.2220446049250313E-15 0 0 1 0 0 1 0 0 + 2 1.738920262895184 -0.8239936510889834E-16 0 0 0 1 0 1 0 0 + 2 0.1267693034948138 -0.5204170427930421E-17 0 0 1 0 0 0 1 0 + 2 -0.1154383242012146E-01 0.1387778780781446E-16 0 0 0 1 0 0 1 0 + 2 -0.2854093682358202 -0.2405301686640682E-34 0 0 1 0 0 0 0 1 + 2 0.2388996613150604 -0.5551115123125783E-16 0 0 0 1 0 0 0 1 + -12 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 530 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 -0.2111389035918150E-06 0.000000000000000 1 0 0 0 0 0 0 0 + 1 -0.5128703059600362E-05 0.000000000000000 0 1 0 0 0 0 0 0 + 1 1.000000000000000 0.000000000000000 0 0 0 0 1 0 0 0 + 1 0.2205173373131304E-27 0.000000000000000 0 0 0 0 0 1 0 0 + 2 -0.3250909097068860E-01 0.000000000000000 2 0 0 0 0 0 0 0 + 2 1.463033638484699 0.000000000000000 1 1 0 0 0 0 0 0 + 2 0.3251392707166587E-01 0.1408687431750624E-37 0 2 0 0 0 0 0 0 + 2 -0.5930226068886798E-02 0.1528459407851808E-35 0 0 2 0 0 0 0 0 + 2 -3.474696356238558 0.000000000000000 0 0 1 1 0 0 0 0 + 2 0.5927276679515037E-02 0.1533723295750092E-35 0 0 0 2 0 0 0 0 + 2 0.5182209068500229E-02 0.000000000000000 1 0 0 0 0 1 0 0 + 2 0.6567274203089519 0.000000000000000 0 1 0 0 0 1 0 0 + 2 -0.1029659397051052 0.000000000000000 1 0 0 0 0 0 1 0 + 2 -0.3255952998481526E-01 0.000000000000000 0 1 0 0 0 0 1 0 + 2 0.1029659338861337 0.000000000000000 1 0 0 0 0 0 0 1 + 2 0.3255879127402354E-01 0.000000000000000 0 1 0 0 0 0 0 1 + 2 -0.5603890957595162E-07 0.000000000000000 0 0 0 0 0 2 0 0 + 2 -0.5212020631774660E-06 0.000000000000000 0 0 0 0 0 1 1 0 + 2 0.5212021892740849E-06 0.000000000000000 0 0 0 0 0 1 0 1 + -19 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 531 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 1.000000000000000 0.000000000000000 0 0 0 0 0 1 0 0 + -1 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + No Spin Matrix + Quaternion + c_quaternion + + 1, NO = 2, NV = 8, INA = 541 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 0 1.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + 2 -0.3656272367183466E-03 0.000000000000000 2 0 0 0 0 0 0 0 + 2 -0.1970644790041756E-01 0.000000000000000 1 1 0 0 0 0 0 0 + 2 -0.2655327953255095 0.000000000000000 0 2 0 0 0 0 0 0 + 2 -17.12576810110231 0.000000000000000 0 0 2 0 0 0 0 0 + 2 -0.2883300980794109E-01 0.000000000000000 0 0 1 1 0 0 0 0 + 2 -6.967620058230136 0.000000000000000 0 0 0 2 0 0 0 0 + 2 0.5169878828456423E-25 0.000000000000000 1 0 0 0 0 1 0 0 + 2 0.9926167350636332E-23 0.000000000000000 0 1 0 0 0 1 0 0 + 2 -0.2465190328815662E-31 0.000000000000000 0 0 0 0 0 2 0 0 + -10 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 542 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 0.6116118172205857E-02 0.000000000000000 0 0 1 0 0 0 0 0 + 1 3.732993311654496 0.000000000000000 0 0 0 1 0 0 0 0 + 2 -0.2443278299575462 0.9947598300641406E-13 1 0 1 0 0 0 0 0 + 2 109.2436427525633 -0.4632405570248467E-13 0 1 1 0 0 0 0 0 + 2 85.72996966047718 -0.5279110482092619E-13 1 0 0 1 0 0 0 0 + 2 -0.7825032607712537 -0.8526512829121201E-13 0 1 0 1 0 0 0 0 + 2 0.1235187017839498E-01 0.7105427357601001E-14 0 0 1 0 0 1 0 0 + 2 61.25133614217027 -0.5023759186428834E-14 0 0 0 1 0 1 0 0 + 2 10.53475422545662 -0.7222237291452134E-33 0 0 1 0 0 0 1 0 + 2 -9.273963311626169 -0.1232595164407831E-31 0 0 0 1 0 0 1 0 + 2 2.790124084553993 -0.8881784197001252E-15 0 0 1 0 0 0 0 1 + 2 -12.49460045688714 0.4440892098500626E-15 0 0 0 1 0 0 0 1 + -12 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 543 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 0.2704171727972763E-01 0.000000000000000 1 0 0 0 0 0 0 0 + 1 0.7287424721059004 0.000000000000000 0 1 0 0 0 0 0 0 + 1 -0.6617444900424221E-23 0.000000000000000 0 0 0 0 0 1 0 0 + 2 -2.011826271918972 0.2220446049250313E-15 2 0 0 0 0 0 0 0 + 2 -5.850441748540714 -0.3081487911019577E-32 1 1 0 0 0 0 0 0 + 2 2.011826271918973 -0.2220446049250314E-15 0 2 0 0 0 0 0 0 + 2 1.082111486715603 0.1783817379350225E-30 0 0 2 0 0 0 0 0 + 2 -23.82368803773213 -0.2220446049250313E-15 0 0 1 1 0 0 0 0 + 2 -1.082111486715608 0.000000000000000 0 0 0 2 0 0 0 0 + 2 0.5617869815516087E-01 0.3078325831215865E-36 1 0 0 0 0 1 0 0 + 2 -2.857948051239356 -0.1175494350822288E-37 0 1 0 0 0 1 0 0 + 2 0.3237737723614194 0.000000000000000 1 0 0 0 0 0 1 0 + 2 -0.8022595317801831E-01 0.000000000000000 0 1 0 0 0 0 1 0 + 2 -0.3237743806053234 0.000000000000000 1 0 0 0 0 0 0 1 + 2 0.8023148227749416E-01 0.000000000000000 0 1 0 0 0 0 0 1 + 2 -0.1387778780781446E-16 0.000000000000000 0 0 0 0 0 1 1 0 + -16 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 544 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 -5.852478004683410 0.000000000000000 0 0 1 0 0 0 0 0 + 1 -0.1025476998462083E-02 0.000000000000000 0 0 0 1 0 0 0 0 + 2 -35.76993929294610 0.6273453978522525E-13 1 0 1 0 0 0 0 0 + 2 2.424757139212243 0.4618527782440651E-13 0 1 1 0 0 0 0 0 + 2 -2.547762920537700 0.3907985046680551E-13 1 0 0 1 0 0 0 0 + 2 34.55658049676584 -0.6895178872312613E-13 0 1 0 1 0 0 0 0 + 2 -62.01949143455906 -0.5023759186428833E-14 0 0 1 0 0 1 0 0 + 2 0.2821363443417674E-01 -0.5916456789157589E-30 0 0 0 1 0 1 0 0 + 2 14.70856647275444 0.2315930758125651E-31 0 0 1 0 0 0 1 0 + 2 12.85892134338535 0.1232595164407831E-31 0 0 0 1 0 0 1 0 + 2 10.79584271946479 -0.1332267629550188E-14 0 0 1 0 0 0 0 1 + 2 3.488645572261075 0.8881784197001252E-15 0 0 0 1 0 0 0 1 + -12 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + Stochastic Radiation + 1 1 (3.50448659610772345E-308,1.75224329805385135E-308) + 1 2 (5.43818303559283309E-310,-7.19911681148759229E-310) + 1 3 (0.0000000000000000,0.0000000000000000) + 1 4 (0.0000000000000000,0.0000000000000000) + 1 5 (-1.68163050528502313E-314,2.60131090882960446E-314) + 1 6 (0.0000000000000000,0.0000000000000000) + 2 1 (5.43818303559298131E-310,7.15078692341526084E-311) + 2 2 (3.66523886152904226E-312,-2.91278374568177670E-311) + 2 3 (0.0000000000000000,0.0000000000000000) + 2 4 (0.0000000000000000,0.0000000000000000) + 2 5 (-9.64765593313416768E-317,6.09257426659024015E-316) + 2 6 (0.0000000000000000,0.0000000000000000) + 3 1 (0.0000000000000000,0.0000000000000000) + 3 2 (0.0000000000000000,0.0000000000000000) + 3 3 (0.0000000000000000,0.0000000000000000) + 3 4 (0.0000000000000000,0.0000000000000000) + 3 5 (0.0000000000000000,0.0000000000000000) + 3 6 (0.0000000000000000,0.0000000000000000) + 4 1 (0.0000000000000000,0.0000000000000000) + 4 2 (0.0000000000000000,0.0000000000000000) + 4 3 (0.0000000000000000,0.0000000000000000) + 4 4 (0.0000000000000000,0.0000000000000000) + 4 5 (0.0000000000000000,0.0000000000000000) + 4 6 (0.0000000000000000,0.0000000000000000) + 5 1 (-1.68163050429689184E-314,-1.45320184530463523E-315) + 5 2 (-9.64765593313416768E-317,5.53934080119998345E-316) + 5 3 (0.0000000000000000,0.0000000000000000) + 5 4 (0.0000000000000000,0.0000000000000000) + 5 5 (2.46538757274782026E-321,-1.50541802287827822E-320) + 5 6 (0.0000000000000000,0.0000000000000000) + 6 1 (0.0000000000000000,0.0000000000000000) + 6 2 (0.0000000000000000,0.0000000000000000) + 6 3 (0.0000000000000000,0.0000000000000000) + 6 4 (0.0000000000000000,0.0000000000000000) + 6 5 (0.0000000000000000,0.0000000000000000) + 6 6 (0.0000000000000000,0.0000000000000000) diff --git a/fpp-ptc-sandbox/code/plot.dat b/fpp-ptc-sandbox/code/plot.dat index e69de29..9d5ee14 100644 --- a/fpp-ptc-sandbox/code/plot.dat +++ b/fpp-ptc-sandbox/code/plot.dat @@ -0,0 +1,100 @@ + 47060.00000000000 0.9999861940741539E-02 2.130070243797641 0.9999998690539408 + 94120.00000000000 0.1999972388148308E-01 2.130270241034320 0.9999999930995638 + 141180.0000000000 0.2999958582222462E-01 2.130470238270999 0.9999999169788700 + 188240.0000000000 0.3999944776296616E-01 2.130670235507678 0.9999999142579082 + 235300.0000000000 0.4999930784106255E-01 2.130870232744357 0.9999999035884019 + 282360.0000000000 0.5999917164444923E-01 2.131070229981036 0.9999999505918226 + 329420.0000000000 0.6999903172254562E-01 2.131270227217716 0.9999999944594228 + 376480.0000000000 0.7999889552593231E-01 2.131470224454394 0.9999999498937586 + 423540.0000000000 0.8999875932931900E-01 2.131670221691074 0.9999999033561825 + 470600.0000000000 0.9999861568212509E-01 2.131870218927753 0.9999999228595640 + 517660.0000000000 0.1099984794855118 2.132070216164432 0.9999998497377904 + 564720.0000000000 0.1199983432888985 2.132270213401111 0.9999999483401635 + 611780.0000000000 0.1299982070922852 2.132470210637790 0.9999997643276962 + 658840.0000000000 0.1399980634450912 2.132670207874469 0.9999999746095959 + 705900.0000000000 0.1499979346990585 2.132870205111148 0.9999999204392338 + 752960.0000000000 0.1599977910518646 2.133070202347827 0.9999997957737058 + 800020.0000000000 0.1699976474046707 2.133270199584506 0.9999997229615297 + 847080.0000000000 0.1799975186586380 2.133470196821185 0.9999999675330148 + 894140.0000000000 0.1899973750114441 2.133670194057864 0.9999997648121702 + 941200.0000000000 0.1999972313642502 2.133870191294543 0.9999999081185219 + 988260.0000000000 0.2099971026182175 2.134070188531222 0.9999996964607345 + 1035320.000000000 0.2199969589710236 2.134270185767901 0.9999996992325101 + 1082380.000000000 0.2299968302249908 2.134470183004580 0.9999998116521314 + 1129440.000000000 0.2399966865777969 2.134670180241260 0.9999996094105336 + 1176500.000000000 0.2499965429306030 2.134870177477938 0.9999997975458368 + 1223560.000000000 0.2599964141845703 2.135070174714618 0.9999995736173370 + 1270620.000000000 0.2699962854385376 2.135270171951297 0.9999995454034117 + 1317680.000000000 0.2799961268901825 2.135470169187976 0.9999998910999252 + 1364740.000000000 0.2899959981441498 2.135670166424655 0.9999997174336765 + 1411800.000000000 0.2999958693981171 2.135870163661334 0.9999998898396736 + 1458860.000000000 0.3099957108497620 2.136070160898013 0.9999995739922685 + 1505920.000000000 0.3199955821037292 2.136270158134692 0.9999996855034164 + 1552980.000000000 0.3299954533576965 2.136470155371371 0.9999994118005752 + 1600040.000000000 0.3399952948093414 2.136670152608050 0.9999995038582928 + 1647100.000000000 0.3499951660633087 2.136870149844729 0.9999992004892030 + 1694160.000000000 0.3599950373172760 2.137070147081408 0.9999992458411114 + 1741220.000000000 0.3699948787689209 2.137270144318087 0.9999989394190700 + 1788280.000000000 0.3799947500228882 2.137470141554766 0.9999995068156391 + 1835340.000000000 0.3899946212768555 2.137670138791445 0.9999993874899621 + 1882400.000000000 0.3999944627285004 2.137870136028124 0.9999983793464552 + 1929460.000000000 0.4099943339824677 2.138070133264804 0.9999986098020546 + 1976520.000000000 0.4199942052364349 2.138270130501482 0.9999976662725276 + 2023580.000000000 0.4299940466880798 2.138470127738162 0.9999982298173526 + 2070640.000000000 0.4399939179420471 2.138670124974841 0.9999973556949662 + 2117700.000000000 0.4499937891960144 2.138870122211520 0.9999944256489829 + 2164760.000000000 0.4599936604499817 2.139070119448199 0.9999915227225489 + 2211820.000000000 0.4699935019016266 2.139270116684878 0.9999860618536174 + 2258880.000000000 0.4799933731555939 2.139470113921557 0.9999669504881820 + 2305940.000000000 0.4899932444095612 2.139670111158236 0.9998517639761982 + 2353000.000000000 0.4999930858612061 2.139870108394915 0.9710040950330100 + 2400060.000000000 0.5099929571151733 2.140070105631594 0.9624724141723591 + 2447120.000000000 0.5199928283691406 2.140270102868273 0.9581016474960029 + 2494180.000000000 0.5299926996231079 2.140470100104952 0.9603139570341057 + 2541240.000000000 0.5399925708770752 2.140670097341631 0.9593282856102132 + 2588300.000000000 0.5499923825263977 2.140870094578310 0.9607426863181403 + 2635360.000000000 0.5599922537803650 2.141070091814989 0.9605947772725700 + 2682420.000000000 0.5699921250343323 2.141270089051668 0.9597700140499462 + 2729480.000000000 0.5799919962882996 2.141470086288348 0.9603651501590875 + 2776540.000000000 0.5899918675422668 2.141670083525026 0.9597629351411721 + 2823600.000000000 0.5999917387962341 2.141870080761706 0.9603948689476888 + 2870660.000000000 0.6099915504455566 2.142070077998385 0.9603651244862944 + 2917720.000000000 0.6199914216995239 2.142270075235063 0.9596818172326196 + 2964780.000000000 0.6299912929534912 2.142470072471743 0.9598600406794281 + 3011840.000000000 0.6399911642074585 2.142670069708422 0.9597771325428934 + 3058900.000000000 0.6499910354614258 2.142870066945101 0.9599254232970938 + 3105960.000000000 0.6599909067153931 2.143070064181780 0.9598052047827081 + 3153020.000000000 0.6699907183647156 2.143270061418459 0.9602005501827096 + 3200080.000000000 0.6799905896186829 2.143470058655138 0.9601122959620750 + 3247140.000000000 0.6899904608726501 2.143670055891817 0.9601119161790992 + 3294200.000000000 0.6999903321266174 2.143870053128496 0.9601804697338535 + 3341260.000000000 0.7099902033805847 2.144070050365175 0.9601783763609725 + 3388320.000000000 0.7199900746345520 2.144270047601854 0.9599016996949961 + 3435380.000000000 0.7299899458885193 2.144470044838533 0.9598963578900317 + 3482440.000000000 0.7399897575378418 2.144670042075212 0.9601630531676477 + 3529500.000000000 0.7499896287918091 2.144870039311892 0.9600841707600341 + 3576560.000000000 0.7599895000457764 2.145070036548570 0.9601492837508474 + 3623620.000000000 0.7699893712997437 2.145270033785249 0.9599452543487355 + 3670680.000000000 0.7799892425537109 2.145470031021929 0.9599887687957550 + 3717740.000000000 0.7899891138076782 2.145670028258607 0.9600568263833108 + 3764800.000000000 0.7999889254570007 2.145870025495287 0.9601275066471210 + 3811860.000000000 0.8099887967109680 2.146070022731966 0.9599437330676168 + 3858920.000000000 0.8199886679649353 2.146270019968645 0.9601225618879408 + 3905980.000000000 0.8299885392189026 2.146470017205324 0.9600926446066826 + 3953040.000000000 0.8399884104728699 2.146670014442003 0.9599271797147144 + 4000100.000000000 0.8499882817268372 2.146870011678682 0.9598987074175923 + 4047160.000000000 0.8599880933761597 2.147070008915361 0.9601013685475019 + 4094220.000000000 0.8699879646301270 2.147270006152040 0.9599517958313661 + 4141280.000000000 0.8799878358840942 2.147470003388719 0.9600487620365100 + 4188340.000000000 0.8899877071380615 2.147670000625398 0.9599350946111253 + 4235400.000000000 0.8999875783920288 2.147869997862077 0.9599430542916887 + 4282460.000000000 0.9099874496459961 2.148069995098756 0.9600152755067196 + 4329520.000000000 0.9199873208999634 2.148269992335435 0.9599208810808922 + 4376580.000000000 0.9299871325492859 2.148469989572114 0.9600410069532400 + 4423640.000000000 0.9399870038032532 2.148669986808793 0.9599234804187090 + 4470700.000000000 0.9499868750572205 2.148869984045473 0.9599255666037199 + 4517760.000000000 0.9599867463111877 2.149069981282151 0.9600780733215439 + 4564820.000000000 0.9699866175651550 2.149269978518831 0.9600293241783109 + 4611880.000000000 0.9799864888191223 2.149469975755510 0.9600792609476374 + 4658940.000000000 0.9899863004684448 2.149669972992189 0.9600207286188650 + 4706000.000000000 0.9999861717224121 2.149869970228868 0.9600693657612053 diff --git a/fpp-ptc-sandbox/code/z_benggston.f90 b/fpp-ptc-sandbox/code/z_benggston.f90 new file mode 100644 index 0000000..57f5117 --- /dev/null +++ b/fpp-ptc-sandbox/code/z_benggston.f90 @@ -0,0 +1,578 @@ +program Guignard_normal_form_average_x + use madx_ptc_module + use pointer_lattice + use c_TPSA + implicit none + + + type(layout), pointer:: ALS + real(dp) prec + real(dp) circ,dtheta,s,total_tunes(4),ts,xx,ph(3),damp(3),misa(6) + type(internal_state),target :: state + logical(lp) :: mis=.false. + type(c_taylor) fonction,fonction_FLOQUET,phase(3) + type(c_damap) one_turn_map, id_s,m01,m12,m23,m3f,a0,a1,a2,as,one_turn_map_linear_floquet + type(c_damap) n01,n12,n23,n3f + type(c_normal_form) normal_form + type(c_vector_field) f_non + integer :: pos =1 + integer i,map_order,mf,mf1,n_s,n_mode,km ,n,i1,i2,i3,np + type(probe) ray_closed + type(probe_8) ray,ray_cs_twiss + type(fibre), pointer :: p + type(integration_node), pointer :: t + character*120 :: dsc + logical :: doit=.true. ,skew=.false. + type(c_vector_field_fourier) H,F,K,f1 + real(dp) DX_AVERAGE_DCS,betax_1,theta + complex(dp) coe, sol(2) + !!!!!!!!!!!!!!!!!!!!! + integer,target :: mybeg, M_res(-1:ndim) + real(dp) , target :: mondelta,closed_orbit(6) + integer selcas + type(c_universal_taylor) h_res + integer , allocatable :: JJ(:) + + use_quaternion=.true. + c_verbose=.false. + prec=1.d-6 ! for printing + use_info = .true. + longprint=.false. + !c_lda_used=c_lda_used + !lda_used=lda_used + n_cai=-i_ + mis=.false. + !state=nocavity0 +time0+spin0 + state=nocavity0 +time0+spin0 + !state=default0 +spin0 +radiation0 + state=nocavity0 +spin0 + skew=.false. + remove_tune_shift=.false. + radfac=100 + call in_bmad_units + + + call ptc_ini_no_append + call append_empty_layout(m_u) + + !call ptc_ini_no_append ! initializes PTC + check_excessive_cutting=.false. + switch_to_drift_kick=.false. + !call read_lattice_append(M_U,'C:\document\my_tex_papers\fpp_handbook\new_pages\tables\g18zd066.txt') + ALS=>m_u%start + call build_lattice(ALS,mis,exact=.false.,thin=.false.,onecell=.false.) + + + closed_orbit=.01d0 + ray_closed=closed_orbit + + p=>als%start + do i=1,200 + if(i==128) then + + write(6,*) i, p%mag%name(1:10) + + write(6,*) p%MAG%usef + write(6,*) p%MAG%useb + write(6,*) p%MAGp%usef + write(6,*) p%MAGp%useb + + endif + + p=>p%next + enddo + + call propagate(als,ray_closed,state,fibre1=128,fibre2=129) + write(6,format4) ray_closed%x(1:4) + + + !goto 1001 + + i1=45 + i2=180 + i3=270 + + np=2 + + map_order=2 !9 + if(np==2) then + p=>als%start + do i=1,als%n + + if(i>=i1.and.i=3) then + if(p%mag%bn(3)/=0) then + write(6,*) p%mag%name + !p%magp%bn(3)=0.d0 + ! call make_it_knob(p%magp%bn(3),1,p%mag%bn(3)) + call make_it_knob(p%magp%bn(2),1) + endif + endif + + endif + + if(i>=i2.and.i=3) then + if(p%mag%bn(3)/=0) then + write(6,*) p%mag%name + ! p%magp%bn(3)=0.d0 + ! call make_it_knob(p%magp%bn(3),2,p%mag%bn(3)) + call make_it_knob(p%magp%bn(2),2) + endif + endif + + endif + + p=>p%next + enddo + endif + + + call init_all(state,map_order,np) + call alloc(one_turn_map, id_s,m01,m12,m23,m3f,a0,a1,a2,one_turn_map_linear_floquet ) + + call alloc(n01,n12,n23,n3f,as ) + call alloc(normal_form); call alloc(ray); + call alloc(f_non); call alloc(fonction,fonction_FLOQUET); + call alloc(phase); + + + + allocate(JJ(1:c_%nd2)) + closed_orbit=0.d0; + call find_orbit_x(als,closed_orbit(1:6),STATE,1.e-8_dp,fibre1=1) + + ray_closed=closed_orbit + id_s=1; + ray=id_s+ray_closed; + call propagate(als,RAY,+state,fibre1=1) + one_turn_map=ray + + call print(one_turn_map) + + + + call c_normal(one_turn_map,normal_form,dospin=state%spin,phase=phase) ! (6) + + !# a_t = ;5 a_s for dir=1 + call c_full_factorise(normal_form%atot,as,a0,a1,a2,dir=1) + m01=as* a0*a1*a2 + m01=m01**(-1)*normal_form%atot + + + call kanalnummer(i,"map.txt") + call print(normal_form%atot,i) + call clean(a0,a0,prec=1.d-10) + call clean(a1,a1,prec=1.d-10) + call clean(a2,a2,prec=1.d-10) + call clean(as,as,prec=1.d-10) + close(i) + pause + call print(a0) + + stop + fonction=normal_form%atot%v(5).par.'1000' + + call print(fonction) + pause 1 + call print(a1) + pause 1 + call print(a2) + pause 1 + call print(as) + pause 1 + fonction=1.d0.cmono.'2' + call average(fonction,a1,fonction_FLOQUET) + + s=(a1%v(1).sub.'1')**2+(a1%v(1).sub.'01')**2 + + write(6,*) s + call print(fonction_FLOQUET) + + stop + + + + call c_fast_canonise(normal_form%A_t ,a0,dospin=state%spin) ! (9a) + + call print(a0) + pause 100 + one_turn_map_linear_floquet=a0**(-1)*one_turn_map*a0 + one_turn_map_linear_floquet=one_turn_map_linear_floquet*(one_turn_map_linear_floquet.cut.2)**(-1) + one_turn_map_linear_floquet=ci_phasor()*one_turn_map_linear_floquet*c_phasor() + f_non=ln(one_turn_map_linear_floquet) + !call c_q0_to_qr(f_non%q,f_non%q) + call print(f_non) + call d_field_for_demin(f_non,H_res) + + call clean(H_res,H_res,prec=1.e-5_dp) + call print(H_res) + jj=0 + jj(2)=1 + jj(3)=2 + fonction=H_res.par.jj + call print(fonction) + stop + pause 666 + jj=0 + jj(2)=5 + fonction= f_non%v(1).par.jj + fonction= fonction /6/i_ + call print(fonction) + call d_field_for_demin(f_non,H_res) + + call clean(H_res,H_res,prec=1.e-5_dp) + + jj=0 + jj(1)=6 + fonction=H_res.par.jj + call print(fonction) + + jj=0 + jj(2)=6 + fonction=H_res.par.jj + call print(fonction) + + + fonction=cgetpb(f_non) + jj=0 + jj(1)=6 + fonction=fonction.par.jj + call print(fonction) + + + + + goto 1001 + !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! + + call propagate(als,RAY,+state,fibre1=1,fibre2=i1) ! (4) + m01=ray + ray_closed=ray + ray=id_s+ray_closed; + call propagate(als,RAY,+state,fibre1=i1,fibre2=i2) ! (4) + m12=ray + ray_closed=ray + ray=id_s+ray_closed; + call propagate(als,RAY,+state,fibre1=i2,fibre2=i3) ! (4) + m23=ray + ray_closed=ray + ray=id_s+ray_closed; + call propagate(als,RAY,+state,fibre1=i3,fibre2=als%n+1) ! (4) + m3f=ray + ray_closed=ray + + one_turn_map=m3f*m23*m12*m01 + + call c_normal(one_turn_map,normal_form,dospin=state%spin,phase=phase) ! (6) + + call c_fast_canonise(normal_form%A_t ,a0) ! (9a) + + one_turn_map_linear_floquet=a0**(-1)*one_turn_map*a0 + one_turn_map_linear_floquet=(one_turn_map_linear_floquet.cut.2)**(-1)*one_turn_map_linear_floquet + one_turn_map_linear_floquet=ci_phasor()*one_turn_map_linear_floquet*c_phasor() + f_non=ln(one_turn_map_linear_floquet) + + jj=0 + jj(2)=5 + fonction= f_non%v(1).par.jj + fonction= fonction /6/i_ + call print(fonction) + call d_field_for_demin(f_non,H_res) + + call clean(H_res,H_res,prec=1.e-5_dp) + + + + fonction=cgetpb(f_non) + jj=0 + jj(1)=6 + fonction=fonction.par.jj + call print(fonction) + + + fonction=H_res.par.jj + call print(fonction) + + + + + stop + + 1001 call ptc_end(graphics_maybe=1,flat_file=.false.) + + contains + + subroutine build_lattice(ALS,mis,error,exact,sl,thin,onecell) + use madx_ptc_module + use pointer_lattice + implicit none + + type(layout), target :: ALS + real(dp),optional :: error(6) + logical, optional :: exact,sl,thin,onecell + real(dp) :: alpha,lbend, cut, ksd, ksf,sig(6) + type(fibre) L1,L2,L3,L4,L5,L6,L7,L8,L9,L10 + type(fibre) L11,L12,L13,L14,L15,L16,L17,L18,L19,L20,CAVM + type(fibre) L21,L22,L23,L24,L25,L26,L27,L27A,L27B,L27C,L27D,DS + type(fibre) QF1,QF2,QD1,QD2,QFA1,QFA2,sf,sd,cav,bend,vc5,bend1 + type(fibre) L27h + type(layout) :: sfline,sdline,sup1,supb + logical(lp) :: mis,thi=.false.,oneperiod + !----------------------------------- + if(present(thin)) thi=thin + + call make_states(.true.) + exact_model = .false.;oneperiod = .false. + if(present(exact)) exact_model=exact + if(present(onecell)) oneperiod=onecell + call update_states + madlength = .false. + + !old_integrator=-100 + call set_mad(energy = 1.5d0, method = 2, step = 1) + + madkind2 = matrix_kick_matrix + !madkind2 = drift_kick_drift + + + L1 = drift("L1 ", 2.832695d0);L2 = drift("L2 ", 0.45698d0); + L3 = drift("L3 ", 0.08902d0);L4 = drift("L4 ", 0.2155d0); + L5 = drift("L5 ", 0.219d0);L6 = drift("L6 ", 0.107078d0); + L7 = drift("L7 ", 0.105716d0);L8 = drift("L8 ", 0.135904d0); + L9 = drift("L9 ", 0.2156993d0);L10 = drift("L10", 0.089084d0); + L11= drift("L11", 0.235416d0);L12= drift("L12", 0.1245d0); + L13= drift("L13", 0.511844d0);L14= drift("L14", 0.1788541d0); + L15= drift("L15", 0.1788483d0);L16= drift("L16", 0.511849d0); + L17= drift("L17", 0.1245d0);L18= drift("L18", 0.235405d0); + L19= drift("L19", 0.089095d0);L20= drift("L20", 0.2157007d0); + L21= drift("L21", 0.177716d0);L22= drift("L22", 0.170981d0); + L23= drift("L23", 0.218997d0);L24 = drift ("L24", 0.215503d0); + L25 = drift ("L25", 0.0890187d0);L26 = drift ("L26", 0.45698d0); + L27 = drift ("L27", 2.832696d0);L27a = drift (" L27a", 0.8596d0); + L27b = drift (" L27b", 0.1524d0);L27c = drift (" L27c", 0.04445d0); + L27d = drift (" L27d", 1.776246d0);ds = drift (" DS ", 0.1015d0); + L27h = drift ("L27", 2.832696d0-0.1d0); + ksd=6.447435260914397D-03 + if(.not.skew)ksd=6.447435260914397D-02-2.5d-2 + !if(skew) ksd=6.447435260914397D-02-2.1d-2 + QF1 = QUADRUPOLE(" QF1 ",0.344D0, K1= 2.2474D0+ksd) + QF2 = QUADRUPOLE(" QF2 ",0.344D0, K1= 2.2474D0) + QD1 = QUADRUPOLE(" QD1 ",0.187D0, K1= -2.3368D0-2.593018157427161D-02); + QD2 = QUADRUPOLE(" QD2 ",0.187D0, K1= -2.3368D0); + QFA1= QUADRUPOLE(" QFA1",0.448D0, K1= 2.8856D0); + QFA2= QUADRUPOLE(" QFA2",0.448D0, K1= 2.8856D0); + + !!! 1/2 mad-x value + ksf=-41.3355516397069748d0; + ksd=56.2564709584745489d0; + + sf=sextupole ("sf",2.d0*0.1015d0, K2= ksf); + sd= sextupole("sd", 2.d0*0.1015d0, K2= ksd); + if(skew) call add(sf,-3,1,300.d0) + VC5=marker("vc5"); + ALPHA=0.17453292519943295769236907684886d0; + + LBEND=0.86621d0; + + BEND = RBEND("BEND", LBEND, ANGLE=ALPHA).q.(-0.778741d0) + BEND1 = RBEND("BEND1", LBEND, ANGLE=ALPHA).q.(-0.778741d0) + + CAVM=MARK("CAVM"); + !CAV=RFCAVITY("CAV",L=0.1000d0,VOLT=-100.0d0,REV_FREQ=500.0d6) + CAV=RFCAVITY("CAV",L=0.0000d0,VOLT=-1000.0d0,REV_FREQ=500.0d6) + + if(thi) then + sf=sextupole ("sf",0.d0, K2= ksf*0.203d0); + sd= sextupole("sd", 0.d0, K2= ksd*0.203d0); + sfline=(ds+sf+ds); + sdline=(ds+sd+ds); + else + sfline=1*sf; + sdline=1*sd; + endif + + if(oneperiod) then + SUP1=L1+L2+L3+QF1+VC5+L4+L5+QD1+L6+L7+L8+VC5+BEND+VC5+L9+sfline+L10+& + L11+QFA1+L12+sdline+L13+ & + L14+BEND+L15+L16+sdline+L17+ & + QFA2+L18+L19+sfline+L20+BEND+L21+& + L22+QD2+L23+L24+QF2+L25+ & + L26+VC5+L27+cav; + else + SUP1=L1+L2+L3+QF1+VC5+L4+L5+QD1+L6+L7+L8+VC5+BEND+VC5+L9+sfline+L10+& + L11+QFA1+L12+sdline+L13+ & + L14+BEND+L15+L16+sdline+L17+ & + QFA2+L18+L19+sfline+L20+BEND+L21+& + L22+QD2+L23+L24+QF2+L25+ & + L26+VC5+L27+cavm; + endif + ! L26+VC5+L27+cavm; + + SUPb=L1+L2+L3+QF1+VC5+L4+L5+QD1+L6+L7+L8+VC5+BEND+VC5+L9+sfline+L10+& + L11+QFA1+L12+sdline+L13+ & + L14+BEND+L15+L16+sdline+L17+ & + QFA2+L18+L19+sfline+L20+BEND1+L21+& + L22+QD2+L23+L24+QF2+L25+ & + L26+VC5+L27+cav; + + if(oneperiod) then + ALS = sup1; !11*sup1+supb; + else + ALS = 11*sup1+supb; + endif + if(present(sl)) then + L1 = drift("L1 ", 2.832695d0); + if( sl ) then + Qf1 = QUADRUPOLE(" QF1 ",L=0.d0, K1= 0.01d0 ); L1 = drift("L1 ",L=0.1d0); + ALS=L1+QF1; + endif + endif + + ALS = .ring.ALS + + call survey(ALS) + + + if(mis) then + sig=1.d-5; cut=4.d0; + if(present(error)) sig=error + call MESS_UP_ALIGNMENT(ALS,SIG,cut); + endif + end subroutine build_lattice + + subroutine phase_advance_n(mf) + implicit none + type(fibre),pointer:: f + type(probe) xs0,xs1 + type(probe_8) xs + type(c_damap) id + type(c_normal_form) n + type(integration_node), pointer :: t + real(dp) phase(3), spin_tune(2),damping(3) + integer i,mff + integer,optional :: mf + type(c_taylor) nu_spin + use_quaternion=.true. + mff=0 + if(present(mf)) mff=mf + + f=>my_ering%start + do i=1,my_start-1 + f=>f%next + enddo + + write(6,*) f%mag%name + my_fix=0 + + my_fix(ndpt_bmad+5)=my_delta + + call find_orbit_x(my_fix,my_estate,1.e-8_dp,fibre1=f) + + if(.not.check_stable) then + write(6,*) "Could not find closed orbit " + write(6,*) " No calculation done " + return + endif + + call init(my_estate,1,0) + + call alloc(id) + call alloc(xs) + call alloc(n) + call alloc(nu_spin) + xs0=my_fix + id=1 + xs=xs0+id + call propagate(xs,my_estate,fibre1=f) + id=xs + + call c_normal(id,n,dospin=my_estate%spin,nu_spin=nu_spin) + + if(mff/=0) then + write(mff,*) " Linear A from c_normal " + call print(n%atot,mff) + endif + + + + + call c_fast_canonise(n%atot,n%atot, dospin=my_estate%spin) + + if(mff/=0) then + write(mff,*) " Linear A canonised " + call print(n%atot,mff) + write(mff,*) " end of Info from phase_advance " + endif + + !if(mff/=0) then + write(mff,*) " Info from map :tunes, damping, spin,quaternion_angle/pi" + write(mff,*) n%tune(1:c_%nd) + write(mff,*) n%damping(1:c_%nd) + write(mff,*) n%spin_tune,n%quaternion_angle/pi + write(mff,*) " end of Info from map " + !endif + phase=0 + spin_tune=0 + damping=0 + t=>f%t1 + + xs=xs0+n%atot + + do i=1,my_ering%t%n + + + call propagate(xs,my_estate,node1=t,node2=t%next) + xs0=xs + n%atot=xs + + f%tm%lf%symplectic=.not.my_estate%radiation + + call c_fast_canonise(n%atot,n%atot,phase=phase,damping=damping,spin_tune=spin_tune,dospin=my_estate%spin) + + xs=xs0+n%atot + t=>t%next + + call compute_lattice_functions(n%atot,t%lf) + t%lf%phase=phase + t%lf%damping=damping + t%lf%spin=spin_tune + t%lf%fix=xs0%x + + if(i==my_ering%t%n) then + t%lf%phase=0 + t%lf%damping=0 + t%lf%spin=0 + endif + + + enddo + + write(6,*) " " + write(6,*) " Phase advance and fractional" + write(6,format3) phase(1:c_%nd) + write(6,format3) n%tune(1:c_%nd) + + write(6,*) " damping advance " + write(6,format3) damping + write(6,*) " spin advance and chromaticity " + write(6,format2) spin_tune + write(6,format1) n%spin_tune + write(6,*) " Closed orbit before and after " + write(6,format6) my_fix + write(6,format6) xs0%x + write(6,*) " " + + write(6,*) "spin from normal form" + call print(nu_spin) + + call kill(id) + call kill(xs) + call kill(n) + call kill(nu_spin) + + end subroutine phase_advance_n + + + + end program Guignard_normal_form_average_x \ No newline at end of file diff --git a/fpp-ptc-sandbox/code/z_spin_res.f90 b/fpp-ptc-sandbox/code/z_spin_res.f90 index 7a5aa9f..9f9c915 100644 --- a/fpp-ptc-sandbox/code/z_spin_res.f90 +++ b/fpp-ptc-sandbox/code/z_spin_res.f90 @@ -164,7 +164,7 @@ program spin_res U_0=u_0*exp(fq) Nc=Ci_phasor()*N*U_0**(-1)*c_phasor() -fq=c_logf_spin(Nc) +fq=ln(Nc) call c_q0_to_qr(fq%q,qf) diff --git a/src/NonlinearNormalForm.jl b/src/NonlinearNormalForm.jl index ab242ec..ab7a6a8 100644 --- a/src/NonlinearNormalForm.jl +++ b/src/NonlinearNormalForm.jl @@ -109,7 +109,8 @@ export TaylorMap, inv_with_log, - equilibrium_moments + equilibrium_moments, + factorize diff --git a/src/normal.jl b/src/normal.jl index 5ed61e4..ef3ae06 100644 --- a/src/normal.jl +++ b/src/normal.jl @@ -148,13 +148,13 @@ function normal(m::DAMap; res=nothing, spin_res=nothing) nrm = sqrt(n0[1]^2+n0[3]^2) Qr = Quaternion(cos(alpha/2), sin(alpha/2)*n0[3]/nrm, 0, -sin(alpha/2)*n0[1]/nrm) # now concatenate m1: - as = DAMap(Q=Qr) + as = DAMap(Q=Qr,idpt=m.idpt) m1 = inv(as)*m1*as # == Quaternion(scalar.(Qr*m.Q*inv(Qr))) nu0 = 2*acos(scalar(m1.Q.q0)) # closed orbit spin tune ( i guess we could have gotten this earlier?) # it is equal to 2*acos(Quaternion(scalar.(m.Q)).q0) (i.e. before transforming) # Now we start killing the spin. The first step is to start with a map # (identity in orbital because we are done with orbital) that does this zero order rotation - Qr_inv = DAMap(Q=inv(Quaternion(scalar.(m1.Q)))) + Qr_inv = DAMap(Q=inv(Quaternion(scalar.(m1.Q))),idpt=m.idpt) # Now store analogous to eg -> egspin egspin = SVector(cos(nu0)+im*sin(nu0), 1, cos(nu0)-im*sin(nu0)) @@ -209,7 +209,7 @@ function normal(m::DAMap; res=nothing, spin_res=nothing) # Exponentiate this part now - Qnr = DAMap(Q=exp(Quaternion(0,na...))) + Qnr = DAMap(Q=exp(Quaternion(0,na...)),idpt=m.idpt) as = as*Qnr # put in normalizing map m1 = inv(Qnr)*m1*Qnr # kill the terms in m1 end @@ -227,6 +227,80 @@ function normal(m::DAMap; res=nothing, spin_res=nothing) end +# to get dbeta/ddelta, first go to fully nonlinear parameter dependent fixed point +# then calculate lattice functions. Lattice functions will be TPSA and then you +# can extract dbeta/ddelta +function factorize(a) + # We get a0*a1*an + if !isnothing(a.idpt) # if coasting, set number of variables executing pseudo-harmonic oscillations + nhv = numvars(a)-2 + ndpt = numvars(a)-1+a.idpt # energy like variable index + sgn = 1-2*a.idpt + nt = ndpt+sgn # timelike variable index + zer = zero(a); zer.x[nt][nt]=1; zer.x[ndpt][ndpt] = 1 + + a0 = a∘zer + + + nhv = numvars(a)-2 + + # ensure poisson bracket does not change + for i=1:Int(nhv/2) + a0.x[nt] += sgn*deriv(a0.x[2*i], ndpt)*mono(2*i-1,use=getdesc(a)) - sgn*deriv(a0.x[2*i-1], ndpt)*mono(2*i,use=getdesc(a)) # first order + #return a0.x[2*i-1] + # dies at second order poisson bracket + #print(sgn*0.5*(deriv(a0.x[2*i], ndpt)*a0.x[2*i-1] - deriv(a0.x[2*i-1], ndpt)*a0.x[2*i])) + #return + a0.x[nt] += sgn*0.5*(deriv(a0.x[2*i], ndpt)*a0.x[2*i-1] - deriv(a0.x[2*i-1], ndpt)*a0.x[2*i]) + end + + eye = zero(a) + setmatrix!(eye, I(nhv)) + if !isnothing(eye.Q) + eye.Q.q0[0] = 1 + end + a0 += eye + else + nhv = numvars(a) + a0 = a∘zero(a)+I + end + + a1 = inv(a0)*a + for i=1:nhv + tmp = zero(a1.x[i]) + ords = zeros(Int, nhv) + for j=1:nhv + ords[j] += 1 + tmp += a1.x[i][(ords..., :)] + ords[j] -= 1 + end + a1.x[i] = tmp + end + + if !isnothing(a.Q) + a1.Q.q0 = 1 + a1.Q.q1 = 0 + a1.Q.q2 = 0 + a1.Q.q3 = 0 + end + + a2 = inv(a1)*inv(a0)*a + + if !isnothing(a.Q) + a2.Q.q0 = 1 + a2.Q.q1 = 0 + a2.Q.q2 = 0 + a2.Q.q3 = 0 + + as = inv(a2)*inv(a1)*inv(a0)*a + return a0, a1, a2, as + end + + return a0, a1, a2 + +end + + """ diff --git a/test/factorize/a.map b/test/factorize/a.map new file mode 100644 index 0000000..cbd2789 --- /dev/null +++ b/test/factorize/a.map @@ -0,0 +1,203 @@ + 6 Dimensional TPSA/DA map + + 1, NO = 2, NV = 8, INA = 526 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 -3.248739248164878 0.000000000000000 1 0 0 0 0 0 0 0 + 1 -0.1666181892190698E-04 0.000000000000000 0 0 0 0 0 1 0 0 + 2 9.481060110485997 0.000000000000000 2 0 0 0 0 0 0 0 + 2 -1.415180496021588 0.000000000000000 1 1 0 0 0 0 0 0 + 2 210.3141529350647 0.000000000000000 0 2 0 0 0 0 0 0 + 2 102.1100983383977 -0.9992007221626407E-15 0 0 2 0 0 0 0 0 + 2 -3.415504135281418 0.4880507905132675E-33 0 0 1 1 0 0 0 0 + 2 -8.850217497579548 -0.9992007221626409E-15 0 0 0 2 0 0 0 0 + 2 -4.752917252847059 0.000000000000000 1 0 0 0 0 1 0 0 + 2 -0.2113546110192288 0.000000000000000 0 1 0 0 0 1 0 0 + 2 -0.6261305504023529 0.000000000000000 1 0 0 0 0 0 1 0 + 2 -0.5772675543170588 0.000000000000000 0 1 0 0 0 0 1 0 + 2 0.6261492747210899 0.000000000000000 1 0 0 0 0 0 0 1 + 2 0.5772697454811320 0.000000000000000 0 1 0 0 0 0 0 1 + 2 1.066755907013854 0.000000000000000 0 0 0 0 0 2 0 0 + 2 -0.1057805123175021 0.000000000000000 0 0 0 0 0 1 1 0 + 2 0.1057781125343333 0.000000000000000 0 0 0 0 0 1 0 1 + -17 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 527 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 -0.1262958087429899E-01 0.000000000000000 1 0 0 0 0 0 0 0 + 1 -0.3078117151337621 0.000000000000000 0 1 0 0 0 0 0 0 + 1 0.2176579745712350E-09 0.000000000000000 0 0 0 0 0 1 0 0 + 2 -1.177767725300229 0.000000000000000 2 0 0 0 0 0 0 0 + 2 -1.802125513034837 0.000000000000000 1 1 0 0 0 0 0 0 + 2 0.8846459970520728 0.000000000000000 0 2 0 0 0 0 0 0 + 2 0.5771097021169346 0.000000000000000 0 0 2 0 0 0 0 0 + 2 -7.959512334345894 -0.3009265538105055E-35 0 0 1 1 0 0 0 0 + 2 -0.4013131322036372 0.000000000000000 0 0 0 2 0 0 0 0 + 2 -0.3850256010038286E-01 0.000000000000000 1 0 0 0 0 1 0 0 + 2 0.4495080028843745 0.000000000000000 0 1 0 0 0 1 0 0 + 2 -0.5712907938500165E-01 0.000000000000000 1 0 0 0 0 0 1 0 + 2 0.5708050329320468E-01 0.000000000000000 0 1 0 0 0 0 1 0 + 2 0.5712935978493520E-01 0.000000000000000 1 0 0 0 0 0 0 1 + 2 -0.5708226886770704E-01 0.000000000000000 0 1 0 0 0 0 0 1 + 2 0.3349377514389771E-02 0.000000000000000 0 0 0 0 0 2 0 0 + 2 0.3128260423381616E-01 0.000000000000000 0 0 0 0 0 1 1 0 + 2 -0.3128261177046555E-01 0.000000000000000 0 0 0 0 0 1 0 1 + -18 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 528 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 1.998246462619670 0.000000000000000 0 0 1 0 0 0 0 0 + 2 -51.58521946741487 -0.1665334536937735E-14 1 0 1 0 0 0 0 0 + 2 2.100820821566119 0.1776356839400250E-14 0 1 1 0 0 0 0 0 + 2 -4.763768354838218 -0.1776356839400250E-14 1 0 0 1 0 0 0 0 + 2 10.88724853368485 -0.1887379141862766E-14 0 1 0 1 0 0 0 0 + 2 -6.943564711368943 0.000000000000000 0 0 1 0 0 1 0 0 + 2 0.2366158864055318E-01 0.000000000000000 0 0 0 1 0 1 0 0 + 2 0.4448179007408107E-01 0.000000000000000 0 0 1 0 0 0 1 0 + 2 0.5063300944771626 0.000000000000000 0 0 0 1 0 0 1 0 + 2 -0.9502844475768885 0.000000000000000 0 0 1 0 0 0 0 1 + 2 -1.142662996648278 0.000000000000000 0 0 0 1 0 0 0 1 + -11 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 529 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 -0.1593841720719695E-02 0.000000000000000 0 0 1 0 0 0 0 0 + 1 0.5004387690440424 0.000000000000000 0 0 0 1 0 0 0 0 + 2 -0.5446393799924619 0.1156344375593035E-33 1 0 1 0 0 0 0 0 + 2 31.45659046770351 -0.3335005882565411E-15 0 1 1 0 0 0 0 0 + 2 12.92274847319755 0.3326332265185527E-15 1 0 0 1 0 0 0 0 + 2 -0.5348112742120309 0.2373376099491244E-33 0 1 0 1 0 0 0 0 + 2 0.1146411108385977E-01 0.000000000000000 0 0 1 0 0 1 0 0 + 2 1.738920262895185 -0.2055531645902959E-36 0 0 0 1 0 1 0 0 + 2 0.1267693034948146 0.000000000000000 0 0 1 0 0 0 1 0 + 2 -0.1154383242012161E-01 0.000000000000000 0 0 0 1 0 0 1 0 + 2 -0.2854093682358197 0.000000000000000 0 0 1 0 0 0 0 1 + 2 0.2388996613150602 0.000000000000000 0 0 0 1 0 0 0 1 + -12 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 530 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 0.2111389035918151E-06 0.000000000000000 1 0 0 0 0 0 0 0 + 1 0.5128703059600361E-05 0.000000000000000 0 1 0 0 0 0 0 0 + 1 1.000000000000000 0.000000000000000 0 0 0 0 1 0 0 0 + 1 0.2019483917365790E-27 0.000000000000000 0 0 0 0 0 1 0 0 + 2 -0.3250909097068332E-01 0.000000000000000 2 0 0 0 0 0 0 0 + 2 1.463033638484707 0.000000000000000 1 1 0 0 0 0 0 0 + 2 0.3251392707166061E-01 0.000000000000000 0 2 0 0 0 0 0 0 + 2 -0.5930226068883956E-02 0.000000000000000 0 0 2 0 0 0 0 0 + 2 -3.474696356238562 0.000000000000000 0 0 1 1 0 0 0 0 + 2 0.5927276679512195E-02 0.2166523348741010E-37 0 0 0 2 0 0 0 0 + 2 -0.5182209068501976E-02 0.000000000000000 1 0 0 0 0 1 0 0 + 2 -0.6567274203089513 0.000000000000000 0 1 0 0 0 1 0 0 + 2 0.1029659397051052 0.000000000000000 1 0 0 0 0 0 1 0 + 2 0.3255952998481525E-01 0.000000000000000 0 1 0 0 0 0 1 0 + 2 -0.1029659338861337 0.000000000000000 1 0 0 0 0 0 0 1 + 2 -0.3255879127402353E-01 0.000000000000000 0 1 0 0 0 0 0 1 + 2 -0.5603890957595181E-07 0.000000000000000 0 0 0 0 0 2 0 0 + 2 -0.5212020631774657E-06 0.000000000000000 0 0 0 0 0 1 1 0 + 2 0.5212021892740850E-06 0.000000000000000 0 0 0 0 0 1 0 1 + -19 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 531 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 1.000000000000000 0.000000000000000 0 0 0 0 0 1 0 0 + -1 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + No Spin Matrix + Quaternion + c_quaternion + + 1, NO = 2, NV = 8, INA = 541 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 0 1.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + 2 -0.3656272367183465E-03 0.000000000000000 2 0 0 0 0 0 0 0 + 2 -0.1970644790041723E-01 0.000000000000000 1 1 0 0 0 0 0 0 + 2 -0.2655327953255107 0.000000000000000 0 2 0 0 0 0 0 0 + 2 -17.12576810110231 0.000000000000000 0 0 2 0 0 0 0 0 + 2 -0.2883300980795104E-01 0.000000000000000 0 0 1 1 0 0 0 0 + 2 -6.967620058230135 0.000000000000000 0 0 0 2 0 0 0 0 + 2 0.2171349107951698E-23 0.000000000000000 1 0 0 0 0 1 0 0 + 2 -0.1985233470127266E-22 0.000000000000000 0 1 0 0 0 1 0 0 + 2 0.9121204216617949E-29 0.000000000000000 0 0 0 0 0 2 0 0 + -10 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 542 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 0.6116118172206402E-02 0.000000000000000 0 0 1 0 0 0 0 0 + 1 3.732993311654496 0.000000000000000 0 0 0 1 0 0 0 0 + 2 0.2443278299565840 -0.7105427357601002E-14 1 0 1 0 0 0 0 0 + 2 -109.2436427525628 -0.1806887972577442E-13 0 1 1 0 0 0 0 0 + 2 -85.72996966047656 -0.1459943277382080E-13 1 0 0 1 0 0 0 0 + 2 0.7825032607721227 0.7105427357601001E-14 0 1 0 1 0 0 0 0 + 2 0.1235187017837736E-01 0.3552713678800500E-14 0 0 1 0 0 1 0 0 + 2 61.25133614217019 -0.7536506141381240E-14 0 0 0 1 0 1 0 0 + 2 10.53475422545662 0.8881784197001252E-15 0 0 1 0 0 0 1 0 + 2 -9.273963311626176 -0.2664535259100376E-14 0 0 0 1 0 0 1 0 + 2 2.790124084553992 0.000000000000000 0 0 1 0 0 0 0 1 + 2 -12.49460045688715 0.000000000000000 0 0 0 1 0 0 0 1 + -12 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 543 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 -0.2704171727972711E-01 0.000000000000000 1 0 0 0 0 0 0 0 + 1 -0.7287424721059020 0.000000000000000 0 1 0 0 0 0 0 0 + 1 0.6617444900424221E-23 0.000000000000000 0 0 0 0 0 1 0 0 + 2 -2.011826271918907 -0.4440892098500626E-15 2 0 0 0 0 0 0 0 + 2 -5.850441748540565 0.1332267629550188E-14 1 1 0 0 0 0 0 0 + 2 2.011826271918907 0.4440892098500626E-15 0 2 0 0 0 0 0 0 + 2 1.082111486715655 0.3721732664138444E-33 0 0 2 0 0 0 0 0 + 2 -23.82368803773236 0.1554312234475219E-14 0 0 1 1 0 0 0 0 + 2 -1.082111486715659 0.1925929944387236E-33 0 0 0 2 0 0 0 0 + 2 -0.5617869815520357E-01 0.2350988701644575E-36 1 0 0 0 0 1 0 0 + 2 2.857948051239394 0.000000000000000 0 1 0 0 0 1 0 0 + 2 -0.3237737723614180 -0.1179611963664229E-15 1 0 0 0 0 0 1 0 + 2 0.8022595317802025E-01 -0.8326672684688675E-16 0 1 0 0 0 0 1 0 + 2 0.3237743806053194 0.000000000000000 1 0 0 0 0 0 0 1 + 2 -0.8023148227749480E-01 0.000000000000000 0 1 0 0 0 0 0 1 + 2 -0.1387778780781446E-16 0.000000000000000 0 0 0 0 0 1 0 1 + -16 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + + 1, NO = 2, NV = 8, INA = 544 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 2 NV = 8 + 1 -5.852478004683411 0.000000000000000 0 0 1 0 0 0 0 0 + 1 -0.1025476998463437E-02 0.000000000000000 0 0 0 1 0 0 0 0 + 2 35.76993929294548 0.2785965902418752E-13 1 0 1 0 0 0 0 0 + 2 -2.424757139212876 -0.1065814103640150E-13 0 1 1 0 0 0 0 0 + 2 2.547762920536604 -0.3552713678800503E-14 1 0 0 1 0 0 0 0 + 2 -34.55658049676505 -0.1986605324688639E-13 0 1 0 1 0 0 0 0 + 2 -62.01949143455905 0.2511879593214417E-14 0 0 1 0 0 1 0 0 + 2 0.2821363443414407E-01 0.7105427357601003E-14 0 0 0 1 0 1 0 0 + 2 14.70856647275445 -0.2664535259100375E-14 0 0 1 0 0 0 1 0 + 2 12.85892134338535 -0.2664535259100376E-14 0 0 0 1 0 0 1 0 + 2 10.79584271946479 0.000000000000000 0 0 1 0 0 0 0 1 + 2 3.488645572261079 0.000000000000000 0 0 0 1 0 0 0 1 + -12 0.000000000000000 0.000000000000000 0 0 0 0 0 0 0 0 + No Stochastic Radiation \ No newline at end of file diff --git a/test/factorize/a4.map b/test/factorize/a4.map new file mode 100644 index 0000000..21e7e1d --- /dev/null +++ b/test/factorize/a4.map @@ -0,0 +1,820 @@ +4 Dimensional TPSA/DA map + + 1, NO = 4, NV = 6, INA = 502 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 4 NV = 6 + 1 -3.248739248164878 0.000000000000000 1 0 0 0 0 0 + 2 9.481060110485997 0.000000000000000 2 0 0 0 0 0 + 2 -1.415180496021588 0.000000000000000 1 1 0 0 0 0 + 2 210.3141529350647 0.000000000000000 0 2 0 0 0 0 + 2 102.1100983383977 -0.9992007221626407E-15 0 0 2 0 0 0 + 2 -3.415504135281418 0.4880507905132675E-33 0 0 1 1 0 0 + 2 -0.6261305504023529 0.000000000000000 1 0 0 0 1 0 + 2 -0.5772675543170588 0.000000000000000 0 1 0 0 1 0 + 2 0.6261492747210899 0.000000000000000 1 0 0 0 0 1 + 2 0.5772697454811320 0.000000000000000 0 1 0 0 0 1 + 2 -8.850217497579548 -0.9992007221626409E-15 0 0 0 2 0 0 + 3 -4685.173279264840 0.3907985046680551E-13 3 0 0 0 0 0 + 3 234.3223107658930 -0.9823870610109868E-31 2 1 0 0 0 0 + 3 19166.48133190933 0.1811883976188255E-11 1 2 0 0 0 0 + 3 301.3283114165062 0.000000000000000 0 3 0 0 0 0 + 3 3678.139976419119 -0.7034373084024992E-12 1 0 2 0 0 0 + 3 1388.671207991414 -0.1818989403545856E-11 0 1 2 0 0 0 + 3 -1626.915756837676 -0.3637978807091714E-11 1 0 1 1 0 0 + 3 -42612.76476344779 0.1691091711109039E-11 0 1 1 1 0 0 + 3 -101.4361857787966 -0.2131628207280300E-13 2 0 0 0 1 0 + 3 232.0833738816372 -0.9947598300641403E-13 1 1 0 0 1 0 + 3 8.832387522478339 -0.7105427357601002E-14 0 2 0 0 1 0 + 3 535.7724175445856 -0.2842170943040401E-13 0 0 2 0 1 0 + 3 101.4321360595929 -0.7105427357601002E-14 2 0 0 0 0 1 + 3 -232.0839803004165 -0.2842170943040401E-13 1 1 0 0 0 1 + 3 -8.825765414943717 -0.7105427357601002E-14 0 2 0 0 0 1 + 3 -346.5190626600597 0.2842170943040401E-13 0 0 2 0 0 1 + 3 655.3786145240102 0.1101341240428155E-11 1 0 0 2 0 0 + 3 -424.6895253606942 -0.1732918607235114E-30 0 1 0 2 0 0 + 3 -430.7799919404032 0.1136868377216160E-12 0 0 1 1 1 0 + 3 -6.408349862891837 0.000000000000000 1 0 0 0 2 0 + 3 -5.713202704273676 0.000000000000000 0 1 0 0 2 0 + 3 660.5579302426089 0.1136868377216160E-12 0 0 1 1 0 1 + 3 -2.069660788587735 0.000000000000000 1 0 0 0 1 1 + 3 4.445557445525357 0.000000000000000 0 1 0 0 1 1 + 3 0.3289068712390379 0.000000000000000 1 0 0 0 0 2 + 3 0.3140163186413550 0.000000000000000 0 1 0 0 0 2 + 3 -404.4587155087173 0.3437400594839027E-31 0 0 0 2 1 0 + 3 116.7870715081967 -0.2842170943040400E-13 0 0 0 2 0 1 + 4 -33667181825.26207 -0.1706182956695557E-05 4 0 0 0 0 0 + 4 -341104325.0654222 -0.1190377424609178E-21 3 1 0 0 0 0 + 4 202004134984.5775 0.1025199890136719E-04 2 2 0 0 0 0 + 4 341161918.6020717 -0.7986926551349102E-22 1 3 0 0 0 0 + 4 -33667110028.01434 -0.1706182956695557E-05 0 4 0 0 0 0 + 4 6213959.888817417 0.1426087692379951E-08 2 0 2 0 0 0 + 4 376011.2968536949 -0.8582410038483518E-25 1 1 2 0 0 0 + 4 -7746383.745573720 0.1455191522836685E-09 0 2 2 0 0 0 + 4 -537220.1464612394 -0.2191882231272756E-09 0 0 4 0 0 0 + 4 186823.6770884130 -0.1862645149230957E-08 2 0 1 1 0 0 + 4 -23659688.37716455 -0.2852175384759903E-08 1 1 1 1 0 0 + 4 -501200.7736398267 0.4253835286682211E-25 0 2 1 1 0 0 + 4 44809.12284060313 0.3492459654808044E-09 0 0 3 1 0 0 + 4 33791.98200378094 -0.9094947017729284E-11 3 0 0 0 1 0 + 4 71045.31625521087 0.1465139121137326E-10 2 1 0 0 1 0 + 4 33541.92126855655 0.4729372449219227E-10 1 2 0 0 1 0 + 4 18141.26952204931 -0.6096456672821660E-11 0 3 0 0 1 0 + 4 126480.6205981716 -0.2000888343900442E-10 1 0 2 0 1 0 + 4 -3854.702078911432 0.1136868377216161E-10 0 1 2 0 1 0 + 4 -23898.97564899129 -0.1818989403545856E-11 3 0 0 0 0 1 + 4 -71045.33902052847 -0.1546140993013978E-10 2 1 0 0 0 1 + 4 -3860.048390878937 0.1091393642127514E-10 1 2 0 0 0 1 + 4 -15825.45127096212 0.2273736754432321E-11 0 3 0 0 0 1 + 4 -15499.11892464682 0.2182787284255027E-10 1 0 2 0 0 1 + 4 -158824.2777224366 -0.2182787284255028E-10 0 1 2 0 0 1 + 4 -5635472.870120240 0.3783497959375381E-09 2 0 0 2 0 0 + 4 -64824.90033233258 0.1862645149230957E-08 1 1 0 2 0 0 + 4 7424152.959663846 -0.4947651177644732E-09 0 2 0 2 0 0 + 4 3269254.300917812 0.6548361852765082E-10 0 0 2 2 0 0 + 4 83253.85471613143 -0.2819433575496078E-10 1 0 1 1 1 0 + 4 99661.25277793621 -0.3637978807091712E-11 0 1 1 1 1 0 + 4 1339.325167986999 0.2842170943040400E-12 2 0 0 0 2 0 + 4 -1138.803123221406 0.6626510555885035E-29 1 1 0 0 2 0 + 4 -2156.824156564607 -0.2842170943040401E-12 0 2 0 0 2 0 + 4 -14160.58516325401 0.6655524071388900E-29 0 0 2 0 2 0 + 4 46481.83646403663 -0.2182787284255027E-10 1 0 1 1 0 1 + 4 187308.4312667271 0.7275957614183428E-11 0 1 1 1 0 1 + 4 310.2005547028029 0.3979039320256561E-12 2 0 0 0 1 1 + 4 1398.218296349253 0.2273736754432320E-12 1 1 0 0 1 1 + 4 -759.7053119568434 -0.2842170943040401E-12 0 2 0 0 1 1 + 4 -2470.665870633805 -0.1818989403545856E-11 0 0 2 0 1 1 + 4 113.0474836289279 0.4263256414560601E-13 2 0 0 0 0 2 + 4 4.514355640863752 0.1065814103640150E-12 1 1 0 0 0 2 + 4 34.40735291421345 0.1421085471520201E-13 0 2 0 0 0 2 + 4 3593.348601322358 -0.9094947017729282E-12 0 0 2 0 0 2 + 4 -45289.65299044129 0.3492459654808044E-09 0 0 1 3 0 0 + 4 42567.91647781737 0.2364686224609613E-10 1 0 0 2 1 0 + 4 120009.0281281247 -0.1136868377216160E-10 0 1 0 2 1 0 + 4 16437.48005618309 -0.9094947017729282E-11 0 0 1 1 2 0 + 4 -34.50208248189761 0.000000000000000 1 0 0 0 3 0 + 4 -28.92352080725891 0.000000000000000 0 1 0 0 3 0 + 4 -74886.12977159879 0.5282129980536478E-27 1 0 0 2 0 1 + 4 41458.04242847222 0.1091393642127514E-10 0 1 0 2 0 1 + 4 -9572.518553154239 -0.5229594535194337E-11 0 0 1 1 1 1 + 4 -39.17198952374213 -0.1776356839400250E-14 1 0 0 0 2 1 + 4 8.108083244749206 0.000000000000000 0 1 0 0 2 1 + 4 -18165.52720298917 0.4657525718287228E-29 0 0 1 1 0 2 + 4 4.207411434908090 -0.8881784197001252E-15 1 0 0 0 1 2 + 4 13.20206318155829 -0.1776356839400250E-14 0 1 0 0 1 2 + 4 -0.6767195444266404E-01 -0.1526556658859590E-15 1 0 0 0 0 3 + 4 0.1333085415517137 0.1387778780781446E-16 0 1 0 0 0 3 + 4 196776.3514213295 0.2082742867060006E-09 0 0 0 4 0 0 + 4 13954.37174124090 -0.4446733036595578E-29 0 0 0 2 2 0 + 4 -386.2652916380869 0.1818989403545856E-11 0 0 0 2 1 1 + 4 -4150.506090117280 0.9094947017729280E-12 0 0 0 2 0 2 + -104 0.000000000000000 0.000000000000000 0 0 0 0 0 0 + + 1, NO = 4, NV = 6, INA = 503 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 4 NV = 6 + 1 -0.1262958087429899E-01 0.000000000000000 1 0 0 0 0 0 + 1 -0.3078117151337621 0.000000000000000 0 1 0 0 0 0 + 2 -1.177767725300229 0.000000000000000 2 0 0 0 0 0 + 2 -1.802125513034837 0.000000000000000 1 1 0 0 0 0 + 2 0.8846459970520728 0.000000000000000 0 2 0 0 0 0 + 2 0.5771097021169346 0.000000000000000 0 0 2 0 0 0 + 2 -7.959512334345894 -0.3009265538105055E-35 0 0 1 1 0 0 + 2 -0.5712907938500165E-01 0.000000000000000 1 0 0 0 1 0 + 2 0.5708050329320468E-01 0.000000000000000 0 1 0 0 1 0 + 2 0.5712935978493520E-01 0.000000000000000 1 0 0 0 0 1 + 2 -0.5708226886770704E-01 0.000000000000000 0 1 0 0 0 1 + 2 -0.4013131322036372 0.000000000000000 0 0 0 2 0 0 + 3 -27.41382149142134 0.2113488546334126E-30 3 0 0 0 0 0 + 3 1321.099387988194 0.4263256414560601E-13 2 1 0 0 0 0 + 3 209.9628897031683 -0.1942565025904414E-29 1 2 0 0 0 0 + 3 -526.6380806981296 -0.7815970093361102E-13 0 3 0 0 0 0 + 3 206.0802895563594 0.1421085471520200E-12 1 0 2 0 0 0 + 3 -842.2304665058783 0.3552713678800500E-13 0 1 2 0 0 0 + 3 -1012.248298644942 0.1136868377216160E-12 1 0 1 1 0 0 + 3 -55.67696376669434 -0.2842170943040400E-12 0 1 1 1 0 0 + 3 1.939997066444821 -0.8881784197001252E-15 2 0 0 0 1 0 + 3 20.36104297299129 0.2664535259100376E-14 1 1 0 0 1 0 + 3 -7.605053312034003 0.2664535259100376E-14 0 2 0 0 1 0 + 3 7.184464324714440 -0.8881784197001249E-14 0 0 2 0 1 0 + 3 -1.940030153044944 -0.8881784197001252E-15 2 0 0 0 0 1 + 3 -20.36029690548624 0.4440892098500627E-14 1 1 0 0 0 1 + 3 7.605095722899843 0.8881784197001254E-15 0 2 0 0 0 1 + 3 -13.09515938889714 -0.5329070518200750E-14 0 0 2 0 0 1 + 3 -77.59010037803536 -0.3126388037344441E-12 1 0 0 2 0 0 + 3 -106.2696317171988 -0.6394884621840902E-13 0 1 0 2 0 0 + 3 62.87687928322872 0.3552713678800501E-14 0 0 1 1 1 0 + 3 -0.5662274980094570 0.000000000000000 1 0 0 0 2 0 + 3 0.5638160086324399 0.000000000000000 0 1 0 0 2 0 + 3 -22.52132587030051 0.7105427357601001E-14 0 0 1 1 0 1 + 3 0.4131620333885505 0.000000000000000 1 0 0 0 1 1 + 3 0.2556841235935167 0.000000000000000 0 1 0 0 1 1 + 3 0.3103107077872048E-01 0.000000000000000 1 0 0 0 0 2 + 3 -0.5109568765498875E-01 0.000000000000000 0 1 0 0 0 2 + 3 -46.85179180996303 0.5329070518200751E-14 0 0 0 2 1 0 + 3 56.94480857526042 0.1776356839400251E-14 0 0 0 2 0 1 + 4 -138965140.2762788 0.4913952974853399E-23 4 0 0 0 0 0 + 4 12758289161.06268 0.6556510925292969E-06 3 1 0 0 0 0 + 4 833793028.1655724 -0.6039369188499941E-23 2 2 0 0 0 0 + 4 -12758452531.17336 -0.6556510925292969E-06 1 3 0 0 0 0 + 4 -138966995.5975010 0.2877667396168854E-24 0 4 0 0 0 0 + 4 44048.98359347768 -0.1164153218269348E-09 2 0 2 0 0 0 + 4 -1262562.107277335 -0.2328306436538696E-09 1 1 2 0 0 0 + 4 -58166.30696803852 0.1364989983916642E-25 0 2 2 0 0 0 + 4 -3486.589790161811 -0.1455191522836686E-10 0 0 4 0 0 0 + 4 -1104985.696353827 0.8731149137020110E-10 2 0 1 1 0 0 + 4 -167932.8295013943 0.2328306436538696E-09 1 1 1 1 0 0 + 4 1160725.209400791 0.1309672370553017E-09 0 2 1 1 0 0 + 4 95251.17035522305 -0.3410605131648481E-10 0 0 3 1 0 0 + 4 5305.341956662323 -0.1136868377216160E-12 3 0 0 0 1 0 + 4 -9439.703503887979 0.1818989403545856E-11 2 1 0 0 1 0 + 4 -6409.978515603184 0.2109994899108679E-28 1 2 0 0 1 0 + 4 -1588.470710488770 -0.4547473508864641E-12 0 3 0 0 1 0 + 4 20.37664973183812 -0.7602807272633072E-12 1 0 2 0 1 0 + 4 -10215.42762659664 0.1477928890381008E-11 0 1 2 0 1 0 + 4 -5047.467168876582 0.6088486430127847E-29 3 0 0 0 0 1 + 4 6627.697841699120 0.9094947017729284E-12 2 1 0 0 0 1 + 4 6525.375596227186 -0.3303641373977116E-28 1 2 0 0 0 1 + 4 660.0038659174996 -0.4547473508864641E-12 0 3 0 0 0 1 + 4 1960.522159424615 -0.2046363078989089E-11 1 0 2 0 0 1 + 4 583.1074072974976 -0.1364242052659392E-11 0 1 2 0 0 1 + 4 -37011.27274590797 0.1164153218269348E-09 2 0 0 2 0 0 + 4 1064583.948954450 0.4365574568510056E-10 1 1 0 2 0 0 + 4 41222.89027577251 0.6542266726571029E-26 0 2 0 2 0 0 + 4 9932.040481477143 0.2910383045673370E-10 0 0 2 2 0 0 + 4 4719.084019776439 0.1591615728102625E-11 1 0 1 1 1 0 + 4 -14237.15739890776 0.1520561454526614E-11 0 1 1 1 1 0 + 4 -58.45193152344128 -0.8446773133120048E-30 2 0 0 0 2 0 + 4 -258.4854886633433 -0.8526512829121205E-13 1 1 0 0 2 0 + 4 84.23737362316231 0.1920359087100780E-30 0 2 0 0 2 0 + 4 -819.6891918475309 0.3410605131648479E-12 0 0 2 0 2 0 + 4 106.7503396705229 -0.1818989403545856E-11 1 0 1 1 0 1 + 4 4375.972678140909 0.2273736754432320E-11 0 1 1 1 0 1 + 4 100.1688910079888 0.8881784197001257E-14 2 0 0 0 1 1 + 4 -41.86646823319943 -0.8526512829121201E-13 1 1 0 0 1 1 + 4 -106.2634964315828 -0.2664535259100376E-13 0 2 0 0 1 1 + 4 488.6408448701696 0.2486899575160351E-12 0 0 2 0 1 1 + 4 -27.31638278565380 0.7993605777301127E-14 2 0 0 0 0 2 + 4 -24.36193535944186 -0.1776356839400251E-14 1 1 0 0 0 2 + 4 3.567638727828086 -0.7993605777301127E-14 0 2 0 0 0 2 + 4 815.6800426683808 -0.4263256414560602E-13 0 0 2 0 0 2 + 4 -175483.3979153981 0.7275957614183426E-11 0 0 1 3 0 0 + 4 6655.252048947346 0.7602807272633075E-12 1 0 0 2 1 0 + 4 -4990.300120438560 -0.1932676241267472E-11 0 1 0 2 1 0 + 4 -2488.784446088441 -0.3410605131648481E-12 0 0 1 1 2 0 + 4 -2.877939923735668 0.000000000000000 1 0 0 0 3 0 + 4 2.734225040434004 0.000000000000000 0 1 0 0 3 0 + 4 -7017.877905581283 0.2273736754432320E-11 1 0 0 2 0 1 + 4 7412.369173851309 0.4547473508864640E-12 0 1 0 2 0 1 + 4 -8.767161701112489 -0.5684341886080799E-12 0 0 1 1 1 1 + 4 0.3797370157412579 0.000000000000000 1 0 0 0 2 1 + 4 4.231284321402064 0.1665334536937735E-15 0 1 0 0 2 1 + 4 463.6298642006412 -0.2371939365409428E-29 0 0 1 1 0 2 + 4 1.510161393575798 -0.1387778780781445E-15 1 0 0 0 1 2 + 4 -0.3866038193896253 0.1110223024625157E-15 0 1 0 0 1 2 + 4 0.9002838628214346E-02 0.7821055249433517E-34 1 0 0 0 0 3 + 4 -0.1973297071025056E-01 0.1474514954580286E-16 0 1 0 0 0 3 + 4 1291.045502378100 -0.1091393642127514E-10 0 0 0 4 0 0 + 4 960.5097072190950 -0.3410605131648481E-12 0 0 0 2 2 0 + 4 -589.0572964623660 -0.2486899575160351E-12 0 0 0 2 1 1 + 4 -775.4429156799174 0.4263256414560600E-13 0 0 0 2 0 2 + -105 0.000000000000000 0.000000000000000 0 0 0 0 0 0 + + 1, NO = 4, NV = 6, INA = 504 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 4 NV = 6 + 1 1.998246462619670 0.000000000000000 0 0 1 0 0 0 + 2 -51.58521946741487 -0.1665334536937735E-14 1 0 1 0 0 0 + 2 2.100820821566119 0.1776356839400250E-14 0 1 1 0 0 0 + 2 -4.763768354838218 -0.1776356839400250E-14 1 0 0 1 0 0 + 2 10.88724853368485 -0.1887379141862766E-14 0 1 0 1 0 0 + 2 0.4448179007408107E-01 0.000000000000000 0 0 1 0 1 0 + 2 -0.9502844475768885 0.000000000000000 0 0 1 0 0 1 + 2 0.5063300944771626 0.000000000000000 0 0 0 1 1 0 + 2 -1.142662996648278 0.000000000000000 0 0 0 1 0 1 + 3 -2437.652456881256 0.2060573933704291E-12 2 0 1 0 0 0 + 3 636.7901138881271 -0.9094947017729282E-12 1 1 1 0 0 0 + 3 16787.52893811051 -0.1492139745096210E-12 0 2 1 0 0 0 + 3 4310.822617948071 -0.3907985046680551E-13 0 0 3 0 0 0 + 3 -463.3741968182840 -0.4547473508864641E-12 2 0 0 1 0 0 + 3 -1020.802893935940 0.3836930773104542E-12 1 1 0 1 0 0 + 3 567.2411346161870 0.4547473508864641E-12 0 2 0 1 0 0 + 3 -85.76398921115168 -0.1612136433382898E-28 0 0 2 1 0 0 + 3 407.7448455822869 -0.7105427357601003E-13 1 0 1 0 1 0 + 3 287.2700519829958 -0.7105427357601002E-14 0 1 1 0 1 0 + 3 -127.4356444509381 0.1278976924368180E-12 1 0 1 0 0 1 + 3 -469.5557197300802 -0.2753353101070388E-13 0 1 1 0 0 1 + 3 -706.1594903792738 0.1421085471520201E-12 0 0 1 2 0 0 + 3 -573.9071279183244 -0.7105427357601002E-14 1 0 0 1 1 0 + 3 494.3171609316425 -0.9947598300641403E-13 0 1 0 1 1 0 + 3 -0.8290718185653297 0.000000000000000 0 0 1 0 2 0 + 3 705.2027400385105 0.2930988785010413E-13 1 0 0 1 0 1 + 3 -144.6990124214067 0.9947598300641404E-13 0 1 0 1 0 1 + 3 -13.18670957317878 0.000000000000000 0 0 1 0 1 1 + 3 -5.292315103875410 0.000000000000000 0 0 1 0 0 2 + 3 52.34624898655001 0.000000000000000 0 0 0 3 0 0 + 3 -10.14479150594820 0.000000000000000 0 0 0 1 2 0 + 3 3.398293482789045 0.000000000000000 0 0 0 1 1 1 + 3 -7.039561137278916 0.000000000000000 0 0 0 1 0 2 + 4 -2375876.821163453 0.7639755494892597E-10 3 0 1 0 0 0 + 4 -186998.5289005467 0.1110029555453186E-27 2 1 1 0 0 0 + 4 7635181.471509907 0.8003553375601768E-09 1 2 1 0 0 0 + 4 109375.0378426104 -0.2328306436538696E-09 0 3 1 0 0 0 + 4 424207.1607729082 -0.1855369191616774E-09 1 0 3 0 0 0 + 4 5072.278290401056 -0.2328306436538696E-09 0 1 3 0 0 0 + 4 -71018.03032064630 -0.2328306436538696E-09 3 0 0 1 0 0 + 4 6889857.590520194 0.5675246939063072E-09 2 1 0 1 0 0 + 4 155689.0582784669 -0.1623889451645243E-25 1 2 0 1 0 0 + 4 -2995702.337719225 0.7639755494892597E-10 0 3 0 1 0 0 + 4 -68253.82837521224 -0.9313225746154784E-09 1 0 2 1 0 0 + 4 -4703456.863359069 -0.6039044819772243E-09 0 1 2 1 0 0 + 4 1505.698831822784 0.1506350599811412E-10 2 0 1 0 1 0 + 4 -89679.34405200719 -0.1953581695597466E-27 1 1 1 0 1 0 + 4 -41941.75362010388 -0.7787548383930698E-11 0 2 1 0 1 0 + 4 10023.43358238364 -0.1648459146963433E-11 0 0 3 0 1 0 + 4 5755.992471105867 0.5456968210637569E-11 2 0 1 0 0 1 + 4 18686.32300165093 0.1818989403545856E-11 1 1 1 0 0 1 + 4 -59179.66916563535 -0.5456968210637569E-11 0 2 1 0 0 1 + 4 -15858.07295985580 0.1818989403545856E-11 0 0 3 0 0 1 + 4 -3380218.487714741 0.1309672370553017E-09 1 0 1 2 0 0 + 4 18993.54510666351 0.9313225746154787E-09 0 1 1 2 0 0 + 4 46352.18025551688 0.4074301947176891E-28 2 0 0 1 1 0 + 4 -59697.27467221003 0.2080469130305574E-10 1 1 0 1 1 0 + 4 -39820.74064639171 0.1024386874536019E-27 0 2 0 1 1 0 + 4 -10822.60900755067 -0.1091393642127514E-10 0 0 2 1 1 0 + 4 -16510.54218610017 0.9322320693172514E-11 1 0 1 0 2 0 + 4 -10653.88440897743 -0.3865352482534945E-11 0 1 1 0 2 0 + 4 -47911.09398561643 0.9094947017729282E-12 2 0 0 1 0 1 + 4 94152.66960436666 -0.3637978807091712E-11 1 1 0 1 0 1 + 4 -64125.29471540351 -0.9094947017729281E-12 0 2 0 1 0 1 + 4 45592.09761052106 -0.1000444171950221E-10 0 0 2 1 0 1 + 4 114.2669403100872 -0.3808509063674137E-11 1 0 1 0 1 1 + 4 5587.723358243999 -0.1421085471520201E-11 0 1 1 0 1 1 + 4 3719.885199363015 -0.6451728040701710E-11 1 0 1 0 0 2 + 4 11099.55909941662 -0.1136868377216160E-11 0 1 1 0 0 2 + 4 23851.57832277831 0.3492459654808044E-09 1 0 0 3 0 0 + 4 -486574.6436311910 -0.2455635694786906E-09 0 1 0 3 0 0 + 4 -39310.40808196925 0.2273736754432321E-11 0 0 1 2 1 0 + 4 11374.62432617859 0.3410605131648481E-11 1 0 0 1 2 0 + 4 -17443.22780515143 0.8867573342286050E-11 0 1 0 1 2 0 + 4 -5.703171374735136 0.000000000000000 0 0 1 0 3 0 + 4 -4192.462142661279 0.6524675676724749E-27 0 0 1 2 0 1 + 4 -6818.207729479502 0.1563194018672220E-11 1 0 0 1 1 1 + 4 764.4337022700502 -0.3467448550509290E-11 0 1 0 1 1 1 + 4 12.64219374034109 0.000000000000000 0 0 1 0 2 1 + 4 -10580.01628785375 0.1591615728102625E-11 1 0 0 1 0 2 + 4 4768.737220177995 -0.6281197784119285E-11 0 1 0 1 0 2 + 4 -104.1247646651883 -0.7105427357601002E-14 0 0 1 0 1 2 + 4 -9.225802567032620 0.000000000000000 0 0 1 0 0 3 + 4 -7009.101249519512 0.4547473508864641E-11 0 0 0 3 1 0 + 4 -35.96482766981881 0.000000000000000 0 0 0 1 3 0 + 4 -27876.14999993544 0.3183231456205249E-11 0 0 0 3 0 1 + 4 4.558536311212326 0.7993605777301127E-14 0 0 0 1 2 1 + 4 -33.47257881814375 0.000000000000000 0 0 0 1 1 2 + 4 -19.35335402181189 0.000000000000000 0 0 0 1 0 3 + -89 0.000000000000000 0.000000000000000 0 0 0 0 0 0 + + 1, NO = 4, NV = 6, INA = 505 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 4 NV = 6 + 1 -0.1593841720719695E-02 0.000000000000000 0 0 1 0 0 0 + 1 0.5004387690440424 0.000000000000000 0 0 0 1 0 0 + 2 -0.5446393799924619 0.1156344375593035E-33 1 0 1 0 0 0 + 2 31.45659046770351 -0.3335005882565411E-15 0 1 1 0 0 0 + 2 12.92274847319755 0.3326332265185527E-15 1 0 0 1 0 0 + 2 -0.5348112742120309 0.2373376099491244E-33 0 1 0 1 0 0 + 2 0.1267693034948146 0.000000000000000 0 0 1 0 1 0 + 2 -0.2854093682358197 0.000000000000000 0 0 1 0 0 1 + 2 -0.1154383242012161E-01 0.000000000000000 0 0 0 1 1 0 + 2 0.2388996613150602 0.000000000000000 0 0 0 1 0 1 + 3 -184.0090065439723 -0.1850411561510876E-30 2 0 1 0 0 0 + 3 2127.506399889023 0.1421085471520200E-13 1 1 1 0 0 0 + 3 231.5924025047428 0.4547473508864640E-12 0 2 1 0 0 0 + 3 -16.61530677996715 -0.4840083374801997E-31 0 0 3 0 0 0 + 3 945.7550544009289 -0.3552713678800502E-13 2 0 0 1 0 0 + 3 -264.0143460238984 0.9094947017729280E-12 1 1 0 1 0 0 + 3 -4032.753366166903 0.7105427357601005E-14 0 2 0 1 0 0 + 3 -2427.238459701021 0.2842170943040400E-13 0 0 2 1 0 0 + 3 -14.69718456096555 -0.3196048724977367E-31 1 0 1 0 1 0 + 3 158.0986057348154 -0.1421085471520200E-13 0 1 1 0 1 0 + 3 39.93347984840327 0.3552713678800500E-14 1 0 1 0 0 1 + 3 -85.55505881842743 0.3730349362740525E-13 0 1 1 0 0 1 + 3 57.94394627124397 0.1433069928553915E-30 0 0 1 2 0 0 + 3 -102.6833277776425 0.1421085471520200E-13 1 0 0 1 1 0 + 3 -63.65246162970041 -0.1882576314932881E-31 0 1 0 1 1 0 + 3 -2.539989765595450 0.000000000000000 0 0 1 0 2 0 + 3 44.65700247606231 -0.3375077994860475E-13 1 0 0 1 0 1 + 3 97.66202742248036 -0.3552713678800501E-14 0 1 0 1 0 1 + 3 0.8615830766926635 0.000000000000000 0 0 1 0 1 1 + 3 -1.758759123702398 0.000000000000000 0 0 1 0 0 2 + 3 81.95264493809327 -0.1065814103640150E-13 0 0 0 3 0 0 + 3 0.2481022663129712 0.000000000000000 0 0 0 1 2 0 + 3 3.144137586901360 0.000000000000000 0 0 0 1 1 1 + 3 1.607834178102175 0.000000000000000 0 0 0 1 0 2 + 4 -16737.41875016026 0.5820766091346741E-10 3 0 1 0 0 0 + 4 1878947.911554564 0.1455191522836685E-09 2 1 1 0 0 0 + 4 84818.68457195300 0.1164153218269348E-09 1 2 1 0 0 0 + 4 -772600.4766642855 0.3819877747446298E-10 0 3 1 0 0 0 + 4 -5371.223602253312 -0.2910383045673370E-10 1 0 3 0 0 0 + 4 -447778.3293925586 -0.5820766091346740E-10 0 1 3 0 0 0 + 4 635813.0424732398 -0.2728484105318786E-11 3 0 0 1 0 0 + 4 15665.03782504595 -0.1164153218269348E-09 2 1 0 1 0 0 + 4 -2129680.094253219 -0.1818989403545856E-09 1 2 0 1 0 0 + 4 -6078.182079045475 0.5820766091346739E-10 0 3 0 1 0 0 + 4 -333702.8009213517 -0.9276845958083867E-10 1 0 2 1 0 0 + 4 -649.5586133701945 0.2910383045673370E-09 0 1 2 1 0 0 + 4 1097.991431426030 0.9094947017729281E-12 2 0 1 0 1 0 + 4 35279.95623999897 -0.2273736754432321E-11 1 1 1 0 1 0 + 4 2523.435088644259 -0.9094947017729281E-12 0 2 1 0 1 0 + 4 -2512.708677225705 0.4547473508864642E-12 0 0 3 0 1 0 + 4 -3345.760085053993 0.2046363078989089E-11 2 0 1 0 0 1 + 4 -2232.155121835779 0.4547473508864646E-12 1 1 1 0 0 1 + 4 -30141.53598417943 -0.1591615728102624E-11 0 2 1 0 0 1 + 4 -894.5213621030287 0.2273736754432321E-12 0 0 3 0 0 1 + 4 35779.13704508558 0.1746229827404022E-09 1 0 1 2 0 0 + 4 1172122.492502885 -0.5456968210637569E-10 0 1 1 2 0 0 + 4 -5621.452430740526 -0.1136868377216160E-11 2 0 0 1 1 0 + 4 10080.57598050411 0.8766323352482859E-29 1 1 0 1 1 0 + 4 19603.80999299334 -0.1136868377216160E-11 0 2 0 1 1 0 + 4 -9842.626419258990 0.1818989403545856E-11 0 0 2 1 1 0 + 4 2671.301735256684 0.1065814103640150E-11 1 0 1 0 2 0 + 4 -4424.717014856576 0.2241762331323116E-11 0 1 1 0 2 0 + 4 1198.845220712692 -0.1591615728102624E-11 2 0 0 1 0 1 + 4 11011.75037699349 0.1364242052659392E-11 1 1 0 1 0 1 + 4 7761.399801271529 0.1591615728102624E-11 0 2 0 1 0 1 + 4 9185.838229065754 -0.4563544265766603E-27 0 0 2 1 0 1 + 4 -1801.379330394780 0.2984279490192421E-12 1 0 1 0 1 1 + 4 -512.3008126806978 -0.8881784197001250E-12 0 1 1 0 1 1 + 4 -2520.023564467173 0.3197442310920451E-12 1 0 1 0 0 2 + 4 1184.075215147580 -0.1563194018672220E-11 0 1 1 0 0 2 + 4 286680.2441106783 -0.1773514668457210E-10 1 0 0 3 0 0 + 4 -4951.650232639699 -0.8731149137020111E-10 0 1 0 3 0 0 + 4 6786.320977680431 0.9094947017729278E-12 0 0 1 2 1 0 + 4 4110.905599310368 -0.2305711177541526E-11 1 0 0 1 2 0 + 4 2582.535172650646 0.7531752999057062E-12 0 1 0 1 2 0 + 4 -9.002760124097602 0.000000000000000 0 0 1 0 3 0 + 4 -15663.52135836102 0.9094947017729284E-12 0 0 1 2 0 1 + 4 186.2091910374405 0.9308109838457310E-12 1 0 0 1 1 1 + 4 -1610.006826937826 0.2842170943040401E-12 0 1 0 1 1 1 + 4 -2.700995546380527 0.2040034807748725E-14 0 0 1 0 2 1 + 4 -928.4083046737449 0.1620037437533028E-11 1 0 0 1 0 2 + 4 -2734.900196835320 0.3055333763768431E-12 0 1 0 1 0 2 + 4 -4.465724969634664 0.1398356251017995E-32 0 0 1 0 1 2 + 4 -4.842168200438589 0.000000000000000 0 0 1 0 0 3 + 4 4069.562187340414 -0.4547473508864640E-12 0 0 0 3 1 0 + 4 0.1594809253216932 0.000000000000000 0 0 0 1 3 0 + 4 330.5129839635260 -0.6821210263296961E-12 0 0 0 3 0 1 + 4 0.2365357869119832 0.7710480085073675E-33 0 0 0 1 2 1 + 4 27.23868941831817 0.1776356839400250E-14 0 0 0 1 1 2 + 4 5.734454405231163 0.000000000000000 0 0 0 1 0 3 + -90 0.000000000000000 0.000000000000000 0 0 0 0 0 0 + No Spin Matrix + Quaternion + c_quaternion + + 1, NO = 4, NV = 6, INA = 515 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 4 NV = 6 + 0 1.000000000000000 0.000000000000000 0 0 0 0 0 0 + 2 -0.3656272367183465E-03 0.000000000000000 2 0 0 0 0 0 + 2 -0.1970644790041723E-01 0.000000000000000 1 1 0 0 0 0 + 2 -0.2655327953255107 0.000000000000000 0 2 0 0 0 0 + 2 -17.12576810110231 0.000000000000000 0 0 2 0 0 0 + 2 -0.2883300980795104E-01 0.000000000000000 0 0 1 1 0 0 + 2 -6.967620058230135 0.000000000000000 0 0 0 2 0 0 + 3 -0.5440323726115835E-01 -0.1047224407260559E-32 3 0 0 0 0 0 + 3 -1.624309242571331 -0.3274080905458301E-31 2 1 0 0 0 0 + 3 -4.209062145481870 0.1110223024625156E-14 1 2 0 0 0 0 + 3 1.466103250845785 0.4440892098500626E-15 0 3 0 0 0 0 + 3 209.3705507558322 0.1731947918415244E-12 1 0 2 0 0 0 + 3 -12.73411019538203 -0.4973799150320701E-13 0 1 2 0 0 0 + 3 13.91543473735720 -0.1421085471520199E-13 1 0 1 1 0 0 + 3 188.1955547313183 -0.4796163466380678E-13 0 1 1 1 0 0 + 3 -0.8755398814788167E-02 0.4092601131822876E-33 2 0 0 0 1 0 + 3 -0.2337782517293764 -0.4163336342344336E-16 1 1 0 0 1 0 + 3 0.5846405944600266E-01 -0.1387778780781446E-16 0 2 0 0 1 0 + 3 86.01712996046119 -0.1065814103640150E-13 0 0 2 0 1 0 + 3 0.8755415262747812E-02 -0.4814824860968090E-34 2 0 0 0 0 1 + 3 0.2337785454661968 -0.6938893903907227E-17 1 1 0 0 0 1 + 3 -0.5846808873562253E-01 0.2775557561562891E-16 0 2 0 0 0 1 + 3 63.16536732907296 0.7041681359165831E-33 0 0 2 0 0 1 + 3 320.0027538702856 0.6528111384795920E-13 1 0 0 2 0 0 + 3 -3.745097017178480 -0.2131628207280301E-13 0 1 0 2 0 0 + 3 36.00219121470011 -0.1065814103640150E-13 0 0 1 1 1 0 + 3 10.08919627276462 0.1776356839400250E-14 0 0 1 1 0 1 + 3 34.63282954289235 0.7105427357601003E-14 0 0 0 2 1 0 + 3 46.64583746314502 -0.8881784197001254E-15 0 0 0 2 0 1 + 4 -4.220180306121361 0.2131628207280300E-13 4 0 0 0 0 0 + 4 -13.33425980098162 -0.1136868377216160E-12 3 1 0 0 0 0 + 4 1550.569246890450 -0.2557953848736360E-12 2 2 0 0 0 0 + 4 260.0536291651096 -0.2273736754432320E-12 1 3 0 0 0 0 + 4 -1034.690262476448 -0.2025046796916286E-12 0 4 0 0 0 0 + 4 -267993.8476833617 -0.7543121682829220E-10 2 0 2 0 0 0 + 4 202121.0735500715 -0.1509761204943061E-09 1 1 2 0 0 0 + 4 1219665.198277253 -0.1460875864722766E-10 0 2 2 0 0 0 + 4 136554.0976086658 -0.3061018105654510E-10 0 0 4 0 0 0 + 4 -168209.8628183886 0.3292370820418000E-09 2 0 1 1 0 0 + 4 -2429794.112376046 -0.8208189683500679E-10 1 1 1 1 0 0 + 4 169527.1467763509 0.4911271389573812E-10 0 2 1 1 0 0 + 4 -166.8553722984571 0.1818989403545856E-11 0 0 3 1 0 0 + 4 0.4512979555731385 0.8881784197001249E-15 3 0 0 0 1 0 + 4 25.81381191579263 -0.2465190328815662E-31 2 1 0 0 1 0 + 4 -58.43800561957799 -0.3552713678800499E-14 1 2 0 0 1 0 + 4 -29.87718239955602 0.5329070518200752E-14 0 3 0 0 1 0 + 4 -26313.14835824511 -0.1659827830735595E-10 1 0 2 0 1 0 + 4 6959.604375357429 0.7730704965069893E-11 0 1 2 0 1 0 + 4 -0.3693996708817561 0.1776356839400250E-14 3 0 0 0 0 1 + 4 -25.00630190021360 -0.8881784197001252E-15 2 1 0 0 0 1 + 4 20.64087604115134 0.1065814103640150E-13 1 2 0 0 0 1 + 4 27.67016406264368 -0.1776356839400250E-14 0 3 0 0 0 1 + 4 11992.81435801683 0.1136868377216160E-11 1 0 2 0 0 1 + 4 -9732.559488254474 0.2387423592153937E-11 0 1 2 0 0 1 + 4 822772.9561142086 -0.4098410499864257E-10 2 0 0 2 0 0 + 4 -133417.2699781793 0.1509761204943061E-09 1 1 0 2 0 0 + 4 -218736.1531302046 -0.5815081749460661E-10 0 2 0 2 0 0 + 4 193149.2511280381 -0.5752553988713772E-10 0 0 2 2 0 0 + 4 14432.77518507465 0.1818989403545856E-11 1 0 1 1 1 0 + 4 -7958.607666140106 0.8185452315956357E-11 0 1 1 1 1 0 + 4 -0.7829139856973288E-01 0.1110223024625156E-15 2 0 0 0 2 0 + 4 -0.6325130963894113 -0.1110223024625157E-15 1 1 0 0 2 0 + 4 1.043950727614690 -0.1665334536937735E-15 0 2 0 0 2 0 + 4 -454.2851238539698 0.1278976924368180E-12 0 0 2 0 2 0 + 4 -18984.69520330591 -0.1136868377216160E-11 1 0 1 1 0 1 + 4 1975.852900865143 -0.7275957614183428E-11 0 1 1 1 0 1 + 4 0.1281712314194305 0.8326672684688673E-16 2 0 0 0 1 1 + 4 0.6085429789579322 -0.5551115123125779E-16 1 1 0 0 1 1 + 4 0.8544310794952416 0.6383782391594650E-15 0 2 0 0 1 1 + 4 -280.8019926069326 0.2629008122312370E-12 0 0 2 0 1 1 + 4 -0.5703823085211572E-01 -0.5204170427930420E-16 2 0 0 0 0 2 + 4 -0.1038681401714410 -0.3469446951953614E-16 1 1 0 0 0 2 + 4 -0.1447639537583246 0.3642919299551295E-15 0 2 0 0 0 2 + 4 -85.80199861133082 0.7370293155926903E-31 0 0 2 0 0 2 + 4 890.3789475009453 -0.1818989403545857E-10 0 0 1 3 0 0 + 4 -11001.84988694094 -0.1250555214937776E-10 1 0 0 2 1 0 + 4 12090.55405205613 0.6821210263296963E-11 0 1 0 2 1 0 + 4 -172.4648577559044 -0.2842170943040400E-13 0 0 1 1 2 0 + 4 5254.881659099812 0.1591615728102625E-11 1 0 0 2 0 1 + 4 -17496.35717054942 0.3069544618483634E-11 0 1 0 2 0 1 + 4 204.1079407754987 0.2140509991477302E-12 0 0 1 1 1 1 + 4 65.01618285851188 -0.7105427357601003E-14 0 0 1 1 0 2 + 4 96990.01973450312 -0.2319211489520968E-10 0 0 0 4 0 0 + 4 -56.13242077714457 -0.8526512829121202E-13 0 0 0 2 2 0 + 4 -290.1698035428681 -0.9237055564881304E-13 0 0 0 2 1 1 + 4 -101.0781143428605 0.000000000000000 0 0 0 2 0 2 + -86 0.000000000000000 0.000000000000000 0 0 0 0 0 0 + + 1, NO = 4, NV = 6, INA = 516 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 4 NV = 6 + 1 0.6116118172206402E-02 0.000000000000000 0 0 1 0 0 0 + 1 3.732993311654496 0.000000000000000 0 0 0 1 0 0 + 2 0.2443278299565911 -0.7105427357601002E-14 1 0 1 0 0 0 + 2 -109.2436427525628 -0.2076117056049043E-13 0 1 1 0 0 0 + 2 -85.72996966047656 -0.1743050148661496E-13 1 0 0 1 0 0 + 2 0.7825032607721139 0.7105427357601001E-14 0 1 0 1 0 0 + 2 10.53475422545664 0.8881784197001252E-15 0 0 1 0 1 0 + 2 2.790124084554003 0.000000000000000 0 0 1 0 0 1 + 2 -9.273963311626169 -0.2664535259100376E-14 0 0 0 1 1 0 + 2 -12.49460045688714 0.2220446049250313E-15 0 0 0 1 0 1 + 3 18154.26970781250 -0.2182787284255027E-10 2 0 1 0 0 0 + 3 248158.5503465072 0.6650680006714539E-11 1 1 1 0 0 0 + 3 -18087.17949499934 0.7275957614183426E-11 0 2 1 0 0 0 + 3 56.21997889920625 0.3410605131648481E-11 0 0 3 0 0 0 + 3 -221397.5897471853 0.1307398633798585E-11 2 0 0 1 0 0 + 3 35660.51805129212 -0.2910383045673370E-10 1 1 0 1 0 0 + 3 58033.59480877340 0.9208633855450898E-11 0 2 0 1 0 0 + 3 -15754.19398872578 0.1421085471520200E-11 0 0 2 1 0 0 + 3 1123.903691641974 0.1818989403545856E-11 1 0 1 0 1 0 + 3 5222.804679981761 0.3183231456205249E-11 0 1 1 0 1 0 + 3 -1860.496719419830 0.6821210263296962E-12 1 0 1 0 0 1 + 3 -2492.333736143510 0.2273736754432321E-12 0 1 1 0 0 1 + 3 -260.0589243783544 0.4547473508864643E-11 0 0 1 2 0 0 + 3 2726.328523004363 0.3183231456205248E-11 1 0 0 1 1 0 + 3 -3113.753895961538 -0.1818989403545856E-11 0 1 0 1 1 0 + 3 -4.243157214318569 0.1598721155460225E-13 0 0 1 0 2 0 + 3 -1697.017081358399 -0.4547473508864641E-12 1 0 0 1 0 1 + 3 4751.797916469847 -0.6821210263296961E-12 0 1 0 1 0 1 + 3 -47.87228700600554 0.1421085471520200E-13 0 0 1 0 1 1 + 3 -4.217779348897849 0.1925929944387236E-33 0 0 1 0 0 2 + 3 -25988.50641002807 0.6608047442568933E-12 0 0 0 3 0 0 + 3 -18.63482852291161 0.1598721155460226E-13 0 0 0 1 2 0 + 3 34.67601365859109 0.1421085471520201E-13 0 0 0 1 1 1 + 3 4.539091951909136 -0.2465190328815662E-31 0 0 0 1 0 2 + 4 2161199.086793333 0.2235174179077148E-07 3 0 1 0 0 0 + 4 -249742573.6430636 0.5308538675308228E-07 2 1 1 0 0 0 + 4 -6621147.742056847 0.1490116119384766E-07 1 2 1 0 0 0 + 4 83397443.42369542 0.1071020960807800E-07 0 3 1 0 0 0 + 4 51845.74070015462 -0.5355104804038998E-08 1 0 3 0 0 0 + 4 3524698.438742683 -0.2182787284255027E-08 0 1 3 0 0 0 + 4 -83389875.74060664 0.2980232238769530E-07 3 0 0 1 0 0 + 4 -6283499.035239162 -0.1490116119384766E-07 2 1 0 1 0 0 + 4 250868960.5927954 -0.1210719347000122E-07 1 2 0 1 0 0 + 4 2097655.653110195 -0.7450580596923830E-08 0 3 0 1 0 0 + 4 7410757.936353696 -0.2153683453798294E-08 1 0 2 1 0 0 + 4 -80923.60273867581 0.8381903171539310E-08 0 1 2 1 0 0 + 4 -1758569.999634604 0.2005720034503611E-08 2 0 1 0 1 0 + 4 -16934620.37242392 0.1587977749295533E-08 1 1 1 0 1 0 + 4 1614789.063963692 -0.1772889390849741E-08 0 2 1 0 1 0 + 4 -85931.35235113505 -0.5093170329928399E-10 0 0 3 0 1 0 + 4 -943570.2452298842 -0.1646185410209000E-09 2 0 1 0 0 1 + 4 -11520061.57085347 -0.2772139851003884E-08 1 1 1 0 0 1 + 4 693877.3798717315 0.1064108801074326E-09 0 2 1 0 0 1 + 4 -24539.74195273014 -0.5275069270282984E-10 0 0 3 0 0 1 + 4 -245332.8632128382 -0.4656612873077393E-08 1 0 1 2 0 0 + 4 -11780497.98234436 -0.3929017111659049E-08 0 1 1 2 0 0 + 4 8738500.882036250 -0.8139977580867705E-09 2 0 0 1 1 0 + 4 -3328100.960672584 0.3206309884262736E-08 1 1 0 1 1 0 + 4 -8194771.464807654 0.8576535037718712E-09 0 2 0 1 1 0 + 4 -28321.19422388317 0.8731149137020111E-10 0 0 2 1 1 0 + 4 -18686.83501247832 -0.7275957614183423E-11 1 0 1 0 2 0 + 4 -103881.3251741582 -0.3274180926382542E-10 0 1 1 0 2 0 + 4 6104587.782064813 0.1473381416872144E-08 2 0 0 1 0 1 + 4 -1778073.870574201 -0.1946318661794066E-09 1 1 0 1 0 1 + 4 -5357270.203180006 -0.1284206518903375E-08 0 2 0 1 0 1 + 4 228255.7705384988 0.2019483917365790E-26 0 0 2 1 0 1 + 4 29065.68964294368 0.9094947017729292E-12 1 0 1 0 1 1 + 4 -13653.84904963208 -0.1455191522836685E-10 0 1 1 0 1 1 + 4 21142.48973467844 0.1818989403545856E-10 1 0 1 0 0 2 + 4 34235.77016243213 -0.2910383045673371E-10 0 1 1 0 0 2 + 4 -3162161.453953072 -0.2211891114711761E-08 1 0 0 3 0 0 + 4 52928.27497186564 0.4656612873077393E-08 0 1 0 3 0 0 + 4 143974.2722831214 0.1891748979687691E-09 0 0 1 2 1 0 + 4 -43714.48833366336 -0.6184563972055911E-10 1 0 0 1 2 0 + 4 56433.61532026789 0.2182787284255028E-10 0 1 0 1 2 0 + 4 -76.38826391854967 -0.2113864638886298E-12 0 0 1 0 3 0 + 4 -176591.9094621278 0.4365574568510056E-10 0 0 1 2 0 1 + 4 -13867.55360832677 0.1009741958682895E-26 1 0 0 1 1 1 + 4 -46920.77599386395 0.6366462912410498E-11 0 1 0 1 1 1 + 4 170.7701102929656 0.2842170943040401E-13 0 0 1 0 2 1 + 4 16385.29352608069 -0.4782469237902384E-27 1 0 0 1 0 2 + 4 -58274.89382991626 -0.5456968210637569E-11 0 1 0 1 0 2 + 4 -17.08640806383911 -0.1030286966852145E-12 0 0 1 0 1 2 + 4 -26.02648711673653 0.1376676550535194E-13 0 0 1 0 0 3 + 4 127399.6670606338 0.4729372449219227E-10 0 0 0 3 1 0 + 4 24.88222272109595 0.1207922650792170E-12 0 0 0 1 3 0 + 4 113762.8098439943 0.2546585164964199E-10 0 0 0 3 0 1 + 4 154.5412960273668 -0.2415845301584341E-12 0 0 0 1 2 1 + 4 -49.44158550382673 0.3907985046680551E-13 0 0 0 1 1 2 + 4 -3.054274603088604 -0.5329070518200752E-14 0 0 0 1 0 3 + -90 0.000000000000000 0.000000000000000 0 0 0 0 0 0 + + 1, NO = 4, NV = 6, INA = 517 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 4 NV = 6 + 1 -0.2704171727972711E-01 0.000000000000000 1 0 0 0 0 0 + 1 -0.7287424721059020 0.000000000000000 0 1 0 0 0 0 + 2 -2.011826271918907 -0.4440892098500626E-15 2 0 0 0 0 0 + 2 -5.850441748540565 0.1332267629550188E-14 1 1 0 0 0 0 + 2 2.011826271918907 0.4440892098500626E-15 0 2 0 0 0 0 + 2 1.082111486715655 0.3721732664138444E-33 0 0 2 0 0 0 + 2 -23.82368803773236 0.1554312234475219E-14 0 0 1 1 0 0 + 2 -0.3237737723614176 -0.6245004513516504E-16 1 0 0 0 1 0 + 2 0.8022595317802000E-01 -0.2775557561562891E-16 0 1 0 0 1 0 + 2 0.3237743806053200 -0.6938893903907227E-17 1 0 0 0 0 1 + 2 -0.8023148227749489E-01 0.2775557561562891E-16 0 1 0 0 0 1 + 2 -1.082111486715659 0.1925929944387236E-33 0 0 0 2 0 0 + 3 -81.22478843955238 -0.3410605131648480E-12 3 0 0 0 0 0 + 3 2131.069278769467 -0.3979039320256560E-12 2 1 0 0 0 0 + 3 393.2898557162020 -0.4547473508864640E-12 1 2 0 0 0 0 + 3 -1417.004395510923 -0.2771116669464390E-12 0 3 0 0 0 0 + 3 412.6196746135161 -0.3410605131648478E-12 1 0 2 0 0 0 + 3 -1887.960602510823 0.4973799150320701E-12 0 1 2 0 0 0 + 3 -4141.800722939841 0.3410605131648481E-12 1 0 1 1 0 0 + 3 -115.1135077438383 -0.4547473508864638E-12 0 1 1 1 0 0 + 3 40.77679407567782 -0.1065814103640150E-13 2 0 0 0 1 0 + 3 -80.21498288808267 -0.6162975822039155E-30 1 1 0 0 1 0 + 3 -40.77679407567783 0.1065814103640150E-13 0 2 0 0 1 0 + 3 122.7660853811677 -0.4263256414560601E-13 0 0 2 0 1 0 + 3 -37.74824895174380 0.8881784197001251E-15 2 0 0 0 0 1 + 3 28.46117356628973 0.2131628207280301E-13 1 1 0 0 0 1 + 3 37.74824895174381 -0.8881784197001252E-15 0 2 0 0 0 1 + 3 -61.09875156714463 -0.2620126338115369E-13 0 0 2 0 0 1 + 3 -175.2370859867859 -0.7958078640513123E-12 1 0 0 2 0 0 + 3 -2083.010752391247 0.1065814103640151E-11 0 1 0 2 0 0 + 3 164.3360210870442 0.1421085471520200E-13 0 0 1 1 1 0 + 3 -0.9569166953224707 0.1155557966632342E-32 1 0 0 0 2 0 + 3 1.436953202919912 -0.5551115123125783E-16 0 1 0 0 2 0 + 3 2.144363684482438 0.1421085471520200E-13 0 0 1 1 0 1 + 3 0.8631692497710721 0.1110223024625156E-15 1 0 0 0 1 1 + 3 1.163640743902481 0.9436895709313829E-15 0 1 0 0 1 1 + 3 -0.1709693965110038 -0.6938893903907226E-16 1 0 0 0 0 2 + 3 -0.1942324124118476 0.4302114220422482E-15 0 1 0 0 0 2 + 3 18.03955507854047 0.4263256414560601E-13 0 0 0 2 1 0 + 3 154.9691785402833 0.2797762022055395E-13 0 0 0 2 0 1 + 4 -299371788.1445808 -0.4768371582031250E-06 4 0 0 0 0 0 + 4 30205366901.07848 0.2026557922363281E-05 3 1 0 0 0 0 + 4 1796241041.678272 0.3814697265625000E-05 2 2 0 0 0 0 + 4 -30205680614.00109 -0.2026557922363281E-05 1 3 0 0 0 0 + 4 -299375226.0653936 -0.4768371582031250E-06 0 4 0 0 0 0 + 4 514703.5159298737 0.1240597988809463E-24 2 0 2 0 0 0 + 4 -3673440.156612800 -0.1105945557355881E-08 1 1 2 0 0 0 + 4 422836.5178208504 -0.4656612873077392E-09 0 2 2 0 0 0 + 4 139318.9082056948 -0.3637978807091712E-10 0 0 4 0 0 0 + 4 -1820987.136308204 -0.5384208634495733E-09 2 0 1 1 0 0 + 4 -548215.2926008935 -0.4656612873077387E-09 1 1 1 1 0 0 + 4 2238014.280409546 0.7421476766467096E-09 0 2 1 1 0 0 + 4 322114.2393925341 0.2364686224609614E-10 0 0 3 1 0 0 + 4 12134.34443118172 -0.9549694368615746E-11 3 0 0 0 1 0 + 4 1912.208326848764 0.5456968210637569E-11 2 1 0 0 1 0 + 4 -24154.03061641012 0.2728484105318784E-11 1 2 0 0 1 0 + 4 -14164.31477805352 -0.5456968210637569E-11 0 3 0 0 1 0 + 4 89682.82107465870 0.2819433575496076E-10 1 0 2 0 1 0 + 4 -178278.2759530711 -0.3183231456205250E-11 0 1 2 0 1 0 + 4 -11869.88517389477 0.7958078640513120E-12 3 0 0 0 0 1 + 4 1613.994287666039 0.1273292582482099E-10 2 1 0 0 0 1 + 4 24809.05192357836 -0.8185452315956353E-11 1 2 0 0 0 1 + 4 7385.046435111839 -0.2728484105318785E-11 0 3 0 0 0 1 + 4 -102597.1432554409 0.3851852170555502E-10 1 0 2 0 0 1 + 4 65553.55932375423 -0.6684786058031021E-10 0 1 2 0 0 1 + 4 428911.4974620307 0.2328306436538695E-09 2 0 0 2 0 0 + 4 2707638.363060511 0.6984919309616089E-09 1 1 0 2 0 0 + 4 542842.6051255808 0.2328306436538697E-09 0 2 0 2 0 0 + 4 311856.2420505609 0.1746229827404022E-09 0 0 2 2 0 0 + 4 241099.4165344072 -0.5002220859751104E-10 1 0 1 1 1 0 + 4 177619.1203758581 0.3092281986027956E-10 0 1 1 1 1 0 + 4 274.7879687966769 0.1136868377216160E-12 2 0 0 0 2 0 + 4 -2229.796745815554 -0.2842170943040400E-12 1 1 0 0 2 0 + 4 -274.7879687966768 -0.1136868377216160E-12 0 2 0 0 2 0 + 4 -2601.158149148785 0.1136868377216160E-11 0 0 2 0 2 0 + 4 -87027.41292879370 0.1427906681783497E-09 1 0 1 1 0 1 + 4 -252743.5341249776 0.1121378545576591E-09 0 1 1 1 0 1 + 4 -118.9998622709010 0.1421085471520201E-13 2 0 0 0 1 1 + 4 -295.0047017981687 -0.2415845301584341E-12 1 1 0 0 1 1 + 4 118.9998622709009 -0.1421085471520200E-13 0 2 0 0 1 1 + 4 692.5871116357573 0.3126388037344443E-12 0 0 2 0 1 1 + 4 -187.4896151356416 0.5684341886080800E-13 2 0 0 0 0 2 + 4 290.4636487662223 0.5684341886080799E-13 1 1 0 0 0 2 + 4 187.4896151356416 -0.5684341886080801E-13 0 2 0 0 0 2 + 4 1498.287351466294 -0.1136868377216161E-12 0 0 2 0 0 2 + 4 -344720.6049563764 -0.2455635694786908E-10 0 0 1 3 0 0 + 4 -97084.77567208475 -0.1364242052659392E-10 1 0 0 2 1 0 + 4 85856.36389682170 -0.4047251422889531E-10 0 1 0 2 1 0 + 4 -5925.202327514604 -0.4547473508864640E-12 0 0 1 1 2 0 + 4 -7.603354265950462 -0.2850376317693109E-31 1 0 0 0 3 0 + 4 6.222968808424137 0.1848892746611746E-31 0 1 0 0 3 0 + 4 125419.4862665349 -0.3851852170555503E-10 1 0 0 2 0 1 + 4 -21239.07812933721 0.7412381819449365E-10 0 1 0 2 0 1 + 4 772.3032880048048 -0.2273736754432320E-11 0 0 1 1 1 1 + 4 0.9740684712787906 0.8881784197001248E-15 1 0 0 0 2 1 + 4 9.330822357796725 0.3330669073875468E-15 0 1 0 0 2 1 + 4 1333.806596515065 -0.6821210263296964E-12 0 0 1 1 0 2 + 4 4.118980604045873 0.2553512956637859E-14 1 0 0 0 1 2 + 4 -1.424573461467301 0.8881784197001251E-15 0 1 0 0 1 2 + 4 0.2439806795698564 0.3139849491518020E-15 1 0 0 0 0 3 + 4 0.1439728766552144E-01 -0.3608224830031759E-15 0 1 0 0 0 3 + 4 148407.1133464848 -0.2182787284255027E-10 0 0 0 4 0 0 + 4 1957.695819663953 -0.1136868377216160E-11 0 0 0 2 2 0 + 4 -1777.465774420309 -0.3126388037344440E-12 0 0 0 2 1 1 + 4 -1934.347675578181 0.1136868377216160E-12 0 0 0 2 0 2 + -105 0.000000000000000 0.000000000000000 0 0 0 0 0 0 + + 1, NO = 4, NV = 6, INA = 518 + ********************************************* + + I COEFFICIENT ORDER EXPONENTS + NO = 4 NV = 6 + 1 -5.852478004683411 0.000000000000000 0 0 1 0 0 0 + 1 -0.1025476998463437E-02 0.000000000000000 0 0 0 1 0 0 + 2 35.76993929294548 0.2967417978005926E-13 1 0 1 0 0 0 + 2 -2.424757139212866 -0.1065814103640150E-13 0 1 1 0 0 0 + 2 2.547762920536613 -0.1065814103640150E-13 1 0 0 1 0 0 + 2 -34.55658049676505 -0.2168057400275814E-13 0 1 0 1 0 0 + 2 14.70856647275444 -0.1776356839400250E-14 0 0 1 0 1 0 + 2 10.79584271946479 0.000000000000000 0 0 1 0 0 1 + 2 12.85892134338537 -0.2664535259100376E-14 0 0 0 1 1 0 + 2 3.488645572261085 0.000000000000000 0 0 0 1 0 1 + 3 -45665.50585706252 -0.1533706495138176E-10 2 0 1 0 0 0 + 3 34732.25915618958 -0.2715694336075103E-10 1 1 1 0 0 0 + 3 209638.9191010812 -0.3243627588744858E-11 0 2 1 0 0 0 + 3 23357.91314319626 -0.8629985615016267E-11 0 0 3 0 0 0 + 3 -17346.02720287200 0.4462918923309189E-10 2 0 0 1 0 0 + 3 -254926.7568701422 -0.9734435479913373E-11 1 1 0 1 0 0 + 3 17459.62334121750 0.1357847168037551E-10 0 2 0 1 0 0 + 3 -17.52764911710167 0.2415845301584343E-12 0 0 2 1 0 0 + 3 -4405.184626316654 -0.2728484105318785E-11 1 0 1 0 1 0 + 3 976.6204991486399 0.1364242052659392E-11 0 1 1 0 1 0 + 3 2113.683511346702 0.4547473508864641E-12 1 0 1 0 0 1 + 3 -1714.546396836433 0.5684341886080801E-12 0 1 1 0 0 1 + 3 23022.60973975051 -0.8796519068710039E-11 0 0 1 2 0 0 + 3 3117.446726500209 0.1364242052659392E-11 1 0 0 1 1 0 + 3 2029.615361155712 0.3296918293926865E-11 0 1 0 1 1 0 + 3 -49.66265368878263 0.1421085471520200E-13 0 0 1 0 2 0 + 3 -4449.430021605681 0.1136868377216160E-12 1 0 0 1 0 1 + 3 -1078.385998462077 -0.1364242052659392E-11 0 1 0 1 0 1 + 3 -15.87535114854447 0.4440892098500627E-13 0 0 1 0 1 1 + 3 -4.042781333651214 0.000000000000000 0 0 1 0 0 2 + 3 -40.49505641309520 -0.3105071755271639E-11 0 0 0 3 0 0 + 3 -16.56225019911051 -0.3197442310920451E-13 0 0 0 1 2 0 + 3 9.955131866807788 -0.2398081733190338E-13 0 0 0 1 1 1 + 3 8.902986283730227 0.000000000000000 0 0 0 1 0 2 + 4 -80430315.80351658 -0.1234002411365509E-07 3 0 1 0 0 0 + 4 -6228924.976235100 0.1135454832538916E-23 2 1 1 0 0 0 + 4 242235563.9549994 -0.4749745130538938E-07 1 2 1 0 0 0 + 4 2026816.825423369 0.2235174179077149E-07 0 3 1 0 0 0 + 4 3201338.627144473 0.5307811079546811E-08 1 0 3 0 0 0 + 4 -44612.23515251315 -0.6752088665962216E-08 0 1 3 0 0 0 + 4 -2079871.960902808 0.7450580596923825E-08 3 0 0 1 0 0 + 4 241646539.5278362 -0.5587935447692870E-08 2 1 0 1 0 0 + 4 6324626.764460130 0.2980232238769531E-07 1 2 0 1 0 0 + 4 -80486754.47939509 0.2887099981307984E-07 0 3 0 1 0 0 + 4 -70194.65217523211 -0.4656612873077396E-08 1 0 2 1 0 0 + 4 -5146184.416501270 -0.4380126483738422E-08 0 1 2 1 0 0 + 4 8009018.427607195 -0.1028183760354295E-08 2 0 1 0 1 0 + 4 -3190294.362352938 0.1993612386286259E-08 1 1 1 0 1 0 + 4 -9075815.334812358 0.1290118234464899E-08 0 2 1 0 1 0 + 4 -269498.1224936843 0.5820766091346743E-10 0 0 3 0 1 0 + 4 5550808.884613156 0.1299895302508958E-08 2 0 1 0 0 1 + 4 -1800517.494292025 0.1095031620934605E-08 1 1 1 0 0 1 + 4 -5753459.009901837 -0.8924416761146859E-09 0 2 1 0 0 1 + 4 -75915.89354105282 0.1018634065985680E-09 0 0 3 0 0 1 + 4 -5946529.188850471 0.5726178642362358E-08 1 0 1 2 0 0 + 4 117770.3914731949 -0.1210719347000122E-07 0 1 1 2 0 0 + 4 1439502.837945917 -0.7494236342608924E-09 2 0 0 1 1 0 + 4 17328633.52131001 -0.2338310878258199E-08 1 1 0 1 1 0 + 4 -1733529.966056132 0.5165929906070232E-09 0 2 0 1 1 0 + 4 137429.9765117652 0.7275957614183425E-10 0 0 2 1 1 0 + 4 122715.2111509288 0.2910383045673372E-10 1 0 1 0 2 0 + 4 1398.967506768271 -0.7275957614183428E-11 0 1 1 0 2 0 + 4 1037249.229686596 -0.6602931534871457E-09 2 0 0 1 0 1 + 4 11178213.66830301 0.2173237589886412E-08 1 1 0 1 0 1 + 4 -767876.7020373951 0.5729816621169449E-09 0 2 0 1 0 1 + 4 -593400.9227058081 -0.5093170329928400E-10 0 0 2 1 0 1 + 4 18187.55517619968 0.2182787284255028E-10 1 0 1 0 1 1 + 4 17864.36322137485 0.4547473508864638E-11 0 1 1 0 1 1 + 4 -31506.65651469816 0.2910383045673369E-10 1 0 1 0 0 2 + 4 1232.864924475040 0.1637090463191271E-10 0 1 1 0 0 2 + 4 13295.73789018139 -0.7334165275096896E-08 1 0 0 3 0 0 + 4 2330589.084725094 -0.4918547347187997E-08 0 1 0 3 0 0 + 4 -8718.243376039998 0.1455191522836685E-09 0 0 1 2 1 0 + 4 -76866.13031311180 -0.2182787284255027E-10 1 0 0 1 2 0 + 4 -18582.50500963759 -0.8003553375601768E-10 0 1 0 1 2 0 + 4 180.7010058020620 -0.7105427357601002E-13 0 0 1 0 3 0 + 4 84765.49263120226 0.1164153218269348E-09 0 0 1 2 0 1 + 4 66588.89117705819 -0.4547473508864642E-11 1 0 0 1 1 1 + 4 -13278.39721702469 0.5143373102040997E-27 0 1 0 1 1 1 + 4 384.1717683157476 0.1607602939657227E-12 0 0 1 0 2 1 + 4 71345.00226059923 0.7275957614183428E-11 1 0 0 1 0 2 + 4 6653.218833778099 -0.1413638742156053E-26 0 1 0 1 0 2 + 4 256.9703941950236 -0.7105427357601003E-13 0 0 1 0 1 2 + 4 73.96957502838900 -0.2131628207280301E-13 0 0 1 0 0 3 + 4 -130979.8174163291 -0.8367351256310941E-10 0 0 0 3 1 0 + 4 253.8950358426778 0.8526512829121202E-13 0 0 0 1 3 0 + 4 248482.1077698862 -0.1091393642127514E-09 0 0 0 3 0 1 + 4 -2.034203798243036 0.1084683744678891E-28 0 0 0 1 2 1 + 4 159.9449552419620 0.1989519660128281E-12 0 0 0 1 1 2 + 4 67.57869690745196 0.2131628207280301E-13 0 0 0 1 0 3 + -90 0.000000000000000 0.000000000000000 0 0 0 0 0 0 + No Stochastic Radiation \ No newline at end of file