How to get the matrix after ordering
Barry Smith
bsmith at mcs.anl.gov
Thu Feb 26 08:25:04 CST 2009
On Feb 26, 2009, at 2:20 AM, Shengli Xu wrote:
>
>
> Another question: Matrix A is known, there are some zero entries in
> A. For example A(2,3)=0.0 and A(4,5)=0.0. I want to get a new matrix
> B, which has the same nonzero structure as A but does not have the
> entries where are zero in A. For example no entries of B(2,3) and
> B(4,5). How to get Matrix B?
There isn't anything like this. It is also not a desirable
operation since the matrix has to be copied into new memory so you
really are not getting a memory savings. Just try to do the best you
can with not having zero values and don't sweat having a few.
MatOption(A, MAT_IGNORE_ZERO_ENTRIES) will cause any call of
MatSetValues() with a zero value to be ignored and not introduce a
"zero" "nonzero location" into the matrix.
Barry
>
>
> I want to use MatGetSeqNonzeroStructure(), But I don't think clearly.
>
> thanks in advance.
>
>
> --
> Shengli Xu
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which
> their experiments lead.
> -- Norbert Wiener
>
>
>
> --
> Shengli Xu
More information about the petsc-users
mailing list