[petsc-users] about MATSOR

[Sun, Hui]

I want to do 3 steps of gauss seidel from a Mat A to another Mat B. Is there a way to do this?

I mean, what I can think of is to get the column vectors of B by

ierr = MatGetColumnVector(B,v,col);CHKERRQ(ierr);


and apply Gauss seidel from A to v:

ierr = MatSOR(A, v, 1, (SOR_FORWARD_SWEEP|SOR_ZERO_INITIAL_GUESS), 0, 3, 0, x);CHKERRQ(ierr);


But then I have to form another matrix C whose columns are composed by v, and I'm not sure how to do that.


Best,

Hui

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150225/f5a2b9f5/attachment.html>