calcSF.c 1.3 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051
  1. t2 = SF(1.0);
  2. t3 = SF(2.0);
  3. t4 = SF(3.0);
  4. t5 = q3*(1.0/2.0);
  5. t6 = q0*t2*(1.0/2.0);
  6. t7 = q2*t4*(1.0/2.0);
  7. t8 = q0*t3*(1.0/2.0);
  8. t9 = q1*t2*(1.0/2.0);
  9. t10 = q3*t4*(1.0/2.0);
  10. t11 = q1*(1.0/2.0);
  11. t12 = q2*t2*(1.0/2.0);
  12. t13 = q3*t3*(1.0/2.0);
  13. t14 = q0*(1.0/2.0);
  14. t15 = q1*t4*(1.0/2.0);
  15. t16 = q2*t3*(1.0/2.0);
  16. t17 = q3*t2*(1.0/2.0);
  17. t18 = q0*t4*(1.0/2.0);
  18. t19 = q2*(1.0/2.0);
  19. t20 = q1*t3*(1.0/2.0);
  20. t21 = q0*q0;
  21. t22 = q1*q1;
  22. t23 = q2*q2;
  23. t24 = q3*q3;
  24. t25 = q0*q1*2.0;
  25. t26 = q0*q3*2.0;
  26. t27 = q0*q2*2.0;
  27. A0[0][0] = daz*(1.0/2.0)-daz_b*(1.0/2.0);
  28. A0[1][0] = day*(1.0/2.0)-day_b*(1.0/2.0);
  29. A0[2][0] = dax*(1.0/2.0)-dax_b*(1.0/2.0);
  30. A0[3][0] = q2*(-1.0/2.0)+t8+t9+t10;
  31. A0[4][0] = t5+t6+t7-q1*t3*(1.0/2.0);
  32. A0[5][0] = t14+t15+t16-q3*t2*(1.0/2.0);
  33. A0[6][0] = t11+t12+t13-q0*t4*(1.0/2.0);
  34. A0[7][0] = t5-t6+t7+t20;
  35. A0[8][0] = t8+t9-t10+t19;
  36. A0[9][0] = -t11+t12+t13+t18;
  37. A0[10][0] = t14+t15-t16+t17;
  38. A0[11][0] = t14-t15+t16+t17;
  39. A0[12][0] = t11-t12+t13+t18;
  40. A0[13][0] = -t8+t9+t10+t19;
  41. A0[14][0] = -t5+t6+t7+t20;
  42. A0[15][0] = -t21-t22-t23-t24;
  43. A0[16][0] = t21-t22-t23+t24;
  44. A0[17][0] = t21-t22+t23-t24;
  45. A0[18][0] = t21+t22-t23-t24;
  46. A0[19][0] = t27-q1*q3*2.0;
  47. A0[20][0] = t25-q2*q3*2.0;
  48. A0[21][0] = t26-q1*q2*2.0;
  49. A0[22][0] = t25+q2*q3*2.0;
  50. A0[23][0] = t26+q1*q2*2.0;
  51. A0[24][0] = t27+q1*q3*2.0;