Example: Then and column Optional D UT L “pointers” index 1 2 3 4 5 6 7 8 9 IJA 6 6 6 7 8 1 2 AM a11 a22 a33 a44 a13 a24 a31 a42 Compute Ax or ATx Procedure MULT( ) N, A, x, y do i = 1:N y(i)=A(i)x(i) enddo Lshift=0 if L is not stored and otherwise IJA(N+1)-IJA(1)
Example: Then D and column “pointers” UT Optional L index 1 2 3 4 5 6 7 8 9 IJA 6 6 6 7 8 1 2 AM a11 a22 a33 a44 a13 a24 a31 a42 Compute Ax or ATx Procedure MULT( N, A, x, y ) do i = 1:N y(i)=A(i)x(i) enddo Lshift=0 if L is not stored and IJA(N+1)-IJA(1) otherwise
2013sw.docx