I know there was a lot of discussion that led up to adding the ShellMatrix object to libMesh.... and looking back through it I see that preconditioning was mentioned... but now that I'm looking at the code I don't see a straight forward way to add a ShellMatrix as a preconditioner.... but maybe I'm missing it.
What I mean is... there needs to be a call to PCSetType(pc,PCSHELL) (which there is) and a call to PCShellSetApply.... which there isn't... or am I just missing it? Currently I'm just using my own hacked up code that doesn't even use ShellMatrix... but it would be nice to use the interfaces (for easier transition to Trilinos). Also... on a related question.... what do you guys think about providing a Preconditioner base class in libMesh? Concrete implementations of this class would be HyprePreconditioner, MLPreconditioner, LinearSolvePreconditioner, etc... I guess they are similar to linear solvers... but not exactly the same.... Derek ------------------------------------------------------------------------------ This SF.net email is sponsored by: SourcForge Community SourceForge wants to tell your story. http://p.sf.net/sfu/sf-spreadtheword _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
