Yes, that works ! Linear firedrake_0_ solve converged due to CONVERGED_ITS iterations 1 Nonlinear firedrake_0_ solve converged due to CONVERGED_ITS iterations 1 David also suggested to use a direct solver but I did not use the option "mat_type": "aij". I'll read the doc about the options you suggested so I (hopefully) don't need to ask next time. Many many thanks, also for answering so fast. I'll do the test of convergence again to make sure it converges to the exact solution, but I expect it still does. Thank you very much, Floriane ________________________________ De : firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk> de la part de Lawrence Mitchell <wence@gmx.li> Envoyé : lundi 5 novembre 2018 09:50 À : firedrake@imperial.ac.uk Objet : Re: [firedrake] local discontinuity Ah, sorry!
On 5 Nov 2018, at 09:30, Floriane Gidel [RPG] <mmfg@leeds.ac.uk> wrote:
Dear Lawrence,
No, I get the following error:
Traceback (most recent call last): File "main.py", line 103, in <module> solve(a==0,u,bcs=[bc_left, bc_right, bc_disk_only,bc_plane_only], solver_parameter=parameters) File "/Users/mmfg/firedrake/src/firedrake/firedrake/solving.py", line 124, in solve _solve_varproblem(*args, **kwargs) File "/Users/mmfg/firedrake/src/firedrake/firedrake/solving.py", line 172, in _solve_varproblem solver.solve() File "/Users/mmfg/firedrake/src/firedrake/firedrake/variational_solver.py", line 220, in solve self.snes.solve(None, v) File "PETSc/SNES.pyx", line 537, in petsc4py.PETSc.SNES.solve (src/petsc4py.PETSc.c:172544) petsc4py.PETSc.Error: error code 56 [0] SNESSolve() line 3967 in /private/var/folders/7h/wbj8xp7n3g5cfbr32ctcmwzcy3jf53/T/pip-yGXqPh-build/src/snes/interface/snes.c [0] SNESSolve_KSPONLY() line 38 in /private/var/folders/7h/wbj8xp7n3g5cfbr32ctcmwzcy3jf53/T/pip-yGXqPh-build/src/snes/impls/ksponly/ksponly.c [0] KSPSolve() line 598 in /private/var/folders/7h/wbj8xp7n3g5cfbr32ctcmwzcy3jf53/T/pip-yGXqPh-build/src/ksp/ksp/interface/itfunc.c [0] KSPSetUp() line 378 in /private/var/folders/7h/wbj8xp7n3g5cfbr32ctcmwzcy3jf53/T/pip-yGXqPh-build/src/ksp/ksp/interface/itfunc.c [0] PCSetUp() line 924 in /private/var/folders/7h/wbj8xp7n3g5cfbr32ctcmwzcy3jf53/T/pip-yGXqPh-build/src/ksp/pc/interface/precon.c [0] PCSetUp_LU() line 84 in /private/var/folders/7h/wbj8xp7n3g5cfbr32ctcmwzcy3jf53/T/pip-yGXqPh-build/src/ksp/pc/impls/factor/lu/lu.c [0] MatGetOrdering() line 250 in /private/var/folders/7h/wbj8xp7n3g5cfbr32ctcmwzcy3jf53/T/pip-yGXqPh-build/src/mat/order/sorder.c [0] MatGetOrdering_ND() line 17 in /private/var/folders/7h/wbj8xp7n3g5cfbr32ctcmwzcy3jf53/T/pip-yGXqPh-build/src/mat/order/spnd.c [0] No support for this operation for this object type [0] Cannot get rows for matrix type nest
I forgot you had a mixed space. Please additionally add the option "mat_type": "aij" to the solver parameters, and try again. Thanks, Lawrence _______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake