Dear all,
Is the an example about resetting the Dirichlet boundary condition after modifying a PETSc matrix?
For example,
Mat A;
A=a(Vh, Vh, tgv=-1);
rhs=a(0,Vh, tgv=-1);
a modification of A and rhs;
We then need to reset the locking conditions in A and rhs.
Best,
Yongxing.