| Version 6 (modified by , 10 years ago) ( diff ) |
|---|
A = x0, x1 b = x3[2]
x1, x2
Step 1: r = b - Ax
Step 2: alpha = (r[i]*r[i]) / (p[i]*A[i][j]*p[j])
alpha = (X3[0]2+X3[1]2) / ((X3[1]2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]2)*X0)
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)
Step 5: beta = (rk[i]*rk[i]) / (r[i]*r[i])
beta = ((X3[1]6)*(X12)+2*((X3[1]5)*X3[0]*X0*X1)-2*((X3[1]5)*X3[0]*X1*X2)+(X3[1]4)*(X3[0]2)*(X02)-2*((X3[1]4)*(X3[0]2)*X0*X2)-1*((X3[1]4)*(X3[0]2)*(X12))+(X3[1]4)*(X3[0]2)*(X22)+(X3[1]2)*(X3[0]4)*(X02)-2*((X3[1]2)*(X3[0]4)*X0*X2)-1*((X3[1]2)*(X3[0]4)*(X12))+(X3[1]2)*(X3[0]4)*(X22)-2*(X3[1]*(X3[0]5)*X0*X1)+2*(X3[1]*(X3[0]5)*X1*X2)+(X3[0]6)*(X12))/((X3[1]6)*(X22)+4*((X3[1]5)*X3[0]*X1*X2)+2*((X3[1]4)*(X3[0]2)*X0*X2)+4*((X3[1]4)*(X3[0]2)*(X12))+(X3[1]4)*(X3[0]2)*(X22)+4*((X3[1]3)*(X3[0]3)*X0*X1)+4*((X3[1]3)*(X3[0]3)*X1*X2)+(X3[1]2)*(X3[0]4)*(X02)+2*((X3[1]2)*(X3[0]4)*X0*X2)+4*((X3[1]2)*(X3[0]4)*(X12))+4*(X3[1]*(X3[0]5)*X0*X1)+(X3[0]6)*(X02))
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]*(X12))+(X3[1]6)*X3[0]*(X22)-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)*(X12))+2*((X3[1]4)*(X3[0]3)*(X22))-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)*(X12)+(X3[1]2)*(X3[0]5)*(X22)-1*(X3[1]*(X3[0]6)*X0*X1)+2*(X3[1]*(X3[0]6)*X1*X2)+(X3[0]7)*(X12))/((X3[1]6)*(X22)+4*((X3[1]5)*X3[0]*X1*X2)+2*((X3[1]4)*(X3[0]2)*X0*X2)+4*((X3[1]4)*(X3[0]2)*(X12))+(X3[1]4)*(X3[0]2)*(X22)+4*((X3[1]3)*(X3[0]3)*X0*X1)+4*((X3[1]3)*(X3[0]3)*X1*X2)+(X3[1]2)*(X3[0]4)*(X02)+2*((X3[1]2)*(X3[0]4)*X0*X2)+4*((X3[1]2)*(X3[0]4)*(X12))+4*(X3[1]*(X3[0]5)*X0*X1)+(X3[0]6)*(X02))
p[1] = ((X3[1]7)*(X12)+2*((X3[1]6)*X3[0]*X0*X1)-1*((X3[1]6)*X3[0]*X1*X2)+(X3[1]5)*(X3[0]2)*(X02)-1*((X3[1]5)*(X3[0]2)*X0*X2)+(X3[1]5)*(X3[0]2)*(X12)+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)*(X02))-2*((X3[1]3)*(X3[0]4)*X0*X2)-1*((X3[1]3)*(X3[0]4)*(X12))-1*((X3[1]2)*(X3[0]5)*X1*X2)+X3[1]*(X3[0]6)*(X02)-1*(X3[1]*(X3[0]6)*X0*X2)-1*(X3[1]*(X3[0]6)*(X12))-1*((X3[0]7)*X0*X1))/((X3[1]6)*(X22)+4*((X3[1]5)*X3[0]*X1*X2)+2*((X3[1]4)*(X3[0]2)*X0*X2)+4*((X3[1]4)*(X3[0]2)*(X12))+(X3[1]4)*(X3[0]2)*(X22)+4*((X3[1]3)*(X3[0]3)*X0*X1)+4*((X3[1]3)*(X3[0]3)*X1*X2)+(X3[1]2)*(X3[0]4)*(X02)+2*((X3[1]2)*(X3[0]4)*X0*X2)+4*((X3[1]2)*(X3[0]4)*(X12))+4*(X3[1]*(X3[0]5)*X0*X1)+(X3[0]6)*(X02))
Step 7: alpha = (r[i]*r[i]) / (p[i]*A[i][j]*p[j])
alpha =
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] =
END
Attachments (4)
-
2x2caseWithoutLoop.cvl
(3.1 KB
) - added by 10 years ago.
CG 2x2 case without loop
-
2x2caseSteps.pdf
(212.0 KB
) - added by 10 years ago.
cg2x2_ intermediateResult
-
2x2caseSimplified.pdf
(113.8 KB
) - added by 10 years ago.
Simplified steps of 2x2 case
- 2x2caseSimplifiedUpdate.pdf (106.8 KB ) - added by 10 years ago.
Download all attachments as: .zip
