<div dir="ltr"><div class="gmail_default" style="font-size:small">​Hi,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm trying to implement modified newton method to solve the nonlinear finite element using petsc.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">As well known,  the advantage of modified newton is the Jacobian matrix is always same during the iteration, which means once the J is factorized at the first iteration, we can store the factors and avoid the factorization for next iteration if we use direct solver, e.g. super_lu. Therefore, the option FACTORED in SUPER_LU is quite useful.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">However, it looks like the option FACTORED is not available in SUPER_LU_DIST in petsc. I tried, and it shows 'unknown option'.<br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Is there alternative way to use the same idea of FACTORED in petsc for super_lu?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Also, I'm wondering whether iterative solver in PETSC is also able to apply the same strategy.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">In other words, in the problem where Jacobian is constant, only residue and solution vectors need to be updated, is there any way to take advantage of such same Jacobian pattern to expedite the computation using iterative solver? </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thank you</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">BTW, though using modified newton will increase the iteration number, however, in the case which is much more expensive to factorize the jacobian, more iterations will probably be worthwhile.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Bests,</div><div class="gmail_default" style="font-size:small">Jinlei</div></div>