[petsc-dev] Question about the pseudo-time stepping

Lulu Liu lulu.liu at kaust.edu.sa
Tue Sep 9 07:07:05 CDT 2014


Hi,

I want to solve a stead-state problem G(u)=0 using the pseudo-transient
continuation method.

I saw the example petsc/src/ts/examples/tutorials/ex1,

  ierr = TSSetRHSJacobian(ts,J,J,*FormJacobian*,&user);CHKERRQ(ierr);

however, the analytical Jacobian in my problem is very complex and it is
not easy to give the function "*FormJacobian*". If I run the code like

 ./ex1 -ts_monitor_pseudo  -snes_mf

Is it correct? Does it involves the time derivative term when using
-snes_mf?

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/20140909/94d1c450/attachment.html>


More information about the petsc-dev mailing list