Changes between Version 55 and Version 56 of PolynomialExpansion


Ignore:
Timestamp:
01/26/16 23:35:16 (10 years ago)
Author:
sili
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • PolynomialExpansion

    v55 v56  
    4040Step 8: r,,2,, = r,,1,, - alpha,,1,,*Ap,,1,,  ('''Need to expand for cancellation''')
    4141
    42   r[0] = (b,,0,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - alpha,,1,,*(a,,0,,p[0]+a,,1,,p[1])
     42  r[0] = (b,,0,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (a,,0,,p[0]+a,,1,,p[1])((b,,0,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,))^2^ + (b,,1,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,))^2^) / (p[0](a,,0,,p[0]+a,,1,,p[1]) + p[1](a,,1,,p[0]+a,,2,,p[1]))(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))^2^ = 0
    4343
    44   r[1] = (b,,1,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - alpha,,1,,*(a,,1,,p[0]+a,,2,,p[1])
     44  r[1] = (b,,1,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (a,,1,,p[0]+a,,2,,p[1])((b,,0,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,))^2^ + (b,,1,,(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)) - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,))^2^) / (p[0](a,,0,,p[0]+a,,1,,p[1]) + p[1](a,,1,,p[0]+a,,2,,p[1]))(b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))^2^ = 0
    4545
    4646Step 9: x,,2,, = x,,1,, + alpha,,1,,*p,,1,,