public inbox for gsl-discuss@sourceware.org
 help / color / mirror / Atom feed
From: Rhys Ulerich <rhys.ulerich@gmail.com>
To: gsl-discuss@sourceware.org
Subject: Implementing hybrid explicit/implicit ODE stepper
Date: Wed, 03 Jun 2009 14:32:00 -0000	[thread overview]
Message-ID: <4a00655d0906030731l8632674o9abb33dca522e689@mail.gmail.com> (raw)

Hi all,

I'm looking into implementing a hybrid explicit, implicit ODE
integrator within the GSL ODE framework.
This one splits the right hand side of u_t = f(u) into an implicit
linear operator L(u) and an explicit nonlinear operator N(u).  The
scheme requires L(u) be available in a matrix or matrix-like form.

What's the "right way" to shoehorn this into the gsl_odeiv_system
struct (http://www.gnu.org/software/gsl/manual/html_node/Defining-the-ODE-System.html)
?
L(u) as gsl_odeiv_system.function with N(u) passed within
gsl_odeiv_system.params ?  The reverse?  Neither?

Thanks,
Rhys

P.S. It's Spalart, Moser, & Rogers 1991 appendix A for anyone interested.

             reply	other threads:[~2009-06-03 14:32 UTC|newest]

Thread overview: 5+ messages / expand[flat|nested]  mbox.gz  Atom feed  top
2009-06-03 14:32 Rhys Ulerich [this message]
2009-06-04 13:38 ` Brian Gough
2009-06-04 17:20   ` Tuomo Keskitalo
2009-06-04 17:40     ` Rhys Ulerich
2009-06-04 18:52       ` Tuomo Keskitalo

Reply instructions:

You may reply publicly to this message via plain-text email
using any one of the following methods:

* Save the following mbox file, import it into your mail client,
  and reply-to-all from there: mbox

  Avoid top-posting and favor interleaved quoting:
  https://en.wikipedia.org/wiki/Posting_style#Interleaved_style

* Reply using the --to, --cc, and --in-reply-to
  switches of git-send-email(1):

  git send-email \
    --in-reply-to=4a00655d0906030731l8632674o9abb33dca522e689@mail.gmail.com \
    --to=rhys.ulerich@gmail.com \
    --cc=gsl-discuss@sourceware.org \
    /path/to/YOUR_REPLY

  https://kernel.org/pub/software/scm/git/docs/git-send-email.html

* If your mail client supports setting the In-Reply-To header
  via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line before the message body. 
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).