On Mon, Mar 23, 2009 at 2:08 PM, Kirk, Benjamin (JSC-EG311) <[email protected]> wrote: > John - since you have at least used SLEPc, would you please give this patch > a look? I think it is great that they did this for us!!
The patch is now checked in. At the moment I could only easily test against Petsc/Slepc 2.3.1, but whenever I get around to installing Petsc 3.0.0 I will take a look at this again. I had to make a few changes to allow backwards-compatibility since EPSKRYLOVSCHUR was not yet available in Slepc 2.3.1. Since I don't know when exactly it showed up, you now only get the Krylov-Schur solver when Petsc/Slepc 3.0.0 is available, otherwise you get ARNOLDI, which was the old default behavior. Perhaps the slepc developers could suggest a different or better default for versions of slepc prior to the introduction of Krylov-Schur? Thanks, -- John > ------ Forwarded Message > From: "Jose E. Roman" <[email protected]> > Date: Mon, 23 Mar 2009 13:54:07 -0500 > To: "Kirk, Benjamin (JSC-EG311)" <[email protected]> > Cc: Francisco Alvaro <[email protected]> > Subject: Upgrade for SLEPc interface in libMesh > > Dear Benjamin, > > I am the lead developer of SLEPc. As you may know, we have recently > released SLEPc 3.0.0, which has some incompatibilities with previous > versions. > > We have done some changes in order to enable libMesh to use SLEPc > 3.0.0. You have done a good job integrating SLEPc in libMesh, so it > was very easy. I am attaching the changes as diff's with respect to > your current svn repository. > > Just a few comments: > - We have set Krylov-Schur as the default eigensolver. This is far > better than the older Arnoldi. > - We have removed the code that set the orthogonalization to MGS, > because the default (iterative CGS) is much better. > - We have added a dummy function solve_adjoint() in order to be able > to compile. > > Please forward this message to the rest of libMesh developers. Let us > know if further interaction is required regarding this. > > Best regards, > Jose E. Roman > > > ------ End of Forwarded Message ------------------------------------------------------------------------------ Apps built with the Adobe(R) Flex(R) framework and Flex Builder(TM) are powering Web 2.0 with engaging, cross-platform capabilities. Quickly and easily build your RIAs with Flex Builder, the Eclipse(TM)based development software that enables intelligent coding and step-through debugging. Download the free 60 day trial. http://p.sf.net/sfu/www-adobe-com _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
