ge_msub.h 1.9 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788
  1. /* qhasm: enter ge_msub */
  2. /* qhasm: fe X1 */
  3. /* qhasm: fe Y1 */
  4. /* qhasm: fe Z1 */
  5. /* qhasm: fe T1 */
  6. /* qhasm: fe ypx2 */
  7. /* qhasm: fe ymx2 */
  8. /* qhasm: fe xy2d2 */
  9. /* qhasm: fe X3 */
  10. /* qhasm: fe Y3 */
  11. /* qhasm: fe Z3 */
  12. /* qhasm: fe T3 */
  13. /* qhasm: fe YpX1 */
  14. /* qhasm: fe YmX1 */
  15. /* qhasm: fe A */
  16. /* qhasm: fe B */
  17. /* qhasm: fe C */
  18. /* qhasm: fe D */
  19. /* qhasm: YpX1 = Y1+X1 */
  20. /* asm 1: fe_add(>YpX1=fe#1,<Y1=fe#12,<X1=fe#11); */
  21. /* asm 2: fe_add(>YpX1=r->X,<Y1=p->Y,<X1=p->X); */
  22. fe_add(r->X,p->Y,p->X);
  23. /* qhasm: YmX1 = Y1-X1 */
  24. /* asm 1: fe_sub(>YmX1=fe#2,<Y1=fe#12,<X1=fe#11); */
  25. /* asm 2: fe_sub(>YmX1=r->Y,<Y1=p->Y,<X1=p->X); */
  26. fe_sub(r->Y,p->Y,p->X);
  27. /* qhasm: A = YpX1*ymx2 */
  28. /* asm 1: fe_mul(>A=fe#3,<YpX1=fe#1,<ymx2=fe#16); */
  29. /* asm 2: fe_mul(>A=r->Z,<YpX1=r->X,<ymx2=q->yminusx); */
  30. fe_mul(r->Z,r->X,q->yminusx);
  31. /* qhasm: B = YmX1*ypx2 */
  32. /* asm 1: fe_mul(>B=fe#2,<YmX1=fe#2,<ypx2=fe#15); */
  33. /* asm 2: fe_mul(>B=r->Y,<YmX1=r->Y,<ypx2=q->yplusx); */
  34. fe_mul(r->Y,r->Y,q->yplusx);
  35. /* qhasm: C = xy2d2*T1 */
  36. /* asm 1: fe_mul(>C=fe#4,<xy2d2=fe#17,<T1=fe#14); */
  37. /* asm 2: fe_mul(>C=r->T,<xy2d2=q->xy2d,<T1=p->T); */
  38. fe_mul(r->T,q->xy2d,p->T);
  39. /* qhasm: D = 2*Z1 */
  40. /* asm 1: fe_add(>D=fe#5,<Z1=fe#13,<Z1=fe#13); */
  41. /* asm 2: fe_add(>D=t0,<Z1=p->Z,<Z1=p->Z); */
  42. fe_add(t0,p->Z,p->Z);
  43. /* qhasm: X3 = A-B */
  44. /* asm 1: fe_sub(>X3=fe#1,<A=fe#3,<B=fe#2); */
  45. /* asm 2: fe_sub(>X3=r->X,<A=r->Z,<B=r->Y); */
  46. fe_sub(r->X,r->Z,r->Y);
  47. /* qhasm: Y3 = A+B */
  48. /* asm 1: fe_add(>Y3=fe#2,<A=fe#3,<B=fe#2); */
  49. /* asm 2: fe_add(>Y3=r->Y,<A=r->Z,<B=r->Y); */
  50. fe_add(r->Y,r->Z,r->Y);
  51. /* qhasm: Z3 = D-C */
  52. /* asm 1: fe_sub(>Z3=fe#3,<D=fe#5,<C=fe#4); */
  53. /* asm 2: fe_sub(>Z3=r->Z,<D=t0,<C=r->T); */
  54. fe_sub(r->Z,t0,r->T);
  55. /* qhasm: T3 = D+C */
  56. /* asm 1: fe_add(>T3=fe#4,<D=fe#5,<C=fe#4); */
  57. /* asm 2: fe_add(>T3=r->T,<D=t0,<C=r->T); */
  58. fe_add(r->T,t0,r->T);
  59. /* qhasm: return */