Changes between Version 3 and Version 4 of PolynomialExpansion

Timestamp:

01/26/16 00:10:36 (10 years ago)

Author:

sili

Comment:

--

PolynomialExpansion

@@ -13 +13 @@
 Step 3: r[i] = r - alpha * A[i][j] * p[j]
   r[0] = (-1*((X3[1]^3)*X1)-1*((X3[1]^2)*X3[0]*X0)+(X3[1]^2)*X3[0]*X2+X3[1]*(X3[0]^2)*X1) / ((X3[1]^2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0)
+
   r[1] = ((X3[1]^2)*X3[0]*X1+X3[1]*(X3[0]^2)*X0-1*(X3[1]*(X3[0]^2)*X2)-1*((X3[0]^3)*X1)) / ((X3[1]^2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0) 
 
 Step 4: x[i] = x + alpha*p[i]
   x[0] = ((X3[1]^2)*X3[0]+X3[0]^3) / ((X3[1]^2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0)
+
   x[1] = (X3[1]^3+X3[1]*(X3[0]^2)) / ((X3[1]^2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0)
 
@@ -24 +26 @@
 Step 6: p[i] = rk[i] +beta * p
   p[0] = (-1*((X3[1]^7)*X1*X2)-1*((X3[1]^6)*X3[0]*X0*X2)-1*((X3[1]^6)*X3[0]*(X1^2))+(X3[1]^6)*X3[0]*(X2^2)-1*((X3[1]^5)*(X3[0]^2)*X0*X1)-2*((X3[1]^4)*(X3[0]^3)*X0*X2)-1*((X3[1]^4)*(X3[0]^3)*(X1^2))+2*((X3[1]^4)*(X3[0]^3)*(X2^2))-2*((X3[1]^3)*(X3[0]^4)*X0*X1)+3*((X3[1]^3)*(X3[0]^4)*X1*X2)-1*((X3[1]^2)*(X3[0]^5)*X0*X2)+(X3[1]^2)*(X3[0]^5)*(X1^2)+(X3[1]^2)*(X3[0]^5)*(X2^2)-1*(X3[1]*(X3[0]^6)*X0*X1)+2*(X3[1]*(X3[0]^6)*X1*X2)+(X3[0]^7)*(X1^2))/((X3[1]^6)*(X2^2)+4*((X3[1]^5)*X3[0]*X1*X2)+2*((X3[1]^4)*(X3[0]^2)*X0*X2)+4*((X3[1]^4)*(X3[0]^2)*(X1^2))+(X3[1]^4)*(X3[0]^2)*(X2^2)+4*((X3[1]^3)*(X3[0]^3)*X0*X1)+4*((X3[1]^3)*(X3[0]^3)*X1*X2)+(X3[1]^2)*(X3[0]^4)*(X0^2)+2*((X3[1]^2)*(X3[0]^4)*X0*X2)+4*((X3[1]^2)*(X3[0]^4)*(X1^2))+4*(X3[1]*(X3[0]^5)*X0*X1)+(X3[0]^6)*(X0^2))
+  
   p[1] = ((X3[1]^7)*(X1^2)+2*((X3[1]^6)*X3[0]*X0*X1)-1*((X3[1]^6)*X3[0]*X1*X2)+(X3[1]^5)*(X3[0]^2)*(X0^2)-1*((X3[1]^5)*(X3[0]^2)*X0*X2)+(X3[1]^5)*(X3[0]^2)*(X1^2)+3*((X3[1]^4)*(X3[0]^3)*X0*X1)-2*((X3[1]^4)*(X3[0]^3)*X1*X2)+2*((X3[1]^3)*(X3[0]^4)*(X0^2))-2*((X3[1]^3)*(X3[0]^4)*X0*X2)-1*((X3[1]^3)*(X3[0]^4)*(X1^2))-1*((X3[1]^2)*(X3[0]^5)*X1*X2)+X3[1]*(X3[0]^6)*(X0^2)-1*(X3[1]*(X3[0]^6)*X0*X2)-1*(X3[1]*(X3[0]^6)*(X1^2))-1*((X3[0]^7)*X0*X1))/((X3[1]^6)*(X2^2)+4*((X3[1]^5)*X3[0]*X1*X2)+2*((X3[1]^4)*(X3[0]^2)*X0*X2)+4*((X3[1]^4)*(X3[0]^2)*(X1^2))+(X3[1]^4)*(X3[0]^2)*(X2^2)+4*((X3[1]^3)*(X3[0]^3)*X0*X1)+4*((X3[1]^3)*(X3[0]^3)*X1*X2)+(X3[1]^2)*(X3[0]^4)*(X0^2)+2*((X3[1]^2)*(X3[0]^4)*X0*X2)+4*((X3[1]^2)*(X3[0]^4)*(X1^2))+4*(X3[1]*(X3[0]^5)*X0*X1)+(X3[0]^6)*(X0^2))
 
@@ -31 +34 @@
 Step 8: r[i] = r - alpha * A[i][j] * p[j]
   r[0] = 0
+
   r[1] = 0
 
 Step 9: x[i] = x + alpha*p[i]
   x[0] =
+
   x[1] =