On Tue, Oct 19, 2010 at 8:41 AM, Roy Stogner <[email protected]> wrote: > > On Tue, 19 Oct 2010, yunfei zhu wrote: > >> I am trying to use arc-length method to track the post-buckling equilibrium >> path. I need to add another parameter and another additional equation to the >> system. >> Has anyone got experience on using arc-length together with libmesh and >> petsc. >> Is there any Class in libmesh that I can use directly? Thanks. > > There's a ContinuationSystem class to do exactly that. I think John > Peterson is the only expert with it; we'll see if he chimes in.
Yes, it's been a while since I've touched it, but the ContinuationSystem class can be used for this purpose. This particular implementation takes a different approach from the SCALAR. It solves two linear systems per Newton iteration using the decomposition method of Keller, rather than a single system with an additional dense row and column. This has the benefit of working with an existing matrix assembly routine for a non-continuation solve, but has the drawback of requiring two solves. H. B. Keller. Numerical solution of bifurcation and nonlinear eigenvalue problems. In Applications of Bifurcation Theory, P. H. Rabinowitz (ed.), pages 359–389. Academic Press, 1977. I have an electronic copy of the above but it's a bit too large to email, let me know if you would like a copy. -- John ------------------------------------------------------------------------------ Download new Adobe(R) Flash(R) Builder(TM) 4 The new Adobe(R) Flex(R) 4 and Flash(R) Builder(TM) 4 (formerly Flex(R) Builder(TM)) enable the development of rich applications that run across multiple browsers and platforms. Download your free trials today! http://p.sf.net/sfu/adobe-dev2dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
