calcSP.c 1.2 KB

12345678910111213141516171819202122232425262728293031323334353637383940
  1. t2 = SF(1.5E1);
  2. t3 = SF(1.4E1);
  3. t4 = SF(1.3E1);
  4. t5 = SF(1.2E1);
  5. t6 = SF(1.1E1);
  6. t7 = SF(1.0E1);
  7. t8 = SF(9.0);
  8. t9 = SF(8.0);
  9. t10 = SF(7.0);
  10. t11 = SF(6.0);
  11. t12 = SF(5.0);
  12. t13 = SF(4.0);
  13. t14 = SF(1.7E1);
  14. t15 = SF(2.3E1);
  15. t16 = SF(2.0E1);
  16. t17 = SF(1.8E1);
  17. t18 = SF(2.1E1);
  18. t19 = SF(2.4E1);
  19. t20 = SF(2.5E1);
  20. t21 = SF(2.2E1);
  21. t22 = SF(1.9E1);
  22. A0[0][0] = q0*t3*-2.0+q1*t2*2.0+q2*t5*2.0+q3*t4*2.0;
  23. A0[1][0] = q0*t7*2.0+q1*t6*2.0+q2*t9*2.0-q3*t8*2.0;
  24. A0[2][0] = q0*t2*2.0+q1*t3*2.0-q2*t4*2.0+q3*t5*2.0;
  25. A0[3][0] = q0*t6*2.0-q1*t7*2.0+q2*t8*2.0+q3*t9*2.0;
  26. A0[4][0] = q0*t5*2.0+q1*t4*2.0+q2*t3*2.0-q3*t2*2.0;
  27. A0[5][0] = q0*t9*-2.0+q1*t8*2.0+q2*t7*2.0+q3*t6*2.0;
  28. A0[6][0] = q0*t11*2.0+q1*t10*2.0-q2*t13*2.0+q3*t12*2.0;
  29. A0[7][0] = q0*t10*-2.0+q1*t11*2.0+q2*t12*2.0+q3*t13*2.0;
  30. A0[8][0] = q0*t13*2.0-q1*t12*2.0+q2*t11*2.0+q3*t10*2.0;
  31. A0[9][0] = dvy*t14-dvz*t15;
  32. A0[10][0] = dvx*t17-dvy*t19;
  33. A0[11][0] = dvx*t20-dvz*t22;
  34. A0[12][0] = dvx*t14+dvz*t16;
  35. A0[13][0] = dvx*t15+dvy*t16;
  36. A0[14][0] = dvy*t18+dvz*t17;
  37. A0[15][0] = dvx*t18+dvz*t19;
  38. A0[16][0] = dvy*t20+dvz*t21;
  39. A0[17][0] = dvx*t21+dvy*t22;
  40. A0[18][0] = SF(1.6E1);