[petsc-users] The exact Jacobian for ASPIN

Lulu Liu lulu.liu at kaust.edu.sa
Sat Nov 6 06:19:59 CDT 2021


Hi,
For ASPIN, the local problem F_{j}(x-T_{j})=0 (j=1,2,...N)  is solved. The
exact Jacobian of the preconditioned Jacobian is
J_exact=\sum_{j=1}^{N}J_{I}(x-T_{i})J(x-T_{j}).

It seems that PETSc uses the approximate Jacobian,
J_exact=\sum_{j=1}^{N}J_{j}(x-T_{i})J(x-\sum_{j=1}^{N}T_{j}).

I want to implement RASPEN, which requires the exact Jacobian of
ASPIN.  Is there
any easy way to compute J(x-T_{j}) (j=1,2,..N)? How can I get the global
vectors like  x-T_{j}  ? PETSc only provides  the vector x-T_{j} on the
subdomain now.

Thanks very much!

-- 
Best wishes,
Lulu Liu

-- 

This message and its contents, including attachments are intended solely 
for the original recipient. If you are not the intended recipient or have 
received this message in error, please notify me immediately and delete 
this message from your computer system. Any unauthorized use or 
distribution is prohibited. Please consider the environment before printing 
this email.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20211106/1aa84edf/attachment.html>


More information about the petsc-users mailing list