[petsc-users] Linesearch algorithm

[Jed Brown]

On Tue, May 15, 2012 at 12:26 AM, behzad baghapour <
behzad.baghapour at gmail.com> wrote:

> Thanks. Actually the jacobian which I obtain has a part of transient as J
> = M/dt + dR/dQ. Then I want to eliminate the time effect before calculating
> R^T*J*dQ as the descent direction. So, which way is correct in my case?


This sounds like a good way to stagnate. Many transient problems _require_
the "steady-state" residual to increase (climbing a hill) before it can
decrease.

Also, I can't tell from your description if you want the search direction
to come from the transient or steady problem.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120515/cd451057/attachment.htm>