Changes between Version 37 and Version 38 of PolynomialExpansion

Timestamp:

01/26/16 19:24:01 (10 years ago)

Author:

sili

Comment:

--

PolynomialExpansion

@@ -14 +14 @@
 
 Step 3: r,,1,, = r,,0,, - alpha,,0,,  Ap,,0,,
-  r[0] = (-b,,1,,(-a,,1,,b,,0,,^2^ + a,,0,,b,,0,,b,,1,, - a,,2,,b,,0,,b,,1,, + a,,1,,b,,1,,^2^)) / (a,,0,,b,,0,,^2^ + 2a,,1,,b,,0,,b,,1,, + a,,2,,b,,1,,^2^)
+  r[0] = (-b,,1,,(-a,,1,,b,,0,,^2^ + a,,0,,b,,0,,b,,1,, - a,,2,,b,,0,,b,,1,, + a,,1,,b,,1,,^2^)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))
 
-  r[1] = (b,,0,,(-a,,1,,b,,0,,^2^ + a,,0,,b,,0,,b,,1,, - a,,2,,b,,0,,b,,1,, + a,,1,,b,,1,,^2^)) / (a,,0,,b,,0,,^2^ + 2a,,1,,b,,0,,b,,1,, + a,,2,,b,,1,,^2^) 
+  r[1] = (b,,0,,(-a,,1,,b,,0,,^2^ + a,,0,,b,,0,,b,,1,, - a,,2,,b,,0,,b,,1,, + a,,1,,b,,1,,^2^)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))
 
 Step 4: x,,1,, = x,,0,, + alphap,,0,,
-  x[0] = (b,,0,,(b,,0,,^2^ + b,,1,,^2^)) / (a,,0,,b,,0,,^2^ + 2a,,1,,b,,0,,b,,1,, + a,,2,,b,,1,,^2^)
+  x[0] = (b,,0,,(b,,0,,^2^ + b,,1,,^2^)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))
 
-  x[1] = (b,,1,,(b,,0,,^2^ + b,,1,,^2^)) / (a,,0,,b,,0,,^2^ + 2a,,1,,b,,0,,b,,1,, + a,,2,,b,,1,,^2^)
+  x[1] = (b,,1,,(b,,0,,^2^ + b,,1,,^2^)) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))
 
 Step 5: beta = rsnew / rsold = <r,,1,,, r,,1,,> / <r,,0,,, r,,0,,>
-  rsnew = ((b,,0,,^2^ + b,,1,,^2^)(-a,,1,,b,,0,,^2^ + a,,0,,b,,0,,b,,1,, - a,,2,,b,,0,,b,,1,, + a,,1,,b,,1,,^2^)^2^) / (a,,0,,b,,0,,^2^ + 2a,,1,,b,,0,,b,,1,, + a,,2,,b,,1,,^2^)^2^
+  rsnew = ((b,,0,,^2^ + b,,1,,^2^)(-a,,1,,b,,0,,^2^ + a,,0,,b,,0,,b,,1,, - a,,2,,b,,0,,b,,1,, + a,,1,,b,,1,,^2^)^2^) / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))^2^
 
-  beta = (a,,1,,b,,0,,^2^ - a,,0,,b,,0,,b,,1,, + a,,2,,b,,0,,b,,1,, - a,,1,,b,,1,,^2^)^2^ / (a,,0,,b,,0,,^2^ + 2a,,1,,b,,0,,b,,1,, + a,,2,,b,,1,,^2^)^2^
+  beta = (a,,1,,b,,0,,^2^ - a,,0,,b,,0,,b,,1,, + a,,2,,b,,0,,b,,1,, - a,,1,,b,,1,,^2^)^2^ / (b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,))^2^
 
 Step 6: p,,1,, = r,,1,, +beta  p,,0,,