[petsc-users] reuse matrix aij structure

[Feng Xing]

Hello everyone,

I want to solve many times A*x=b, where the non-zero structure of A (aij) is always same, but its values change. I would like to ask if it is possible to call assembly_begin/end only one time, then change the values of A each time. 

I find it is possible for a vector(Vec) through VecGetArray. It returns a pointer. Is there similar way for matrix? Thank you very much.

Regards,

Feng Xing
Postdoc at Inria