S11
FF (a, b, c, d, x[ 0], S11); /* 1 */
FF (a, b, c, d, x[ 4], S11); /* 5 */
FF (a, b, c, d, x[ 8], S11); /* 9 */
FF (a, b, c, d, x[12], S11); /* 13 */
FF (a, b, c, d, x[ 0], S11, 0xd76aa478); /* 1 */
FF (a, b, c, d, x[ 4], S11, 0xf57c0faf); /* 5 */
FF (a, b, c, d, x[ 8], S11, 0x698098d8); /* 9 */
FF (a, b, c, d, x[12], S11, 0x6b901122); /* 13 */
S11 = -0.38680813379701966970673724299207480965452616911420e-22,
z*(S9+z*(S10+z*(S11+z*S12)))))))));
RETURNI(((S17*x2 + S15)*x4 + (S13*x2 + S11))*(x2*x*x4*x4) +
S17)*x2 + S15)*x2 + S13)*x2 + S11)*x2 + S9)*x2 + S7)*x2 +
S11 = 2.5052108406704084e-8, /* 0x1ae6456857530f.0p-78 */
S11 = 2.50521083854417187749675637460977997e-8L, /* 0x1ae64567f544e38fe59b3380d7413.0p-138L */
S11 = -0.38680813379701966970673724299207480965452616911420e-22,
z * (S9 + z * (S10 + z * (S11 + z * S12)))))))));
M = (S11 * S22 - S12 * S12);
B = (S2Y * S11 - S1Y * S12) * scale / M;
int64_t S11, S22, S12;
S11 = X1X1s - n * X1a * X1a;
printf("S11=%d S22=%d S12=%d\n", (int)S11, (int)S22, (int)S12);