Hi Tony, Sorry for the delay in responding to this.
Have you checked which solver library is being used in each of the different cases? You can actually import the solver from the library directly so you can do "from fipy.solvers.pysaprse.linearLUSolver import LinearLUSolver" for example. FiPy may select the Trilinos solver suite if it thinks it is running in parallel. Definitely make sure you know which solver suite is being used in each case. One of the solvers may not be converging properly or not to enough accuracy. Cheers, Daniel On Fri, Mar 13, 2015 at 4:16 PM, <[email protected]> wrote: > # Evolve system > solver = fipy.solvers.LinearPCGSolver(tolerance=10**-20, \ > iterations=1000, precon=None) Import from fipy.solvers.pysparse here or trilinos and then you'll see, which one is working. It may be that PCG doesn't work without the default preconditioning in one of the cases. -- Daniel Wheeler _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
