diff options
-rw-r--r-- | src/lib/netlist/generated/static_solvers.cpp | 4216 | ||||
-rw-r--r-- | src/mame/audio/nl_brdrline.cpp | 42 |
2 files changed, 3454 insertions, 804 deletions
diff --git a/src/lib/netlist/generated/static_solvers.cpp b/src/lib/netlist/generated/static_solvers.cpp index 605ce8a7886..ce58a075fa7 100644 --- a/src/lib/netlist/generated/static_solvers.cpp +++ b/src/lib/netlist/generated/static_solvers.cpp @@ -3039,6 +3039,246 @@ static void nl_gcr_1fad5cda2646cf42_30_double_double(double * __restrict V, cons V[0] = (RHS0 - tmp0) / m_A0; } +// brdrline +static void nl_gcr_1fec893028a7b809_34_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) + +{ + + plib::unused_var(cnV); + double m_A0(0.0); + double m_A1(0.0); + double m_A2(0.0); + double m_A3(0.0); + double m_A4(0.0); + double m_A5(0.0); + double m_A6(0.0); + double m_A7(0.0); + double m_A8(0.0); + double m_A9(0.0); + double m_A10(0.0); + double m_A11(0.0); + double m_A12(0.0); + double m_A13(0.0); + double m_A14(0.0); + double m_A15(0.0); + double m_A16(0.0); + double m_A17(0.0); + double m_A18(0.0); + double m_A19(0.0); + double m_A20(0.0); + double m_A21(0.0); + double m_A22(0.0); + double m_A23(0.0); + double m_A24(0.0); + double m_A25(0.0); + double m_A26(0.0); + double m_A27(0.0); + double m_A28(0.0); + double m_A29(0.0); + double m_A30(0.0); + double m_A31(0.0); + double m_A32(0.0); + double m_A33(0.0); + m_A0 += gt[0]; + m_A0 += gt[1]; + m_A1 += go[0]; + double RHS0 = Idr[0]; + RHS0 += Idr[1]; + RHS0 -= go[1] * *cnV[1]; + m_A2 += gt[2]; + m_A2 += gt[3]; + m_A3 += go[2]; + double RHS1 = Idr[2]; + RHS1 += Idr[3]; + RHS1 -= go[3] * *cnV[3]; + m_A4 += gt[4]; + m_A4 += gt[5]; + m_A5 += go[4]; + double RHS2 = Idr[4]; + RHS2 += Idr[5]; + RHS2 -= go[5] * *cnV[5]; + m_A6 += gt[6]; + m_A6 += gt[7]; + m_A7 += go[6]; + double RHS3 = Idr[6]; + RHS3 += Idr[7]; + RHS3 -= go[7] * *cnV[7]; + m_A8 += gt[8]; + m_A8 += gt[9]; + m_A9 += go[8]; + double RHS4 = Idr[8]; + RHS4 += Idr[9]; + RHS4 -= go[9] * *cnV[9]; + m_A10 += gt[10]; + m_A10 += gt[11]; + m_A11 += go[10]; + double RHS5 = Idr[10]; + RHS5 += Idr[11]; + RHS5 -= go[11] * *cnV[11]; + m_A12 += gt[12]; + m_A12 += gt[13]; + m_A13 += go[12]; + double RHS6 = Idr[12]; + RHS6 += Idr[13]; + RHS6 -= go[13] * *cnV[13]; + m_A14 += gt[14]; + m_A14 += gt[15]; + m_A15 += go[14]; + double RHS7 = Idr[14]; + RHS7 += Idr[15]; + RHS7 -= go[15] * *cnV[15]; + m_A16 += gt[16]; + m_A16 += gt[17]; + m_A17 += go[16]; + double RHS8 = Idr[16]; + RHS8 += Idr[17]; + RHS8 -= go[17] * *cnV[17]; + m_A18 += gt[18]; + m_A18 += gt[19]; + m_A18 += gt[20]; + m_A18 += gt[21]; + m_A18 += gt[22]; + m_A18 += gt[23]; + m_A18 += gt[24]; + m_A19 += go[18]; + double RHS9 = Idr[18]; + RHS9 += Idr[19]; + RHS9 += Idr[20]; + RHS9 += Idr[21]; + RHS9 += Idr[22]; + RHS9 += Idr[23]; + RHS9 += Idr[24]; + RHS9 -= go[19] * *cnV[19]; + RHS9 -= go[20] * *cnV[20]; + RHS9 -= go[21] * *cnV[21]; + RHS9 -= go[22] * *cnV[22]; + RHS9 -= go[23] * *cnV[23]; + RHS9 -= go[24] * *cnV[24]; + m_A28 += gt[25]; + m_A28 += gt[26]; + m_A28 += gt[27]; + m_A28 += gt[28]; + m_A28 += gt[29]; + m_A28 += gt[30]; + m_A28 += gt[31]; + m_A28 += gt[32]; + m_A28 += gt[33]; + m_A28 += gt[34]; + m_A29 += go[25]; + m_A27 += go[26]; + m_A26 += go[27]; + m_A25 += go[28]; + m_A24 += go[29]; + m_A23 += go[30]; + m_A22 += go[31]; + m_A21 += go[32]; + m_A20 += go[33]; + double RHS10 = Idr[25]; + RHS10 += Idr[26]; + RHS10 += Idr[27]; + RHS10 += Idr[28]; + RHS10 += Idr[29]; + RHS10 += Idr[30]; + RHS10 += Idr[31]; + RHS10 += Idr[32]; + RHS10 += Idr[33]; + RHS10 += Idr[34]; + RHS10 -= go[34] * *cnV[34]; + m_A33 += gt[35]; + m_A33 += gt[36]; + m_A33 += gt[37]; + m_A33 += gt[38]; + m_A33 += gt[39]; + m_A31 += go[35]; + m_A30 += go[36]; + m_A32 += go[37]; + double RHS11 = Idr[35]; + RHS11 += Idr[36]; + RHS11 += Idr[37]; + RHS11 += Idr[38]; + RHS11 += Idr[39]; + RHS11 -= go[38] * *cnV[38]; + RHS11 -= go[39] * *cnV[39]; + const double f0 = 1.0 / m_A0; + const double f0_10 = -f0 * m_A20; + m_A28 += m_A1 * f0_10; + RHS10 += f0_10 * RHS0; + const double f1 = 1.0 / m_A2; + const double f1_10 = -f1 * m_A21; + m_A28 += m_A3 * f1_10; + RHS10 += f1_10 * RHS1; + const double f2 = 1.0 / m_A4; + const double f2_11 = -f2 * m_A30; + m_A33 += m_A5 * f2_11; + RHS11 += f2_11 * RHS2; + const double f3 = 1.0 / m_A6; + const double f3_10 = -f3 * m_A22; + m_A28 += m_A7 * f3_10; + RHS10 += f3_10 * RHS3; + const double f4 = 1.0 / m_A8; + const double f4_10 = -f4 * m_A23; + m_A28 += m_A9 * f4_10; + RHS10 += f4_10 * RHS4; + const double f5 = 1.0 / m_A10; + const double f5_10 = -f5 * m_A24; + m_A28 += m_A11 * f5_10; + RHS10 += f5_10 * RHS5; + const double f6 = 1.0 / m_A12; + const double f6_10 = -f6 * m_A25; + m_A28 += m_A13 * f6_10; + RHS10 += f6_10 * RHS6; + const double f7 = 1.0 / m_A14; + const double f7_10 = -f7 * m_A26; + m_A28 += m_A15 * f7_10; + RHS10 += f7_10 * RHS7; + const double f8 = 1.0 / m_A16; + const double f8_10 = -f8 * m_A27; + m_A28 += m_A17 * f8_10; + RHS10 += f8_10 * RHS8; + const double f9 = 1.0 / m_A18; + const double f9_11 = -f9 * m_A31; + m_A32 += m_A19 * f9_11; + RHS11 += f9_11 * RHS9; + const double f10 = 1.0 / m_A28; + const double f10_11 = -f10 * m_A32; + m_A33 += m_A29 * f10_11; + RHS11 += f10_11 * RHS10; + V[11] = RHS11 / m_A33; + double tmp10 = 0.0; + tmp10 += m_A29 * V[11]; + V[10] = (RHS10 - tmp10) / m_A28; + double tmp9 = 0.0; + tmp9 += m_A19 * V[10]; + V[9] = (RHS9 - tmp9) / m_A18; + double tmp8 = 0.0; + tmp8 += m_A17 * V[10]; + V[8] = (RHS8 - tmp8) / m_A16; + double tmp7 = 0.0; + tmp7 += m_A15 * V[10]; + V[7] = (RHS7 - tmp7) / m_A14; + double tmp6 = 0.0; + tmp6 += m_A13 * V[10]; + V[6] = (RHS6 - tmp6) / m_A12; + double tmp5 = 0.0; + tmp5 += m_A11 * V[10]; + V[5] = (RHS5 - tmp5) / m_A10; + double tmp4 = 0.0; + tmp4 += m_A9 * V[10]; + V[4] = (RHS4 - tmp4) / m_A8; + double tmp3 = 0.0; + tmp3 += m_A7 * V[10]; + V[3] = (RHS3 - tmp3) / m_A6; + double tmp2 = 0.0; + tmp2 += m_A5 * V[11]; + V[2] = (RHS2 - tmp2) / m_A4; + double tmp1 = 0.0; + tmp1 += m_A3 * V[10]; + V[1] = (RHS1 - tmp1) / m_A2; + double tmp0 = 0.0; + tmp0 += m_A1 * V[10]; + V[0] = (RHS0 - tmp0) / m_A0; +} + // sspeedr static void nl_gcr_2294220d3c91e762_176_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) @@ -4145,6 +4385,718 @@ static void nl_gcr_2294220d3c91e762_176_double_double(double * __restrict V, con V[0] = (RHS0 - tmp0) / m_A0; } +// headon +static void nl_gcr_24023c8529f6bdb3_97_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) + +{ + + plib::unused_var(cnV); + double m_A0(0.0); + double m_A1(0.0); + double m_A2(0.0); + double m_A3(0.0); + double m_A4(0.0); + double m_A5(0.0); + double m_A6(0.0); + double m_A7(0.0); + double m_A8(0.0); + double m_A9(0.0); + double m_A10(0.0); + double m_A11(0.0); + double m_A12(0.0); + double m_A13(0.0); + double m_A14(0.0); + double m_A15(0.0); + double m_A16(0.0); + double m_A17(0.0); + double m_A18(0.0); + double m_A19(0.0); + double m_A20(0.0); + double m_A21(0.0); + double m_A22(0.0); + double m_A23(0.0); + double m_A24(0.0); + double m_A25(0.0); + double m_A26(0.0); + double m_A27(0.0); + double m_A28(0.0); + double m_A29(0.0); + double m_A30(0.0); + double m_A31(0.0); + double m_A32(0.0); + double m_A33(0.0); + double m_A34(0.0); + double m_A35(0.0); + double m_A36(0.0); + double m_A37(0.0); + double m_A38(0.0); + double m_A39(0.0); + double m_A40(0.0); + double m_A41(0.0); + double m_A42(0.0); + double m_A43(0.0); + double m_A44(0.0); + double m_A45(0.0); + double m_A46(0.0); + double m_A47(0.0); + double m_A48(0.0); + double m_A49(0.0); + double m_A50(0.0); + double m_A51(0.0); + double m_A52(0.0); + double m_A53(0.0); + double m_A54(0.0); + double m_A55(0.0); + double m_A56(0.0); + double m_A57(0.0); + double m_A58(0.0); + double m_A59(0.0); + double m_A60(0.0); + double m_A61(0.0); + double m_A62(0.0); + double m_A63(0.0); + double m_A64(0.0); + double m_A65(0.0); + double m_A66(0.0); + double m_A67(0.0); + double m_A68(0.0); + double m_A69(0.0); + double m_A70(0.0); + double m_A71(0.0); + double m_A72(0.0); + double m_A73(0.0); + double m_A74(0.0); + double m_A75(0.0); + double m_A76(0.0); + double m_A77(0.0); + double m_A78(0.0); + double m_A79(0.0); + double m_A80(0.0); + double m_A81(0.0); + double m_A82(0.0); + double m_A83(0.0); + double m_A84(0.0); + double m_A85(0.0); + double m_A86(0.0); + double m_A87(0.0); + double m_A88(0.0); + double m_A89(0.0); + double m_A90(0.0); + double m_A91(0.0); + double m_A92(0.0); + double m_A93(0.0); + double m_A94(0.0); + double m_A95(0.0); + double m_A96(0.0); + m_A0 += gt[0]; + m_A0 += gt[1]; + m_A0 += gt[2]; + m_A1 += go[0]; + double RHS0 = Idr[0]; + RHS0 += Idr[1]; + RHS0 += Idr[2]; + RHS0 -= go[1] * *cnV[1]; + RHS0 -= go[2] * *cnV[2]; + m_A2 += gt[3]; + m_A2 += gt[4]; + m_A2 += gt[5]; + m_A2 += gt[6]; + m_A2 += gt[7]; + m_A2 += gt[8]; + m_A2 += gt[9]; + m_A3 += go[3]; + m_A3 += go[4]; + m_A3 += go[5]; + m_A4 += go[6]; + m_A4 += go[7]; + double RHS1 = Idr[3]; + RHS1 += Idr[4]; + RHS1 += Idr[5]; + RHS1 += Idr[6]; + RHS1 += Idr[7]; + RHS1 += Idr[8]; + RHS1 += Idr[9]; + RHS1 -= go[8] * *cnV[8]; + RHS1 -= go[9] * *cnV[9]; + m_A5 += gt[10]; + m_A5 += gt[11]; + m_A6 += go[10]; + double RHS2 = Idr[10]; + RHS2 += Idr[11]; + RHS2 -= go[11] * *cnV[11]; + m_A7 += gt[12]; + m_A7 += gt[13]; + m_A8 += go[12]; + double RHS3 = Idr[12]; + RHS3 += Idr[13]; + RHS3 -= go[13] * *cnV[13]; + m_A9 += gt[14]; + m_A9 += gt[15]; + m_A9 += gt[16]; + m_A9 += gt[17]; + m_A9 += gt[18]; + m_A11 += go[14]; + m_A10 += go[15]; + m_A12 += go[16]; + m_A12 += go[17]; + double RHS4 = Idr[14]; + RHS4 += Idr[15]; + RHS4 += Idr[16]; + RHS4 += Idr[17]; + RHS4 += Idr[18]; + RHS4 -= go[18] * *cnV[18]; + m_A13 += gt[19]; + m_A13 += gt[20]; + m_A13 += gt[21]; + m_A14 += go[19]; + double RHS5 = Idr[19]; + RHS5 += Idr[20]; + RHS5 += Idr[21]; + RHS5 -= go[20] * *cnV[20]; + RHS5 -= go[21] * *cnV[21]; + m_A15 += gt[22]; + m_A15 += gt[23]; + m_A15 += gt[24]; + m_A15 += gt[25]; + m_A15 += gt[26]; + m_A15 += gt[27]; + m_A15 += gt[28]; + m_A16 += go[22]; + m_A16 += go[23]; + m_A16 += go[24]; + m_A17 += go[25]; + m_A17 += go[26]; + double RHS6 = Idr[22]; + RHS6 += Idr[23]; + RHS6 += Idr[24]; + RHS6 += Idr[25]; + RHS6 += Idr[26]; + RHS6 += Idr[27]; + RHS6 += Idr[28]; + RHS6 -= go[27] * *cnV[27]; + RHS6 -= go[28] * *cnV[28]; + m_A18 += gt[29]; + m_A18 += gt[30]; + m_A19 += go[29]; + double RHS7 = Idr[29]; + RHS7 += Idr[30]; + RHS7 -= go[30] * *cnV[30]; + m_A20 += gt[31]; + m_A20 += gt[32]; + m_A21 += go[31]; + double RHS8 = Idr[31]; + RHS8 += Idr[32]; + RHS8 -= go[32] * *cnV[32]; + m_A22 += gt[33]; + m_A22 += gt[34]; + m_A22 += gt[35]; + m_A23 += go[33]; + double RHS9 = Idr[33]; + RHS9 += Idr[34]; + RHS9 += Idr[35]; + RHS9 -= go[34] * *cnV[34]; + RHS9 -= go[35] * *cnV[35]; + m_A24 += gt[36]; + m_A24 += gt[37]; + m_A24 += gt[38]; + m_A24 += gt[39]; + m_A24 += gt[40]; + m_A25 += go[36]; + m_A26 += go[37]; + m_A27 += go[38]; + m_A27 += go[39]; + double RHS10 = Idr[36]; + RHS10 += Idr[37]; + RHS10 += Idr[38]; + RHS10 += Idr[39]; + RHS10 += Idr[40]; + RHS10 -= go[40] * *cnV[40]; + m_A29 += gt[41]; + m_A29 += gt[42]; + m_A29 += gt[43]; + m_A29 += gt[44]; + m_A29 += gt[45]; + m_A29 += gt[46]; + m_A29 += gt[47]; + m_A30 += go[41]; + m_A30 += go[42]; + m_A30 += go[43]; + m_A28 += go[44]; + double RHS11 = Idr[41]; + RHS11 += Idr[42]; + RHS11 += Idr[43]; + RHS11 += Idr[44]; + RHS11 += Idr[45]; + RHS11 += Idr[46]; + RHS11 += Idr[47]; + RHS11 -= go[45] * *cnV[45]; + RHS11 -= go[46] * *cnV[46]; + RHS11 -= go[47] * *cnV[47]; + m_A31 += gt[48]; + m_A31 += gt[49]; + m_A31 += gt[50]; + m_A31 += gt[51]; + m_A31 += gt[52]; + m_A31 += gt[53]; + m_A32 += go[48]; + m_A32 += go[49]; + m_A32 += go[50]; + m_A33 += go[51]; + m_A34 += go[52]; + double RHS12 = Idr[48]; + RHS12 += Idr[49]; + RHS12 += Idr[50]; + RHS12 += Idr[51]; + RHS12 += Idr[52]; + RHS12 += Idr[53]; + RHS12 -= go[53] * *cnV[53]; + m_A37 += gt[54]; + m_A37 += gt[55]; + m_A37 += gt[56]; + m_A37 += gt[57]; + m_A37 += gt[58]; + m_A36 += go[54]; + m_A35 += go[55]; + m_A35 += go[56]; + m_A35 += go[57]; + m_A39 += go[58]; + double RHS13 = Idr[54]; + RHS13 += Idr[55]; + RHS13 += Idr[56]; + RHS13 += Idr[57]; + RHS13 += Idr[58]; + m_A42 += gt[59]; + m_A42 += gt[60]; + m_A42 += gt[61]; + m_A42 += gt[62]; + m_A42 += gt[63]; + m_A42 += gt[64]; + m_A42 += gt[65]; + m_A43 += go[59]; + m_A43 += go[60]; + m_A43 += go[61]; + m_A41 += go[62]; + double RHS14 = Idr[59]; + RHS14 += Idr[60]; + RHS14 += Idr[61]; + RHS14 += Idr[62]; + RHS14 += Idr[63]; + RHS14 += Idr[64]; + RHS14 += Idr[65]; + RHS14 -= go[63] * *cnV[63]; + RHS14 -= go[64] * *cnV[64]; + RHS14 -= go[65] * *cnV[65]; + m_A47 += gt[66]; + m_A47 += gt[67]; + m_A47 += gt[68]; + m_A47 += gt[69]; + m_A47 += gt[70]; + m_A44 += go[66]; + m_A45 += go[67]; + m_A49 += go[68]; + m_A49 += go[69]; + double RHS15 = Idr[66]; + RHS15 += Idr[67]; + RHS15 += Idr[68]; + RHS15 += Idr[69]; + RHS15 += Idr[70]; + RHS15 -= go[70] * *cnV[70]; + m_A52 += gt[71]; + m_A52 += gt[72]; + m_A52 += gt[73]; + m_A52 += gt[74]; + m_A52 += gt[75]; + m_A52 += gt[76]; + m_A52 += gt[77]; + m_A51 += go[71]; + m_A51 += go[72]; + m_A51 += go[73]; + m_A50 += go[74]; + double RHS16 = Idr[71]; + RHS16 += Idr[72]; + RHS16 += Idr[73]; + RHS16 += Idr[74]; + RHS16 += Idr[75]; + RHS16 += Idr[76]; + RHS16 += Idr[77]; + RHS16 -= go[75] * *cnV[75]; + RHS16 -= go[76] * *cnV[76]; + RHS16 -= go[77] * *cnV[77]; + m_A57 += gt[78]; + m_A57 += gt[79]; + m_A57 += gt[80]; + m_A57 += gt[81]; + m_A57 += gt[82]; + m_A55 += go[78]; + m_A56 += go[79]; + m_A59 += go[80]; + m_A59 += go[81]; + double RHS17 = Idr[78]; + RHS17 += Idr[79]; + RHS17 += Idr[80]; + RHS17 += Idr[81]; + RHS17 += Idr[82]; + RHS17 -= go[82] * *cnV[82]; + m_A64 += gt[83]; + m_A64 += gt[84]; + m_A64 += gt[85]; + m_A64 += gt[86]; + m_A64 += gt[87]; + m_A64 += gt[88]; + m_A64 += gt[89]; + m_A64 += gt[90]; + m_A62 += go[83]; + m_A65 += go[84]; + m_A60 += go[85]; + m_A60 += go[86]; + m_A61 += go[87]; + m_A61 += go[88]; + m_A61 += go[89]; + double RHS18 = Idr[83]; + RHS18 += Idr[84]; + RHS18 += Idr[85]; + RHS18 += Idr[86]; + RHS18 += Idr[87]; + RHS18 += Idr[88]; + RHS18 += Idr[89]; + RHS18 += Idr[90]; + RHS18 -= go[90] * *cnV[90]; + m_A69 += gt[91]; + m_A69 += gt[92]; + m_A69 += gt[93]; + m_A69 += gt[94]; + m_A69 += gt[95]; + m_A67 += go[91]; + m_A66 += go[92]; + m_A66 += go[93]; + m_A66 += go[94]; + m_A71 += go[95]; + double RHS19 = Idr[91]; + RHS19 += Idr[92]; + RHS19 += Idr[93]; + RHS19 += Idr[94]; + RHS19 += Idr[95]; + m_A78 += gt[96]; + m_A78 += gt[97]; + m_A78 += gt[98]; + m_A78 += gt[99]; + m_A78 += gt[100]; + m_A78 += gt[101]; + m_A78 += gt[102]; + m_A74 += go[96]; + m_A76 += go[97]; + m_A76 += go[98]; + m_A73 += go[99]; + m_A73 += go[100]; + m_A72 += go[101]; + m_A80 += go[102]; + double RHS20 = Idr[96]; + RHS20 += Idr[97]; + RHS20 += Idr[98]; + RHS20 += Idr[99]; + RHS20 += Idr[100]; + RHS20 += Idr[101]; + RHS20 += Idr[102]; + m_A89 += gt[103]; + m_A89 += gt[104]; + m_A89 += gt[105]; + m_A89 += gt[106]; + m_A89 += gt[107]; + m_A89 += gt[108]; + m_A89 += gt[109]; + m_A83 += go[103]; + m_A85 += go[104]; + m_A85 += go[105]; + m_A82 += go[106]; + m_A82 += go[107]; + m_A81 += go[108]; + m_A87 += go[109]; + double RHS21 = Idr[103]; + RHS21 += Idr[104]; + RHS21 += Idr[105]; + RHS21 += Idr[106]; + RHS21 += Idr[107]; + RHS21 += Idr[108]; + RHS21 += Idr[109]; + m_A96 += gt[110]; + m_A96 += gt[111]; + m_A96 += gt[112]; + m_A96 += gt[113]; + m_A96 += gt[114]; + m_A96 += gt[115]; + m_A96 += gt[116]; + m_A96 += gt[117]; + m_A93 += go[110]; + m_A94 += go[111]; + m_A91 += go[112]; + m_A91 += go[113]; + m_A92 += go[114]; + m_A92 += go[115]; + m_A92 += go[116]; + double RHS22 = Idr[110]; + RHS22 += Idr[111]; + RHS22 += Idr[112]; + RHS22 += Idr[113]; + RHS22 += Idr[114]; + RHS22 += Idr[115]; + RHS22 += Idr[116]; + RHS22 += Idr[117]; + RHS22 -= go[117] * *cnV[117]; + const double f0 = 1.0 / m_A0; + const double f0_11 = -f0 * m_A28; + m_A29 += m_A1 * f0_11; + RHS11 += f0_11 * RHS0; + const double f1 = 1.0 / m_A2; + const double f1_13 = -f1 * m_A35; + m_A37 += m_A3 * f1_13; + m_A39 += m_A4 * f1_13; + RHS13 += f1_13 * RHS1; + const double f1_18 = -f1 * m_A60; + m_A62 += m_A3 * f1_18; + m_A64 += m_A4 * f1_18; + RHS18 += f1_18 * RHS1; + const double f2 = 1.0 / m_A5; + const double f2_21 = -f2 * m_A81; + m_A89 += m_A6 * f2_21; + RHS21 += f2_21 * RHS2; + const double f3 = 1.0 / m_A7; + const double f3_15 = -f3 * m_A44; + m_A47 += m_A8 * f3_15; + RHS15 += f3_15 * RHS3; + const double f4 = 1.0 / m_A9; + const double f4_13 = -f4 * m_A36; + m_A37 += m_A10 * f4_13; + m_A38 += m_A11 * f4_13; + m_A40 += m_A12 * f4_13; + RHS13 += f4_13 * RHS4; + const double f4_15 = -f4 * m_A45; + m_A46 += m_A10 * f4_15; + m_A47 += m_A11 * f4_15; + m_A49 += m_A12 * f4_15; + RHS15 += f4_15 * RHS4; + const double f4_21 = -f4 * m_A82; + m_A84 += m_A10 * f4_21; + m_A85 += m_A11 * f4_21; + m_A89 += m_A12 * f4_21; + RHS21 += f4_21 * RHS4; + const double f5 = 1.0 / m_A13; + const double f5_14 = -f5 * m_A41; + m_A42 += m_A14 * f5_14; + RHS14 += f5_14 * RHS5; + const double f6 = 1.0 / m_A15; + const double f6_19 = -f6 * m_A66; + m_A69 += m_A16 * f6_19; + m_A71 += m_A17 * f6_19; + RHS19 += f6_19 * RHS6; + const double f6_22 = -f6 * m_A91; + m_A93 += m_A16 * f6_22; + m_A96 += m_A17 * f6_22; + RHS22 += f6_22 * RHS6; + const double f7 = 1.0 / m_A18; + const double f7_17 = -f7 * m_A55; + m_A57 += m_A19 * f7_17; + RHS17 += f7_17 * RHS7; + const double f8 = 1.0 / m_A20; + const double f8_20 = -f8 * m_A72; + m_A78 += m_A21 * f8_20; + RHS20 += f8_20 * RHS8; + const double f9 = 1.0 / m_A22; + const double f9_16 = -f9 * m_A50; + m_A52 += m_A23 * f9_16; + RHS16 += f9_16 * RHS9; + const double f10 = 1.0 / m_A24; + const double f10_17 = -f10 * m_A56; + m_A57 += m_A25 * f10_17; + m_A58 += m_A26 * f10_17; + m_A59 += m_A27 * f10_17; + RHS17 += f10_17 * RHS10; + const double f10_19 = -f10 * m_A67; + m_A68 += m_A25 * f10_19; + m_A69 += m_A26 * f10_19; + m_A70 += m_A27 * f10_19; + RHS19 += f10_19 * RHS10; + const double f10_20 = -f10 * m_A73; + m_A76 += m_A25 * f10_20; + m_A77 += m_A26 * f10_20; + m_A78 += m_A27 * f10_20; + RHS20 += f10_20 * RHS10; + const double f11 = 1.0 / m_A29; + const double f11_18 = -f11 * m_A61; + m_A64 += m_A30 * f11_18; + RHS18 += f11_18 * RHS11; + const double f12 = 1.0 / m_A31; + const double f12_16 = -f12 * m_A51; + m_A52 += m_A32 * f12_16; + m_A53 += m_A33 * f12_16; + m_A54 += m_A34 * f12_16; + RHS16 += f12_16 * RHS12; + const double f12_20 = -f12 * m_A74; + m_A75 += m_A32 * f12_20; + m_A78 += m_A33 * f12_20; + m_A79 += m_A34 * f12_20; + RHS20 += f12_20 * RHS12; + const double f12_21 = -f12 * m_A83; + m_A86 += m_A32 * f12_21; + m_A88 += m_A33 * f12_21; + m_A89 += m_A34 * f12_21; + RHS21 += f12_21 * RHS12; + const double f13 = 1.0 / m_A37; + const double f13_15 = -f13 * m_A46; + m_A47 += m_A38 * f13_15; + m_A48 += m_A39 * f13_15; + m_A49 += m_A40 * f13_15; + RHS15 += f13_15 * RHS13; + const double f13_18 = -f13 * m_A62; + m_A63 += m_A38 * f13_18; + m_A64 += m_A39 * f13_18; + m_A65 += m_A40 * f13_18; + RHS18 += f13_18 * RHS13; + const double f13_21 = -f13 * m_A84; + m_A85 += m_A38 * f13_21; + m_A87 += m_A39 * f13_21; + m_A89 += m_A40 * f13_21; + RHS21 += f13_21 * RHS13; + const double f14 = 1.0 / m_A42; + const double f14_22 = -f14 * m_A92; + m_A96 += m_A43 * f14_22; + RHS22 += f14_22 * RHS14; + const double f15 = 1.0 / m_A47; + const double f15_18 = -f15 * m_A63; + m_A64 += m_A48 * f15_18; + m_A65 += m_A49 * f15_18; + RHS18 += f15_18 * RHS15; + const double f15_21 = -f15 * m_A85; + m_A87 += m_A48 * f15_21; + m_A89 += m_A49 * f15_21; + RHS21 += f15_21 * RHS15; + const double f16 = 1.0 / m_A52; + const double f16_20 = -f16 * m_A75; + m_A78 += m_A53 * f16_20; + m_A79 += m_A54 * f16_20; + RHS20 += f16_20 * RHS16; + const double f16_21 = -f16 * m_A86; + m_A88 += m_A53 * f16_21; + m_A89 += m_A54 * f16_21; + RHS21 += f16_21 * RHS16; + const double f17 = 1.0 / m_A57; + const double f17_19 = -f17 * m_A68; + m_A69 += m_A58 * f17_19; + m_A70 += m_A59 * f17_19; + RHS19 += f17_19 * RHS17; + const double f17_20 = -f17 * m_A76; + m_A77 += m_A58 * f17_20; + m_A78 += m_A59 * f17_20; + RHS20 += f17_20 * RHS17; + const double f18 = 1.0 / m_A64; + const double f18_21 = -f18 * m_A87; + m_A89 += m_A65 * f18_21; + RHS21 += f18_21 * RHS18; + const double f19 = 1.0 / m_A69; + const double f19_20 = -f19 * m_A77; + m_A78 += m_A70 * f19_20; + m_A80 += m_A71 * f19_20; + RHS20 += f19_20 * RHS19; + const double f19_22 = -f19 * m_A93; + m_A94 += m_A70 * f19_22; + m_A96 += m_A71 * f19_22; + RHS22 += f19_22 * RHS19; + const double f20 = 1.0 / m_A78; + const double f20_21 = -f20 * m_A88; + m_A89 += m_A79 * f20_21; + m_A90 += m_A80 * f20_21; + RHS21 += f20_21 * RHS20; + const double f20_22 = -f20 * m_A94; + m_A95 += m_A79 * f20_22; + m_A96 += m_A80 * f20_22; + RHS22 += f20_22 * RHS20; + const double f21 = 1.0 / m_A89; + const double f21_22 = -f21 * m_A95; + m_A96 += m_A90 * f21_22; + RHS22 += f21_22 * RHS21; + V[22] = RHS22 / m_A96; + double tmp21 = 0.0; + tmp21 += m_A90 * V[22]; + V[21] = (RHS21 - tmp21) / m_A89; + double tmp20 = 0.0; + tmp20 += m_A79 * V[21]; + tmp20 += m_A80 * V[22]; + V[20] = (RHS20 - tmp20) / m_A78; + double tmp19 = 0.0; + tmp19 += m_A70 * V[20]; + tmp19 += m_A71 * V[22]; + V[19] = (RHS19 - tmp19) / m_A69; + double tmp18 = 0.0; + tmp18 += m_A65 * V[21]; + V[18] = (RHS18 - tmp18) / m_A64; + double tmp17 = 0.0; + tmp17 += m_A58 * V[19]; + tmp17 += m_A59 * V[20]; + V[17] = (RHS17 - tmp17) / m_A57; + double tmp16 = 0.0; + tmp16 += m_A53 * V[20]; + tmp16 += m_A54 * V[21]; + V[16] = (RHS16 - tmp16) / m_A52; + double tmp15 = 0.0; + tmp15 += m_A48 * V[18]; + tmp15 += m_A49 * V[21]; + V[15] = (RHS15 - tmp15) / m_A47; + double tmp14 = 0.0; + tmp14 += m_A43 * V[22]; + V[14] = (RHS14 - tmp14) / m_A42; + double tmp13 = 0.0; + tmp13 += m_A38 * V[15]; + tmp13 += m_A39 * V[18]; + tmp13 += m_A40 * V[21]; + V[13] = (RHS13 - tmp13) / m_A37; + double tmp12 = 0.0; + tmp12 += m_A32 * V[16]; + tmp12 += m_A33 * V[20]; + tmp12 += m_A34 * V[21]; + V[12] = (RHS12 - tmp12) / m_A31; + double tmp11 = 0.0; + tmp11 += m_A30 * V[18]; + V[11] = (RHS11 - tmp11) / m_A29; + double tmp10 = 0.0; + tmp10 += m_A25 * V[17]; + tmp10 += m_A26 * V[19]; + tmp10 += m_A27 * V[20]; + V[10] = (RHS10 - tmp10) / m_A24; + double tmp9 = 0.0; + tmp9 += m_A23 * V[16]; + V[9] = (RHS9 - tmp9) / m_A22; + double tmp8 = 0.0; + tmp8 += m_A21 * V[20]; + V[8] = (RHS8 - tmp8) / m_A20; + double tmp7 = 0.0; + tmp7 += m_A19 * V[17]; + V[7] = (RHS7 - tmp7) / m_A18; + double tmp6 = 0.0; + tmp6 += m_A16 * V[19]; + tmp6 += m_A17 * V[22]; + V[6] = (RHS6 - tmp6) / m_A15; + double tmp5 = 0.0; + tmp5 += m_A14 * V[14]; + V[5] = (RHS5 - tmp5) / m_A13; + double tmp4 = 0.0; + tmp4 += m_A10 * V[13]; + tmp4 += m_A11 * V[15]; + tmp4 += m_A12 * V[21]; + V[4] = (RHS4 - tmp4) / m_A9; + double tmp3 = 0.0; + tmp3 += m_A8 * V[15]; + V[3] = (RHS3 - tmp3) / m_A7; + double tmp2 = 0.0; + tmp2 += m_A6 * V[21]; + V[2] = (RHS2 - tmp2) / m_A5; + double tmp1 = 0.0; + tmp1 += m_A3 * V[13]; + tmp1 += m_A4 * V[18]; + V[1] = (RHS1 - tmp1) / m_A2; + double tmp0 = 0.0; + tmp0 += m_A1 * V[11]; + V[0] = (RHS0 - tmp0) / m_A0; +} + // 280zzzap static void nl_gcr_24643c159711f292_95_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) @@ -20170,6 +21122,101 @@ static void nl_gcr_546396f65ce48700_12_double_double(double * __restrict V, cons } // brdrline +static void nl_gcr_576f391d0418cb76_13_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) + +{ + + plib::unused_var(cnV); + double m_A0(0.0); + double m_A1(0.0); + double m_A2(0.0); + double m_A3(0.0); + double m_A4(0.0); + double m_A5(0.0); + double m_A6(0.0); + double m_A7(0.0); + double m_A8(0.0); + double m_A9(0.0); + double m_A10(0.0); + double m_A11(0.0); + double m_A12(0.0); + m_A0 += gt[0]; + m_A0 += gt[1]; + m_A0 += gt[2]; + m_A1 += go[0]; + double RHS0 = Idr[0]; + RHS0 += Idr[1]; + RHS0 += Idr[2]; + RHS0 -= go[1] * *cnV[1]; + RHS0 -= go[2] * *cnV[2]; + m_A2 += gt[3]; + m_A2 += gt[4]; + m_A2 += gt[5]; + m_A3 += go[3]; + double RHS1 = Idr[3]; + RHS1 += Idr[4]; + RHS1 += Idr[5]; + RHS1 -= go[4] * *cnV[4]; + RHS1 -= go[5] * *cnV[5]; + m_A4 += gt[6]; + m_A4 += gt[7]; + m_A4 += gt[8]; + m_A5 += go[6]; + double RHS2 = Idr[6]; + RHS2 += Idr[7]; + RHS2 += Idr[8]; + RHS2 -= go[7] * *cnV[7]; + RHS2 -= go[8] * *cnV[8]; + m_A6 += gt[9]; + m_A7 += go[9]; + double RHS3 = Idr[9]; + m_A12 += gt[10]; + m_A12 += gt[11]; + m_A12 += gt[12]; + m_A12 += gt[13]; + m_A12 += gt[14]; + m_A11 += go[10]; + m_A10 += go[11]; + m_A9 += go[12]; + m_A8 += go[13]; + double RHS4 = Idr[10]; + RHS4 += Idr[11]; + RHS4 += Idr[12]; + RHS4 += Idr[13]; + RHS4 += Idr[14]; + RHS4 -= go[14] * *cnV[14]; + const double f0 = 1.0 / m_A0; + const double f0_4 = -f0 * m_A8; + m_A12 += m_A1 * f0_4; + RHS4 += f0_4 * RHS0; + const double f1 = 1.0 / m_A2; + const double f1_4 = -f1 * m_A9; + m_A12 += m_A3 * f1_4; + RHS4 += f1_4 * RHS1; + const double f2 = 1.0 / m_A4; + const double f2_4 = -f2 * m_A10; + m_A12 += m_A5 * f2_4; + RHS4 += f2_4 * RHS2; + const double f3 = 1.0 / m_A6; + const double f3_4 = -f3 * m_A11; + m_A12 += m_A7 * f3_4; + RHS4 += f3_4 * RHS3; + V[4] = RHS4 / m_A12; + double tmp3 = 0.0; + tmp3 += m_A7 * V[4]; + V[3] = (RHS3 - tmp3) / m_A6; + double tmp2 = 0.0; + tmp2 += m_A5 * V[4]; + V[2] = (RHS2 - tmp2) / m_A4; + double tmp1 = 0.0; + tmp1 += m_A3 * V[4]; + V[1] = (RHS1 - tmp1) / m_A2; + double tmp0 = 0.0; + tmp0 += m_A1 * V[4]; + V[0] = (RHS0 - tmp0) / m_A0; +} + +// brdrline static void nl_gcr_59cb6bf7cb9d17dc_7_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) { @@ -55697,6 +56744,153 @@ static void nl_gcr_c4cec7aed23b7b94_23_double_double(double * __restrict V, cons V[0] = (RHS0 - tmp0) / m_A0; } +// brdrline +static void nl_gcr_c5b523ed6a9677e6_21_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) + +{ + + plib::unused_var(cnV); + double m_A0(0.0); + double m_A1(0.0); + double m_A2(0.0); + double m_A3(0.0); + double m_A4(0.0); + double m_A5(0.0); + double m_A6(0.0); + double m_A7(0.0); + double m_A8(0.0); + double m_A9(0.0); + double m_A10(0.0); + double m_A11(0.0); + double m_A12(0.0); + double m_A13(0.0); + double m_A14(0.0); + double m_A15(0.0); + double m_A16(0.0); + double m_A17(0.0); + double m_A18(0.0); + double m_A19(0.0); + double m_A20(0.0); + m_A0 += gt[0]; + m_A0 += gt[1]; + m_A1 += go[0]; + double RHS0 = Idr[0]; + RHS0 += Idr[1]; + RHS0 -= go[1] * *cnV[1]; + m_A2 += gt[2]; + m_A2 += gt[3]; + m_A4 += go[2]; + m_A3 += go[3]; + double RHS1 = Idr[2]; + RHS1 += Idr[3]; + m_A5 += gt[4]; + m_A5 += gt[5]; + m_A5 += gt[6]; + m_A5 += gt[7]; + m_A5 += gt[8]; + m_A5 += gt[9]; + m_A5 += gt[10]; + m_A6 += go[4]; + double RHS2 = Idr[4]; + RHS2 += Idr[5]; + RHS2 += Idr[6]; + RHS2 += Idr[7]; + RHS2 += Idr[8]; + RHS2 += Idr[9]; + RHS2 += Idr[10]; + RHS2 -= go[5] * *cnV[5]; + RHS2 -= go[6] * *cnV[6]; + RHS2 -= go[7] * *cnV[7]; + RHS2 -= go[8] * *cnV[8]; + RHS2 -= go[9] * *cnV[9]; + RHS2 -= go[10] * *cnV[10]; + m_A7 += gt[11]; + m_A8 += go[11]; + double RHS3 = Idr[11]; + m_A11 += gt[12]; + m_A11 += gt[13]; + m_A11 += gt[14]; + m_A10 += go[12]; + m_A9 += go[13]; + double RHS4 = Idr[12]; + RHS4 += Idr[13]; + RHS4 += Idr[14]; + RHS4 -= go[14] * *cnV[14]; + m_A15 += gt[15]; + m_A15 += gt[16]; + m_A15 += gt[17]; + m_A16 += go[15]; + m_A13 += go[16]; + double RHS5 = Idr[15]; + RHS5 += Idr[16]; + RHS5 += Idr[17]; + RHS5 -= go[17] * *cnV[17]; + m_A20 += gt[18]; + m_A20 += gt[19]; + m_A20 += gt[20]; + m_A20 += gt[21]; + m_A20 += gt[22]; + m_A17 += go[18]; + m_A18 += go[19]; + m_A19 += go[20]; + double RHS6 = Idr[18]; + RHS6 += Idr[19]; + RHS6 += Idr[20]; + RHS6 += Idr[21]; + RHS6 += Idr[22]; + RHS6 -= go[21] * *cnV[21]; + RHS6 -= go[22] * *cnV[22]; + const double f0 = 1.0 / m_A0; + const double f0_4 = -f0 * m_A9; + m_A11 += m_A1 * f0_4; + RHS4 += f0_4 * RHS0; + const double f1 = 1.0 / m_A2; + const double f1_4 = -f1 * m_A10; + m_A11 += m_A3 * f1_4; + m_A12 += m_A4 * f1_4; + RHS4 += f1_4 * RHS1; + const double f1_5 = -f1 * m_A13; + m_A14 += m_A3 * f1_5; + m_A15 += m_A4 * f1_5; + RHS5 += f1_5 * RHS1; + const double f2 = 1.0 / m_A5; + const double f2_6 = -f2 * m_A17; + m_A19 += m_A6 * f2_6; + RHS6 += f2_6 * RHS2; + const double f3 = 1.0 / m_A7; + const double f3_6 = -f3 * m_A18; + m_A20 += m_A8 * f3_6; + RHS6 += f3_6 * RHS3; + const double f4 = 1.0 / m_A11; + const double f4_5 = -f4 * m_A14; + m_A15 += m_A12 * f4_5; + RHS5 += f4_5 * RHS4; + const double f5 = 1.0 / m_A15; + const double f5_6 = -f5 * m_A19; + m_A20 += m_A16 * f5_6; + RHS6 += f5_6 * RHS5; + V[6] = RHS6 / m_A20; + double tmp5 = 0.0; + tmp5 += m_A16 * V[6]; + V[5] = (RHS5 - tmp5) / m_A15; + double tmp4 = 0.0; + tmp4 += m_A12 * V[5]; + V[4] = (RHS4 - tmp4) / m_A11; + double tmp3 = 0.0; + tmp3 += m_A8 * V[6]; + V[3] = (RHS3 - tmp3) / m_A7; + double tmp2 = 0.0; + tmp2 += m_A6 * V[5]; + V[2] = (RHS2 - tmp2) / m_A5; + double tmp1 = 0.0; + tmp1 += m_A3 * V[4]; + tmp1 += m_A4 * V[5]; + V[1] = (RHS1 - tmp1) / m_A2; + double tmp0 = 0.0; + tmp0 += m_A1 * V[4]; + V[0] = (RHS0 - tmp0) / m_A0; +} + // segausb static void nl_gcr_c61e08cf5e35918_84_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) @@ -56635,6 +57829,2216 @@ static void nl_gcr_c74b1a65978d7121_7_double_double(double * __restrict V, const V[0] = (RHS0 - tmp0) / m_A0; } +// headon +static void nl_gcr_c90daff82be2d052_350_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) + +{ + + plib::unused_var(cnV); + double m_A0(0.0); + double m_A1(0.0); + double m_A2(0.0); + double m_A3(0.0); + double m_A4(0.0); + double m_A5(0.0); + double m_A6(0.0); + double m_A7(0.0); + double m_A8(0.0); + double m_A9(0.0); + double m_A10(0.0); + double m_A11(0.0); + double m_A12(0.0); + double m_A13(0.0); + double m_A14(0.0); + double m_A15(0.0); + double m_A16(0.0); + double m_A17(0.0); + double m_A18(0.0); + double m_A19(0.0); + double m_A20(0.0); + double m_A21(0.0); + double m_A22(0.0); + double m_A23(0.0); + double m_A24(0.0); + double m_A25(0.0); + double m_A26(0.0); + double m_A27(0.0); + double m_A28(0.0); + double m_A29(0.0); + double m_A30(0.0); + double m_A31(0.0); + double m_A32(0.0); + double m_A33(0.0); + double m_A34(0.0); + double m_A35(0.0); + double m_A36(0.0); + double m_A37(0.0); + double m_A38(0.0); + double m_A39(0.0); + double m_A40(0.0); + double m_A41(0.0); + double m_A42(0.0); + double m_A43(0.0); + double m_A44(0.0); + double m_A45(0.0); + double m_A46(0.0); + double m_A47(0.0); + double m_A48(0.0); + double m_A49(0.0); + double m_A50(0.0); + double m_A51(0.0); + double m_A52(0.0); + double m_A53(0.0); + double m_A54(0.0); + double m_A55(0.0); + double m_A56(0.0); + double m_A57(0.0); + double m_A58(0.0); + double m_A59(0.0); + double m_A60(0.0); + double m_A61(0.0); + double m_A62(0.0); + double m_A63(0.0); + double m_A64(0.0); + double m_A65(0.0); + double m_A66(0.0); + double m_A67(0.0); + double m_A68(0.0); + double m_A69(0.0); + double m_A70(0.0); + double m_A71(0.0); + double m_A72(0.0); + double m_A73(0.0); + double m_A74(0.0); + double m_A75(0.0); + double m_A76(0.0); + double m_A77(0.0); + double m_A78(0.0); + double m_A79(0.0); + double m_A80(0.0); + double m_A81(0.0); + double m_A82(0.0); + double m_A83(0.0); + double m_A84(0.0); + double m_A85(0.0); + double m_A86(0.0); + double m_A87(0.0); + double m_A88(0.0); + double m_A89(0.0); + double m_A90(0.0); + double m_A91(0.0); + double m_A92(0.0); + double m_A93(0.0); + double m_A94(0.0); + double m_A95(0.0); + double m_A96(0.0); + double m_A97(0.0); + double m_A98(0.0); + double m_A99(0.0); + double m_A100(0.0); + double m_A101(0.0); + double m_A102(0.0); + double m_A103(0.0); + double m_A104(0.0); + double m_A105(0.0); + double m_A106(0.0); + double m_A107(0.0); + double m_A108(0.0); + double m_A109(0.0); + double m_A110(0.0); + double m_A111(0.0); + double m_A112(0.0); + double m_A113(0.0); + double m_A114(0.0); + double m_A115(0.0); + double m_A116(0.0); + double m_A117(0.0); + double m_A118(0.0); + double m_A119(0.0); + double m_A120(0.0); + double m_A121(0.0); + double m_A122(0.0); + double m_A123(0.0); + double m_A124(0.0); + double m_A125(0.0); + double m_A126(0.0); + double m_A127(0.0); + double m_A128(0.0); + double m_A129(0.0); + double m_A130(0.0); + double m_A131(0.0); + double m_A132(0.0); + double m_A133(0.0); + double m_A134(0.0); + double m_A135(0.0); + double m_A136(0.0); + double m_A137(0.0); + double m_A138(0.0); + double m_A139(0.0); + double m_A140(0.0); + double m_A141(0.0); + double m_A142(0.0); + double m_A143(0.0); + double m_A144(0.0); + double m_A145(0.0); + double m_A146(0.0); + double m_A147(0.0); + double m_A148(0.0); + double m_A149(0.0); + double m_A150(0.0); + double m_A151(0.0); + double m_A152(0.0); + double m_A153(0.0); + double m_A154(0.0); + double m_A155(0.0); + double m_A156(0.0); + double m_A157(0.0); + double m_A158(0.0); + double m_A159(0.0); + double m_A160(0.0); + double m_A161(0.0); + double m_A162(0.0); + double m_A163(0.0); + double m_A164(0.0); + double m_A165(0.0); + double m_A166(0.0); + double m_A167(0.0); + double m_A168(0.0); + double m_A169(0.0); + double m_A170(0.0); + double m_A171(0.0); + double m_A172(0.0); + double m_A173(0.0); + double m_A174(0.0); + double m_A175(0.0); + double m_A176(0.0); + double m_A177(0.0); + double m_A178(0.0); + double m_A179(0.0); + double m_A180(0.0); + double m_A181(0.0); + double m_A182(0.0); + double m_A183(0.0); + double m_A184(0.0); + double m_A185(0.0); + double m_A186(0.0); + double m_A187(0.0); + double m_A188(0.0); + double m_A189(0.0); + double m_A190(0.0); + double m_A191(0.0); + double m_A192(0.0); + double m_A193(0.0); + double m_A194(0.0); + double m_A195(0.0); + double m_A196(0.0); + double m_A197(0.0); + double m_A198(0.0); + double m_A199(0.0); + double m_A200(0.0); + double m_A201(0.0); + double m_A202(0.0); + double m_A203(0.0); + double m_A204(0.0); + double m_A205(0.0); + double m_A206(0.0); + double m_A207(0.0); + double m_A208(0.0); + double m_A209(0.0); + double m_A210(0.0); + double m_A211(0.0); + double m_A212(0.0); + double m_A213(0.0); + double m_A214(0.0); + double m_A215(0.0); + double m_A216(0.0); + double m_A217(0.0); + double m_A218(0.0); + double m_A219(0.0); + double m_A220(0.0); + double m_A221(0.0); + double m_A222(0.0); + double m_A223(0.0); + double m_A224(0.0); + double m_A225(0.0); + double m_A226(0.0); + double m_A227(0.0); + double m_A228(0.0); + double m_A229(0.0); + double m_A230(0.0); + double m_A231(0.0); + double m_A232(0.0); + double m_A233(0.0); + double m_A234(0.0); + double m_A235(0.0); + double m_A236(0.0); + double m_A237(0.0); + double m_A238(0.0); + double m_A239(0.0); + double m_A240(0.0); + double m_A241(0.0); + double m_A242(0.0); + double m_A243(0.0); + double m_A244(0.0); + double m_A245(0.0); + double m_A246(0.0); + double m_A247(0.0); + double m_A248(0.0); + double m_A249(0.0); + double m_A250(0.0); + double m_A251(0.0); + double m_A252(0.0); + double m_A253(0.0); + double m_A254(0.0); + double m_A255(0.0); + double m_A256(0.0); + double m_A257(0.0); + double m_A258(0.0); + double m_A259(0.0); + double m_A260(0.0); + double m_A261(0.0); + double m_A262(0.0); + double m_A263(0.0); + double m_A264(0.0); + double m_A265(0.0); + double m_A266(0.0); + double m_A267(0.0); + double m_A268(0.0); + double m_A269(0.0); + double m_A270(0.0); + double m_A271(0.0); + double m_A272(0.0); + double m_A273(0.0); + double m_A274(0.0); + double m_A275(0.0); + double m_A276(0.0); + double m_A277(0.0); + double m_A278(0.0); + double m_A279(0.0); + double m_A280(0.0); + double m_A281(0.0); + double m_A282(0.0); + double m_A283(0.0); + double m_A284(0.0); + double m_A285(0.0); + double m_A286(0.0); + double m_A287(0.0); + double m_A288(0.0); + double m_A289(0.0); + double m_A290(0.0); + double m_A291(0.0); + double m_A292(0.0); + double m_A293(0.0); + double m_A294(0.0); + double m_A295(0.0); + double m_A296(0.0); + double m_A297(0.0); + double m_A298(0.0); + double m_A299(0.0); + double m_A300(0.0); + double m_A301(0.0); + double m_A302(0.0); + double m_A303(0.0); + double m_A304(0.0); + double m_A305(0.0); + double m_A306(0.0); + double m_A307(0.0); + double m_A308(0.0); + double m_A309(0.0); + double m_A310(0.0); + double m_A311(0.0); + double m_A312(0.0); + double m_A313(0.0); + double m_A314(0.0); + double m_A315(0.0); + double m_A316(0.0); + double m_A317(0.0); + double m_A318(0.0); + double m_A319(0.0); + double m_A320(0.0); + double m_A321(0.0); + double m_A322(0.0); + double m_A323(0.0); + double m_A324(0.0); + double m_A325(0.0); + double m_A326(0.0); + double m_A327(0.0); + double m_A328(0.0); + double m_A329(0.0); + double m_A330(0.0); + double m_A331(0.0); + double m_A332(0.0); + double m_A333(0.0); + double m_A334(0.0); + double m_A335(0.0); + double m_A336(0.0); + double m_A337(0.0); + double m_A338(0.0); + double m_A339(0.0); + double m_A340(0.0); + double m_A341(0.0); + double m_A342(0.0); + double m_A343(0.0); + double m_A344(0.0); + double m_A345(0.0); + double m_A346(0.0); + double m_A347(0.0); + double m_A348(0.0); + double m_A349(0.0); + m_A0 += gt[0]; + m_A0 += gt[1]; + m_A0 += gt[2]; + m_A1 += go[0]; + double RHS0 = Idr[0]; + RHS0 += Idr[1]; + RHS0 += Idr[2]; + RHS0 -= go[1] * *cnV[1]; + RHS0 -= go[2] * *cnV[2]; + m_A2 += gt[3]; + m_A2 += gt[4]; + m_A3 += go[3]; + double RHS1 = Idr[3]; + RHS1 += Idr[4]; + RHS1 -= go[4] * *cnV[4]; + m_A4 += gt[5]; + m_A4 += gt[6]; + m_A4 += gt[7]; + m_A4 += gt[8]; + m_A6 += go[5]; + m_A5 += go[6]; + m_A5 += go[7]; + double RHS2 = Idr[5]; + RHS2 += Idr[6]; + RHS2 += Idr[7]; + RHS2 += Idr[8]; + RHS2 -= go[8] * *cnV[8]; + m_A7 += gt[9]; + m_A7 += gt[10]; + m_A8 += go[9]; + double RHS3 = Idr[9]; + RHS3 += Idr[10]; + RHS3 -= go[10] * *cnV[10]; + m_A9 += gt[11]; + m_A9 += gt[12]; + m_A9 += gt[13]; + m_A10 += go[11]; + m_A11 += go[12]; + double RHS4 = Idr[11]; + RHS4 += Idr[12]; + RHS4 += Idr[13]; + RHS4 -= go[13] * *cnV[13]; + m_A12 += gt[14]; + m_A12 += gt[15]; + m_A13 += go[14]; + double RHS5 = Idr[14]; + RHS5 += Idr[15]; + RHS5 -= go[15] * *cnV[15]; + m_A14 += gt[16]; + m_A14 += gt[17]; + m_A15 += go[16]; + double RHS6 = Idr[16]; + RHS6 += Idr[17]; + RHS6 -= go[17] * *cnV[17]; + m_A16 += gt[18]; + m_A16 += gt[19]; + m_A17 += go[18]; + double RHS7 = Idr[18]; + RHS7 += Idr[19]; + RHS7 -= go[19] * *cnV[19]; + m_A18 += gt[20]; + m_A18 += gt[21]; + m_A19 += go[20]; + double RHS8 = Idr[20]; + RHS8 += Idr[21]; + RHS8 -= go[21] * *cnV[21]; + m_A20 += gt[22]; + m_A20 += gt[23]; + m_A20 += gt[24]; + m_A20 += gt[25]; + m_A21 += go[22]; + m_A22 += go[23]; + m_A23 += go[24]; + double RHS9 = Idr[22]; + RHS9 += Idr[23]; + RHS9 += Idr[24]; + RHS9 += Idr[25]; + RHS9 -= go[25] * *cnV[25]; + m_A24 += gt[26]; + m_A24 += gt[27]; + m_A24 += gt[28]; + m_A24 += gt[29]; + m_A24 += gt[30]; + m_A25 += go[26]; + m_A27 += go[27]; + m_A26 += go[28]; + double RHS10 = Idr[26]; + RHS10 += Idr[27]; + RHS10 += Idr[28]; + RHS10 += Idr[29]; + RHS10 += Idr[30]; + RHS10 -= go[29] * *cnV[29]; + RHS10 -= go[30] * *cnV[30]; + m_A28 += gt[31]; + m_A28 += gt[32]; + m_A28 += gt[33]; + m_A28 += gt[34]; + m_A28 += gt[35]; + m_A28 += gt[36]; + m_A28 += gt[37]; + m_A29 += go[31]; + m_A30 += go[32]; + double RHS11 = Idr[31]; + RHS11 += Idr[32]; + RHS11 += Idr[33]; + RHS11 += Idr[34]; + RHS11 += Idr[35]; + RHS11 += Idr[36]; + RHS11 += Idr[37]; + RHS11 -= go[33] * *cnV[33]; + RHS11 -= go[34] * *cnV[34]; + RHS11 -= go[35] * *cnV[35]; + RHS11 -= go[36] * *cnV[36]; + RHS11 -= go[37] * *cnV[37]; + m_A31 += gt[38]; + m_A31 += gt[39]; + m_A31 += gt[40]; + m_A31 += gt[41]; + m_A33 += go[38]; + m_A34 += go[39]; + m_A32 += go[40]; + double RHS12 = Idr[38]; + RHS12 += Idr[39]; + RHS12 += Idr[40]; + RHS12 += Idr[41]; + RHS12 -= go[41] * *cnV[41]; + m_A35 += gt[42]; + m_A35 += gt[43]; + m_A35 += gt[44]; + m_A37 += go[42]; + m_A36 += go[43]; + double RHS13 = Idr[42]; + RHS13 += Idr[43]; + RHS13 += Idr[44]; + RHS13 -= go[44] * *cnV[44]; + m_A38 += gt[45]; + m_A38 += gt[46]; + m_A38 += gt[47]; + m_A38 += gt[48]; + m_A38 += gt[49]; + m_A38 += gt[50]; + m_A38 += gt[51]; + m_A39 += go[45]; + double RHS14 = Idr[45]; + RHS14 += Idr[46]; + RHS14 += Idr[47]; + RHS14 += Idr[48]; + RHS14 += Idr[49]; + RHS14 += Idr[50]; + RHS14 += Idr[51]; + RHS14 -= go[46] * *cnV[46]; + RHS14 -= go[47] * *cnV[47]; + RHS14 -= go[48] * *cnV[48]; + RHS14 -= go[49] * *cnV[49]; + RHS14 -= go[50] * *cnV[50]; + RHS14 -= go[51] * *cnV[51]; + m_A40 += gt[52]; + m_A40 += gt[53]; + m_A40 += gt[54]; + m_A40 += gt[55]; + m_A43 += go[52]; + m_A42 += go[53]; + m_A41 += go[54]; + double RHS15 = Idr[52]; + RHS15 += Idr[53]; + RHS15 += Idr[54]; + RHS15 += Idr[55]; + RHS15 -= go[55] * *cnV[55]; + m_A44 += gt[56]; + m_A44 += gt[57]; + m_A45 += go[56]; + m_A46 += go[57]; + double RHS16 = Idr[56]; + RHS16 += Idr[57]; + m_A47 += gt[58]; + m_A47 += gt[59]; + m_A48 += go[58]; + m_A49 += go[59]; + double RHS17 = Idr[58]; + RHS17 += Idr[59]; + m_A50 += gt[60]; + m_A50 += gt[61]; + m_A50 += gt[62]; + m_A50 += gt[63]; + m_A50 += gt[64]; + m_A50 += gt[65]; + m_A50 += gt[66]; + m_A51 += go[60]; + double RHS18 = Idr[60]; + RHS18 += Idr[61]; + RHS18 += Idr[62]; + RHS18 += Idr[63]; + RHS18 += Idr[64]; + RHS18 += Idr[65]; + RHS18 += Idr[66]; + RHS18 -= go[61] * *cnV[61]; + RHS18 -= go[62] * *cnV[62]; + RHS18 -= go[63] * *cnV[63]; + RHS18 -= go[64] * *cnV[64]; + RHS18 -= go[65] * *cnV[65]; + RHS18 -= go[66] * *cnV[66]; + m_A52 += gt[67]; + m_A52 += gt[68]; + m_A52 += gt[69]; + m_A53 += go[67]; + double RHS19 = Idr[67]; + RHS19 += Idr[68]; + RHS19 += Idr[69]; + RHS19 -= go[68] * *cnV[68]; + RHS19 -= go[69] * *cnV[69]; + m_A54 += gt[70]; + m_A54 += gt[71]; + m_A54 += gt[72]; + m_A54 += gt[73]; + m_A54 += gt[74]; + m_A56 += go[70]; + m_A55 += go[71]; + m_A55 += go[72]; + double RHS20 = Idr[70]; + RHS20 += Idr[71]; + RHS20 += Idr[72]; + RHS20 += Idr[73]; + RHS20 += Idr[74]; + RHS20 -= go[73] * *cnV[73]; + RHS20 -= go[74] * *cnV[74]; + m_A57 += gt[75]; + m_A57 += gt[76]; + m_A57 += gt[77]; + m_A57 += gt[78]; + m_A57 += gt[79]; + m_A59 += go[75]; + m_A58 += go[76]; + m_A58 += go[77]; + double RHS21 = Idr[75]; + RHS21 += Idr[76]; + RHS21 += Idr[77]; + RHS21 += Idr[78]; + RHS21 += Idr[79]; + RHS21 -= go[78] * *cnV[78]; + RHS21 -= go[79] * *cnV[79]; + m_A60 += gt[80]; + m_A60 += gt[81]; + m_A60 += gt[82]; + m_A60 += gt[83]; + m_A62 += go[80]; + m_A63 += go[81]; + m_A61 += go[82]; + double RHS22 = Idr[80]; + RHS22 += Idr[81]; + RHS22 += Idr[82]; + RHS22 += Idr[83]; + RHS22 -= go[83] * *cnV[83]; + m_A64 += gt[84]; + m_A64 += gt[85]; + m_A64 += gt[86]; + m_A65 += go[84]; + double RHS23 = Idr[84]; + RHS23 += Idr[85]; + RHS23 += Idr[86]; + RHS23 -= go[85] * *cnV[85]; + RHS23 -= go[86] * *cnV[86]; + m_A66 += gt[87]; + m_A66 += gt[88]; + m_A66 += gt[89]; + m_A66 += gt[90]; + m_A66 += gt[91]; + m_A67 += go[87]; + double RHS24 = Idr[87]; + RHS24 += Idr[88]; + RHS24 += Idr[89]; + RHS24 += Idr[90]; + RHS24 += Idr[91]; + RHS24 -= go[88] * *cnV[88]; + RHS24 -= go[89] * *cnV[89]; + RHS24 -= go[90] * *cnV[90]; + RHS24 -= go[91] * *cnV[91]; + m_A68 += gt[92]; + m_A69 += go[92]; + double RHS25 = Idr[92]; + m_A70 += gt[93]; + m_A70 += gt[94]; + m_A70 += gt[95]; + m_A70 += gt[96]; + m_A72 += go[93]; + m_A71 += go[94]; + m_A71 += go[95]; + double RHS26 = Idr[93]; + RHS26 += Idr[94]; + RHS26 += Idr[95]; + RHS26 += Idr[96]; + RHS26 -= go[96] * *cnV[96]; + m_A73 += gt[97]; + m_A73 += gt[98]; + m_A73 += gt[99]; + m_A73 += gt[100]; + m_A75 += go[97]; + m_A76 += go[98]; + m_A74 += go[99]; + double RHS27 = Idr[97]; + RHS27 += Idr[98]; + RHS27 += Idr[99]; + RHS27 += Idr[100]; + RHS27 -= go[100] * *cnV[100]; + m_A77 += gt[101]; + m_A77 += gt[102]; + m_A78 += go[101]; + double RHS28 = Idr[101]; + RHS28 += Idr[102]; + RHS28 -= go[102] * *cnV[102]; + m_A79 += gt[103]; + m_A79 += gt[104]; + m_A79 += gt[105]; + m_A80 += go[103]; + double RHS29 = Idr[103]; + RHS29 += Idr[104]; + RHS29 += Idr[105]; + RHS29 -= go[104] * *cnV[104]; + RHS29 -= go[105] * *cnV[105]; + m_A83 += gt[106]; + m_A83 += gt[107]; + m_A83 += gt[108]; + m_A83 += gt[109]; + m_A83 += gt[110]; + m_A83 += gt[111]; + m_A82 += go[106]; + m_A82 += go[107]; + m_A81 += go[108]; + double RHS30 = Idr[106]; + RHS30 += Idr[107]; + RHS30 += Idr[108]; + RHS30 += Idr[109]; + RHS30 += Idr[110]; + RHS30 += Idr[111]; + RHS30 -= go[109] * *cnV[109]; + RHS30 -= go[110] * *cnV[110]; + RHS30 -= go[111] * *cnV[111]; + m_A86 += gt[112]; + m_A86 += gt[113]; + m_A85 += go[112]; + double RHS31 = Idr[112]; + RHS31 += Idr[113]; + RHS31 -= go[113] * *cnV[113]; + m_A89 += gt[114]; + m_A89 += gt[115]; + m_A89 += gt[116]; + m_A88 += go[114]; + m_A90 += go[115]; + double RHS32 = Idr[114]; + RHS32 += Idr[115]; + RHS32 += Idr[116]; + RHS32 -= go[116] * *cnV[116]; + m_A93 += gt[117]; + m_A93 += gt[118]; + m_A92 += go[117]; + m_A91 += go[118]; + double RHS33 = Idr[117]; + RHS33 += Idr[118]; + m_A98 += gt[119]; + m_A98 += gt[120]; + m_A98 += gt[121]; + m_A97 += go[119]; + m_A99 += go[120]; + double RHS34 = Idr[119]; + RHS34 += Idr[120]; + RHS34 += Idr[121]; + RHS34 -= go[121] * *cnV[121]; + m_A101 += gt[122]; + m_A101 += gt[123]; + m_A101 += gt[124]; + m_A101 += gt[125]; + m_A101 += gt[126]; + m_A103 += go[122]; + m_A100 += go[123]; + m_A102 += go[124]; + double RHS35 = Idr[122]; + RHS35 += Idr[123]; + RHS35 += Idr[124]; + RHS35 += Idr[125]; + RHS35 += Idr[126]; + RHS35 -= go[125] * *cnV[125]; + RHS35 -= go[126] * *cnV[126]; + m_A105 += gt[127]; + m_A105 += gt[128]; + m_A105 += gt[129]; + m_A105 += gt[130]; + m_A105 += gt[131]; + m_A105 += gt[132]; + m_A105 += gt[133]; + m_A106 += go[127]; + m_A104 += go[128]; + double RHS36 = Idr[127]; + RHS36 += Idr[128]; + RHS36 += Idr[129]; + RHS36 += Idr[130]; + RHS36 += Idr[131]; + RHS36 += Idr[132]; + RHS36 += Idr[133]; + RHS36 -= go[129] * *cnV[129]; + RHS36 -= go[130] * *cnV[130]; + RHS36 -= go[131] * *cnV[131]; + RHS36 -= go[132] * *cnV[132]; + RHS36 -= go[133] * *cnV[133]; + m_A112 += gt[134]; + m_A112 += gt[135]; + m_A112 += gt[136]; + m_A110 += go[134]; + m_A109 += go[135]; + m_A115 += go[136]; + double RHS37 = Idr[134]; + RHS37 += Idr[135]; + RHS37 += Idr[136]; + m_A119 += gt[137]; + m_A119 += gt[138]; + m_A119 += gt[139]; + m_A119 += gt[140]; + m_A118 += go[137]; + m_A121 += go[138]; + m_A117 += go[139]; + m_A116 += go[140]; + double RHS38 = Idr[137]; + RHS38 += Idr[138]; + RHS38 += Idr[139]; + RHS38 += Idr[140]; + m_A123 += gt[141]; + m_A123 += gt[142]; + m_A123 += gt[143]; + m_A123 += gt[144]; + m_A123 += gt[145]; + m_A123 += gt[146]; + m_A123 += gt[147]; + m_A125 += go[141]; + m_A122 += go[142]; + double RHS39 = Idr[141]; + RHS39 += Idr[142]; + RHS39 += Idr[143]; + RHS39 += Idr[144]; + RHS39 += Idr[145]; + RHS39 += Idr[146]; + RHS39 += Idr[147]; + RHS39 -= go[143] * *cnV[143]; + RHS39 -= go[144] * *cnV[144]; + RHS39 -= go[145] * *cnV[145]; + RHS39 -= go[146] * *cnV[146]; + RHS39 -= go[147] * *cnV[147]; + m_A130 += gt[148]; + m_A130 += gt[149]; + m_A130 += gt[150]; + m_A133 += go[148]; + m_A127 += go[149]; + m_A126 += go[150]; + double RHS40 = Idr[148]; + RHS40 += Idr[149]; + RHS40 += Idr[150]; + m_A138 += gt[151]; + m_A138 += gt[152]; + m_A138 += gt[153]; + m_A138 += gt[154]; + m_A138 += gt[155]; + m_A134 += go[151]; + m_A136 += go[152]; + m_A141 += go[153]; + double RHS41 = Idr[151]; + RHS41 += Idr[152]; + RHS41 += Idr[153]; + RHS41 += Idr[154]; + RHS41 += Idr[155]; + RHS41 -= go[154] * *cnV[154]; + RHS41 -= go[155] * *cnV[155]; + m_A145 += gt[156]; + m_A145 += gt[157]; + m_A145 += gt[158]; + m_A148 += go[156]; + m_A143 += go[157]; + double RHS42 = Idr[156]; + RHS42 += Idr[157]; + RHS42 += Idr[158]; + RHS42 -= go[158] * *cnV[158]; + m_A150 += gt[159]; + m_A150 += gt[160]; + m_A150 += gt[161]; + m_A150 += gt[162]; + m_A150 += gt[163]; + m_A150 += gt[164]; + m_A151 += go[159]; + m_A151 += go[160]; + m_A149 += go[161]; + m_A149 += go[162]; + double RHS43 = Idr[159]; + RHS43 += Idr[160]; + RHS43 += Idr[161]; + RHS43 += Idr[162]; + RHS43 += Idr[163]; + RHS43 += Idr[164]; + RHS43 -= go[163] * *cnV[163]; + RHS43 -= go[164] * *cnV[164]; + m_A153 += gt[165]; + m_A153 += gt[166]; + m_A153 += gt[167]; + m_A153 += gt[168]; + m_A153 += gt[169]; + m_A153 += gt[170]; + m_A154 += go[165]; + m_A154 += go[166]; + m_A152 += go[167]; + m_A152 += go[168]; + double RHS44 = Idr[165]; + RHS44 += Idr[166]; + RHS44 += Idr[167]; + RHS44 += Idr[168]; + RHS44 += Idr[169]; + RHS44 += Idr[170]; + RHS44 -= go[169] * *cnV[169]; + RHS44 -= go[170] * *cnV[170]; + m_A158 += gt[171]; + m_A158 += gt[172]; + m_A155 += go[171]; + m_A156 += go[172]; + double RHS45 = Idr[171]; + RHS45 += Idr[172]; + m_A164 += gt[173]; + m_A164 += gt[174]; + m_A164 += gt[175]; + m_A164 += gt[176]; + m_A164 += gt[177]; + m_A164 += gt[178]; + m_A163 += go[173]; + m_A163 += go[174]; + m_A162 += go[175]; + double RHS46 = Idr[173]; + RHS46 += Idr[174]; + RHS46 += Idr[175]; + RHS46 += Idr[176]; + RHS46 += Idr[177]; + RHS46 += Idr[178]; + RHS46 -= go[176] * *cnV[176]; + RHS46 -= go[177] * *cnV[177]; + RHS46 -= go[178] * *cnV[178]; + m_A172 += gt[179]; + m_A172 += gt[180]; + m_A172 += gt[181]; + m_A172 += gt[182]; + m_A172 += gt[183]; + m_A166 += go[179]; + m_A167 += go[180]; + m_A168 += go[181]; + double RHS47 = Idr[179]; + RHS47 += Idr[180]; + RHS47 += Idr[181]; + RHS47 += Idr[182]; + RHS47 += Idr[183]; + RHS47 -= go[182] * *cnV[182]; + RHS47 -= go[183] * *cnV[183]; + m_A181 += gt[184]; + m_A181 += gt[185]; + m_A181 += gt[186]; + m_A181 += gt[187]; + m_A181 += gt[188]; + m_A182 += go[184]; + m_A178 += go[185]; + m_A179 += go[186]; + double RHS48 = Idr[184]; + RHS48 += Idr[185]; + RHS48 += Idr[186]; + RHS48 += Idr[187]; + RHS48 += Idr[188]; + RHS48 -= go[187] * *cnV[187]; + RHS48 -= go[188] * *cnV[188]; + m_A183 += gt[189]; + m_A183 += gt[190]; + m_A183 += gt[191]; + m_A184 += go[189]; + double RHS49 = Idr[189]; + RHS49 += Idr[190]; + RHS49 += Idr[191]; + RHS49 -= go[190] * *cnV[190]; + RHS49 -= go[191] * *cnV[191]; + m_A192 += gt[192]; + m_A192 += gt[193]; + m_A192 += gt[194]; + m_A192 += gt[195]; + m_A192 += gt[196]; + m_A191 += go[192]; + m_A186 += go[193]; + m_A190 += go[194]; + m_A189 += go[195]; + m_A185 += go[196]; + double RHS50 = Idr[192]; + RHS50 += Idr[193]; + RHS50 += Idr[194]; + RHS50 += Idr[195]; + RHS50 += Idr[196]; + m_A201 += gt[197]; + m_A201 += gt[198]; + m_A201 += gt[199]; + m_A201 += gt[200]; + m_A197 += go[197]; + m_A195 += go[198]; + m_A194 += go[199]; + m_A196 += go[200]; + double RHS51 = Idr[197]; + RHS51 += Idr[198]; + RHS51 += Idr[199]; + RHS51 += Idr[200]; + m_A214 += gt[201]; + m_A214 += gt[202]; + m_A214 += gt[203]; + m_A214 += gt[204]; + m_A214 += gt[205]; + m_A214 += gt[206]; + m_A209 += go[201]; + m_A207 += go[202]; + m_A206 += go[203]; + m_A210 += go[204]; + double RHS52 = Idr[201]; + RHS52 += Idr[202]; + RHS52 += Idr[203]; + RHS52 += Idr[204]; + RHS52 += Idr[205]; + RHS52 += Idr[206]; + RHS52 -= go[205] * *cnV[205]; + RHS52 -= go[206] * *cnV[206]; + m_A225 += gt[207]; + m_A225 += gt[208]; + m_A225 += gt[209]; + m_A225 += gt[210]; + m_A225 += gt[211]; + m_A231 += go[207]; + m_A231 += go[208]; + m_A220 += go[209]; + m_A229 += go[210]; + m_A221 += go[211]; + double RHS53 = Idr[207]; + RHS53 += Idr[208]; + RHS53 += Idr[209]; + RHS53 += Idr[210]; + RHS53 += Idr[211]; + m_A240 += gt[212]; + m_A240 += gt[213]; + m_A240 += gt[214]; + m_A234 += go[212]; + m_A232 += go[213]; + m_A233 += go[214]; + double RHS54 = Idr[212]; + RHS54 += Idr[213]; + RHS54 += Idr[214]; + m_A257 += gt[215]; + m_A257 += gt[216]; + m_A257 += gt[217]; + m_A257 += gt[218]; + m_A257 += gt[219]; + m_A257 += gt[220]; + m_A246 += go[215]; + m_A248 += go[216]; + m_A250 += go[217]; + m_A247 += go[218]; + double RHS55 = Idr[215]; + RHS55 += Idr[216]; + RHS55 += Idr[217]; + RHS55 += Idr[218]; + RHS55 += Idr[219]; + RHS55 += Idr[220]; + RHS55 -= go[219] * *cnV[219]; + RHS55 -= go[220] * *cnV[220]; + m_A270 += gt[221]; + m_A270 += gt[222]; + m_A270 += gt[223]; + m_A270 += gt[224]; + m_A275 += go[221]; + m_A262 += go[222]; + m_A272 += go[223]; + m_A263 += go[224]; + double RHS56 = Idr[221]; + RHS56 += Idr[222]; + RHS56 += Idr[223]; + RHS56 += Idr[224]; + m_A282 += gt[225]; + m_A282 += gt[226]; + m_A282 += gt[227]; + m_A282 += gt[228]; + m_A282 += gt[229]; + m_A278 += go[225]; + m_A277 += go[226]; + m_A277 += go[227]; + m_A276 += go[228]; + double RHS57 = Idr[225]; + RHS57 += Idr[226]; + RHS57 += Idr[227]; + RHS57 += Idr[228]; + RHS57 += Idr[229]; + RHS57 -= go[229] * *cnV[229]; + m_A287 += gt[230]; + m_A287 += gt[231]; + m_A287 += gt[232]; + m_A288 += go[230]; + double RHS58 = Idr[230]; + RHS58 += Idr[231]; + RHS58 += Idr[232]; + RHS58 -= go[231] * *cnV[231]; + RHS58 -= go[232] * *cnV[232]; + m_A293 += gt[233]; + m_A293 += gt[234]; + m_A293 += gt[235]; + m_A293 += gt[236]; + m_A293 += gt[237]; + m_A291 += go[233]; + m_A290 += go[234]; + m_A290 += go[235]; + m_A289 += go[236]; + double RHS59 = Idr[233]; + RHS59 += Idr[234]; + RHS59 += Idr[235]; + RHS59 += Idr[236]; + RHS59 += Idr[237]; + RHS59 -= go[237] * *cnV[237]; + m_A297 += gt[238]; + m_A297 += gt[239]; + m_A297 += gt[240]; + m_A298 += go[238]; + double RHS60 = Idr[238]; + RHS60 += Idr[239]; + RHS60 += Idr[240]; + RHS60 -= go[239] * *cnV[239]; + RHS60 -= go[240] * *cnV[240]; + m_A299 += gt[241]; + m_A299 += gt[242]; + m_A299 += gt[243]; + m_A300 += go[241]; + double RHS61 = Idr[241]; + RHS61 += Idr[242]; + RHS61 += Idr[243]; + RHS61 -= go[242] * *cnV[242]; + RHS61 -= go[243] * *cnV[243]; + m_A301 += gt[244]; + m_A301 += gt[245]; + m_A301 += gt[246]; + m_A302 += go[244]; + double RHS62 = Idr[244]; + RHS62 += Idr[245]; + RHS62 += Idr[246]; + RHS62 -= go[245] * *cnV[245]; + RHS62 -= go[246] * *cnV[246]; + m_A314 += gt[247]; + m_A314 += gt[248]; + m_A314 += gt[249]; + m_A314 += gt[250]; + m_A314 += gt[251]; + m_A303 += go[247]; + m_A317 += go[248]; + m_A304 += go[249]; + double RHS63 = Idr[247]; + RHS63 += Idr[248]; + RHS63 += Idr[249]; + RHS63 += Idr[250]; + RHS63 += Idr[251]; + RHS63 -= go[250] * *cnV[250]; + RHS63 -= go[251] * *cnV[251]; + m_A320 += gt[252]; + m_A320 += gt[253]; + m_A320 += gt[254]; + m_A321 += go[252]; + m_A318 += go[253]; + m_A319 += go[254]; + double RHS64 = Idr[252]; + RHS64 += Idr[253]; + RHS64 += Idr[254]; + m_A329 += gt[255]; + m_A329 += gt[256]; + m_A329 += gt[257]; + m_A331 += go[255]; + m_A322 += go[256]; + m_A322 += go[257]; + double RHS65 = Idr[255]; + RHS65 += Idr[256]; + RHS65 += Idr[257]; + m_A337 += gt[258]; + m_A337 += gt[259]; + m_A337 += gt[260]; + m_A337 += gt[261]; + m_A338 += go[258]; + m_A332 += go[259]; + double RHS66 = Idr[258]; + RHS66 += Idr[259]; + RHS66 += Idr[260]; + RHS66 += Idr[261]; + RHS66 -= go[260] * *cnV[260]; + RHS66 -= go[261] * *cnV[261]; + m_A349 += gt[262]; + m_A349 += gt[263]; + m_A349 += gt[264]; + m_A349 += gt[265]; + m_A349 += gt[266]; + m_A349 += gt[267]; + m_A349 += gt[268]; + m_A349 += gt[269]; + m_A349 += gt[270]; + m_A349 += gt[271]; + m_A346 += go[262]; + m_A347 += go[263]; + m_A348 += go[264]; + m_A345 += go[265]; + m_A344 += go[266]; + m_A343 += go[267]; + m_A342 += go[268]; + m_A341 += go[269]; + m_A340 += go[270]; + m_A339 += go[271]; + double RHS67 = Idr[262]; + RHS67 += Idr[263]; + RHS67 += Idr[264]; + RHS67 += Idr[265]; + RHS67 += Idr[266]; + RHS67 += Idr[267]; + RHS67 += Idr[268]; + RHS67 += Idr[269]; + RHS67 += Idr[270]; + RHS67 += Idr[271]; + const double f0 = 1.0 / m_A0; + const double f0_30 = -f0 * m_A81; + m_A83 += m_A1 * f0_30; + RHS30 += f0_30 * RHS0; + const double f1 = 1.0 / m_A2; + const double f1_32 = -f1 * m_A88; + m_A89 += m_A3 * f1_32; + RHS32 += f1_32 * RHS1; + const double f2 = 1.0 / m_A4; + const double f2_30 = -f2 * m_A82; + m_A83 += m_A5 * f2_30; + m_A84 += m_A6 * f2_30; + RHS30 += f2_30 * RHS2; + const double f2_50 = -f2 * m_A185; + m_A187 += m_A5 * f2_50; + m_A192 += m_A6 * f2_50; + RHS50 += f2_50 * RHS2; + const double f3 = 1.0 / m_A7; + const double f3_34 = -f3 * m_A97; + m_A98 += m_A8 * f3_34; + RHS34 += f3_34 * RHS3; + const double f4 = 1.0 / m_A9; + const double f4_31 = -f4 * m_A85; + m_A86 += m_A10 * f4_31; + m_A87 += m_A11 * f4_31; + RHS31 += f4_31 * RHS4; + const double f4_50 = -f4 * m_A186; + m_A188 += m_A10 * f4_50; + m_A192 += m_A11 * f4_50; + RHS50 += f4_50 * RHS4; + const double f5 = 1.0 / m_A12; + const double f5_35 = -f5 * m_A100; + m_A101 += m_A13 * f5_35; + RHS35 += f5_35 * RHS5; + const double f6 = 1.0 / m_A14; + const double f6_38 = -f6 * m_A116; + m_A119 += m_A15 * f6_38; + RHS38 += f6_38 * RHS6; + const double f7 = 1.0 / m_A16; + const double f7_51 = -f7 * m_A194; + m_A201 += m_A17 * f7_51; + RHS51 += f7_51 * RHS7; + const double f8 = 1.0 / m_A18; + const double f8_38 = -f8 * m_A117; + m_A119 += m_A19 * f8_38; + RHS38 += f8_38 * RHS8; + const double f9 = 1.0 / m_A20; + const double f9_36 = -f9 * m_A104; + m_A106 += m_A21 * f9_36; + m_A107 += m_A22 * f9_36; + m_A108 += m_A23 * f9_36; + RHS36 += f9_36 * RHS9; + const double f9_37 = -f9 * m_A109; + m_A112 += m_A21 * f9_37; + m_A114 += m_A22 * f9_37; + m_A115 += m_A23 * f9_37; + RHS37 += f9_37 * RHS9; + const double f9_47 = -f9 * m_A166; + m_A169 += m_A21 * f9_47; + m_A172 += m_A22 * f9_47; + m_A173 += m_A23 * f9_47; + RHS47 += f9_47 * RHS9; + const double f9_51 = -f9 * m_A195; + m_A197 += m_A21 * f9_51; + m_A199 += m_A22 * f9_51; + m_A201 += m_A23 * f9_51; + RHS51 += f9_51 * RHS9; + const double f10 = 1.0 / m_A24; + const double f10_37 = -f10 * m_A110; + m_A111 += m_A25 * f10_37; + m_A112 += m_A26 * f10_37; + m_A113 += m_A27 * f10_37; + RHS37 += f10_37 * RHS10; + const double f10_40 = -f10 * m_A126; + m_A128 += m_A25 * f10_40; + m_A129 += m_A26 * f10_40; + m_A130 += m_A27 * f10_40; + RHS40 += f10_40 * RHS10; + const double f11 = 1.0 / m_A28; + const double f11_41 = -f11 * m_A134; + m_A135 += m_A29 * f11_41; + m_A141 += m_A30 * f11_41; + RHS41 += f11_41 * RHS11; + const double f12 = 1.0 / m_A31; + const double f12_38 = -f12 * m_A118; + m_A119 += m_A32 * f12_38; + m_A120 += m_A33 * f12_38; + m_A121 += m_A34 * f12_38; + RHS38 += f12_38 * RHS12; + const double f12_41 = -f12 * m_A135; + m_A137 += m_A32 * f12_41; + m_A140 += m_A33 * f12_41; + m_A141 += m_A34 * f12_41; + RHS41 += f12_41 * RHS12; + const double f12_47 = -f12 * m_A167; + m_A170 += m_A32 * f12_47; + m_A172 += m_A33 * f12_47; + m_A176 += m_A34 * f12_47; + RHS47 += f12_47 * RHS12; + const double f12_54 = -f12 * m_A232; + m_A233 += m_A32 * f12_54; + m_A236 += m_A33 * f12_54; + m_A240 += m_A34 * f12_54; + RHS54 += f12_54 * RHS12; + const double f13 = 1.0 / m_A35; + const double f13_39 = -f13 * m_A122; + m_A124 += m_A36 * f13_39; + m_A125 += m_A37 * f13_39; + RHS39 += f13_39 * RHS13; + const double f13_40 = -f13 * m_A127; + m_A130 += m_A36 * f13_40; + m_A133 += m_A37 * f13_40; + RHS40 += f13_40 * RHS13; + const double f13_52 = -f13 * m_A206; + m_A210 += m_A36 * f13_52; + m_A214 += m_A37 * f13_52; + RHS52 += f13_52 * RHS13; + const double f14 = 1.0 / m_A38; + const double f14_55 = -f14 * m_A246; + m_A250 += m_A39 * f14_55; + RHS55 += f14_55 * RHS14; + const double f15 = 1.0 / m_A40; + const double f15_41 = -f15 * m_A136; + m_A138 += m_A41 * f15_41; + m_A139 += m_A42 * f15_41; + m_A142 += m_A43 * f15_41; + RHS41 += f15_41 * RHS15; + const double f15_42 = -f15 * m_A143; + m_A144 += m_A41 * f15_42; + m_A145 += m_A42 * f15_42; + m_A148 += m_A43 * f15_42; + RHS42 += f15_42 * RHS15; + const double f15_55 = -f15 * m_A247; + m_A249 += m_A41 * f15_55; + m_A250 += m_A42 * f15_55; + m_A257 += m_A43 * f15_55; + RHS55 += f15_55 * RHS15; + const double f16 = 1.0 / m_A44; + const double f16_33 = -f16 * m_A91; + m_A93 += m_A45 * f16_33; + m_A95 += m_A46 * f16_33; + RHS33 += f16_33 * RHS16; + const double f16_52 = -f16 * m_A207; + m_A208 += m_A45 * f16_52; + m_A214 += m_A46 * f16_52; + RHS52 += f16_52 * RHS16; + const double f17 = 1.0 / m_A47; + const double f17_45 = -f17 * m_A155; + m_A158 += m_A48 * f17_45; + m_A160 += m_A49 * f17_45; + RHS45 += f17_45 * RHS17; + const double f17_55 = -f17 * m_A248; + m_A251 += m_A48 * f17_55; + m_A257 += m_A49 * f17_55; + RHS55 += f17_55 * RHS17; + const double f18 = 1.0 / m_A50; + const double f18_63 = -f18 * m_A303; + m_A304 += m_A51 * f18_63; + RHS63 += f18_63 * RHS18; + const double f19 = 1.0 / m_A52; + const double f19_46 = -f19 * m_A162; + m_A164 += m_A53 * f19_46; + RHS46 += f19_46 * RHS19; + const double f20 = 1.0 / m_A54; + const double f20_43 = -f20 * m_A149; + m_A150 += m_A55 * f20_43; + m_A151 += m_A56 * f20_43; + RHS43 += f20_43 * RHS20; + const double f20_59 = -f20 * m_A289; + m_A290 += m_A55 * f20_59; + m_A293 += m_A56 * f20_59; + RHS59 += f20_59 * RHS20; + const double f21 = 1.0 / m_A57; + const double f21_44 = -f21 * m_A152; + m_A153 += m_A58 * f21_44; + m_A154 += m_A59 * f21_44; + RHS44 += f21_44 * RHS21; + const double f21_57 = -f21 * m_A276; + m_A277 += m_A58 * f21_57; + m_A282 += m_A59 * f21_57; + RHS57 += f21_57 * RHS21; + const double f22 = 1.0 / m_A60; + const double f22_33 = -f22 * m_A92; + m_A93 += m_A61 * f22_33; + m_A94 += m_A62 * f22_33; + m_A96 += m_A63 * f22_33; + RHS33 += f22_33 * RHS22; + const double f22_45 = -f22 * m_A156; + m_A157 += m_A61 * f22_45; + m_A158 += m_A62 * f22_45; + m_A161 += m_A63 * f22_45; + RHS45 += f22_45 * RHS22; + const double f22_63 = -f22 * m_A304; + m_A305 += m_A61 * f22_63; + m_A306 += m_A62 * f22_63; + m_A314 += m_A63 * f22_63; + RHS63 += f22_63 * RHS22; + const double f23 = 1.0 / m_A64; + const double f23_56 = -f23 * m_A262; + m_A270 += m_A65 * f23_56; + RHS56 += f23_56 * RHS23; + const double f24 = 1.0 / m_A66; + const double f24_53 = -f24 * m_A220; + m_A225 += m_A67 * f24_53; + RHS53 += f24_53 * RHS24; + const double f25 = 1.0 / m_A68; + const double f25_48 = -f25 * m_A178; + m_A181 += m_A69 * f25_48; + RHS48 += f25_48 * RHS25; + const double f26 = 1.0 / m_A70; + const double f26_46 = -f26 * m_A163; + m_A164 += m_A71 * f26_46; + m_A165 += m_A72 * f26_46; + RHS46 += f26_46 * RHS26; + const double f26_48 = -f26 * m_A179; + m_A180 += m_A71 * f26_48; + m_A181 += m_A72 * f26_48; + RHS48 += f26_48 * RHS26; + const double f27 = 1.0 / m_A73; + const double f27_47 = -f27 * m_A168; + m_A172 += m_A74 * f27_47; + m_A175 += m_A75 * f27_47; + m_A177 += m_A76 * f27_47; + RHS47 += f27_47 * RHS27; + const double f27_53 = -f27 * m_A221; + m_A222 += m_A74 * f27_53; + m_A225 += m_A75 * f27_53; + m_A228 += m_A76 * f27_53; + RHS53 += f27_53 * RHS27; + const double f27_56 = -f27 * m_A263; + m_A264 += m_A74 * f27_56; + m_A267 += m_A75 * f27_56; + m_A270 += m_A76 * f27_56; + RHS56 += f27_56 * RHS27; + const double f28 = 1.0 / m_A77; + const double f28_64 = -f28 * m_A318; + m_A320 += m_A78 * f28_64; + RHS64 += f28_64 * RHS28; + const double f29 = 1.0 / m_A79; + const double f29_67 = -f29 * m_A339; + m_A349 += m_A80 * f29_67; + RHS67 += f29_67 * RHS29; + const double f30 = 1.0 / m_A83; + const double f30_50 = -f30 * m_A187; + m_A192 += m_A84 * f30_50; + RHS50 += f30_50 * RHS30; + const double f31 = 1.0 / m_A86; + const double f31_50 = -f31 * m_A188; + m_A192 += m_A87 * f31_50; + RHS50 += f31_50 * RHS31; + const double f32 = 1.0 / m_A89; + const double f32_50 = -f32 * m_A189; + m_A192 += m_A90 * f32_50; + RHS50 += f32_50 * RHS32; + const double f33 = 1.0 / m_A93; + const double f33_45 = -f33 * m_A157; + m_A158 += m_A94 * f33_45; + m_A159 += m_A95 * f33_45; + m_A161 += m_A96 * f33_45; + RHS45 += f33_45 * RHS33; + const double f33_52 = -f33 * m_A208; + m_A211 += m_A94 * f33_52; + m_A214 += m_A95 * f33_52; + m_A219 += m_A96 * f33_52; + RHS52 += f33_52 * RHS33; + const double f33_63 = -f33 * m_A305; + m_A306 += m_A94 * f33_63; + m_A307 += m_A95 * f33_63; + m_A314 += m_A96 * f33_63; + RHS63 += f33_63 * RHS33; + const double f34 = 1.0 / m_A98; + const double f34_50 = -f34 * m_A190; + m_A192 += m_A99 * f34_50; + RHS50 += f34_50 * RHS34; + const double f35 = 1.0 / m_A101; + const double f35_50 = -f35 * m_A191; + m_A192 += m_A102 * f35_50; + m_A193 += m_A103 * f35_50; + RHS50 += f35_50 * RHS35; + const double f35_51 = -f35 * m_A196; + m_A200 += m_A102 * f35_51; + m_A201 += m_A103 * f35_51; + RHS51 += f35_51 * RHS35; + const double f36 = 1.0 / m_A105; + const double f36_37 = -f36 * m_A111; + m_A112 += m_A106 * f36_37; + m_A114 += m_A107 * f36_37; + m_A115 += m_A108 * f36_37; + RHS37 += f36_37 * RHS36; + const double f36_40 = -f36 * m_A128; + m_A129 += m_A106 * f36_40; + m_A131 += m_A107 * f36_40; + m_A132 += m_A108 * f36_40; + RHS40 += f36_40 * RHS36; + const double f37 = 1.0 / m_A112; + const double f37_40 = -f37 * m_A129; + m_A130 += m_A113 * f37_40; + m_A131 += m_A114 * f37_40; + m_A132 += m_A115 * f37_40; + RHS40 += f37_40 * RHS37; + const double f37_47 = -f37 * m_A169; + m_A171 += m_A113 * f37_47; + m_A172 += m_A114 * f37_47; + m_A173 += m_A115 * f37_47; + RHS47 += f37_47 * RHS37; + const double f37_51 = -f37 * m_A197; + m_A198 += m_A113 * f37_51; + m_A199 += m_A114 * f37_51; + m_A201 += m_A115 * f37_51; + RHS51 += f37_51 * RHS37; + const double f38 = 1.0 / m_A119; + const double f38_41 = -f38 * m_A137; + m_A140 += m_A120 * f38_41; + m_A141 += m_A121 * f38_41; + RHS41 += f38_41 * RHS38; + const double f38_47 = -f38 * m_A170; + m_A172 += m_A120 * f38_47; + m_A176 += m_A121 * f38_47; + RHS47 += f38_47 * RHS38; + const double f38_54 = -f38 * m_A233; + m_A236 += m_A120 * f38_54; + m_A240 += m_A121 * f38_54; + RHS54 += f38_54 * RHS38; + const double f39 = 1.0 / m_A123; + const double f39_52 = -f39 * m_A209; + m_A210 += m_A124 * f39_52; + m_A214 += m_A125 * f39_52; + RHS52 += f39_52 * RHS39; + const double f40 = 1.0 / m_A130; + const double f40_47 = -f40 * m_A171; + m_A172 += m_A131 * f40_47; + m_A173 += m_A132 * f40_47; + m_A174 += m_A133 * f40_47; + RHS47 += f40_47 * RHS40; + const double f40_51 = -f40 * m_A198; + m_A199 += m_A131 * f40_51; + m_A201 += m_A132 * f40_51; + m_A202 += m_A133 * f40_51; + RHS51 += f40_51 * RHS40; + const double f40_52 = -f40 * m_A210; + m_A212 += m_A131 * f40_52; + m_A213 += m_A132 * f40_52; + m_A214 += m_A133 * f40_52; + RHS52 += f40_52 * RHS40; + const double f41 = 1.0 / m_A138; + const double f41_42 = -f41 * m_A144; + m_A145 += m_A139 * f41_42; + m_A146 += m_A140 * f41_42; + m_A147 += m_A141 * f41_42; + m_A148 += m_A142 * f41_42; + RHS42 += f41_42 * RHS41; + const double f41_54 = -f41 * m_A234; + m_A235 += m_A139 * f41_54; + m_A236 += m_A140 * f41_54; + m_A240 += m_A141 * f41_54; + m_A241 += m_A142 * f41_54; + RHS54 += f41_54 * RHS41; + const double f41_55 = -f41 * m_A249; + m_A250 += m_A139 * f41_55; + m_A252 += m_A140 * f41_55; + m_A256 += m_A141 * f41_55; + m_A257 += m_A142 * f41_55; + RHS55 += f41_55 * RHS41; + const double f42 = 1.0 / m_A145; + const double f42_54 = -f42 * m_A235; + m_A236 += m_A146 * f42_54; + m_A240 += m_A147 * f42_54; + m_A241 += m_A148 * f42_54; + RHS54 += f42_54 * RHS42; + const double f42_55 = -f42 * m_A250; + m_A252 += m_A146 * f42_55; + m_A256 += m_A147 * f42_55; + m_A257 += m_A148 * f42_55; + RHS55 += f42_55 * RHS42; + const double f43 = 1.0 / m_A150; + const double f43_59 = -f43 * m_A290; + m_A293 += m_A151 * f43_59; + RHS59 += f43_59 * RHS43; + const double f44 = 1.0 / m_A153; + const double f44_57 = -f44 * m_A277; + m_A282 += m_A154 * f44_57; + RHS57 += f44_57 * RHS44; + const double f45 = 1.0 / m_A158; + const double f45_52 = -f45 * m_A211; + m_A214 += m_A159 * f45_52; + m_A217 += m_A160 * f45_52; + m_A219 += m_A161 * f45_52; + RHS52 += f45_52 * RHS45; + const double f45_55 = -f45 * m_A251; + m_A254 += m_A159 * f45_55; + m_A257 += m_A160 * f45_55; + m_A260 += m_A161 * f45_55; + RHS55 += f45_55 * RHS45; + const double f45_63 = -f45 * m_A306; + m_A307 += m_A159 * f45_63; + m_A310 += m_A160 * f45_63; + m_A314 += m_A161 * f45_63; + RHS63 += f45_63 * RHS45; + const double f46 = 1.0 / m_A164; + const double f46_48 = -f46 * m_A180; + m_A181 += m_A165 * f46_48; + RHS48 += f46_48 * RHS46; + const double f47 = 1.0 / m_A172; + const double f47_51 = -f47 * m_A199; + m_A201 += m_A173 * f47_51; + m_A202 += m_A174 * f47_51; + m_A203 += m_A175 * f47_51; + m_A204 += m_A176 * f47_51; + m_A205 += m_A177 * f47_51; + RHS51 += f47_51 * RHS47; + const double f47_52 = -f47 * m_A212; + m_A213 += m_A173 * f47_52; + m_A214 += m_A174 * f47_52; + m_A215 += m_A175 * f47_52; + m_A216 += m_A176 * f47_52; + m_A218 += m_A177 * f47_52; + RHS52 += f47_52 * RHS47; + const double f47_53 = -f47 * m_A222; + m_A223 += m_A173 * f47_53; + m_A224 += m_A174 * f47_53; + m_A225 += m_A175 * f47_53; + m_A226 += m_A176 * f47_53; + m_A228 += m_A177 * f47_53; + RHS53 += f47_53 * RHS47; + const double f47_54 = -f47 * m_A236; + m_A237 += m_A173 * f47_54; + m_A238 += m_A174 * f47_54; + m_A239 += m_A175 * f47_54; + m_A240 += m_A176 * f47_54; + m_A242 += m_A177 * f47_54; + RHS54 += f47_54 * RHS47; + const double f47_55 = -f47 * m_A252; + m_A253 += m_A173 * f47_55; + m_A254 += m_A174 * f47_55; + m_A255 += m_A175 * f47_55; + m_A256 += m_A176 * f47_55; + m_A258 += m_A177 * f47_55; + RHS55 += f47_55 * RHS47; + const double f47_56 = -f47 * m_A264; + m_A265 += m_A173 * f47_56; + m_A266 += m_A174 * f47_56; + m_A267 += m_A175 * f47_56; + m_A268 += m_A176 * f47_56; + m_A270 += m_A177 * f47_56; + RHS56 += f47_56 * RHS47; + const double f48 = 1.0 / m_A181; + const double f48_64 = -f48 * m_A319; + m_A320 += m_A182 * f48_64; + RHS64 += f48_64 * RHS48; + const double f49 = 1.0 / m_A183; + const double f49_67 = -f49 * m_A340; + m_A349 += m_A184 * f49_67; + RHS67 += f49_67 * RHS49; + const double f50 = 1.0 / m_A192; + const double f50_51 = -f50 * m_A200; + m_A201 += m_A193 * f50_51; + RHS51 += f50_51 * RHS50; + const double f51 = 1.0 / m_A201; + const double f51_52 = -f51 * m_A213; + m_A214 += m_A202 * f51_52; + m_A215 += m_A203 * f51_52; + m_A216 += m_A204 * f51_52; + m_A218 += m_A205 * f51_52; + RHS52 += f51_52 * RHS51; + const double f51_53 = -f51 * m_A223; + m_A224 += m_A202 * f51_53; + m_A225 += m_A203 * f51_53; + m_A226 += m_A204 * f51_53; + m_A228 += m_A205 * f51_53; + RHS53 += f51_53 * RHS51; + const double f51_54 = -f51 * m_A237; + m_A238 += m_A202 * f51_54; + m_A239 += m_A203 * f51_54; + m_A240 += m_A204 * f51_54; + m_A242 += m_A205 * f51_54; + RHS54 += f51_54 * RHS51; + const double f51_55 = -f51 * m_A253; + m_A254 += m_A202 * f51_55; + m_A255 += m_A203 * f51_55; + m_A256 += m_A204 * f51_55; + m_A258 += m_A205 * f51_55; + RHS55 += f51_55 * RHS51; + const double f51_56 = -f51 * m_A265; + m_A266 += m_A202 * f51_56; + m_A267 += m_A203 * f51_56; + m_A268 += m_A204 * f51_56; + m_A270 += m_A205 * f51_56; + RHS56 += f51_56 * RHS51; + const double f52 = 1.0 / m_A214; + const double f52_53 = -f52 * m_A224; + m_A225 += m_A215 * f52_53; + m_A226 += m_A216 * f52_53; + m_A227 += m_A217 * f52_53; + m_A228 += m_A218 * f52_53; + m_A230 += m_A219 * f52_53; + RHS53 += f52_53 * RHS52; + const double f52_54 = -f52 * m_A238; + m_A239 += m_A215 * f52_54; + m_A240 += m_A216 * f52_54; + m_A241 += m_A217 * f52_54; + m_A242 += m_A218 * f52_54; + m_A244 += m_A219 * f52_54; + RHS54 += f52_54 * RHS52; + const double f52_55 = -f52 * m_A254; + m_A255 += m_A215 * f52_55; + m_A256 += m_A216 * f52_55; + m_A257 += m_A217 * f52_55; + m_A258 += m_A218 * f52_55; + m_A260 += m_A219 * f52_55; + RHS55 += f52_55 * RHS52; + const double f52_56 = -f52 * m_A266; + m_A267 += m_A215 * f52_56; + m_A268 += m_A216 * f52_56; + m_A269 += m_A217 * f52_56; + m_A270 += m_A218 * f52_56; + m_A273 += m_A219 * f52_56; + RHS56 += f52_56 * RHS52; + const double f52_63 = -f52 * m_A307; + m_A308 += m_A215 * f52_63; + m_A309 += m_A216 * f52_63; + m_A310 += m_A217 * f52_63; + m_A311 += m_A218 * f52_63; + m_A314 += m_A219 * f52_63; + RHS63 += f52_63 * RHS52; + const double f53 = 1.0 / m_A225; + const double f53_54 = -f53 * m_A239; + m_A240 += m_A226 * f53_54; + m_A241 += m_A227 * f53_54; + m_A242 += m_A228 * f53_54; + m_A243 += m_A229 * f53_54; + m_A244 += m_A230 * f53_54; + m_A245 += m_A231 * f53_54; + RHS54 += f53_54 * RHS53; + const double f53_55 = -f53 * m_A255; + m_A256 += m_A226 * f53_55; + m_A257 += m_A227 * f53_55; + m_A258 += m_A228 * f53_55; + m_A259 += m_A229 * f53_55; + m_A260 += m_A230 * f53_55; + m_A261 += m_A231 * f53_55; + RHS55 += f53_55 * RHS53; + const double f53_56 = -f53 * m_A267; + m_A268 += m_A226 * f53_56; + m_A269 += m_A227 * f53_56; + m_A270 += m_A228 * f53_56; + m_A271 += m_A229 * f53_56; + m_A273 += m_A230 * f53_56; + m_A274 += m_A231 * f53_56; + RHS56 += f53_56 * RHS53; + const double f53_57 = -f53 * m_A278; + m_A279 += m_A226 * f53_57; + m_A280 += m_A227 * f53_57; + m_A281 += m_A228 * f53_57; + m_A282 += m_A229 * f53_57; + m_A284 += m_A230 * f53_57; + m_A285 += m_A231 * f53_57; + RHS57 += f53_57 * RHS53; + const double f53_63 = -f53 * m_A308; + m_A309 += m_A226 * f53_63; + m_A310 += m_A227 * f53_63; + m_A311 += m_A228 * f53_63; + m_A312 += m_A229 * f53_63; + m_A314 += m_A230 * f53_63; + m_A315 += m_A231 * f53_63; + RHS63 += f53_63 * RHS53; + const double f53_65 = -f53 * m_A322; + m_A323 += m_A226 * f53_65; + m_A324 += m_A227 * f53_65; + m_A325 += m_A228 * f53_65; + m_A326 += m_A229 * f53_65; + m_A328 += m_A230 * f53_65; + m_A329 += m_A231 * f53_65; + RHS65 += f53_65 * RHS53; + const double f54 = 1.0 / m_A240; + const double f54_55 = -f54 * m_A256; + m_A257 += m_A241 * f54_55; + m_A258 += m_A242 * f54_55; + m_A259 += m_A243 * f54_55; + m_A260 += m_A244 * f54_55; + m_A261 += m_A245 * f54_55; + RHS55 += f54_55 * RHS54; + const double f54_56 = -f54 * m_A268; + m_A269 += m_A241 * f54_56; + m_A270 += m_A242 * f54_56; + m_A271 += m_A243 * f54_56; + m_A273 += m_A244 * f54_56; + m_A274 += m_A245 * f54_56; + RHS56 += f54_56 * RHS54; + const double f54_57 = -f54 * m_A279; + m_A280 += m_A241 * f54_57; + m_A281 += m_A242 * f54_57; + m_A282 += m_A243 * f54_57; + m_A284 += m_A244 * f54_57; + m_A285 += m_A245 * f54_57; + RHS57 += f54_57 * RHS54; + const double f54_63 = -f54 * m_A309; + m_A310 += m_A241 * f54_63; + m_A311 += m_A242 * f54_63; + m_A312 += m_A243 * f54_63; + m_A314 += m_A244 * f54_63; + m_A315 += m_A245 * f54_63; + RHS63 += f54_63 * RHS54; + const double f54_65 = -f54 * m_A323; + m_A324 += m_A241 * f54_65; + m_A325 += m_A242 * f54_65; + m_A326 += m_A243 * f54_65; + m_A328 += m_A244 * f54_65; + m_A329 += m_A245 * f54_65; + RHS65 += f54_65 * RHS54; + const double f55 = 1.0 / m_A257; + const double f55_56 = -f55 * m_A269; + m_A270 += m_A258 * f55_56; + m_A271 += m_A259 * f55_56; + m_A273 += m_A260 * f55_56; + m_A274 += m_A261 * f55_56; + RHS56 += f55_56 * RHS55; + const double f55_57 = -f55 * m_A280; + m_A281 += m_A258 * f55_57; + m_A282 += m_A259 * f55_57; + m_A284 += m_A260 * f55_57; + m_A285 += m_A261 * f55_57; + RHS57 += f55_57 * RHS55; + const double f55_63 = -f55 * m_A310; + m_A311 += m_A258 * f55_63; + m_A312 += m_A259 * f55_63; + m_A314 += m_A260 * f55_63; + m_A315 += m_A261 * f55_63; + RHS63 += f55_63 * RHS55; + const double f55_65 = -f55 * m_A324; + m_A325 += m_A258 * f55_65; + m_A326 += m_A259 * f55_65; + m_A328 += m_A260 * f55_65; + m_A329 += m_A261 * f55_65; + RHS65 += f55_65 * RHS55; + const double f56 = 1.0 / m_A270; + const double f56_57 = -f56 * m_A281; + m_A282 += m_A271 * f56_57; + m_A283 += m_A272 * f56_57; + m_A284 += m_A273 * f56_57; + m_A285 += m_A274 * f56_57; + m_A286 += m_A275 * f56_57; + RHS57 += f56_57 * RHS56; + const double f56_59 = -f56 * m_A291; + m_A292 += m_A271 * f56_59; + m_A293 += m_A272 * f56_59; + m_A294 += m_A273 * f56_59; + m_A295 += m_A274 * f56_59; + m_A296 += m_A275 * f56_59; + RHS59 += f56_59 * RHS56; + const double f56_63 = -f56 * m_A311; + m_A312 += m_A271 * f56_63; + m_A313 += m_A272 * f56_63; + m_A314 += m_A273 * f56_63; + m_A315 += m_A274 * f56_63; + m_A316 += m_A275 * f56_63; + RHS63 += f56_63 * RHS56; + const double f56_65 = -f56 * m_A325; + m_A326 += m_A271 * f56_65; + m_A327 += m_A272 * f56_65; + m_A328 += m_A273 * f56_65; + m_A329 += m_A274 * f56_65; + m_A330 += m_A275 * f56_65; + RHS65 += f56_65 * RHS56; + const double f56_66 = -f56 * m_A332; + m_A333 += m_A271 * f56_66; + m_A334 += m_A272 * f56_66; + m_A335 += m_A273 * f56_66; + m_A336 += m_A274 * f56_66; + m_A337 += m_A275 * f56_66; + RHS66 += f56_66 * RHS56; + const double f57 = 1.0 / m_A282; + const double f57_59 = -f57 * m_A292; + m_A293 += m_A283 * f57_59; + m_A294 += m_A284 * f57_59; + m_A295 += m_A285 * f57_59; + m_A296 += m_A286 * f57_59; + RHS59 += f57_59 * RHS57; + const double f57_63 = -f57 * m_A312; + m_A313 += m_A283 * f57_63; + m_A314 += m_A284 * f57_63; + m_A315 += m_A285 * f57_63; + m_A316 += m_A286 * f57_63; + RHS63 += f57_63 * RHS57; + const double f57_65 = -f57 * m_A326; + m_A327 += m_A283 * f57_65; + m_A328 += m_A284 * f57_65; + m_A329 += m_A285 * f57_65; + m_A330 += m_A286 * f57_65; + RHS65 += f57_65 * RHS57; + const double f57_66 = -f57 * m_A333; + m_A334 += m_A283 * f57_66; + m_A335 += m_A284 * f57_66; + m_A336 += m_A285 * f57_66; + m_A337 += m_A286 * f57_66; + RHS66 += f57_66 * RHS57; + const double f58 = 1.0 / m_A287; + const double f58_67 = -f58 * m_A341; + m_A349 += m_A288 * f58_67; + RHS67 += f58_67 * RHS58; + const double f59 = 1.0 / m_A293; + const double f59_63 = -f59 * m_A313; + m_A314 += m_A294 * f59_63; + m_A315 += m_A295 * f59_63; + m_A316 += m_A296 * f59_63; + RHS63 += f59_63 * RHS59; + const double f59_65 = -f59 * m_A327; + m_A328 += m_A294 * f59_65; + m_A329 += m_A295 * f59_65; + m_A330 += m_A296 * f59_65; + RHS65 += f59_65 * RHS59; + const double f59_66 = -f59 * m_A334; + m_A335 += m_A294 * f59_66; + m_A336 += m_A295 * f59_66; + m_A337 += m_A296 * f59_66; + RHS66 += f59_66 * RHS59; + const double f60 = 1.0 / m_A297; + const double f60_67 = -f60 * m_A342; + m_A349 += m_A298 * f60_67; + RHS67 += f60_67 * RHS60; + const double f61 = 1.0 / m_A299; + const double f61_67 = -f61 * m_A343; + m_A349 += m_A300 * f61_67; + RHS67 += f61_67 * RHS61; + const double f62 = 1.0 / m_A301; + const double f62_67 = -f62 * m_A344; + m_A349 += m_A302 * f62_67; + RHS67 += f62_67 * RHS62; + const double f63 = 1.0 / m_A314; + const double f63_65 = -f63 * m_A328; + m_A329 += m_A315 * f63_65; + m_A330 += m_A316 * f63_65; + m_A331 += m_A317 * f63_65; + RHS65 += f63_65 * RHS63; + const double f63_66 = -f63 * m_A335; + m_A336 += m_A315 * f63_66; + m_A337 += m_A316 * f63_66; + m_A338 += m_A317 * f63_66; + RHS66 += f63_66 * RHS63; + const double f63_67 = -f63 * m_A345; + m_A347 += m_A315 * f63_67; + m_A348 += m_A316 * f63_67; + m_A349 += m_A317 * f63_67; + RHS67 += f63_67 * RHS63; + const double f64 = 1.0 / m_A320; + const double f64_67 = -f64 * m_A346; + m_A349 += m_A321 * f64_67; + RHS67 += f64_67 * RHS64; + const double f65 = 1.0 / m_A329; + const double f65_66 = -f65 * m_A336; + m_A337 += m_A330 * f65_66; + m_A338 += m_A331 * f65_66; + RHS66 += f65_66 * RHS65; + const double f65_67 = -f65 * m_A347; + m_A348 += m_A330 * f65_67; + m_A349 += m_A331 * f65_67; + RHS67 += f65_67 * RHS65; + const double f66 = 1.0 / m_A337; + const double f66_67 = -f66 * m_A348; + m_A349 += m_A338 * f66_67; + RHS67 += f66_67 * RHS66; + V[67] = RHS67 / m_A349; + double tmp66 = 0.0; + tmp66 += m_A338 * V[67]; + V[66] = (RHS66 - tmp66) / m_A337; + double tmp65 = 0.0; + tmp65 += m_A330 * V[66]; + tmp65 += m_A331 * V[67]; + V[65] = (RHS65 - tmp65) / m_A329; + double tmp64 = 0.0; + tmp64 += m_A321 * V[67]; + V[64] = (RHS64 - tmp64) / m_A320; + double tmp63 = 0.0; + tmp63 += m_A315 * V[65]; + tmp63 += m_A316 * V[66]; + tmp63 += m_A317 * V[67]; + V[63] = (RHS63 - tmp63) / m_A314; + double tmp62 = 0.0; + tmp62 += m_A302 * V[67]; + V[62] = (RHS62 - tmp62) / m_A301; + double tmp61 = 0.0; + tmp61 += m_A300 * V[67]; + V[61] = (RHS61 - tmp61) / m_A299; + double tmp60 = 0.0; + tmp60 += m_A298 * V[67]; + V[60] = (RHS60 - tmp60) / m_A297; + double tmp59 = 0.0; + tmp59 += m_A294 * V[63]; + tmp59 += m_A295 * V[65]; + tmp59 += m_A296 * V[66]; + V[59] = (RHS59 - tmp59) / m_A293; + double tmp58 = 0.0; + tmp58 += m_A288 * V[67]; + V[58] = (RHS58 - tmp58) / m_A287; + double tmp57 = 0.0; + tmp57 += m_A283 * V[59]; + tmp57 += m_A284 * V[63]; + tmp57 += m_A285 * V[65]; + tmp57 += m_A286 * V[66]; + V[57] = (RHS57 - tmp57) / m_A282; + double tmp56 = 0.0; + tmp56 += m_A271 * V[57]; + tmp56 += m_A272 * V[59]; + tmp56 += m_A273 * V[63]; + tmp56 += m_A274 * V[65]; + tmp56 += m_A275 * V[66]; + V[56] = (RHS56 - tmp56) / m_A270; + double tmp55 = 0.0; + tmp55 += m_A258 * V[56]; + tmp55 += m_A259 * V[57]; + tmp55 += m_A260 * V[63]; + tmp55 += m_A261 * V[65]; + V[55] = (RHS55 - tmp55) / m_A257; + double tmp54 = 0.0; + tmp54 += m_A241 * V[55]; + tmp54 += m_A242 * V[56]; + tmp54 += m_A243 * V[57]; + tmp54 += m_A244 * V[63]; + tmp54 += m_A245 * V[65]; + V[54] = (RHS54 - tmp54) / m_A240; + double tmp53 = 0.0; + tmp53 += m_A226 * V[54]; + tmp53 += m_A227 * V[55]; + tmp53 += m_A228 * V[56]; + tmp53 += m_A229 * V[57]; + tmp53 += m_A230 * V[63]; + tmp53 += m_A231 * V[65]; + V[53] = (RHS53 - tmp53) / m_A225; + double tmp52 = 0.0; + tmp52 += m_A215 * V[53]; + tmp52 += m_A216 * V[54]; + tmp52 += m_A217 * V[55]; + tmp52 += m_A218 * V[56]; + tmp52 += m_A219 * V[63]; + V[52] = (RHS52 - tmp52) / m_A214; + double tmp51 = 0.0; + tmp51 += m_A202 * V[52]; + tmp51 += m_A203 * V[53]; + tmp51 += m_A204 * V[54]; + tmp51 += m_A205 * V[56]; + V[51] = (RHS51 - tmp51) / m_A201; + double tmp50 = 0.0; + tmp50 += m_A193 * V[51]; + V[50] = (RHS50 - tmp50) / m_A192; + double tmp49 = 0.0; + tmp49 += m_A184 * V[67]; + V[49] = (RHS49 - tmp49) / m_A183; + double tmp48 = 0.0; + tmp48 += m_A182 * V[64]; + V[48] = (RHS48 - tmp48) / m_A181; + double tmp47 = 0.0; + tmp47 += m_A173 * V[51]; + tmp47 += m_A174 * V[52]; + tmp47 += m_A175 * V[53]; + tmp47 += m_A176 * V[54]; + tmp47 += m_A177 * V[56]; + V[47] = (RHS47 - tmp47) / m_A172; + double tmp46 = 0.0; + tmp46 += m_A165 * V[48]; + V[46] = (RHS46 - tmp46) / m_A164; + double tmp45 = 0.0; + tmp45 += m_A159 * V[52]; + tmp45 += m_A160 * V[55]; + tmp45 += m_A161 * V[63]; + V[45] = (RHS45 - tmp45) / m_A158; + double tmp44 = 0.0; + tmp44 += m_A154 * V[57]; + V[44] = (RHS44 - tmp44) / m_A153; + double tmp43 = 0.0; + tmp43 += m_A151 * V[59]; + V[43] = (RHS43 - tmp43) / m_A150; + double tmp42 = 0.0; + tmp42 += m_A146 * V[47]; + tmp42 += m_A147 * V[54]; + tmp42 += m_A148 * V[55]; + V[42] = (RHS42 - tmp42) / m_A145; + double tmp41 = 0.0; + tmp41 += m_A139 * V[42]; + tmp41 += m_A140 * V[47]; + tmp41 += m_A141 * V[54]; + tmp41 += m_A142 * V[55]; + V[41] = (RHS41 - tmp41) / m_A138; + double tmp40 = 0.0; + tmp40 += m_A131 * V[47]; + tmp40 += m_A132 * V[51]; + tmp40 += m_A133 * V[52]; + V[40] = (RHS40 - tmp40) / m_A130; + double tmp39 = 0.0; + tmp39 += m_A124 * V[40]; + tmp39 += m_A125 * V[52]; + V[39] = (RHS39 - tmp39) / m_A123; + double tmp38 = 0.0; + tmp38 += m_A120 * V[47]; + tmp38 += m_A121 * V[54]; + V[38] = (RHS38 - tmp38) / m_A119; + double tmp37 = 0.0; + tmp37 += m_A113 * V[40]; + tmp37 += m_A114 * V[47]; + tmp37 += m_A115 * V[51]; + V[37] = (RHS37 - tmp37) / m_A112; + double tmp36 = 0.0; + tmp36 += m_A106 * V[37]; + tmp36 += m_A107 * V[47]; + tmp36 += m_A108 * V[51]; + V[36] = (RHS36 - tmp36) / m_A105; + double tmp35 = 0.0; + tmp35 += m_A102 * V[50]; + tmp35 += m_A103 * V[51]; + V[35] = (RHS35 - tmp35) / m_A101; + double tmp34 = 0.0; + tmp34 += m_A99 * V[50]; + V[34] = (RHS34 - tmp34) / m_A98; + double tmp33 = 0.0; + tmp33 += m_A94 * V[45]; + tmp33 += m_A95 * V[52]; + tmp33 += m_A96 * V[63]; + V[33] = (RHS33 - tmp33) / m_A93; + double tmp32 = 0.0; + tmp32 += m_A90 * V[50]; + V[32] = (RHS32 - tmp32) / m_A89; + double tmp31 = 0.0; + tmp31 += m_A87 * V[50]; + V[31] = (RHS31 - tmp31) / m_A86; + double tmp30 = 0.0; + tmp30 += m_A84 * V[50]; + V[30] = (RHS30 - tmp30) / m_A83; + double tmp29 = 0.0; + tmp29 += m_A80 * V[67]; + V[29] = (RHS29 - tmp29) / m_A79; + double tmp28 = 0.0; + tmp28 += m_A78 * V[64]; + V[28] = (RHS28 - tmp28) / m_A77; + double tmp27 = 0.0; + tmp27 += m_A74 * V[47]; + tmp27 += m_A75 * V[53]; + tmp27 += m_A76 * V[56]; + V[27] = (RHS27 - tmp27) / m_A73; + double tmp26 = 0.0; + tmp26 += m_A71 * V[46]; + tmp26 += m_A72 * V[48]; + V[26] = (RHS26 - tmp26) / m_A70; + double tmp25 = 0.0; + tmp25 += m_A69 * V[48]; + V[25] = (RHS25 - tmp25) / m_A68; + double tmp24 = 0.0; + tmp24 += m_A67 * V[53]; + V[24] = (RHS24 - tmp24) / m_A66; + double tmp23 = 0.0; + tmp23 += m_A65 * V[56]; + V[23] = (RHS23 - tmp23) / m_A64; + double tmp22 = 0.0; + tmp22 += m_A61 * V[33]; + tmp22 += m_A62 * V[45]; + tmp22 += m_A63 * V[63]; + V[22] = (RHS22 - tmp22) / m_A60; + double tmp21 = 0.0; + tmp21 += m_A58 * V[44]; + tmp21 += m_A59 * V[57]; + V[21] = (RHS21 - tmp21) / m_A57; + double tmp20 = 0.0; + tmp20 += m_A55 * V[43]; + tmp20 += m_A56 * V[59]; + V[20] = (RHS20 - tmp20) / m_A54; + double tmp19 = 0.0; + tmp19 += m_A53 * V[46]; + V[19] = (RHS19 - tmp19) / m_A52; + double tmp18 = 0.0; + tmp18 += m_A51 * V[22]; + V[18] = (RHS18 - tmp18) / m_A50; + double tmp17 = 0.0; + tmp17 += m_A48 * V[45]; + tmp17 += m_A49 * V[55]; + V[17] = (RHS17 - tmp17) / m_A47; + double tmp16 = 0.0; + tmp16 += m_A45 * V[33]; + tmp16 += m_A46 * V[52]; + V[16] = (RHS16 - tmp16) / m_A44; + double tmp15 = 0.0; + tmp15 += m_A41 * V[41]; + tmp15 += m_A42 * V[42]; + tmp15 += m_A43 * V[55]; + V[15] = (RHS15 - tmp15) / m_A40; + double tmp14 = 0.0; + tmp14 += m_A39 * V[42]; + V[14] = (RHS14 - tmp14) / m_A38; + double tmp13 = 0.0; + tmp13 += m_A36 * V[40]; + tmp13 += m_A37 * V[52]; + V[13] = (RHS13 - tmp13) / m_A35; + double tmp12 = 0.0; + tmp12 += m_A32 * V[38]; + tmp12 += m_A33 * V[47]; + tmp12 += m_A34 * V[54]; + V[12] = (RHS12 - tmp12) / m_A31; + double tmp11 = 0.0; + tmp11 += m_A29 * V[12]; + tmp11 += m_A30 * V[54]; + V[11] = (RHS11 - tmp11) / m_A28; + double tmp10 = 0.0; + tmp10 += m_A25 * V[36]; + tmp10 += m_A26 * V[37]; + tmp10 += m_A27 * V[40]; + V[10] = (RHS10 - tmp10) / m_A24; + double tmp9 = 0.0; + tmp9 += m_A21 * V[37]; + tmp9 += m_A22 * V[47]; + tmp9 += m_A23 * V[51]; + V[9] = (RHS9 - tmp9) / m_A20; + double tmp8 = 0.0; + tmp8 += m_A19 * V[38]; + V[8] = (RHS8 - tmp8) / m_A18; + double tmp7 = 0.0; + tmp7 += m_A17 * V[51]; + V[7] = (RHS7 - tmp7) / m_A16; + double tmp6 = 0.0; + tmp6 += m_A15 * V[38]; + V[6] = (RHS6 - tmp6) / m_A14; + double tmp5 = 0.0; + tmp5 += m_A13 * V[35]; + V[5] = (RHS5 - tmp5) / m_A12; + double tmp4 = 0.0; + tmp4 += m_A10 * V[31]; + tmp4 += m_A11 * V[50]; + V[4] = (RHS4 - tmp4) / m_A9; + double tmp3 = 0.0; + tmp3 += m_A8 * V[34]; + V[3] = (RHS3 - tmp3) / m_A7; + double tmp2 = 0.0; + tmp2 += m_A5 * V[30]; + tmp2 += m_A6 * V[50]; + V[2] = (RHS2 - tmp2) / m_A4; + double tmp1 = 0.0; + tmp1 += m_A3 * V[32]; + V[1] = (RHS1 - tmp1) / m_A2; + double tmp0 = 0.0; + tmp0 += m_A1 * V[30]; + V[0] = (RHS0 - tmp0) / m_A0; +} + // stuntcyc static void nl_gcr_c924fe5960b1479e_20_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) @@ -63829,806 +67233,6 @@ static void nl_gcr_d5e1c37cfa2d2853_30_double_double(double * __restrict V, cons V[0] = (RHS0 - tmp0) / m_A0; } -// brdrline -static void nl_gcr_d74d89ba31f2cb81_126_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) - -{ - - plib::unused_var(cnV); - double m_A0(0.0); - double m_A1(0.0); - double m_A2(0.0); - double m_A3(0.0); - double m_A4(0.0); - double m_A5(0.0); - double m_A6(0.0); - double m_A7(0.0); - double m_A8(0.0); - double m_A9(0.0); - double m_A10(0.0); - double m_A11(0.0); - double m_A12(0.0); - double m_A13(0.0); - double m_A14(0.0); - double m_A15(0.0); - double m_A16(0.0); - double m_A17(0.0); - double m_A18(0.0); - double m_A19(0.0); - double m_A20(0.0); - double m_A21(0.0); - double m_A22(0.0); - double m_A23(0.0); - double m_A24(0.0); - double m_A25(0.0); - double m_A26(0.0); - double m_A27(0.0); - double m_A28(0.0); - double m_A29(0.0); - double m_A30(0.0); - double m_A31(0.0); - double m_A32(0.0); - double m_A33(0.0); - double m_A34(0.0); - double m_A35(0.0); - double m_A36(0.0); - double m_A37(0.0); - double m_A38(0.0); - double m_A39(0.0); - double m_A40(0.0); - double m_A41(0.0); - double m_A42(0.0); - double m_A43(0.0); - double m_A44(0.0); - double m_A45(0.0); - double m_A46(0.0); - double m_A47(0.0); - double m_A48(0.0); - double m_A49(0.0); - double m_A50(0.0); - double m_A51(0.0); - double m_A52(0.0); - double m_A53(0.0); - double m_A54(0.0); - double m_A55(0.0); - double m_A56(0.0); - double m_A57(0.0); - double m_A58(0.0); - double m_A59(0.0); - double m_A60(0.0); - double m_A61(0.0); - double m_A62(0.0); - double m_A63(0.0); - double m_A64(0.0); - double m_A65(0.0); - double m_A66(0.0); - double m_A67(0.0); - double m_A68(0.0); - double m_A69(0.0); - double m_A70(0.0); - double m_A71(0.0); - double m_A72(0.0); - double m_A73(0.0); - double m_A74(0.0); - double m_A75(0.0); - double m_A76(0.0); - double m_A77(0.0); - double m_A78(0.0); - double m_A79(0.0); - double m_A80(0.0); - double m_A81(0.0); - double m_A82(0.0); - double m_A83(0.0); - double m_A84(0.0); - double m_A85(0.0); - double m_A86(0.0); - double m_A87(0.0); - double m_A88(0.0); - double m_A89(0.0); - double m_A90(0.0); - double m_A91(0.0); - double m_A92(0.0); - double m_A93(0.0); - double m_A94(0.0); - double m_A95(0.0); - double m_A96(0.0); - double m_A97(0.0); - double m_A98(0.0); - double m_A99(0.0); - double m_A100(0.0); - double m_A101(0.0); - double m_A102(0.0); - double m_A103(0.0); - double m_A104(0.0); - double m_A105(0.0); - double m_A106(0.0); - double m_A107(0.0); - double m_A108(0.0); - double m_A109(0.0); - double m_A110(0.0); - double m_A111(0.0); - double m_A112(0.0); - double m_A113(0.0); - double m_A114(0.0); - double m_A115(0.0); - double m_A116(0.0); - double m_A117(0.0); - double m_A118(0.0); - double m_A119(0.0); - double m_A120(0.0); - double m_A121(0.0); - double m_A122(0.0); - double m_A123(0.0); - double m_A124(0.0); - double m_A125(0.0); - m_A0 += gt[0]; - m_A0 += gt[1]; - m_A1 += go[0]; - double RHS0 = Idr[0]; - RHS0 += Idr[1]; - RHS0 -= go[1] * *cnV[1]; - m_A2 += gt[2]; - m_A2 += gt[3]; - m_A2 += gt[4]; - m_A3 += go[2]; - double RHS1 = Idr[2]; - RHS1 += Idr[3]; - RHS1 += Idr[4]; - RHS1 -= go[3] * *cnV[3]; - RHS1 -= go[4] * *cnV[4]; - m_A4 += gt[5]; - m_A4 += gt[6]; - m_A4 += gt[7]; - m_A5 += go[5]; - double RHS2 = Idr[5]; - RHS2 += Idr[6]; - RHS2 += Idr[7]; - RHS2 -= go[6] * *cnV[6]; - RHS2 -= go[7] * *cnV[7]; - m_A6 += gt[8]; - m_A6 += gt[9]; - m_A6 += gt[10]; - m_A7 += go[8]; - double RHS3 = Idr[8]; - RHS3 += Idr[9]; - RHS3 += Idr[10]; - RHS3 -= go[9] * *cnV[9]; - RHS3 -= go[10] * *cnV[10]; - m_A8 += gt[11]; - m_A8 += gt[12]; - m_A9 += go[11]; - double RHS4 = Idr[11]; - RHS4 += Idr[12]; - RHS4 -= go[12] * *cnV[12]; - m_A10 += gt[13]; - m_A10 += gt[14]; - m_A11 += go[13]; - double RHS5 = Idr[13]; - RHS5 += Idr[14]; - RHS5 -= go[14] * *cnV[14]; - m_A12 += gt[15]; - m_A12 += gt[16]; - m_A14 += go[15]; - m_A13 += go[16]; - double RHS6 = Idr[15]; - RHS6 += Idr[16]; - m_A15 += gt[17]; - m_A15 += gt[18]; - m_A15 += gt[19]; - m_A15 += gt[20]; - m_A15 += gt[21]; - m_A15 += gt[22]; - m_A15 += gt[23]; - m_A16 += go[17]; - double RHS7 = Idr[17]; - RHS7 += Idr[18]; - RHS7 += Idr[19]; - RHS7 += Idr[20]; - RHS7 += Idr[21]; - RHS7 += Idr[22]; - RHS7 += Idr[23]; - RHS7 -= go[18] * *cnV[18]; - RHS7 -= go[19] * *cnV[19]; - RHS7 -= go[20] * *cnV[20]; - RHS7 -= go[21] * *cnV[21]; - RHS7 -= go[22] * *cnV[22]; - RHS7 -= go[23] * *cnV[23]; - m_A17 += gt[24]; - m_A17 += gt[25]; - m_A18 += go[24]; - double RHS8 = Idr[24]; - RHS8 += Idr[25]; - RHS8 -= go[25] * *cnV[25]; - m_A19 += gt[26]; - m_A19 += gt[27]; - m_A20 += go[26]; - double RHS9 = Idr[26]; - RHS9 += Idr[27]; - RHS9 -= go[27] * *cnV[27]; - m_A21 += gt[28]; - m_A21 += gt[29]; - m_A22 += go[28]; - double RHS10 = Idr[28]; - RHS10 += Idr[29]; - RHS10 -= go[29] * *cnV[29]; - m_A23 += gt[30]; - m_A23 += gt[31]; - m_A24 += go[30]; - double RHS11 = Idr[30]; - RHS11 += Idr[31]; - RHS11 -= go[31] * *cnV[31]; - m_A25 += gt[32]; - m_A25 += gt[33]; - m_A25 += gt[34]; - m_A25 += gt[35]; - m_A27 += go[32]; - m_A25 += go[33]; - m_A25 += go[34]; - m_A26 += go[35]; - double RHS12 = Idr[32]; - RHS12 += Idr[33]; - RHS12 += Idr[34]; - RHS12 += Idr[35]; - m_A28 += gt[36]; - m_A28 += gt[37]; - m_A28 += gt[38]; - m_A28 += gt[39]; - m_A29 += go[36]; - m_A28 += go[37]; - m_A28 += go[38]; - m_A30 += go[39]; - double RHS13 = Idr[36]; - RHS13 += Idr[37]; - RHS13 += Idr[38]; - RHS13 += Idr[39]; - m_A31 += gt[40]; - m_A31 += gt[41]; - m_A33 += go[40]; - m_A32 += go[41]; - double RHS14 = Idr[40]; - RHS14 += Idr[41]; - m_A34 += gt[42]; - m_A34 += gt[43]; - m_A36 += go[42]; - m_A35 += go[43]; - double RHS15 = Idr[42]; - RHS15 += Idr[43]; - m_A37 += gt[44]; - m_A37 += gt[45]; - m_A38 += go[44]; - double RHS16 = Idr[44]; - RHS16 += Idr[45]; - RHS16 -= go[45] * *cnV[45]; - m_A39 += gt[46]; - m_A39 += gt[47]; - m_A41 += go[46]; - m_A40 += go[47]; - double RHS17 = Idr[46]; - RHS17 += Idr[47]; - m_A42 += gt[48]; - m_A42 += gt[49]; - m_A44 += go[48]; - m_A43 += go[49]; - double RHS18 = Idr[48]; - RHS18 += Idr[49]; - m_A45 += gt[50]; - m_A45 += gt[51]; - m_A47 += go[50]; - m_A46 += go[51]; - double RHS19 = Idr[50]; - RHS19 += Idr[51]; - m_A48 += gt[52]; - m_A48 += gt[53]; - m_A50 += go[52]; - m_A49 += go[53]; - double RHS20 = Idr[52]; - RHS20 += Idr[53]; - m_A53 += gt[54]; - m_A53 += gt[55]; - m_A53 += gt[56]; - m_A52 += go[54]; - m_A51 += go[55]; - double RHS21 = Idr[54]; - RHS21 += Idr[55]; - RHS21 += Idr[56]; - RHS21 -= go[56] * *cnV[56]; - m_A59 += gt[57]; - m_A59 += gt[58]; - m_A59 += gt[59]; - m_A59 += gt[60]; - m_A59 += gt[61]; - m_A58 += go[57]; - m_A57 += go[58]; - m_A56 += go[59]; - m_A55 += go[60]; - double RHS22 = Idr[57]; - RHS22 += Idr[58]; - RHS22 += Idr[59]; - RHS22 += Idr[60]; - RHS22 += Idr[61]; - RHS22 -= go[61] * *cnV[61]; - m_A63 += gt[62]; - m_A63 += gt[63]; - m_A63 += gt[64]; - m_A63 += gt[65]; - m_A62 += go[62]; - m_A63 += go[63]; - m_A63 += go[64]; - m_A61 += go[65]; - double RHS23 = Idr[62]; - RHS23 += Idr[63]; - RHS23 += Idr[64]; - RHS23 += Idr[65]; - m_A67 += gt[66]; - m_A67 += gt[67]; - m_A67 += gt[68]; - m_A67 += gt[69]; - m_A66 += go[66]; - m_A67 += go[67]; - m_A67 += go[68]; - m_A65 += go[69]; - double RHS24 = Idr[66]; - RHS24 += Idr[67]; - RHS24 += Idr[68]; - RHS24 += Idr[69]; - m_A69 += gt[70]; - m_A69 += gt[71]; - m_A69 += gt[72]; - m_A69 += gt[73]; - m_A69 += gt[74]; - m_A69 += gt[75]; - m_A69 += gt[76]; - m_A70 += go[70]; - double RHS25 = Idr[70]; - RHS25 += Idr[71]; - RHS25 += Idr[72]; - RHS25 += Idr[73]; - RHS25 += Idr[74]; - RHS25 += Idr[75]; - RHS25 += Idr[76]; - RHS25 -= go[71] * *cnV[71]; - RHS25 -= go[72] * *cnV[72]; - RHS25 -= go[73] * *cnV[73]; - RHS25 -= go[74] * *cnV[74]; - RHS25 -= go[75] * *cnV[75]; - RHS25 -= go[76] * *cnV[76]; - m_A73 += gt[77]; - m_A73 += gt[78]; - m_A73 += gt[79]; - m_A73 += gt[80]; - m_A72 += go[77]; - m_A73 += go[78]; - m_A73 += go[79]; - m_A71 += go[80]; - double RHS26 = Idr[77]; - RHS26 += Idr[78]; - RHS26 += Idr[79]; - RHS26 += Idr[80]; - m_A77 += gt[81]; - m_A77 += gt[82]; - m_A77 += gt[83]; - m_A77 += gt[84]; - m_A76 += go[81]; - m_A77 += go[82]; - m_A77 += go[83]; - m_A75 += go[84]; - double RHS27 = Idr[81]; - RHS27 += Idr[82]; - RHS27 += Idr[83]; - RHS27 += Idr[84]; - m_A81 += gt[85]; - m_A81 += gt[86]; - m_A81 += gt[87]; - m_A81 += gt[88]; - m_A80 += go[85]; - m_A81 += go[86]; - m_A81 += go[87]; - m_A79 += go[88]; - double RHS28 = Idr[85]; - RHS28 += Idr[86]; - RHS28 += Idr[87]; - RHS28 += Idr[88]; - m_A85 += gt[89]; - m_A85 += gt[90]; - m_A85 += gt[91]; - m_A86 += go[89]; - m_A83 += go[90]; - double RHS29 = Idr[89]; - RHS29 += Idr[90]; - RHS29 += Idr[91]; - RHS29 -= go[91] * *cnV[91]; - m_A89 += gt[92]; - m_A89 += gt[93]; - m_A90 += go[92]; - m_A87 += go[93]; - double RHS30 = Idr[92]; - RHS30 += Idr[93]; - m_A92 += gt[94]; - m_A92 += gt[95]; - m_A94 += go[94]; - m_A91 += go[95]; - double RHS31 = Idr[94]; - RHS31 += Idr[95]; - m_A97 += gt[96]; - m_A97 += gt[97]; - m_A97 += gt[98]; - m_A97 += gt[99]; - m_A96 += go[96]; - m_A97 += go[97]; - m_A97 += go[98]; - m_A95 += go[99]; - double RHS32 = Idr[96]; - RHS32 += Idr[97]; - RHS32 += Idr[98]; - RHS32 += Idr[99]; - m_A101 += gt[100]; - m_A101 += gt[101]; - m_A101 += gt[102]; - m_A101 += gt[103]; - m_A101 += gt[104]; - m_A100 += go[100]; - m_A99 += go[101]; - m_A102 += go[102]; - double RHS33 = Idr[100]; - RHS33 += Idr[101]; - RHS33 += Idr[102]; - RHS33 += Idr[103]; - RHS33 += Idr[104]; - RHS33 -= go[103] * *cnV[103]; - RHS33 -= go[104] * *cnV[104]; - m_A107 += gt[105]; - m_A107 += gt[106]; - m_A107 += gt[107]; - m_A107 += gt[108]; - m_A107 += gt[109]; - m_A103 += go[105]; - m_A104 += go[106]; - m_A105 += go[107]; - double RHS34 = Idr[105]; - RHS34 += Idr[106]; - RHS34 += Idr[107]; - RHS34 += Idr[108]; - RHS34 += Idr[109]; - RHS34 -= go[108] * *cnV[108]; - RHS34 -= go[109] * *cnV[109]; - m_A125 += gt[110]; - m_A125 += gt[111]; - m_A125 += gt[112]; - m_A125 += gt[113]; - m_A125 += gt[114]; - m_A125 += gt[115]; - m_A125 += gt[116]; - m_A125 += gt[117]; - m_A125 += gt[118]; - m_A125 += gt[119]; - m_A123 += go[110]; - m_A114 += go[111]; - m_A121 += go[112]; - m_A113 += go[113]; - m_A112 += go[114]; - m_A111 += go[115]; - m_A110 += go[116]; - m_A109 += go[117]; - m_A120 += go[118]; - double RHS35 = Idr[110]; - RHS35 += Idr[111]; - RHS35 += Idr[112]; - RHS35 += Idr[113]; - RHS35 += Idr[114]; - RHS35 += Idr[115]; - RHS35 += Idr[116]; - RHS35 += Idr[117]; - RHS35 += Idr[118]; - RHS35 += Idr[119]; - RHS35 -= go[119] * *cnV[119]; - const double f0 = 1.0 / m_A0; - const double f0_23 = -f0 * m_A61; - m_A63 += m_A1 * f0_23; - RHS23 += f0_23 * RHS0; - const double f1 = 1.0 / m_A2; - const double f1_22 = -f1 * m_A55; - m_A59 += m_A3 * f1_22; - RHS22 += f1_22 * RHS1; - const double f2 = 1.0 / m_A4; - const double f2_22 = -f2 * m_A56; - m_A59 += m_A5 * f2_22; - RHS22 += f2_22 * RHS2; - const double f3 = 1.0 / m_A6; - const double f3_22 = -f3 * m_A57; - m_A59 += m_A7 * f3_22; - RHS22 += f3_22 * RHS3; - const double f4 = 1.0 / m_A8; - const double f4_21 = -f4 * m_A51; - m_A53 += m_A9 * f4_21; - RHS21 += f4_21 * RHS4; - const double f5 = 1.0 / m_A10; - const double f5_24 = -f5 * m_A65; - m_A67 += m_A11 * f5_24; - RHS24 += f5_24 * RHS5; - const double f6 = 1.0 / m_A12; - const double f6_21 = -f6 * m_A52; - m_A53 += m_A13 * f6_21; - m_A54 += m_A14 * f6_21; - RHS21 += f6_21 * RHS6; - const double f6_29 = -f6 * m_A83; - m_A84 += m_A13 * f6_29; - m_A85 += m_A14 * f6_29; - RHS29 += f6_29 * RHS6; - const double f7 = 1.0 / m_A15; - const double f7_34 = -f7 * m_A103; - m_A105 += m_A16 * f7_34; - RHS34 += f7_34 * RHS7; - const double f8 = 1.0 / m_A17; - const double f8_26 = -f8 * m_A71; - m_A73 += m_A18 * f8_26; - RHS26 += f8_26 * RHS8; - const double f9 = 1.0 / m_A19; - const double f9_27 = -f9 * m_A75; - m_A77 += m_A20 * f9_27; - RHS27 += f9_27 * RHS9; - const double f10 = 1.0 / m_A21; - const double f10_28 = -f10 * m_A79; - m_A81 += m_A22 * f10_28; - RHS28 += f10_28 * RHS10; - const double f11 = 1.0 / m_A23; - const double f11_32 = -f11 * m_A95; - m_A97 += m_A24 * f11_32; - RHS32 += f11_32 * RHS11; - const double f12 = 1.0 / m_A25; - const double f12_22 = -f12 * m_A58; - m_A59 += m_A26 * f12_22; - m_A60 += m_A27 * f12_22; - RHS22 += f12_22 * RHS12; - const double f12_30 = -f12 * m_A87; - m_A88 += m_A26 * f12_30; - m_A89 += m_A27 * f12_30; - RHS30 += f12_30 * RHS12; - const double f13 = 1.0 / m_A28; - const double f13_31 = -f13 * m_A91; - m_A92 += m_A29 * f13_31; - m_A93 += m_A30 * f13_31; - RHS31 += f13_31 * RHS13; - const double f13_34 = -f13 * m_A104; - m_A106 += m_A29 * f13_34; - m_A107 += m_A30 * f13_34; - RHS34 += f13_34 * RHS13; - const double f14 = 1.0 / m_A31; - const double f14_23 = -f14 * m_A62; - m_A63 += m_A32 * f14_23; - m_A64 += m_A33 * f14_23; - RHS23 += f14_23 * RHS14; - const double f14_35 = -f14 * m_A109; - m_A115 += m_A32 * f14_35; - m_A125 += m_A33 * f14_35; - RHS35 += f14_35 * RHS14; - const double f15 = 1.0 / m_A34; - const double f15_24 = -f15 * m_A66; - m_A67 += m_A35 * f15_24; - m_A68 += m_A36 * f15_24; - RHS24 += f15_24 * RHS15; - const double f15_35 = -f15 * m_A110; - m_A116 += m_A35 * f15_35; - m_A125 += m_A36 * f15_35; - RHS35 += f15_35 * RHS15; - const double f16 = 1.0 / m_A37; - const double f16_33 = -f16 * m_A99; - m_A101 += m_A38 * f16_33; - RHS33 += f16_33 * RHS16; - const double f17 = 1.0 / m_A39; - const double f17_26 = -f17 * m_A72; - m_A73 += m_A40 * f17_26; - m_A74 += m_A41 * f17_26; - RHS26 += f17_26 * RHS17; - const double f17_35 = -f17 * m_A111; - m_A117 += m_A40 * f17_35; - m_A125 += m_A41 * f17_35; - RHS35 += f17_35 * RHS17; - const double f18 = 1.0 / m_A42; - const double f18_27 = -f18 * m_A76; - m_A77 += m_A43 * f18_27; - m_A78 += m_A44 * f18_27; - RHS27 += f18_27 * RHS18; - const double f18_35 = -f18 * m_A112; - m_A118 += m_A43 * f18_35; - m_A125 += m_A44 * f18_35; - RHS35 += f18_35 * RHS18; - const double f19 = 1.0 / m_A45; - const double f19_28 = -f19 * m_A80; - m_A81 += m_A46 * f19_28; - m_A82 += m_A47 * f19_28; - RHS28 += f19_28 * RHS19; - const double f19_35 = -f19 * m_A113; - m_A119 += m_A46 * f19_35; - m_A125 += m_A47 * f19_35; - RHS35 += f19_35 * RHS19; - const double f20 = 1.0 / m_A48; - const double f20_32 = -f20 * m_A96; - m_A97 += m_A49 * f20_32; - m_A98 += m_A50 * f20_32; - RHS32 += f20_32 * RHS20; - const double f20_35 = -f20 * m_A114; - m_A122 += m_A49 * f20_35; - m_A125 += m_A50 * f20_35; - RHS35 += f20_35 * RHS20; - const double f21 = 1.0 / m_A53; - const double f21_29 = -f21 * m_A84; - m_A85 += m_A54 * f21_29; - RHS29 += f21_29 * RHS21; - const double f22 = 1.0 / m_A59; - const double f22_30 = -f22 * m_A88; - m_A89 += m_A60 * f22_30; - RHS30 += f22_30 * RHS22; - const double f23 = 1.0 / m_A63; - const double f23_35 = -f23 * m_A115; - m_A125 += m_A64 * f23_35; - RHS35 += f23_35 * RHS23; - const double f24 = 1.0 / m_A67; - const double f24_35 = -f24 * m_A116; - m_A125 += m_A68 * f24_35; - RHS35 += f24_35 * RHS24; - const double f25 = 1.0 / m_A69; - const double f25_33 = -f25 * m_A100; - m_A102 += m_A70 * f25_33; - RHS33 += f25_33 * RHS25; - const double f26 = 1.0 / m_A73; - const double f26_35 = -f26 * m_A117; - m_A125 += m_A74 * f26_35; - RHS35 += f26_35 * RHS26; - const double f27 = 1.0 / m_A77; - const double f27_35 = -f27 * m_A118; - m_A125 += m_A78 * f27_35; - RHS35 += f27_35 * RHS27; - const double f28 = 1.0 / m_A81; - const double f28_35 = -f28 * m_A119; - m_A125 += m_A82 * f28_35; - RHS35 += f28_35 * RHS28; - const double f29 = 1.0 / m_A85; - const double f29_34 = -f29 * m_A105; - m_A107 += m_A86 * f29_34; - RHS34 += f29_34 * RHS29; - const double f30 = 1.0 / m_A89; - const double f30_35 = -f30 * m_A120; - m_A125 += m_A90 * f30_35; - RHS35 += f30_35 * RHS30; - const double f31 = 1.0 / m_A92; - const double f31_34 = -f31 * m_A106; - m_A107 += m_A93 * f31_34; - m_A108 += m_A94 * f31_34; - RHS34 += f31_34 * RHS31; - const double f31_35 = -f31 * m_A121; - m_A124 += m_A93 * f31_35; - m_A125 += m_A94 * f31_35; - RHS35 += f31_35 * RHS31; - const double f32 = 1.0 / m_A97; - const double f32_35 = -f32 * m_A122; - m_A125 += m_A98 * f32_35; - RHS35 += f32_35 * RHS32; - const double f33 = 1.0 / m_A101; - const double f33_35 = -f33 * m_A123; - m_A125 += m_A102 * f33_35; - RHS35 += f33_35 * RHS33; - const double f34 = 1.0 / m_A107; - const double f34_35 = -f34 * m_A124; - m_A125 += m_A108 * f34_35; - RHS35 += f34_35 * RHS34; - V[35] = RHS35 / m_A125; - double tmp34 = 0.0; - tmp34 += m_A108 * V[35]; - V[34] = (RHS34 - tmp34) / m_A107; - double tmp33 = 0.0; - tmp33 += m_A102 * V[35]; - V[33] = (RHS33 - tmp33) / m_A101; - double tmp32 = 0.0; - tmp32 += m_A98 * V[35]; - V[32] = (RHS32 - tmp32) / m_A97; - double tmp31 = 0.0; - tmp31 += m_A93 * V[34]; - tmp31 += m_A94 * V[35]; - V[31] = (RHS31 - tmp31) / m_A92; - double tmp30 = 0.0; - tmp30 += m_A90 * V[35]; - V[30] = (RHS30 - tmp30) / m_A89; - double tmp29 = 0.0; - tmp29 += m_A86 * V[34]; - V[29] = (RHS29 - tmp29) / m_A85; - double tmp28 = 0.0; - tmp28 += m_A82 * V[35]; - V[28] = (RHS28 - tmp28) / m_A81; - double tmp27 = 0.0; - tmp27 += m_A78 * V[35]; - V[27] = (RHS27 - tmp27) / m_A77; - double tmp26 = 0.0; - tmp26 += m_A74 * V[35]; - V[26] = (RHS26 - tmp26) / m_A73; - double tmp25 = 0.0; - tmp25 += m_A70 * V[35]; - V[25] = (RHS25 - tmp25) / m_A69; - double tmp24 = 0.0; - tmp24 += m_A68 * V[35]; - V[24] = (RHS24 - tmp24) / m_A67; - double tmp23 = 0.0; - tmp23 += m_A64 * V[35]; - V[23] = (RHS23 - tmp23) / m_A63; - double tmp22 = 0.0; - tmp22 += m_A60 * V[30]; - V[22] = (RHS22 - tmp22) / m_A59; - double tmp21 = 0.0; - tmp21 += m_A54 * V[29]; - V[21] = (RHS21 - tmp21) / m_A53; - double tmp20 = 0.0; - tmp20 += m_A49 * V[32]; - tmp20 += m_A50 * V[35]; - V[20] = (RHS20 - tmp20) / m_A48; - double tmp19 = 0.0; - tmp19 += m_A46 * V[28]; - tmp19 += m_A47 * V[35]; - V[19] = (RHS19 - tmp19) / m_A45; - double tmp18 = 0.0; - tmp18 += m_A43 * V[27]; - tmp18 += m_A44 * V[35]; - V[18] = (RHS18 - tmp18) / m_A42; - double tmp17 = 0.0; - tmp17 += m_A40 * V[26]; - tmp17 += m_A41 * V[35]; - V[17] = (RHS17 - tmp17) / m_A39; - double tmp16 = 0.0; - tmp16 += m_A38 * V[33]; - V[16] = (RHS16 - tmp16) / m_A37; - double tmp15 = 0.0; - tmp15 += m_A35 * V[24]; - tmp15 += m_A36 * V[35]; - V[15] = (RHS15 - tmp15) / m_A34; - double tmp14 = 0.0; - tmp14 += m_A32 * V[23]; - tmp14 += m_A33 * V[35]; - V[14] = (RHS14 - tmp14) / m_A31; - double tmp13 = 0.0; - tmp13 += m_A29 * V[31]; - tmp13 += m_A30 * V[34]; - V[13] = (RHS13 - tmp13) / m_A28; - double tmp12 = 0.0; - tmp12 += m_A26 * V[22]; - tmp12 += m_A27 * V[30]; - V[12] = (RHS12 - tmp12) / m_A25; - double tmp11 = 0.0; - tmp11 += m_A24 * V[32]; - V[11] = (RHS11 - tmp11) / m_A23; - double tmp10 = 0.0; - tmp10 += m_A22 * V[28]; - V[10] = (RHS10 - tmp10) / m_A21; - double tmp9 = 0.0; - tmp9 += m_A20 * V[27]; - V[9] = (RHS9 - tmp9) / m_A19; - double tmp8 = 0.0; - tmp8 += m_A18 * V[26]; - V[8] = (RHS8 - tmp8) / m_A17; - double tmp7 = 0.0; - tmp7 += m_A16 * V[29]; - V[7] = (RHS7 - tmp7) / m_A15; - double tmp6 = 0.0; - tmp6 += m_A13 * V[21]; - tmp6 += m_A14 * V[29]; - V[6] = (RHS6 - tmp6) / m_A12; - double tmp5 = 0.0; - tmp5 += m_A11 * V[24]; - V[5] = (RHS5 - tmp5) / m_A10; - double tmp4 = 0.0; - tmp4 += m_A9 * V[21]; - V[4] = (RHS4 - tmp4) / m_A8; - double tmp3 = 0.0; - tmp3 += m_A7 * V[22]; - V[3] = (RHS3 - tmp3) / m_A6; - double tmp2 = 0.0; - tmp2 += m_A5 * V[22]; - V[2] = (RHS2 - tmp2) / m_A4; - double tmp1 = 0.0; - tmp1 += m_A3 * V[22]; - V[1] = (RHS1 - tmp1) / m_A2; - double tmp0 = 0.0; - tmp0 += m_A1 * V[23]; - V[0] = (RHS0 - tmp0) / m_A0; -} - // boxingb static void nl_gcr_d7d45dc58b08cab9_10_double_double(double * __restrict V, const double * __restrict go, const double * __restrict gt, const double * __restrict Idr, const double * const * __restrict cnV) @@ -80455,8 +83059,12 @@ const plib::dynlib_static_sym nl_static_solver_syms[] = { {"nl_gcr_1f1086787c94f97c_40_double_double", reinterpret_cast<void *>(&nl_gcr_1f1086787c94f97c_40_double_double)}, // NOLINT // zac1b11142 {"nl_gcr_1fad5cda2646cf42_30_double_double", reinterpret_cast<void *>(&nl_gcr_1fad5cda2646cf42_30_double_double)}, // NOLINT +// brdrline + {"nl_gcr_1fec893028a7b809_34_double_double", reinterpret_cast<void *>(&nl_gcr_1fec893028a7b809_34_double_double)}, // NOLINT // sspeedr {"nl_gcr_2294220d3c91e762_176_double_double", reinterpret_cast<void *>(&nl_gcr_2294220d3c91e762_176_double_double)}, // NOLINT +// headon + {"nl_gcr_24023c8529f6bdb3_97_double_double", reinterpret_cast<void *>(&nl_gcr_24023c8529f6bdb3_97_double_double)}, // NOLINT // 280zzzap {"nl_gcr_24643c159711f292_95_double_double", reinterpret_cast<void *>(&nl_gcr_24643c159711f292_95_double_double)}, // NOLINT // frogs @@ -80530,6 +83138,8 @@ const plib::dynlib_static_sym nl_static_solver_syms[] = { // cheekyms {"nl_gcr_546396f65ce48700_12_double_double", reinterpret_cast<void *>(&nl_gcr_546396f65ce48700_12_double_double)}, // NOLINT // brdrline + {"nl_gcr_576f391d0418cb76_13_double_double", reinterpret_cast<void *>(&nl_gcr_576f391d0418cb76_13_double_double)}, // NOLINT +// brdrline {"nl_gcr_59cb6bf7cb9d17dc_7_double_double", reinterpret_cast<void *>(&nl_gcr_59cb6bf7cb9d17dc_7_double_double)}, // NOLINT // breakout {"nl_gcr_5a3419e2809520de_13_double_double", reinterpret_cast<void *>(&nl_gcr_5a3419e2809520de_13_double_double)}, // NOLINT @@ -80737,12 +83347,16 @@ const plib::dynlib_static_sym nl_static_solver_syms[] = { {"nl_gcr_c2e616f3de30f15b_31_double_double", reinterpret_cast<void *>(&nl_gcr_c2e616f3de30f15b_31_double_double)}, // NOLINT // tailg {"nl_gcr_c4cec7aed23b7b94_23_double_double", reinterpret_cast<void *>(&nl_gcr_c4cec7aed23b7b94_23_double_double)}, // NOLINT +// brdrline + {"nl_gcr_c5b523ed6a9677e6_21_double_double", reinterpret_cast<void *>(&nl_gcr_c5b523ed6a9677e6_21_double_double)}, // NOLINT // segausb {"nl_gcr_c61e08cf5e35918_84_double_double", reinterpret_cast<void *>(&nl_gcr_c61e08cf5e35918_84_double_double)}, // NOLINT // popeye {"nl_gcr_c6f25bb06e161d1c_50_double_double", reinterpret_cast<void *>(&nl_gcr_c6f25bb06e161d1c_50_double_double)}, // NOLINT // flyball {"nl_gcr_c74b1a65978d7121_7_double_double", reinterpret_cast<void *>(&nl_gcr_c74b1a65978d7121_7_double_double)}, // NOLINT +// headon + {"nl_gcr_c90daff82be2d052_350_double_double", reinterpret_cast<void *>(&nl_gcr_c90daff82be2d052_350_double_double)}, // NOLINT // stuntcyc {"nl_gcr_c924fe5960b1479e_20_double_double", reinterpret_cast<void *>(&nl_gcr_c924fe5960b1479e_20_double_double)}, // NOLINT // dpatrol @@ -80779,8 +83393,6 @@ const plib::dynlib_static_sym nl_static_solver_syms[] = { {"nl_gcr_d4c34516ff6aa139_46_double_double", reinterpret_cast<void *>(&nl_gcr_d4c34516ff6aa139_46_double_double)}, // NOLINT // brdrline {"nl_gcr_d5e1c37cfa2d2853_30_double_double", reinterpret_cast<void *>(&nl_gcr_d5e1c37cfa2d2853_30_double_double)}, // NOLINT -// brdrline - {"nl_gcr_d74d89ba31f2cb81_126_double_double", reinterpret_cast<void *>(&nl_gcr_d74d89ba31f2cb81_126_double_double)}, // NOLINT // boxingb {"nl_gcr_d7d45dc58b08cab9_10_double_double", reinterpret_cast<void *>(&nl_gcr_d7d45dc58b08cab9_10_double_double)}, // NOLINT // kidniki diff --git a/src/mame/audio/nl_brdrline.cpp b/src/mame/audio/nl_brdrline.cpp index 61a1b93ed3b..e08121f90fa 100644 --- a/src/mame/audio/nl_brdrline.cpp +++ b/src/mame/audio/nl_brdrline.cpp @@ -19,7 +19,7 @@ // * it takes up to 30seconds to charge C24 to its working point // in ANIMAL_SOUND // -// * noise generation with MM5837_DIP eats a lot of CPU power +// * noise generation with MM5837_DIP eats a lot of CPU power // @@ -63,6 +63,8 @@ #include "netlist/devices/net_lib.h" +#define ENABLE_FRONTIERS (1) + /* ---------------------------------------------------------------------------- * Library section content START @@ -1422,6 +1424,41 @@ static NETLIST_START(brdrline_sound_out) // -------------------------------- // INPUT +#if (ENABLE_FRONTIERS) + // using AFUNCs here tends to remove the DC bias, so add + // it back in manually + AFUNC(JEEP_F, 1, "5+A0") + ALIAS(JEEP_SOUND, JEEP_F.A0) + NET_C(JEEP_F.Q, VR8.1) + + AFUNC(GUN_F, 1, "6+A0") + ALIAS(GUN_SOUND, GUN_F.A0) + NET_C(GUN_F.Q, VR1.1) + + AFUNC(POINT_F, 1, "6+A0") + ALIAS(POINT_SOUND, POINT_F.A0) + NET_C(POINT_F.Q, VR6.1) + + AFUNC(HIT_F, 1, "6+A0") + ALIAS(HIT_SOUND, HIT_F.A0) + NET_C(HIT_F.Q, VR7.1) + + AFUNC(ANIMAL_F, 1, "6+A0") + ALIAS(ANIMAL_SOUND, ANIMAL_F.A0) + NET_C(ANIMAL_F.Q, VR5.1) + + AFUNC(EMAR_F, 1, "6+A0") + ALIAS(EMAR_SOUND, EMAR_F.A0) + NET_C(EMAR_F.Q, VR4.1) + + AFUNC(WALK_F, 1, "A0") + ALIAS(WALK_SOUND, WALK_F.A0) + NET_C(WALK_F.Q, VR3.1) + + AFUNC(CRY_F, 1, "6+A0") + ALIAS(CRY_SOUND, CRY_F.A0) + NET_C(CRY_F.Q, VR2.1) +#else ALIAS(JEEP_SOUND, VR8.1) ALIAS(GUN_SOUND, VR1.1) ALIAS(POINT_SOUND, VR6.1) @@ -1430,6 +1467,7 @@ static NETLIST_START(brdrline_sound_out) ALIAS(EMAR_SOUND, VR4.1) ALIAS(WALK_SOUND, VR3.1) ALIAS(CRY_SOUND, VR2.1) +#endif // -------------------------------- // OUTPUT @@ -1454,7 +1492,7 @@ NETLIST_START(brdrline) #if 1 SOLVER(Solver, 1000) PARAM(Solver.DYNAMIC_TS, 1) - PARAM(Solver.DYNAMIC_MIN_TIMESTEP, 2e-5) + PARAM(Solver.DYNAMIC_MIN_TIMESTEP, 4e-5) #else SOLVER(solver, 48000) #endif |