On Thu, 14 Jan 2010 14:28:30 -0700, Derek Gaston <[email protected]> wrote: > BTW we do everything using JFNK (Jacobian Free Newton Krylov)... and > heavily use threads... and see some good gains (as you say you > should). And of course we're using libMesh for everything.
Cool. Dana told me a little about this MOOSE thing. I think moose are tasty. I'd be very interested to see how general linear methods (TSGL in PETSc-dev) perform on your stiff problems. These are A- and L-stable integrators with diagonally implicit computational structure and stage-order 1 through 5. My implementation is adaptive order, adaptive step although the controllers are probably weak (but it has a plugin architecture). > We're actually in the middle of putting together our talk for next > month's SIAM Parallel Processing conference where we're on the hook to > give a talk entitled: "Hybrid Parallelism for Preconditioned > Jacobian-Free Newton-Krylov Multiphysics Applications". After the > conference I can put the slides somewhere for interested parties. Neat. Jed ------------------------------------------------------------------------------ Throughout its 18-year history, RSA Conference consistently attracts the world's best and brightest in the field, creating opportunities for Conference attendees to learn about information security's most important issues through interactions with peers, luminaries and emerging and established companies. http://p.sf.net/sfu/rsaconf-dev2dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
