First time here? Checkout the FAQ!
x
0 votes
by (160 points)
Hi,

I have a script that works for monodomain simulation on a 3D block with a planar stimulus. For some reason when I try to run bidomain simulations instead of monodomain the following error occurs during prepacing:

L5 : elliptic PDE (PETSc) solver diverged. Reason: indefinite preconditioning.

What is the problem here and do you have any suggestions on what I could try to fix it?

Thank you a lot in advance!

Nils
by (17.4k points)
Hi Nils,

Please provide a minimal working example that allows reproducing the error.
by (160 points)
I found out what caused the error, while trying to provide a minimal working example. The block contains a conduction block area that has a low conductivity.  
The selected conductivity value (g=0.00000001) of this area caused the error. Changing the conductivity to a higher value that is still sufficiently low (g=0.0001) to create a conduction block, fixed the problem.

2 Answers

+1 vote
by (17.4k points)

Copying the comment as answer:

I found out what caused the error, while trying to provide a minimal working example. The block contains a conduction block area that has a low conductivity.  
The selected conductivity value (g=0.00000001) of this area caused the error. Changing the conductivity to a higher value that is still sufficiently low (g=0.0001) to create a conduction block, fixed the problem.

+1 vote
by (7.8k points)
Hi!

Extreme conductivity heterogeneities will cause problems when solving the elliptic PDE. Even if the solver doesnt crash, convergence will be reduced. Conduction blocks can be more easily created by turning of Na in the ionic model parameters (GNa*0).

Cheers, Aurel
Welcome to openCARP Q&A. Ask questions and receive answers from other members of the community. For best support, please use appropriate TAGS!
architecture, carputils, documentation, experiments, installation-containers-packages, limpet, slimfem, website, governance
...