multibolzmatvectorrectlast.fp

PDE simulator on GPU.

!!FP1.0
#Date:2003/12/31
#try to compute sparse Matrix - Vector Multiplication Just For TEXRECT
# Vector X  ----TEX0 
# Diagonal Element   ----TEX1 -- AiiTexID  // record the diagonal elements
# First NonZero Postion--TEX1 -- RTexID;   // record the first non zero & non diagonal element position J
# NonZero NonDiagonal----TEX2 -- AijTexID; // record the nondiagonal & nonzero elements
# For reference in X ----TEX2 -- CTexID;   // record the according position for Vector X which multiply with Aij
# Previous  ====TEX5 //////record the previous result

# fetch Ai[i] from TEX1
TEX R0, f[WPOS], TEX1, RECT;

# fetch X[i] from TEX0
TEX R1, f[WPOS], TEX0, RECT;

# fetch previous result
TEX R2, f[WPOS], TEX3, RECT;

# Ai[i]*X[i] + PreviousResult
MAD o[COLR], R0, R1, R2;

END