I did some more testing, the trilinos LinearGMRESsolver works even better than the Pysparse LinearPCG and allows me to continue working, but I still don't know why the trilinos LinearPCG is failing for me in that one case.
I am willing to debug it if the developers are interested, otherwise I'll go on my merry way. Thanks, -mike waters On 3/8/16 2:09 PM, Michael Waters wrote: > Hi, I am trying to solve the Poisson equation like this: > > > > from fipy.solvers.trilinos import LinearPCGSolver as mysolver > from fipy.solvers.trilinos.preconditioners import JacobiPreconditioner > as myprecon > > rho_h = CellVariable(name = 'Charge Density', mesh=mesh, value = 0.0) > phi_h = CellVariable(name = 'Potential', mesh=mesh, value = 0.0) > epsilon = CellVariable(name = 'Dielectric Permitivity', mesh=mesh, > value = epsilon0) > > > phi_h.equation = (0.0 == DiffusionTerm(coeff = epsilon) + rho_h) > > > phi_h_res = phi_h.equation.sweep(var = phi_h, solver = mysolver > (iterations = phi_solver_iterations_per_step, precon = myprecon() ) ) > > What is strange is that if I use the PySparse equivalent > LinearPCGSolver and JacobiPreconditioner, the problem solves normally. > Also, in my code I solve a similar equation using the same Trilinos > solvers without any problems. > > Thoughts? > -mike waters _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
