Results 1 to 2 of 2

Math Help - Implementation of psuedo-arclength continuation

  1. #1
    Newbie
    Joined
    Jul 2009
    Posts
    5

    Question Implementation of psuedo-arclength continuation

    Hi there,

    I have a quick question regarding pseudo-arclength continuation. As some background, I am a chemical engineer, not a mathematician, applied or otherwise, so my knowledge of numerical methods is limited to the "standard" stuff.

    I've been reading up as extensively as I can on pseudo-arclength continuation, but unfortunately, it's all from second-hand sources; I don't have access to Keller's original paper. Here's what I understand at this point:

    We want to solve a problem F(x,\lambda)=0. We assume that the solution is known at x^0 and \lambda^0. To avoid the singularity of the Jacobian, and therefore the breakdown of Newton's method, at turning points, x and \lambda both become parameterised by arclength ( s), and we end up with an augmented system of equations to solve:

    F(x,\lambda)=0
    \left(u-u^{0}\right)\mathrm{d}u^{0}/\mathrm{d}s+\left(\lambda-\lambda^{0}\right)\mathrm{d}\lambda^{0}/\mathrm{d}s-\Delta S=0

    While this seems simple enough, how does one obtain the derivatives w.r.t s? Not a single text seems to mention this. Ideas that spring to mind are forward differences using, say, cubic splines; however, that seems horrendously inefficient to me. There must be a better way!

    Thanks,
    NChE
    Follow Math Help Forum on Facebook and Google+

  2. #2
    Newbie
    Joined
    Jul 2009
    Posts
    5

    Re: Implementation of psuedo-arclength continuation

    I found the answer: a procedure is described in Kubicek's Algorithm 502 in ACM Trans. Math. Software. dl.acm.org/citation.cfm?id=355675
    Follow Math Help Forum on Facebook and Google+

Similar Math Help Forum Discussions

  1. Replies: 3
    Last Post: September 11th 2011, 08:17 AM
  2. LaTex implementation on MHF?
    Posted in the LaTeX Help Forum
    Replies: 0
    Last Post: July 7th 2011, 04:13 PM
  3. Discrete Fourier Transform implementation
    Posted in the Advanced Math Topics Forum
    Replies: 3
    Last Post: March 8th 2011, 05:05 PM
  4. Continuation
    Posted in the Geometry Forum
    Replies: 2
    Last Post: January 23rd 2011, 01:37 AM
  5. Probability Transition Matrix Implementation Algorithm
    Posted in the Advanced Statistics Forum
    Replies: 2
    Last Post: February 12th 2009, 04:41 AM

/mathhelpforum @mathhelpforum