[petsc-users] about finite difference jacobian approximations

[陈乐平（Leping Chen）]

hello, petsc teams,

If I want to use finite difference jacobian approximations,how should I define 

the function FormJacobian()? I don't get analytic jacobians.

thanks,
leping

2010-05-02 



  



发件人： Jed Brown 
发送时间： 2010-05-02  20:10:20 
收件人： chenleping; PETSc users list 
抄送： 
主题： Re: [petsc-users] about finite difference jacobian approximations 
 
On Sat, 1 May 2010 15:35:25 +0800, "=?gb2312?B?s8LA1sa9o6hMZXBpbmcgQ2hlbqOp?=" <chenleping at yahoo.cn> wrote:
> hello, petsc teams,
> 
> When I use finite difference jacobian approximations, I only find the
> computation methods and results of jacobian matrix cannot influence
> iterative process of snes.
I don't know what you're asking.  The accuracy of the Jacobian certainly
can affect the convergence rate; for simple or well-conditioned
problems, finite difference Jacobians are usually reliable.  Analytic
Jacobians, or Jacobians from automatic differentiation (AD) are usually
better when available and not too expensive to compute.  Another common
thing is to assemble a "cheap" preconditioning matrix (dropping terms,
or with a lower order discretization) for use with a matrix-free
(usually finite difference or AD) operator.
Jed
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100501/fcc599d2/attachment.htm>