Hi, The demo in nls/nonlinearpoisson/python (for me) produces a seg fault during the application of the DirichletBC. The problem is related to the specification of the Vector x and the Function u.
the code has x = Vector() u = Function(element, mesh, x) (then later) nonlinear_solver.solve(nonlinear_problem, x) But this doesn't work as expected. What does work however is nonlinear_solver.solve(nonlinear_problem, u.vector()) Now in the cpp demo, we assign x as GenericVector& x = u.vector(); which is why it works, so my question is what is the equivalent way to initialize the function u and then declare x as the vector of u? Jake **************** Jake Ostien 1524 - Solid Mechanics Sandia National Laboratories 505.844.2230 **************** _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
