Changes between Version 30 and Version 31 of PolynomialExpansion


Ignore:
Timestamp:
01/26/16 16:13:19 (10 years ago)
Author:
sili
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • PolynomialExpansion

    v30 v31  
    11
    2 A = {{a,,0,,, a,,1,,}, {a,,1,,, a,,2,,}},  b = {b,,0,,, b,,1,,}
     2A = {{a,,0,,, a,,1,,}, {a,,1,,, a,,2,,}},  b = {b,,0,,, b,,1,,}, x = {0, 0}
    33
    4 x[0] = 0, x[1] = 0
    5 
    6 k = 0
    74
    85Step 1: r,,0,, = b - Ax,,0,,
    … …  
    1310  p,,0,, = r,,0,,
    1411
    15 Step 2: alpha = <r, r> / <p, Ap>
     12Step 2: alpha,,0,, = <r,,0,,, r,,0,,> / <p,,0,,, Ap,,0,,>
    1613  <p, Ap> = a,,0,,*b,,0,,^2^ + 2*a,,1,,*b,,0,,*b,,1,, + a,,2,,*b,,1,,^2^
    1714
    18   alpha = (b,,0,,^2^+b,,1,,^2^) / ((b,,1,,^2^)*a,,2,,+2*(b,,1,,*b,,0,,*a,,1,,)+(b,,0,,^2^)*a,,0,,)
     15  alpha,,0,, = (b,,0,,^2^+b,,1,,^2^) / ((b,,1,,^2^)*a,,2,,+2*(b,,1,,*b,,0,,*a,,1,,)+(b,,0,,^2^)*a,,0,,)
    1916
    20 Step 3: r[i] = r - alpha * Ap
     17Step 3: r,,1,, = r,,0,, - alpha,,0,, * Ap,,0,,
    2118  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^ + 2*a,,1,,*b,,0,,*b,,1,, + a,,2,,*b,,1,,^2^)
    2219