[petsc-dev] The question about FD Jacobian

[Lulu Liu]

The Jacobian of the nonlinear system (unknowns are p,s) has the form of J=[
Jpp,Jps;Jsp,Jss].
However, I want to compute each block of Jacobian at different points.
Jpp at x=(p1,s1); Jps at x(p2,s2);
Is it possible to approximate each block (evaluated at different points)
using FD? I have the residual expression for the nonlinear system:
f[j][i].p=....
f[j][i].S=....

Thanks

-- 

------------------------------
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-dev/attachments/20131206/634d07a4/attachment.html>