Hello:
So I'm a bit of a newbie, but I am trying to adapt example 13 to
isotropic/nematic liquid crystal phase transition model that I have been
using Comsol to solve for a while now. I think I've made some decent
progress, but I have a bit of a problem when I run the code I get:
*** Solving time step 0, time = 1e-07 ***
ERROR: LASPACK Error:
in ILUFactor for Mat:
Factorization produces zero pivot elements.
[0] src/numerics/laspack_linear_solver.C, line 276, compiled Apr 21 2008
at 14:01:19
Aborted
Basically I added my expressions for the residual and Jacobian and
removed all of the boundary condition code because I need to impose
Neuman BCs (which from what I understand is implicitely taken care of by
the FEM algorithm). I initialized t=0 with an acceptable solution to
the equations, so I just expected some smooth sailing...
I'll keep debugging of course, but I just wanted to through the error
out there in case it has to do with my understanding of how to implement
Neuman conditions.
I am also wondering about DiffSystem, should I just scrap my work with
ex13 and start with DiffSystem?
Nasser
-------------------------------------------------------------------------
This SF.net email is sponsored by the 2008 JavaOne(SM) Conference
Don't miss this year's exciting event. There's still time to save $100.
Use priority code J8TL2D2.
http://ad.doubleclick.net/clk;198757673;13503038;p?http://java.sun.com/javaone
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
