Howdy, Just curious - I’m talking to the PETSC folks, and shared this error message *** ------------------------------------------------------------------------- *** DOLFIN encountered an error. If you are not able to resolve this issue *** using the information listed below, you can ask for hel…

Do you still get the error when solving a much smaller problem? E.g: a 8^3 cube?

Hi Nate, For 32^2 and 64^2 everything is going great.

How many degrees of freedom arise in the 128^3 problem (E.g. V.dim()), and do you think using a direct solver such as mumps is appropriate with the hardware you have available? If not I’d recommend an iterative solver. It’s likely you’re running out of memory.

Hey Nate, The DOF is 4293378. What would you suggest as a solver (I really have no clue which is the correct solver to call)? Here is a piece of my code: state=inner(grad(u_r),grad(v_r))dx + inner(grad(u_i),grad(v_i))dx +kou_iv_rds -kou_rv_ids -v_r*(kosquared*(1+qw)u_r +kosquaredqwcos_f)dx -v_i(…

Looks sesquilinear. Is this the discretisation of the complex valued modified helmholtz equation? Preconditioning may be challenging. My usual naive route is to try gmres with hypre.

Hi Nate Yes it is! To do it do I how would I pass it correctly to solve(state == 0, u, J=Jac,solver_parameters={‘newton_solver’: {‘linear_solver’: ‘mumps’}})?

Thanks a lot Nate! Because my BC are in the weak form do I just ignore implementing/doing things with boundary conditions?

yes, or just an empty list. List of KSP types List of PC types

Dolfin error message converted to PETSC form?

nate January 8, 2020, 5:28pm 9

I recommend taking control and setting the options directly rather than using the black-box method solve(). See here for example.

Topic		Replies	Views
PETSc Krylov Solver RuntimeError with fine mesh I/O	21	5706	October 28, 2022
Error: Unable to solve linear system using PETSc Krylov solver Errors	1	560	August 10, 2022
PETSc error when running code on a fine mesh Linear Algebra	0	488	October 7, 2022
Solve linear system in Mixed finite element method General	6	647	March 1, 2023
Error with PETSc Linear Algebra	4	379	October 4, 2023

Dolfin error message converted to PETSC form?

Related topics