From 00c2543c7b5edf476cfeb461d6583f9b6574a078 Mon Sep 17 00:00:00 2001 From: Don Cross Date: Wed, 3 Jul 2019 20:36:20 -0400 Subject: [PATCH] Python: AddSol translation gets time down to 49 seconds. --- generate/fixaddsol.py | 30 +- generate/template/astronomy.py | 1082 +++++++++++++++++-- generate/template/{astronomy.old => old.py} | 10 +- generate/test.py | 0 source/python/astronomy.py | 1082 +++++++++++++++++-- 5 files changed, 1991 insertions(+), 213 deletions(-) rename generate/template/{astronomy.old => old.py} (99%) mode change 100644 => 100755 generate/test.py diff --git a/generate/fixaddsol.py b/generate/fixaddsol.py index 4ba1e865..cf90d6a9 100755 --- a/generate/fixaddsol.py +++ b/generate/fixaddsol.py @@ -5,7 +5,33 @@ import re def Translate(line): m = re.match(r'^ AddSol\(\s*([0-9\+\-\.]+)\s*,\s*([0-9\+\-\.]+)\s*,\s*([0-9\+\-\.]+)\s*,\s*([0-9\+\-\.]+)\s*,\s*([\+\-]?\d+)\s*,\s*([\+\-]?\d+)\s*,\s*([\+\-]?\d+)\s*,\s*([\+\-]?\d+)\s*\)\s*$', line) if m: - return '' + cl = float(m.group(1)) + cs = float(m.group(2)) + cg = float(m.group(3)) + cp = float(m.group(4)) + p = int(m.group(5)) + q = int(m.group(6)) + r = int(m.group(7)) + s = int(m.group(8)) + + text = '\n' + text += ' # AddSol({}, {}, {}, {}, {}, {}, {}, {})\n'.format(cl, cs, cg, cp, p, q, r, s) + text += ' (a, b) = (1, 0)\n' + if p != 0: + text += ' (a, b) = AddThe(a, b, co[{0:d}][1], si[{0:d}][1])\n'.format(p) + if q != 0: + text += ' (a, b) = AddThe(a, b, co[{0:d}][2], si[{0:d}][2])\n'.format(q) + if r != 0: + text += ' (a, b) = AddThe(a, b, co[{0:d}][3], si[{0:d}][3])\n'.format(r) + if s != 0: + text += ' (a, b) = AddThe(a, b, co[{0:d}][4], si[{0:d}][4])\n'.format(s) + + text += ' DLAM += {} * b\n'.format(cl) + text += ' DS += {} * b\n'.format(cs) + text += ' GAM1C += {} * a\n'.format(cg) + text += ' SINPI += {} * a\n'.format(cp) + + return text return line def FixAddSol(inFileName, outFileName): @@ -18,5 +44,5 @@ def FixAddSol(inFileName, outFileName): outfile.write(xlat) if __name__ == '__main__': - FixAddSol('template/astronomy.old', 'template/astronomy.py') + FixAddSol('template/old.py', 'template/astronomy.py') sys.exit(0) diff --git a/generate/template/astronomy.py b/generate/template/astronomy.py index 720771e3..7d6f1c2b 100644 --- a/generate/template/astronomy.py +++ b/generate/template/astronomy.py @@ -735,118 +735,990 @@ def _CalcMoon(time): def AddSol(coeffl, coeffs, coeffg, coeffp, p, q, r, s): nonlocal DLAM, DS, GAM1C, SINPI - (a, b) = Term(p, q, r, s) + (a, b) = (1, 0) + if p != 0: + (a, b) = AddThe(a, b, co[p][1], si[p][1]) + if q != 0: + (a, b) = AddThe(a, b, co[q][2], si[q][2]) + if r != 0: + (a, b) = AddThe(a, b, co[r][3], si[r][3]) + if s != 0: + (a, b) = AddThe(a, b, co[s][4], si[s][4]) DLAM += coeffl * b DS += coeffs * b GAM1C += coeffg * a SINPI += coeffp * a - AddSol( 13.902, 14.06,-0.001, 0.2607,0, 0, 0, 4) - AddSol( 0.403, -4.01,+0.394, 0.0023,0, 0, 0, 3) - AddSol( 2369.912, 2373.36,+0.601, 28.2333,0, 0, 0, 2) - AddSol( -125.154, -112.79,-0.725, -0.9781,0, 0, 0, 1) - AddSol( 1.979, 6.98,-0.445, 0.0433,1, 0, 0, 4) - AddSol( 191.953, 192.72,+0.029, 3.0861,1, 0, 0, 2) - AddSol( -8.466, -13.51,+0.455, -0.1093,1, 0, 0, 1) - AddSol( 22639.500,22609.07,+0.079, 186.5398,1, 0, 0, 0) - AddSol( 18.609, 3.59,-0.094, 0.0118,1, 0, 0,-1) - AddSol( -4586.465,-4578.13,-0.077, 34.3117,1, 0, 0,-2) - AddSol( +3.215, 5.44,+0.192, -0.0386,1, 0, 0,-3) - AddSol( -38.428, -38.64,+0.001, 0.6008,1, 0, 0,-4) - AddSol( -0.393, -1.43,-0.092, 0.0086,1, 0, 0,-6) - AddSol( -0.289, -1.59,+0.123, -0.0053,0, 1, 0, 4) - AddSol( -24.420, -25.10,+0.040, -0.3000,0, 1, 0, 2) - AddSol( 18.023, 17.93,+0.007, 0.1494,0, 1, 0, 1) - AddSol( -668.146, -126.98,-1.302, -0.3997,0, 1, 0, 0) - AddSol( 0.560, 0.32,-0.001, -0.0037,0, 1, 0,-1) - AddSol( -165.145, -165.06,+0.054, 1.9178,0, 1, 0,-2) - AddSol( -1.877, -6.46,-0.416, 0.0339,0, 1, 0,-4) - AddSol( 0.213, 1.02,-0.074, 0.0054,2, 0, 0, 4) - AddSol( 14.387, 14.78,-0.017, 0.2833,2, 0, 0, 2) - AddSol( -0.586, -1.20,+0.054, -0.0100,2, 0, 0, 1) - AddSol( 769.016, 767.96,+0.107, 10.1657,2, 0, 0, 0) - AddSol( +1.750, 2.01,-0.018, 0.0155,2, 0, 0,-1) - AddSol( -211.656, -152.53,+5.679, -0.3039,2, 0, 0,-2) - AddSol( +1.225, 0.91,-0.030, -0.0088,2, 0, 0,-3) - AddSol( -30.773, -34.07,-0.308, 0.3722,2, 0, 0,-4) - AddSol( -0.570, -1.40,-0.074, 0.0109,2, 0, 0,-6) - AddSol( -2.921, -11.75,+0.787, -0.0484,1, 1, 0, 2) - AddSol( +1.267, 1.52,-0.022, 0.0164,1, 1, 0, 1) - AddSol( -109.673, -115.18,+0.461, -0.9490,1, 1, 0, 0) - AddSol( -205.962, -182.36,+2.056, +1.4437,1, 1, 0,-2) - AddSol( 0.233, 0.36, 0.012, -0.0025,1, 1, 0,-3) - AddSol( -4.391, -9.66,-0.471, 0.0673,1, 1, 0,-4) - AddSol( 0.283, 1.53,-0.111, +0.0060,1,-1, 0,+4) - AddSol( 14.577, 31.70,-1.540, +0.2302,1,-1, 0, 2) - AddSol( 147.687, 138.76,+0.679, +1.1528,1,-1, 0, 0) - AddSol( -1.089, 0.55,+0.021, 0.0 ,1,-1, 0,-1) - AddSol( 28.475, 23.59,-0.443, -0.2257,1,-1, 0,-2) - AddSol( -0.276, -0.38,-0.006, -0.0036,1,-1, 0,-3) - AddSol( 0.636, 2.27,+0.146, -0.0102,1,-1, 0,-4) - AddSol( -0.189, -1.68,+0.131, -0.0028,0, 2, 0, 2) - AddSol( -7.486, -0.66,-0.037, -0.0086,0, 2, 0, 0) - AddSol( -8.096, -16.35,-0.740, 0.0918,0, 2, 0,-2) - AddSol( -5.741, -0.04, 0.0 , -0.0009,0, 0, 2, 2) - AddSol( 0.255, 0.0 , 0.0 , 0.0 ,0, 0, 2, 1) - AddSol( -411.608, -0.20, 0.0 , -0.0124,0, 0, 2, 0) - AddSol( 0.584, 0.84, 0.0 , +0.0071,0, 0, 2,-1) - AddSol( -55.173, -52.14, 0.0 , -0.1052,0, 0, 2,-2) - AddSol( 0.254, 0.25, 0.0 , -0.0017,0, 0, 2,-3) - AddSol( +0.025, -1.67, 0.0 , +0.0031,0, 0, 2,-4) - AddSol( 1.060, 2.96,-0.166, 0.0243,3, 0, 0,+2) - AddSol( 36.124, 50.64,-1.300, 0.6215,3, 0, 0, 0) - AddSol( -13.193, -16.40,+0.258, -0.1187,3, 0, 0,-2) - AddSol( -1.187, -0.74,+0.042, 0.0074,3, 0, 0,-4) - AddSol( -0.293, -0.31,-0.002, 0.0046,3, 0, 0,-6) - AddSol( -0.290, -1.45,+0.116, -0.0051,2, 1, 0, 2) - AddSol( -7.649, -10.56,+0.259, -0.1038,2, 1, 0, 0) - AddSol( -8.627, -7.59,+0.078, -0.0192,2, 1, 0,-2) - AddSol( -2.740, -2.54,+0.022, 0.0324,2, 1, 0,-4) - AddSol( 1.181, 3.32,-0.212, 0.0213,2,-1, 0,+2) - AddSol( 9.703, 11.67,-0.151, 0.1268,2,-1, 0, 0) - AddSol( -0.352, -0.37,+0.001, -0.0028,2,-1, 0,-1) - AddSol( -2.494, -1.17,-0.003, -0.0017,2,-1, 0,-2) - AddSol( 0.360, 0.20,-0.012, -0.0043,2,-1, 0,-4) - AddSol( -1.167, -1.25,+0.008, -0.0106,1, 2, 0, 0) - AddSol( -7.412, -6.12,+0.117, 0.0484,1, 2, 0,-2) - AddSol( -0.311, -0.65,-0.032, 0.0044,1, 2, 0,-4) - AddSol( +0.757, 1.82,-0.105, 0.0112,1,-2, 0, 2) - AddSol( +2.580, 2.32,+0.027, 0.0196,1,-2, 0, 0) - AddSol( +2.533, 2.40,-0.014, -0.0212,1,-2, 0,-2) - AddSol( -0.344, -0.57,-0.025, +0.0036,0, 3, 0,-2) - AddSol( -0.992, -0.02, 0.0 , 0.0 ,1, 0, 2, 2) - AddSol( -45.099, -0.02, 0.0 , -0.0010,1, 0, 2, 0) - AddSol( -0.179, -9.52, 0.0 , -0.0833,1, 0, 2,-2) - AddSol( -0.301, -0.33, 0.0 , 0.0014,1, 0, 2,-4) - AddSol( -6.382, -3.37, 0.0 , -0.0481,1, 0,-2, 2) - AddSol( 39.528, 85.13, 0.0 , -0.7136,1, 0,-2, 0) - AddSol( 9.366, 0.71, 0.0 , -0.0112,1, 0,-2,-2) - AddSol( 0.202, 0.02, 0.0 , 0.0 ,1, 0,-2,-4) + # AddSol(13.902, 14.06, -0.001, 0.2607, 0, 0, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += 13.902 * b + DS += 14.06 * b + GAM1C += -0.001 * a + SINPI += 0.2607 * a - AddSol( 0.415, 0.10, 0.0 , 0.0013,0, 1, 2, 0) - AddSol( -2.152, -2.26, 0.0 , -0.0066,0, 1, 2,-2) - AddSol( -1.440, -1.30, 0.0 , +0.0014,0, 1,-2, 2) - AddSol( 0.384, -0.04, 0.0 , 0.0 ,0, 1,-2,-2) - AddSol( +1.938, +3.60,-0.145, +0.0401,4, 0, 0, 0) - AddSol( -0.952, -1.58,+0.052, -0.0130,4, 0, 0,-2) - AddSol( -0.551, -0.94,+0.032, -0.0097,3, 1, 0, 0) - AddSol( -0.482, -0.57,+0.005, -0.0045,3, 1, 0,-2) - AddSol( 0.681, 0.96,-0.026, 0.0115,3,-1, 0, 0) - AddSol( -0.297, -0.27, 0.002, -0.0009,2, 2, 0,-2) - AddSol( 0.254, +0.21,-0.003, 0.0 ,2,-2, 0,-2) - AddSol( -0.250, -0.22, 0.004, 0.0014,1, 3, 0,-2) - AddSol( -3.996, 0.0 , 0.0 , +0.0004,2, 0, 2, 0) - AddSol( 0.557, -0.75, 0.0 , -0.0090,2, 0, 2,-2) - AddSol( -0.459, -0.38, 0.0 , -0.0053,2, 0,-2, 2) - AddSol( -1.298, 0.74, 0.0 , +0.0004,2, 0,-2, 0) - AddSol( 0.538, 1.14, 0.0 , -0.0141,2, 0,-2,-2) - AddSol( 0.263, 0.02, 0.0 , 0.0 ,1, 1, 2, 0) - AddSol( 0.426, +0.07, 0.0 , -0.0006,1, 1,-2,-2) - AddSol( -0.304, +0.03, 0.0 , +0.0003,1,-1, 2, 0) - AddSol( -0.372, -0.19, 0.0 , -0.0027,1,-1,-2, 2) - AddSol( +0.418, 0.0 , 0.0 , 0.0 ,0, 0, 4, 0) - AddSol( -0.330, -0.04, 0.0 , 0.0 ,3, 0, 2, 0) + # AddSol(0.403, -4.01, 0.394, 0.0023, 0, 0, 0, 3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][4], si[3][4]) + DLAM += 0.403 * b + DS += -4.01 * b + GAM1C += 0.394 * a + SINPI += 0.0023 * a + + # AddSol(2369.912, 2373.36, 0.601, 28.2333, 0, 0, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 2369.912 * b + DS += 2373.36 * b + GAM1C += 0.601 * a + SINPI += 28.2333 * a + + # AddSol(-125.154, -112.79, -0.725, -0.9781, 0, 0, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += -125.154 * b + DS += -112.79 * b + GAM1C += -0.725 * a + SINPI += -0.9781 * a + + # AddSol(1.979, 6.98, -0.445, 0.0433, 1, 0, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += 1.979 * b + DS += 6.98 * b + GAM1C += -0.445 * a + SINPI += 0.0433 * a + + # AddSol(191.953, 192.72, 0.029, 3.0861, 1, 0, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 191.953 * b + DS += 192.72 * b + GAM1C += 0.029 * a + SINPI += 3.0861 * a + + # AddSol(-8.466, -13.51, 0.455, -0.1093, 1, 0, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += -8.466 * b + DS += -13.51 * b + GAM1C += 0.455 * a + SINPI += -0.1093 * a + + # AddSol(22639.5, 22609.07, 0.079, 186.5398, 1, 0, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + DLAM += 22639.5 * b + DS += 22609.07 * b + GAM1C += 0.079 * a + SINPI += 186.5398 * a + + # AddSol(18.609, 3.59, -0.094, 0.0118, 1, 0, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += 18.609 * b + DS += 3.59 * b + GAM1C += -0.094 * a + SINPI += 0.0118 * a + + # AddSol(-4586.465, -4578.13, -0.077, 34.3117, 1, 0, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -4586.465 * b + DS += -4578.13 * b + GAM1C += -0.077 * a + SINPI += 34.3117 * a + + # AddSol(3.215, 5.44, 0.192, -0.0386, 1, 0, 0, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += 3.215 * b + DS += 5.44 * b + GAM1C += 0.192 * a + SINPI += -0.0386 * a + + # AddSol(-38.428, -38.64, 0.001, 0.6008, 1, 0, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -38.428 * b + DS += -38.64 * b + GAM1C += 0.001 * a + SINPI += 0.6008 * a + + # AddSol(-0.393, -1.43, -0.092, 0.0086, 1, 0, 0, -6) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-6][4], si[-6][4]) + DLAM += -0.393 * b + DS += -1.43 * b + GAM1C += -0.092 * a + SINPI += 0.0086 * a + + # AddSol(-0.289, -1.59, 0.123, -0.0053, 0, 1, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += -0.289 * b + DS += -1.59 * b + GAM1C += 0.123 * a + SINPI += -0.0053 * a + + # AddSol(-24.42, -25.1, 0.04, -0.3, 0, 1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -24.42 * b + DS += -25.1 * b + GAM1C += 0.04 * a + SINPI += -0.3 * a + + # AddSol(18.023, 17.93, 0.007, 0.1494, 0, 1, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += 18.023 * b + DS += 17.93 * b + GAM1C += 0.007 * a + SINPI += 0.1494 * a + + # AddSol(-668.146, -126.98, -1.302, -0.3997, 0, 1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + DLAM += -668.146 * b + DS += -126.98 * b + GAM1C += -1.302 * a + SINPI += -0.3997 * a + + # AddSol(0.56, 0.32, -0.001, -0.0037, 0, 1, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += 0.56 * b + DS += 0.32 * b + GAM1C += -0.001 * a + SINPI += -0.0037 * a + + # AddSol(-165.145, -165.06, 0.054, 1.9178, 0, 1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -165.145 * b + DS += -165.06 * b + GAM1C += 0.054 * a + SINPI += 1.9178 * a + + # AddSol(-1.877, -6.46, -0.416, 0.0339, 0, 1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -1.877 * b + DS += -6.46 * b + GAM1C += -0.416 * a + SINPI += 0.0339 * a + + # AddSol(0.213, 1.02, -0.074, 0.0054, 2, 0, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += 0.213 * b + DS += 1.02 * b + GAM1C += -0.074 * a + SINPI += 0.0054 * a + + # AddSol(14.387, 14.78, -0.017, 0.2833, 2, 0, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 14.387 * b + DS += 14.78 * b + GAM1C += -0.017 * a + SINPI += 0.2833 * a + + # AddSol(-0.586, -1.2, 0.054, -0.01, 2, 0, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += -0.586 * b + DS += -1.2 * b + GAM1C += 0.054 * a + SINPI += -0.01 * a + + # AddSol(769.016, 767.96, 0.107, 10.1657, 2, 0, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + DLAM += 769.016 * b + DS += 767.96 * b + GAM1C += 0.107 * a + SINPI += 10.1657 * a + + # AddSol(1.75, 2.01, -0.018, 0.0155, 2, 0, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += 1.75 * b + DS += 2.01 * b + GAM1C += -0.018 * a + SINPI += 0.0155 * a + + # AddSol(-211.656, -152.53, 5.679, -0.3039, 2, 0, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -211.656 * b + DS += -152.53 * b + GAM1C += 5.679 * a + SINPI += -0.3039 * a + + # AddSol(1.225, 0.91, -0.03, -0.0088, 2, 0, 0, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += 1.225 * b + DS += 0.91 * b + GAM1C += -0.03 * a + SINPI += -0.0088 * a + + # AddSol(-30.773, -34.07, -0.308, 0.3722, 2, 0, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -30.773 * b + DS += -34.07 * b + GAM1C += -0.308 * a + SINPI += 0.3722 * a + + # AddSol(-0.57, -1.4, -0.074, 0.0109, 2, 0, 0, -6) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-6][4], si[-6][4]) + DLAM += -0.57 * b + DS += -1.4 * b + GAM1C += -0.074 * a + SINPI += 0.0109 * a + + # AddSol(-2.921, -11.75, 0.787, -0.0484, 1, 1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -2.921 * b + DS += -11.75 * b + GAM1C += 0.787 * a + SINPI += -0.0484 * a + + # AddSol(1.267, 1.52, -0.022, 0.0164, 1, 1, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += 1.267 * b + DS += 1.52 * b + GAM1C += -0.022 * a + SINPI += 0.0164 * a + + # AddSol(-109.673, -115.18, 0.461, -0.949, 1, 1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + DLAM += -109.673 * b + DS += -115.18 * b + GAM1C += 0.461 * a + SINPI += -0.949 * a + + # AddSol(-205.962, -182.36, 2.056, 1.4437, 1, 1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -205.962 * b + DS += -182.36 * b + GAM1C += 2.056 * a + SINPI += 1.4437 * a + + # AddSol(0.233, 0.36, 0.012, -0.0025, 1, 1, 0, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += 0.233 * b + DS += 0.36 * b + GAM1C += 0.012 * a + SINPI += -0.0025 * a + + # AddSol(-4.391, -9.66, -0.471, 0.0673, 1, 1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -4.391 * b + DS += -9.66 * b + GAM1C += -0.471 * a + SINPI += 0.0673 * a + + + # AddSol(0.283, 1.53, -0.111, 0.006, 1, -1, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += 0.283 * b + DS += 1.53 * b + GAM1C += -0.111 * a + SINPI += 0.006 * a + + # AddSol(14.577, 31.7, -1.54, 0.2302, 1, -1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 14.577 * b + DS += 31.7 * b + GAM1C += -1.54 * a + SINPI += 0.2302 * a + + # AddSol(147.687, 138.76, 0.679, 1.1528, 1, -1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + DLAM += 147.687 * b + DS += 138.76 * b + GAM1C += 0.679 * a + SINPI += 1.1528 * a + + # AddSol(-1.089, 0.55, 0.021, 0.0, 1, -1, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += -1.089 * b + DS += 0.55 * b + GAM1C += 0.021 * a + SINPI += 0.0 * a + + # AddSol(28.475, 23.59, -0.443, -0.2257, 1, -1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 28.475 * b + DS += 23.59 * b + GAM1C += -0.443 * a + SINPI += -0.2257 * a + + # AddSol(-0.276, -0.38, -0.006, -0.0036, 1, -1, 0, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += -0.276 * b + DS += -0.38 * b + GAM1C += -0.006 * a + SINPI += -0.0036 * a + + # AddSol(0.636, 2.27, 0.146, -0.0102, 1, -1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += 0.636 * b + DS += 2.27 * b + GAM1C += 0.146 * a + SINPI += -0.0102 * a + + # AddSol(-0.189, -1.68, 0.131, -0.0028, 0, 2, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.189 * b + DS += -1.68 * b + GAM1C += 0.131 * a + SINPI += -0.0028 * a + + # AddSol(-7.486, -0.66, -0.037, -0.0086, 0, 2, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + DLAM += -7.486 * b + DS += -0.66 * b + GAM1C += -0.037 * a + SINPI += -0.0086 * a + + # AddSol(-8.096, -16.35, -0.74, 0.0918, 0, 2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -8.096 * b + DS += -16.35 * b + GAM1C += -0.74 * a + SINPI += 0.0918 * a + + # AddSol(-5.741, -0.04, 0.0, -0.0009, 0, 0, 2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -5.741 * b + DS += -0.04 * b + GAM1C += 0.0 * a + SINPI += -0.0009 * a + + # AddSol(0.255, 0.0, 0.0, 0.0, 0, 0, 2, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += 0.255 * b + DS += 0.0 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(-411.608, -0.2, 0.0, -0.0124, 0, 0, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -411.608 * b + DS += -0.2 * b + GAM1C += 0.0 * a + SINPI += -0.0124 * a + + # AddSol(0.584, 0.84, 0.0, 0.0071, 0, 0, 2, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += 0.584 * b + DS += 0.84 * b + GAM1C += 0.0 * a + SINPI += 0.0071 * a + + # AddSol(-55.173, -52.14, 0.0, -0.1052, 0, 0, 2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -55.173 * b + DS += -52.14 * b + GAM1C += 0.0 * a + SINPI += -0.1052 * a + + # AddSol(0.254, 0.25, 0.0, -0.0017, 0, 0, 2, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += 0.254 * b + DS += 0.25 * b + GAM1C += 0.0 * a + SINPI += -0.0017 * a + + # AddSol(0.025, -1.67, 0.0, 0.0031, 0, 0, 2, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += 0.025 * b + DS += -1.67 * b + GAM1C += 0.0 * a + SINPI += 0.0031 * a + + # AddSol(1.06, 2.96, -0.166, 0.0243, 3, 0, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 1.06 * b + DS += 2.96 * b + GAM1C += -0.166 * a + SINPI += 0.0243 * a + + # AddSol(36.124, 50.64, -1.3, 0.6215, 3, 0, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + DLAM += 36.124 * b + DS += 50.64 * b + GAM1C += -1.3 * a + SINPI += 0.6215 * a + + # AddSol(-13.193, -16.4, 0.258, -0.1187, 3, 0, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -13.193 * b + DS += -16.4 * b + GAM1C += 0.258 * a + SINPI += -0.1187 * a + + # AddSol(-1.187, -0.74, 0.042, 0.0074, 3, 0, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -1.187 * b + DS += -0.74 * b + GAM1C += 0.042 * a + SINPI += 0.0074 * a + + # AddSol(-0.293, -0.31, -0.002, 0.0046, 3, 0, 0, -6) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[-6][4], si[-6][4]) + DLAM += -0.293 * b + DS += -0.31 * b + GAM1C += -0.002 * a + SINPI += 0.0046 * a + + # AddSol(-0.29, -1.45, 0.116, -0.0051, 2, 1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.29 * b + DS += -1.45 * b + GAM1C += 0.116 * a + SINPI += -0.0051 * a + + # AddSol(-7.649, -10.56, 0.259, -0.1038, 2, 1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + DLAM += -7.649 * b + DS += -10.56 * b + GAM1C += 0.259 * a + SINPI += -0.1038 * a + + # AddSol(-8.627, -7.59, 0.078, -0.0192, 2, 1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -8.627 * b + DS += -7.59 * b + GAM1C += 0.078 * a + SINPI += -0.0192 * a + + # AddSol(-2.74, -2.54, 0.022, 0.0324, 2, 1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -2.74 * b + DS += -2.54 * b + GAM1C += 0.022 * a + SINPI += 0.0324 * a + + # AddSol(1.181, 3.32, -0.212, 0.0213, 2, -1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 1.181 * b + DS += 3.32 * b + GAM1C += -0.212 * a + SINPI += 0.0213 * a + + # AddSol(9.703, 11.67, -0.151, 0.1268, 2, -1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + DLAM += 9.703 * b + DS += 11.67 * b + GAM1C += -0.151 * a + SINPI += 0.1268 * a + + # AddSol(-0.352, -0.37, 0.001, -0.0028, 2, -1, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += -0.352 * b + DS += -0.37 * b + GAM1C += 0.001 * a + SINPI += -0.0028 * a + + # AddSol(-2.494, -1.17, -0.003, -0.0017, 2, -1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -2.494 * b + DS += -1.17 * b + GAM1C += -0.003 * a + SINPI += -0.0017 * a + + # AddSol(0.36, 0.2, -0.012, -0.0043, 2, -1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += 0.36 * b + DS += 0.2 * b + GAM1C += -0.012 * a + SINPI += -0.0043 * a + + # AddSol(-1.167, -1.25, 0.008, -0.0106, 1, 2, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + DLAM += -1.167 * b + DS += -1.25 * b + GAM1C += 0.008 * a + SINPI += -0.0106 * a + + # AddSol(-7.412, -6.12, 0.117, 0.0484, 1, 2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -7.412 * b + DS += -6.12 * b + GAM1C += 0.117 * a + SINPI += 0.0484 * a + + # AddSol(-0.311, -0.65, -0.032, 0.0044, 1, 2, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -0.311 * b + DS += -0.65 * b + GAM1C += -0.032 * a + SINPI += 0.0044 * a + + # AddSol(0.757, 1.82, -0.105, 0.0112, 1, -2, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][2], si[-2][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 0.757 * b + DS += 1.82 * b + GAM1C += -0.105 * a + SINPI += 0.0112 * a + + # AddSol(2.58, 2.32, 0.027, 0.0196, 1, -2, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][2], si[-2][2]) + DLAM += 2.58 * b + DS += 2.32 * b + GAM1C += 0.027 * a + SINPI += 0.0196 * a + + # AddSol(2.533, 2.4, -0.014, -0.0212, 1, -2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][2], si[-2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 2.533 * b + DS += 2.4 * b + GAM1C += -0.014 * a + SINPI += -0.0212 * a + + # AddSol(-0.344, -0.57, -0.025, 0.0036, 0, 3, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][2], si[3][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.344 * b + DS += -0.57 * b + GAM1C += -0.025 * a + SINPI += 0.0036 * a + + # AddSol(-0.992, -0.02, 0.0, 0.0, 1, 0, 2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.992 * b + DS += -0.02 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(-45.099, -0.02, 0.0, -0.001, 1, 0, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -45.099 * b + DS += -0.02 * b + GAM1C += 0.0 * a + SINPI += -0.001 * a + + # AddSol(-0.179, -9.52, 0.0, -0.0833, 1, 0, 2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.179 * b + DS += -9.52 * b + GAM1C += 0.0 * a + SINPI += -0.0833 * a + + # AddSol(-0.301, -0.33, 0.0, 0.0014, 1, 0, 2, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -0.301 * b + DS += -0.33 * b + GAM1C += 0.0 * a + SINPI += 0.0014 * a + + # AddSol(-6.382, -3.37, 0.0, -0.0481, 1, 0, -2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -6.382 * b + DS += -3.37 * b + GAM1C += 0.0 * a + SINPI += -0.0481 * a + + # AddSol(39.528, 85.13, 0.0, -0.7136, 1, 0, -2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + DLAM += 39.528 * b + DS += 85.13 * b + GAM1C += 0.0 * a + SINPI += -0.7136 * a + + # AddSol(9.366, 0.71, 0.0, -0.0112, 1, 0, -2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 9.366 * b + DS += 0.71 * b + GAM1C += 0.0 * a + SINPI += -0.0112 * a + + # AddSol(0.202, 0.02, 0.0, 0.0, 1, 0, -2, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += 0.202 * b + DS += 0.02 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + + # AddSol(0.415, 0.1, 0.0, 0.0013, 0, 1, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += 0.415 * b + DS += 0.1 * b + GAM1C += 0.0 * a + SINPI += 0.0013 * a + + # AddSol(-2.152, -2.26, 0.0, -0.0066, 0, 1, 2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -2.152 * b + DS += -2.26 * b + GAM1C += 0.0 * a + SINPI += -0.0066 * a + + # AddSol(-1.44, -1.3, 0.0, 0.0014, 0, 1, -2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -1.44 * b + DS += -1.3 * b + GAM1C += 0.0 * a + SINPI += 0.0014 * a + + # AddSol(0.384, -0.04, 0.0, 0.0, 0, 1, -2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.384 * b + DS += -0.04 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(1.938, 3.6, -0.145, 0.0401, 4, 0, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[4][1], si[4][1]) + DLAM += 1.938 * b + DS += 3.6 * b + GAM1C += -0.145 * a + SINPI += 0.0401 * a + + # AddSol(-0.952, -1.58, 0.052, -0.013, 4, 0, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[4][1], si[4][1]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.952 * b + DS += -1.58 * b + GAM1C += 0.052 * a + SINPI += -0.013 * a + + # AddSol(-0.551, -0.94, 0.032, -0.0097, 3, 1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + DLAM += -0.551 * b + DS += -0.94 * b + GAM1C += 0.032 * a + SINPI += -0.0097 * a + + # AddSol(-0.482, -0.57, 0.005, -0.0045, 3, 1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.482 * b + DS += -0.57 * b + GAM1C += 0.005 * a + SINPI += -0.0045 * a + + # AddSol(0.681, 0.96, -0.026, 0.0115, 3, -1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + DLAM += 0.681 * b + DS += 0.96 * b + GAM1C += -0.026 * a + SINPI += 0.0115 * a + + # AddSol(-0.297, -0.27, 0.002, -0.0009, 2, 2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.297 * b + DS += -0.27 * b + GAM1C += 0.002 * a + SINPI += -0.0009 * a + + # AddSol(0.254, 0.21, -0.003, 0.0, 2, -2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][2], si[-2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.254 * b + DS += 0.21 * b + GAM1C += -0.003 * a + SINPI += 0.0 * a + + # AddSol(-0.25, -0.22, 0.004, 0.0014, 1, 3, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[3][2], si[3][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.25 * b + DS += -0.22 * b + GAM1C += 0.004 * a + SINPI += 0.0014 * a + + # AddSol(-3.996, 0.0, 0.0, 0.0004, 2, 0, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -3.996 * b + DS += 0.0 * b + GAM1C += 0.0 * a + SINPI += 0.0004 * a + + # AddSol(0.557, -0.75, 0.0, -0.009, 2, 0, 2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.557 * b + DS += -0.75 * b + GAM1C += 0.0 * a + SINPI += -0.009 * a + + # AddSol(-0.459, -0.38, 0.0, -0.0053, 2, 0, -2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.459 * b + DS += -0.38 * b + GAM1C += 0.0 * a + SINPI += -0.0053 * a + + # AddSol(-1.298, 0.74, 0.0, 0.0004, 2, 0, -2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + DLAM += -1.298 * b + DS += 0.74 * b + GAM1C += 0.0 * a + SINPI += 0.0004 * a + + # AddSol(0.538, 1.14, 0.0, -0.0141, 2, 0, -2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.538 * b + DS += 1.14 * b + GAM1C += 0.0 * a + SINPI += -0.0141 * a + + # AddSol(0.263, 0.02, 0.0, 0.0, 1, 1, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += 0.263 * b + DS += 0.02 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(0.426, 0.07, 0.0, -0.0006, 1, 1, -2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.426 * b + DS += 0.07 * b + GAM1C += 0.0 * a + SINPI += -0.0006 * a + + # AddSol(-0.304, 0.03, 0.0, 0.0003, 1, -1, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -0.304 * b + DS += 0.03 * b + GAM1C += 0.0 * a + SINPI += 0.0003 * a + + # AddSol(-0.372, -0.19, 0.0, -0.0027, 1, -1, -2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.372 * b + DS += -0.19 * b + GAM1C += 0.0 * a + SINPI += -0.0027 * a + + # AddSol(0.418, 0.0, 0.0, 0.0, 0, 0, 4, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[4][3], si[4][3]) + DLAM += 0.418 * b + DS += 0.0 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(-0.33, -0.04, 0.0, 0.0, 3, 0, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -0.33 * b + DS += -0.04 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a def ADDN(coeffn, p, q, r, s): return coeffn * Term(p, q, r, s)[1] diff --git a/generate/template/astronomy.old b/generate/template/old.py similarity index 99% rename from generate/template/astronomy.old rename to generate/template/old.py index 720771e3..99caf667 100644 --- a/generate/template/astronomy.old +++ b/generate/template/old.py @@ -735,7 +735,15 @@ def _CalcMoon(time): def AddSol(coeffl, coeffs, coeffg, coeffp, p, q, r, s): nonlocal DLAM, DS, GAM1C, SINPI - (a, b) = Term(p, q, r, s) + (a, b) = (1, 0) + if p != 0: + (a, b) = AddThe(a, b, co[p][1], si[p][1]) + if q != 0: + (a, b) = AddThe(a, b, co[q][2], si[q][2]) + if r != 0: + (a, b) = AddThe(a, b, co[r][3], si[r][3]) + if s != 0: + (a, b) = AddThe(a, b, co[s][4], si[s][4]) DLAM += coeffl * b DS += coeffs * b GAM1C += coeffg * a diff --git a/generate/test.py b/generate/test.py old mode 100644 new mode 100755 diff --git a/source/python/astronomy.py b/source/python/astronomy.py index 112fd745..0789b7bd 100644 --- a/source/python/astronomy.py +++ b/source/python/astronomy.py @@ -826,118 +826,990 @@ def _CalcMoon(time): def AddSol(coeffl, coeffs, coeffg, coeffp, p, q, r, s): nonlocal DLAM, DS, GAM1C, SINPI - (a, b) = Term(p, q, r, s) + (a, b) = (1, 0) + if p != 0: + (a, b) = AddThe(a, b, co[p][1], si[p][1]) + if q != 0: + (a, b) = AddThe(a, b, co[q][2], si[q][2]) + if r != 0: + (a, b) = AddThe(a, b, co[r][3], si[r][3]) + if s != 0: + (a, b) = AddThe(a, b, co[s][4], si[s][4]) DLAM += coeffl * b DS += coeffs * b GAM1C += coeffg * a SINPI += coeffp * a - AddSol( 13.902, 14.06,-0.001, 0.2607,0, 0, 0, 4) - AddSol( 0.403, -4.01,+0.394, 0.0023,0, 0, 0, 3) - AddSol( 2369.912, 2373.36,+0.601, 28.2333,0, 0, 0, 2) - AddSol( -125.154, -112.79,-0.725, -0.9781,0, 0, 0, 1) - AddSol( 1.979, 6.98,-0.445, 0.0433,1, 0, 0, 4) - AddSol( 191.953, 192.72,+0.029, 3.0861,1, 0, 0, 2) - AddSol( -8.466, -13.51,+0.455, -0.1093,1, 0, 0, 1) - AddSol( 22639.500,22609.07,+0.079, 186.5398,1, 0, 0, 0) - AddSol( 18.609, 3.59,-0.094, 0.0118,1, 0, 0,-1) - AddSol( -4586.465,-4578.13,-0.077, 34.3117,1, 0, 0,-2) - AddSol( +3.215, 5.44,+0.192, -0.0386,1, 0, 0,-3) - AddSol( -38.428, -38.64,+0.001, 0.6008,1, 0, 0,-4) - AddSol( -0.393, -1.43,-0.092, 0.0086,1, 0, 0,-6) - AddSol( -0.289, -1.59,+0.123, -0.0053,0, 1, 0, 4) - AddSol( -24.420, -25.10,+0.040, -0.3000,0, 1, 0, 2) - AddSol( 18.023, 17.93,+0.007, 0.1494,0, 1, 0, 1) - AddSol( -668.146, -126.98,-1.302, -0.3997,0, 1, 0, 0) - AddSol( 0.560, 0.32,-0.001, -0.0037,0, 1, 0,-1) - AddSol( -165.145, -165.06,+0.054, 1.9178,0, 1, 0,-2) - AddSol( -1.877, -6.46,-0.416, 0.0339,0, 1, 0,-4) - AddSol( 0.213, 1.02,-0.074, 0.0054,2, 0, 0, 4) - AddSol( 14.387, 14.78,-0.017, 0.2833,2, 0, 0, 2) - AddSol( -0.586, -1.20,+0.054, -0.0100,2, 0, 0, 1) - AddSol( 769.016, 767.96,+0.107, 10.1657,2, 0, 0, 0) - AddSol( +1.750, 2.01,-0.018, 0.0155,2, 0, 0,-1) - AddSol( -211.656, -152.53,+5.679, -0.3039,2, 0, 0,-2) - AddSol( +1.225, 0.91,-0.030, -0.0088,2, 0, 0,-3) - AddSol( -30.773, -34.07,-0.308, 0.3722,2, 0, 0,-4) - AddSol( -0.570, -1.40,-0.074, 0.0109,2, 0, 0,-6) - AddSol( -2.921, -11.75,+0.787, -0.0484,1, 1, 0, 2) - AddSol( +1.267, 1.52,-0.022, 0.0164,1, 1, 0, 1) - AddSol( -109.673, -115.18,+0.461, -0.9490,1, 1, 0, 0) - AddSol( -205.962, -182.36,+2.056, +1.4437,1, 1, 0,-2) - AddSol( 0.233, 0.36, 0.012, -0.0025,1, 1, 0,-3) - AddSol( -4.391, -9.66,-0.471, 0.0673,1, 1, 0,-4) - AddSol( 0.283, 1.53,-0.111, +0.0060,1,-1, 0,+4) - AddSol( 14.577, 31.70,-1.540, +0.2302,1,-1, 0, 2) - AddSol( 147.687, 138.76,+0.679, +1.1528,1,-1, 0, 0) - AddSol( -1.089, 0.55,+0.021, 0.0 ,1,-1, 0,-1) - AddSol( 28.475, 23.59,-0.443, -0.2257,1,-1, 0,-2) - AddSol( -0.276, -0.38,-0.006, -0.0036,1,-1, 0,-3) - AddSol( 0.636, 2.27,+0.146, -0.0102,1,-1, 0,-4) - AddSol( -0.189, -1.68,+0.131, -0.0028,0, 2, 0, 2) - AddSol( -7.486, -0.66,-0.037, -0.0086,0, 2, 0, 0) - AddSol( -8.096, -16.35,-0.740, 0.0918,0, 2, 0,-2) - AddSol( -5.741, -0.04, 0.0 , -0.0009,0, 0, 2, 2) - AddSol( 0.255, 0.0 , 0.0 , 0.0 ,0, 0, 2, 1) - AddSol( -411.608, -0.20, 0.0 , -0.0124,0, 0, 2, 0) - AddSol( 0.584, 0.84, 0.0 , +0.0071,0, 0, 2,-1) - AddSol( -55.173, -52.14, 0.0 , -0.1052,0, 0, 2,-2) - AddSol( 0.254, 0.25, 0.0 , -0.0017,0, 0, 2,-3) - AddSol( +0.025, -1.67, 0.0 , +0.0031,0, 0, 2,-4) - AddSol( 1.060, 2.96,-0.166, 0.0243,3, 0, 0,+2) - AddSol( 36.124, 50.64,-1.300, 0.6215,3, 0, 0, 0) - AddSol( -13.193, -16.40,+0.258, -0.1187,3, 0, 0,-2) - AddSol( -1.187, -0.74,+0.042, 0.0074,3, 0, 0,-4) - AddSol( -0.293, -0.31,-0.002, 0.0046,3, 0, 0,-6) - AddSol( -0.290, -1.45,+0.116, -0.0051,2, 1, 0, 2) - AddSol( -7.649, -10.56,+0.259, -0.1038,2, 1, 0, 0) - AddSol( -8.627, -7.59,+0.078, -0.0192,2, 1, 0,-2) - AddSol( -2.740, -2.54,+0.022, 0.0324,2, 1, 0,-4) - AddSol( 1.181, 3.32,-0.212, 0.0213,2,-1, 0,+2) - AddSol( 9.703, 11.67,-0.151, 0.1268,2,-1, 0, 0) - AddSol( -0.352, -0.37,+0.001, -0.0028,2,-1, 0,-1) - AddSol( -2.494, -1.17,-0.003, -0.0017,2,-1, 0,-2) - AddSol( 0.360, 0.20,-0.012, -0.0043,2,-1, 0,-4) - AddSol( -1.167, -1.25,+0.008, -0.0106,1, 2, 0, 0) - AddSol( -7.412, -6.12,+0.117, 0.0484,1, 2, 0,-2) - AddSol( -0.311, -0.65,-0.032, 0.0044,1, 2, 0,-4) - AddSol( +0.757, 1.82,-0.105, 0.0112,1,-2, 0, 2) - AddSol( +2.580, 2.32,+0.027, 0.0196,1,-2, 0, 0) - AddSol( +2.533, 2.40,-0.014, -0.0212,1,-2, 0,-2) - AddSol( -0.344, -0.57,-0.025, +0.0036,0, 3, 0,-2) - AddSol( -0.992, -0.02, 0.0 , 0.0 ,1, 0, 2, 2) - AddSol( -45.099, -0.02, 0.0 , -0.0010,1, 0, 2, 0) - AddSol( -0.179, -9.52, 0.0 , -0.0833,1, 0, 2,-2) - AddSol( -0.301, -0.33, 0.0 , 0.0014,1, 0, 2,-4) - AddSol( -6.382, -3.37, 0.0 , -0.0481,1, 0,-2, 2) - AddSol( 39.528, 85.13, 0.0 , -0.7136,1, 0,-2, 0) - AddSol( 9.366, 0.71, 0.0 , -0.0112,1, 0,-2,-2) - AddSol( 0.202, 0.02, 0.0 , 0.0 ,1, 0,-2,-4) + # AddSol(13.902, 14.06, -0.001, 0.2607, 0, 0, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += 13.902 * b + DS += 14.06 * b + GAM1C += -0.001 * a + SINPI += 0.2607 * a - AddSol( 0.415, 0.10, 0.0 , 0.0013,0, 1, 2, 0) - AddSol( -2.152, -2.26, 0.0 , -0.0066,0, 1, 2,-2) - AddSol( -1.440, -1.30, 0.0 , +0.0014,0, 1,-2, 2) - AddSol( 0.384, -0.04, 0.0 , 0.0 ,0, 1,-2,-2) - AddSol( +1.938, +3.60,-0.145, +0.0401,4, 0, 0, 0) - AddSol( -0.952, -1.58,+0.052, -0.0130,4, 0, 0,-2) - AddSol( -0.551, -0.94,+0.032, -0.0097,3, 1, 0, 0) - AddSol( -0.482, -0.57,+0.005, -0.0045,3, 1, 0,-2) - AddSol( 0.681, 0.96,-0.026, 0.0115,3,-1, 0, 0) - AddSol( -0.297, -0.27, 0.002, -0.0009,2, 2, 0,-2) - AddSol( 0.254, +0.21,-0.003, 0.0 ,2,-2, 0,-2) - AddSol( -0.250, -0.22, 0.004, 0.0014,1, 3, 0,-2) - AddSol( -3.996, 0.0 , 0.0 , +0.0004,2, 0, 2, 0) - AddSol( 0.557, -0.75, 0.0 , -0.0090,2, 0, 2,-2) - AddSol( -0.459, -0.38, 0.0 , -0.0053,2, 0,-2, 2) - AddSol( -1.298, 0.74, 0.0 , +0.0004,2, 0,-2, 0) - AddSol( 0.538, 1.14, 0.0 , -0.0141,2, 0,-2,-2) - AddSol( 0.263, 0.02, 0.0 , 0.0 ,1, 1, 2, 0) - AddSol( 0.426, +0.07, 0.0 , -0.0006,1, 1,-2,-2) - AddSol( -0.304, +0.03, 0.0 , +0.0003,1,-1, 2, 0) - AddSol( -0.372, -0.19, 0.0 , -0.0027,1,-1,-2, 2) - AddSol( +0.418, 0.0 , 0.0 , 0.0 ,0, 0, 4, 0) - AddSol( -0.330, -0.04, 0.0 , 0.0 ,3, 0, 2, 0) + # AddSol(0.403, -4.01, 0.394, 0.0023, 0, 0, 0, 3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][4], si[3][4]) + DLAM += 0.403 * b + DS += -4.01 * b + GAM1C += 0.394 * a + SINPI += 0.0023 * a + + # AddSol(2369.912, 2373.36, 0.601, 28.2333, 0, 0, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 2369.912 * b + DS += 2373.36 * b + GAM1C += 0.601 * a + SINPI += 28.2333 * a + + # AddSol(-125.154, -112.79, -0.725, -0.9781, 0, 0, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += -125.154 * b + DS += -112.79 * b + GAM1C += -0.725 * a + SINPI += -0.9781 * a + + # AddSol(1.979, 6.98, -0.445, 0.0433, 1, 0, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += 1.979 * b + DS += 6.98 * b + GAM1C += -0.445 * a + SINPI += 0.0433 * a + + # AddSol(191.953, 192.72, 0.029, 3.0861, 1, 0, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 191.953 * b + DS += 192.72 * b + GAM1C += 0.029 * a + SINPI += 3.0861 * a + + # AddSol(-8.466, -13.51, 0.455, -0.1093, 1, 0, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += -8.466 * b + DS += -13.51 * b + GAM1C += 0.455 * a + SINPI += -0.1093 * a + + # AddSol(22639.5, 22609.07, 0.079, 186.5398, 1, 0, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + DLAM += 22639.5 * b + DS += 22609.07 * b + GAM1C += 0.079 * a + SINPI += 186.5398 * a + + # AddSol(18.609, 3.59, -0.094, 0.0118, 1, 0, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += 18.609 * b + DS += 3.59 * b + GAM1C += -0.094 * a + SINPI += 0.0118 * a + + # AddSol(-4586.465, -4578.13, -0.077, 34.3117, 1, 0, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -4586.465 * b + DS += -4578.13 * b + GAM1C += -0.077 * a + SINPI += 34.3117 * a + + # AddSol(3.215, 5.44, 0.192, -0.0386, 1, 0, 0, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += 3.215 * b + DS += 5.44 * b + GAM1C += 0.192 * a + SINPI += -0.0386 * a + + # AddSol(-38.428, -38.64, 0.001, 0.6008, 1, 0, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -38.428 * b + DS += -38.64 * b + GAM1C += 0.001 * a + SINPI += 0.6008 * a + + # AddSol(-0.393, -1.43, -0.092, 0.0086, 1, 0, 0, -6) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-6][4], si[-6][4]) + DLAM += -0.393 * b + DS += -1.43 * b + GAM1C += -0.092 * a + SINPI += 0.0086 * a + + # AddSol(-0.289, -1.59, 0.123, -0.0053, 0, 1, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += -0.289 * b + DS += -1.59 * b + GAM1C += 0.123 * a + SINPI += -0.0053 * a + + # AddSol(-24.42, -25.1, 0.04, -0.3, 0, 1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -24.42 * b + DS += -25.1 * b + GAM1C += 0.04 * a + SINPI += -0.3 * a + + # AddSol(18.023, 17.93, 0.007, 0.1494, 0, 1, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += 18.023 * b + DS += 17.93 * b + GAM1C += 0.007 * a + SINPI += 0.1494 * a + + # AddSol(-668.146, -126.98, -1.302, -0.3997, 0, 1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + DLAM += -668.146 * b + DS += -126.98 * b + GAM1C += -1.302 * a + SINPI += -0.3997 * a + + # AddSol(0.56, 0.32, -0.001, -0.0037, 0, 1, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += 0.56 * b + DS += 0.32 * b + GAM1C += -0.001 * a + SINPI += -0.0037 * a + + # AddSol(-165.145, -165.06, 0.054, 1.9178, 0, 1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -165.145 * b + DS += -165.06 * b + GAM1C += 0.054 * a + SINPI += 1.9178 * a + + # AddSol(-1.877, -6.46, -0.416, 0.0339, 0, 1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -1.877 * b + DS += -6.46 * b + GAM1C += -0.416 * a + SINPI += 0.0339 * a + + # AddSol(0.213, 1.02, -0.074, 0.0054, 2, 0, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += 0.213 * b + DS += 1.02 * b + GAM1C += -0.074 * a + SINPI += 0.0054 * a + + # AddSol(14.387, 14.78, -0.017, 0.2833, 2, 0, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 14.387 * b + DS += 14.78 * b + GAM1C += -0.017 * a + SINPI += 0.2833 * a + + # AddSol(-0.586, -1.2, 0.054, -0.01, 2, 0, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += -0.586 * b + DS += -1.2 * b + GAM1C += 0.054 * a + SINPI += -0.01 * a + + # AddSol(769.016, 767.96, 0.107, 10.1657, 2, 0, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + DLAM += 769.016 * b + DS += 767.96 * b + GAM1C += 0.107 * a + SINPI += 10.1657 * a + + # AddSol(1.75, 2.01, -0.018, 0.0155, 2, 0, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += 1.75 * b + DS += 2.01 * b + GAM1C += -0.018 * a + SINPI += 0.0155 * a + + # AddSol(-211.656, -152.53, 5.679, -0.3039, 2, 0, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -211.656 * b + DS += -152.53 * b + GAM1C += 5.679 * a + SINPI += -0.3039 * a + + # AddSol(1.225, 0.91, -0.03, -0.0088, 2, 0, 0, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += 1.225 * b + DS += 0.91 * b + GAM1C += -0.03 * a + SINPI += -0.0088 * a + + # AddSol(-30.773, -34.07, -0.308, 0.3722, 2, 0, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -30.773 * b + DS += -34.07 * b + GAM1C += -0.308 * a + SINPI += 0.3722 * a + + # AddSol(-0.57, -1.4, -0.074, 0.0109, 2, 0, 0, -6) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-6][4], si[-6][4]) + DLAM += -0.57 * b + DS += -1.4 * b + GAM1C += -0.074 * a + SINPI += 0.0109 * a + + # AddSol(-2.921, -11.75, 0.787, -0.0484, 1, 1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -2.921 * b + DS += -11.75 * b + GAM1C += 0.787 * a + SINPI += -0.0484 * a + + # AddSol(1.267, 1.52, -0.022, 0.0164, 1, 1, 0, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += 1.267 * b + DS += 1.52 * b + GAM1C += -0.022 * a + SINPI += 0.0164 * a + + # AddSol(-109.673, -115.18, 0.461, -0.949, 1, 1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + DLAM += -109.673 * b + DS += -115.18 * b + GAM1C += 0.461 * a + SINPI += -0.949 * a + + # AddSol(-205.962, -182.36, 2.056, 1.4437, 1, 1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -205.962 * b + DS += -182.36 * b + GAM1C += 2.056 * a + SINPI += 1.4437 * a + + # AddSol(0.233, 0.36, 0.012, -0.0025, 1, 1, 0, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += 0.233 * b + DS += 0.36 * b + GAM1C += 0.012 * a + SINPI += -0.0025 * a + + # AddSol(-4.391, -9.66, -0.471, 0.0673, 1, 1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -4.391 * b + DS += -9.66 * b + GAM1C += -0.471 * a + SINPI += 0.0673 * a + + + # AddSol(0.283, 1.53, -0.111, 0.006, 1, -1, 0, 4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[4][4], si[4][4]) + DLAM += 0.283 * b + DS += 1.53 * b + GAM1C += -0.111 * a + SINPI += 0.006 * a + + # AddSol(14.577, 31.7, -1.54, 0.2302, 1, -1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 14.577 * b + DS += 31.7 * b + GAM1C += -1.54 * a + SINPI += 0.2302 * a + + # AddSol(147.687, 138.76, 0.679, 1.1528, 1, -1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + DLAM += 147.687 * b + DS += 138.76 * b + GAM1C += 0.679 * a + SINPI += 1.1528 * a + + # AddSol(-1.089, 0.55, 0.021, 0.0, 1, -1, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += -1.089 * b + DS += 0.55 * b + GAM1C += 0.021 * a + SINPI += 0.0 * a + + # AddSol(28.475, 23.59, -0.443, -0.2257, 1, -1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 28.475 * b + DS += 23.59 * b + GAM1C += -0.443 * a + SINPI += -0.2257 * a + + # AddSol(-0.276, -0.38, -0.006, -0.0036, 1, -1, 0, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += -0.276 * b + DS += -0.38 * b + GAM1C += -0.006 * a + SINPI += -0.0036 * a + + # AddSol(0.636, 2.27, 0.146, -0.0102, 1, -1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += 0.636 * b + DS += 2.27 * b + GAM1C += 0.146 * a + SINPI += -0.0102 * a + + # AddSol(-0.189, -1.68, 0.131, -0.0028, 0, 2, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.189 * b + DS += -1.68 * b + GAM1C += 0.131 * a + SINPI += -0.0028 * a + + # AddSol(-7.486, -0.66, -0.037, -0.0086, 0, 2, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + DLAM += -7.486 * b + DS += -0.66 * b + GAM1C += -0.037 * a + SINPI += -0.0086 * a + + # AddSol(-8.096, -16.35, -0.74, 0.0918, 0, 2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -8.096 * b + DS += -16.35 * b + GAM1C += -0.74 * a + SINPI += 0.0918 * a + + # AddSol(-5.741, -0.04, 0.0, -0.0009, 0, 0, 2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -5.741 * b + DS += -0.04 * b + GAM1C += 0.0 * a + SINPI += -0.0009 * a + + # AddSol(0.255, 0.0, 0.0, 0.0, 0, 0, 2, 1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[1][4], si[1][4]) + DLAM += 0.255 * b + DS += 0.0 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(-411.608, -0.2, 0.0, -0.0124, 0, 0, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -411.608 * b + DS += -0.2 * b + GAM1C += 0.0 * a + SINPI += -0.0124 * a + + # AddSol(0.584, 0.84, 0.0, 0.0071, 0, 0, 2, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += 0.584 * b + DS += 0.84 * b + GAM1C += 0.0 * a + SINPI += 0.0071 * a + + # AddSol(-55.173, -52.14, 0.0, -0.1052, 0, 0, 2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -55.173 * b + DS += -52.14 * b + GAM1C += 0.0 * a + SINPI += -0.1052 * a + + # AddSol(0.254, 0.25, 0.0, -0.0017, 0, 0, 2, -3) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-3][4], si[-3][4]) + DLAM += 0.254 * b + DS += 0.25 * b + GAM1C += 0.0 * a + SINPI += -0.0017 * a + + # AddSol(0.025, -1.67, 0.0, 0.0031, 0, 0, 2, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += 0.025 * b + DS += -1.67 * b + GAM1C += 0.0 * a + SINPI += 0.0031 * a + + # AddSol(1.06, 2.96, -0.166, 0.0243, 3, 0, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 1.06 * b + DS += 2.96 * b + GAM1C += -0.166 * a + SINPI += 0.0243 * a + + # AddSol(36.124, 50.64, -1.3, 0.6215, 3, 0, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + DLAM += 36.124 * b + DS += 50.64 * b + GAM1C += -1.3 * a + SINPI += 0.6215 * a + + # AddSol(-13.193, -16.4, 0.258, -0.1187, 3, 0, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -13.193 * b + DS += -16.4 * b + GAM1C += 0.258 * a + SINPI += -0.1187 * a + + # AddSol(-1.187, -0.74, 0.042, 0.0074, 3, 0, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -1.187 * b + DS += -0.74 * b + GAM1C += 0.042 * a + SINPI += 0.0074 * a + + # AddSol(-0.293, -0.31, -0.002, 0.0046, 3, 0, 0, -6) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[-6][4], si[-6][4]) + DLAM += -0.293 * b + DS += -0.31 * b + GAM1C += -0.002 * a + SINPI += 0.0046 * a + + # AddSol(-0.29, -1.45, 0.116, -0.0051, 2, 1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.29 * b + DS += -1.45 * b + GAM1C += 0.116 * a + SINPI += -0.0051 * a + + # AddSol(-7.649, -10.56, 0.259, -0.1038, 2, 1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + DLAM += -7.649 * b + DS += -10.56 * b + GAM1C += 0.259 * a + SINPI += -0.1038 * a + + # AddSol(-8.627, -7.59, 0.078, -0.0192, 2, 1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -8.627 * b + DS += -7.59 * b + GAM1C += 0.078 * a + SINPI += -0.0192 * a + + # AddSol(-2.74, -2.54, 0.022, 0.0324, 2, 1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -2.74 * b + DS += -2.54 * b + GAM1C += 0.022 * a + SINPI += 0.0324 * a + + # AddSol(1.181, 3.32, -0.212, 0.0213, 2, -1, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 1.181 * b + DS += 3.32 * b + GAM1C += -0.212 * a + SINPI += 0.0213 * a + + # AddSol(9.703, 11.67, -0.151, 0.1268, 2, -1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + DLAM += 9.703 * b + DS += 11.67 * b + GAM1C += -0.151 * a + SINPI += 0.1268 * a + + # AddSol(-0.352, -0.37, 0.001, -0.0028, 2, -1, 0, -1) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-1][4], si[-1][4]) + DLAM += -0.352 * b + DS += -0.37 * b + GAM1C += 0.001 * a + SINPI += -0.0028 * a + + # AddSol(-2.494, -1.17, -0.003, -0.0017, 2, -1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -2.494 * b + DS += -1.17 * b + GAM1C += -0.003 * a + SINPI += -0.0017 * a + + # AddSol(0.36, 0.2, -0.012, -0.0043, 2, -1, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += 0.36 * b + DS += 0.2 * b + GAM1C += -0.012 * a + SINPI += -0.0043 * a + + # AddSol(-1.167, -1.25, 0.008, -0.0106, 1, 2, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + DLAM += -1.167 * b + DS += -1.25 * b + GAM1C += 0.008 * a + SINPI += -0.0106 * a + + # AddSol(-7.412, -6.12, 0.117, 0.0484, 1, 2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -7.412 * b + DS += -6.12 * b + GAM1C += 0.117 * a + SINPI += 0.0484 * a + + # AddSol(-0.311, -0.65, -0.032, 0.0044, 1, 2, 0, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -0.311 * b + DS += -0.65 * b + GAM1C += -0.032 * a + SINPI += 0.0044 * a + + # AddSol(0.757, 1.82, -0.105, 0.0112, 1, -2, 0, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][2], si[-2][2]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += 0.757 * b + DS += 1.82 * b + GAM1C += -0.105 * a + SINPI += 0.0112 * a + + # AddSol(2.58, 2.32, 0.027, 0.0196, 1, -2, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][2], si[-2][2]) + DLAM += 2.58 * b + DS += 2.32 * b + GAM1C += 0.027 * a + SINPI += 0.0196 * a + + # AddSol(2.533, 2.4, -0.014, -0.0212, 1, -2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][2], si[-2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 2.533 * b + DS += 2.4 * b + GAM1C += -0.014 * a + SINPI += -0.0212 * a + + # AddSol(-0.344, -0.57, -0.025, 0.0036, 0, 3, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][2], si[3][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.344 * b + DS += -0.57 * b + GAM1C += -0.025 * a + SINPI += 0.0036 * a + + # AddSol(-0.992, -0.02, 0.0, 0.0, 1, 0, 2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.992 * b + DS += -0.02 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(-45.099, -0.02, 0.0, -0.001, 1, 0, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -45.099 * b + DS += -0.02 * b + GAM1C += 0.0 * a + SINPI += -0.001 * a + + # AddSol(-0.179, -9.52, 0.0, -0.0833, 1, 0, 2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.179 * b + DS += -9.52 * b + GAM1C += 0.0 * a + SINPI += -0.0833 * a + + # AddSol(-0.301, -0.33, 0.0, 0.0014, 1, 0, 2, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += -0.301 * b + DS += -0.33 * b + GAM1C += 0.0 * a + SINPI += 0.0014 * a + + # AddSol(-6.382, -3.37, 0.0, -0.0481, 1, 0, -2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -6.382 * b + DS += -3.37 * b + GAM1C += 0.0 * a + SINPI += -0.0481 * a + + # AddSol(39.528, 85.13, 0.0, -0.7136, 1, 0, -2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + DLAM += 39.528 * b + DS += 85.13 * b + GAM1C += 0.0 * a + SINPI += -0.7136 * a + + # AddSol(9.366, 0.71, 0.0, -0.0112, 1, 0, -2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 9.366 * b + DS += 0.71 * b + GAM1C += 0.0 * a + SINPI += -0.0112 * a + + # AddSol(0.202, 0.02, 0.0, 0.0, 1, 0, -2, -4) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-4][4], si[-4][4]) + DLAM += 0.202 * b + DS += 0.02 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + + # AddSol(0.415, 0.1, 0.0, 0.0013, 0, 1, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += 0.415 * b + DS += 0.1 * b + GAM1C += 0.0 * a + SINPI += 0.0013 * a + + # AddSol(-2.152, -2.26, 0.0, -0.0066, 0, 1, 2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -2.152 * b + DS += -2.26 * b + GAM1C += 0.0 * a + SINPI += -0.0066 * a + + # AddSol(-1.44, -1.3, 0.0, 0.0014, 0, 1, -2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -1.44 * b + DS += -1.3 * b + GAM1C += 0.0 * a + SINPI += 0.0014 * a + + # AddSol(0.384, -0.04, 0.0, 0.0, 0, 1, -2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.384 * b + DS += -0.04 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(1.938, 3.6, -0.145, 0.0401, 4, 0, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[4][1], si[4][1]) + DLAM += 1.938 * b + DS += 3.6 * b + GAM1C += -0.145 * a + SINPI += 0.0401 * a + + # AddSol(-0.952, -1.58, 0.052, -0.013, 4, 0, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[4][1], si[4][1]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.952 * b + DS += -1.58 * b + GAM1C += 0.052 * a + SINPI += -0.013 * a + + # AddSol(-0.551, -0.94, 0.032, -0.0097, 3, 1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + DLAM += -0.551 * b + DS += -0.94 * b + GAM1C += 0.032 * a + SINPI += -0.0097 * a + + # AddSol(-0.482, -0.57, 0.005, -0.0045, 3, 1, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.482 * b + DS += -0.57 * b + GAM1C += 0.005 * a + SINPI += -0.0045 * a + + # AddSol(0.681, 0.96, -0.026, 0.0115, 3, -1, 0, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + DLAM += 0.681 * b + DS += 0.96 * b + GAM1C += -0.026 * a + SINPI += 0.0115 * a + + # AddSol(-0.297, -0.27, 0.002, -0.0009, 2, 2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[2][2], si[2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.297 * b + DS += -0.27 * b + GAM1C += 0.002 * a + SINPI += -0.0009 * a + + # AddSol(0.254, 0.21, -0.003, 0.0, 2, -2, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][2], si[-2][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.254 * b + DS += 0.21 * b + GAM1C += -0.003 * a + SINPI += 0.0 * a + + # AddSol(-0.25, -0.22, 0.004, 0.0014, 1, 3, 0, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[3][2], si[3][2]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += -0.25 * b + DS += -0.22 * b + GAM1C += 0.004 * a + SINPI += 0.0014 * a + + # AddSol(-3.996, 0.0, 0.0, 0.0004, 2, 0, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -3.996 * b + DS += 0.0 * b + GAM1C += 0.0 * a + SINPI += 0.0004 * a + + # AddSol(0.557, -0.75, 0.0, -0.009, 2, 0, 2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.557 * b + DS += -0.75 * b + GAM1C += 0.0 * a + SINPI += -0.009 * a + + # AddSol(-0.459, -0.38, 0.0, -0.0053, 2, 0, -2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.459 * b + DS += -0.38 * b + GAM1C += 0.0 * a + SINPI += -0.0053 * a + + # AddSol(-1.298, 0.74, 0.0, 0.0004, 2, 0, -2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + DLAM += -1.298 * b + DS += 0.74 * b + GAM1C += 0.0 * a + SINPI += 0.0004 * a + + # AddSol(0.538, 1.14, 0.0, -0.0141, 2, 0, -2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[2][1], si[2][1]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.538 * b + DS += 1.14 * b + GAM1C += 0.0 * a + SINPI += -0.0141 * a + + # AddSol(0.263, 0.02, 0.0, 0.0, 1, 1, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += 0.263 * b + DS += 0.02 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(0.426, 0.07, 0.0, -0.0006, 1, 1, -2, -2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[1][2], si[1][2]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[-2][4], si[-2][4]) + DLAM += 0.426 * b + DS += 0.07 * b + GAM1C += 0.0 * a + SINPI += -0.0006 * a + + # AddSol(-0.304, 0.03, 0.0, 0.0003, 1, -1, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -0.304 * b + DS += 0.03 * b + GAM1C += 0.0 * a + SINPI += 0.0003 * a + + # AddSol(-0.372, -0.19, 0.0, -0.0027, 1, -1, -2, 2) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[1][1], si[1][1]) + (a, b) = AddThe(a, b, co[-1][2], si[-1][2]) + (a, b) = AddThe(a, b, co[-2][3], si[-2][3]) + (a, b) = AddThe(a, b, co[2][4], si[2][4]) + DLAM += -0.372 * b + DS += -0.19 * b + GAM1C += 0.0 * a + SINPI += -0.0027 * a + + # AddSol(0.418, 0.0, 0.0, 0.0, 0, 0, 4, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[4][3], si[4][3]) + DLAM += 0.418 * b + DS += 0.0 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a + + # AddSol(-0.33, -0.04, 0.0, 0.0, 3, 0, 2, 0) + (a, b) = (1, 0) + (a, b) = AddThe(a, b, co[3][1], si[3][1]) + (a, b) = AddThe(a, b, co[2][3], si[2][3]) + DLAM += -0.33 * b + DS += -0.04 * b + GAM1C += 0.0 * a + SINPI += 0.0 * a def ADDN(coeffn, p, q, r, s): return coeffn * Term(p, q, r, s)[1]