ge_sub.q 773 B

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849
  1. :name:fe:r->X:r->Y:r->Z:r->T:t0:t1:t2:t3:t4:t5:p->X:p->Y:p->Z:p->T:q->YplusX:q->YminusX:q->Z:q->T2d:
  2. fe r:var/r=fe:
  3. enter f:enter/f:>X1=fe#11:>Y1=fe#12:>Z1=fe#13:>T1=fe#14:>YpX2=fe#15:>YmX2=fe#16:>Z2=fe#17:>T2d2=fe#18:
  4. return:nofallthrough:<X3=fe#1:<Y3=fe#2:<Z3=fe#3:<T3=fe#4:leave:
  5. h=f+g:<f=fe:<g=fe:>h=fe:asm/fe_add(>h,<f,<g);:
  6. h=f-g:<f=fe:<g=fe:>h=fe:asm/fe_sub(>h,<f,<g);:
  7. h=f*g:<f=fe:<g=fe:>h=fe:asm/fe_mul(>h,<f,<g);:
  8. h=f^2:<f=fe:>h=fe:asm/fe_sq(>h,<f);:
  9. h=2*g:<g=fe:>h=fe:asm/fe_add(>h,<g,<g);:
  10. :
  11. enter ge_sub
  12. fe X1
  13. fe Y1
  14. fe Z1
  15. fe Z2
  16. fe T1
  17. fe ZZ
  18. fe YpX2
  19. fe YmX2
  20. fe T2d2
  21. fe X3
  22. fe Y3
  23. fe Z3
  24. fe T3
  25. fe YpX1
  26. fe YmX1
  27. fe A
  28. fe B
  29. fe C
  30. fe D
  31. YpX1 = Y1+X1
  32. YmX1 = Y1-X1
  33. A = YpX1*YmX2
  34. B = YmX1*YpX2
  35. C = T2d2*T1
  36. ZZ = Z1*Z2
  37. D = 2*ZZ
  38. X3 = A-B
  39. Y3 = A+B
  40. Z3 = D-C
  41. T3 = D+C
  42. return