Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

User-Supplied Jacobian in NDSolve

  • To: mathgroup at smc.vnet.net
  • Subject: [mg91152] User-Supplied Jacobian in NDSolve
  • From: Eric Brown <eric.c.brown at mac.com>
  • Date: Thu, 7 Aug 2008 04:42:30 -0400 (EDT)

Dear Mathgroup,

  I am using compiled functions (values computed with Compile or
Mathlink) and so I wrap my functions to ensure only numerical values
are supplied to these functions:

f[x_Real]:=fCompiled[x];

Therefore, if NDSolve could derive a symbolic Jacobian, then I assume
that NDSolve can not elucidate it for f[x_Real] since the functions are
numeric.

I guess I have two questions:

1) Does NDSolve use an analytic Jacobian an any case
and
2) Is it possible to supply a User-Defined Jacobian to NDSolve

I think that it can be done in general, given experience with GSL, LSODA,=
 etc.

Many thanks for your advice,

Eric








  • Prev by Date: Multiplying a vector over multiple vectors
  • Next by Date: Re: Find count of binary number pattern within concatenated number
  • Previous by thread: Re: Multiplying a vector over multiple vectors
  • Next by thread: NDSolve::mconly in algebraic differential equations