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