dolfin team mailing list archive

[Alexander Jarosch <alexanj@xxxxx>]

I got a bit lost in the current updates and now I just have a short 
question. I set up my solver like this:

 // Discretize equation
 Matrix A;
 Vector x, b;
 FEM::assemble(*a, *L, A, b, mesh, bc);

 // Solve the linear system
 KrylovSolver solver(Preconditioner::ilu);
 KSP ksp;
 ksp = solver.solver();

 PC pc;
 KSPGetPC(ksp, &pc);
 PCFactorSetShiftNonzero(pc, PETSC_DECIDE);
 //PCFactorSetShiftPd(pc, PETSC_TRUE);

 solver.set("Krylov relative tolerance", 1e-14);

 solver.solve(A, x, b);

but Petsc still complains about zero Pivot. should the line 
"PCFactorSetShiftNonzero(pc, PETSC_DECIDE);" not take care of that.

Alex

--
Alexander H. Jarosch

Jarðvísindastofnun Háskólans
Institute of Earth Sciences, University of Iceland
Náttúrufræðahús, Askja
Building of Natural Sciences, Askja
Sturlugata 7
IS - 101 Reykjavík
Iceland

Tel.: +354 525 4906
http://raunvis.hi.is/~jarosch/