Difference between a Newton solver and a linear solver for a linear problem.

Linear Algebra
1

Hi,

My question is: is there a difference between a linear problem solved with LinearProblem() and the same linear problem solved with NonlinearProblem()?
If I pass

petsc_options = {"ksp_type": "preonly",
                 "pc_type": "lu",
                 "pc_factor_mat_solver_type": "mumps",
}

then both methods are equivalent. Am I correct?

Thanks.

2

Solving it with nonlinear-problem will do an extra evaluation of the residual (ie the initial residual).
If you choose a Newton solver with no linesearch it should converge in one iteration,

1 Like