Changes between Version 74 and Version 75 of PolynomialExpansion

Timestamp:

01/27/16 16:56:50 (10 years ago)

Author:

sili

Comment:

--

PolynomialExpansion

@@ -24 +24 @@
 
 Step 5: beta = rsnew / rsold = <r,,1,,, r,,1,,> / <r,,0,,, r,,0,,>  ('''No expansion''')
-  to make it simple to read, let's devote m = b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)
+  **to make it simple to read, let's devote m = b,,0,,(a,,0,,b,,0,, + a,,1,,b,,1,,) + b,,1,,(a,,1,,b,,0,, + a,,2,,b,,1,,)
   
   rsnew =
@@ -37 +37 @@
 
 Step 6: p,,1,, = r,,1,, +beta*p,,0,,  ('''No expansion''')
-  p,,1,,[0] = m(b,,0,,^2^+b,,1,,^2^)(b,,0,,m - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,)) + b,,0,,[(b,,0,,m - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,))^2^ + (b,,1,,m - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,))^2^]\\
-   ----------------------------------------------------------------------\\ 
+  p,,1,,[0] = 
+   m(b,,0,,^2^+b,,1,,^2^)[b,,0,,m - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,)] + b,,0,,{[b,,0,,m - (b,,0,,^2^+b,,1,,^2^)(a,,0,,b,,0,,+a,,1,,b,,1,,)]^2^ + [b,,1,,m - (b,,0,,^2^+b,,1,,^2^)(a,,1,,b,,0,,+a,,2,,b,,1,,)]^2^}\\
+   --------------------------------------------------------------------------------------------------------------------\\ 
    (b,,0,,^2^+b,,1,,^2^)m^2^