New question #138742 on DOLFIN: https://answers.launchpad.net/dolfin/+question/138742
hi, where can I find what VariationalProblem(...nonlinear=True...) is really doing, I can follow until the cpp.py where the classes have been defined, but cannot find the .cpp data where the functions are written, my aim to find out is why, if I use my own Newton iteration procedure, like ... problem = VariationalProblem(G, F, bcs=bc, ...) dv = problem.solve() differ = dv.vector().array() eps = numpy.linalg.norm(differ, ord=2) print 'K=%g : after iteration=%d , the L2-Norm of perturbation: %g ' % (mat['K'],iter, eps) v.vector()[:] += dv.vector() ... vs. the automatized one ... problem = VariationalProblem(G, F, bcs=bc, nonlinear=True, ...) problem.solve(v) ... they do not have the same convergence tendency, neither the same solution (if both would converge to one)? -- You received this question notification because you are a member of DOLFIN Team, which is an answer contact for DOLFIN. _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
