[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Help with ODE solver
From: |
Marco Antoniotti |
Subject: |
Re: Help with ODE solver |
Date: |
Thu, 12 Sep 2002 10:08:48 -0400 |
> X-Authentication-Warning: bevo.che.wisc.edu: list set sender to
> address@hidden using -f
> From: "John W. Eaton" <address@hidden>
> Date: Thu, 12 Sep 2002 02:22:51 -0500
> Cc: Marco Antoniotti <address@hidden>, address@hidden
> Resent-From: address@hidden
> X-Mailing-List: <address@hidden>
> X-Loop: address@hidden
> Resent-Sender: address@hidden
>
> On 12-Sep-2002, Douglas Eck <address@hidden> wrote:
>
> | My experience is that it's better to break down the integration into
> | intervals. Also, if you later need to speed up your code some, it can
> | be much faster to write the xdot = f(x,t) as a compiled .oct file.
>
> Lsode may integrate past the final output point and interpolate back
> to find the state values at the final output point, so simply breaking
> up the integration into subintervals may not be enough to avoid
> numerical problems if parameters change in your RHS function at the
> final output point of each interval. Using the optional tcrit input
> for lsode may be needed anyway.
Thanks to all. I hear you.
Could anybody please explain into a little more detains the usage of
the TCRIT argument? The manual is not that illuminating and I am not
a numerical analyst (despite the fact that I could just bother a few
across the street :) )
Cheers
--
Marco Antoniotti ========================================================
NYU Courant Bioinformatics Group tel. +1 - 212 - 998 3488
715 Broadway 10th Floor fax +1 - 212 - 995 4122
New York, NY 10003, USA http://bioinformatics.cat.nyu.edu
"Hello New York! We'll do what we can!"
Bill Murray in `Ghostbusters'.
-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------
- Help with ODE solver, Marco Antoniotti, 2002/09/11
- Re: Help with ODE solver, Michael Hanke, 2002/09/12
- Re: Help with ODE solver, Douglas Eck, 2002/09/12
- Re: Help with ODE solver, John W. Eaton, 2002/09/12
- Re: Help with ODE solver,
Marco Antoniotti <=
- Message not available
- Re: Help with ODE solver, John W. Eaton, 2002/09/12
- Re: Help with ODE solver, Marco Antoniotti, 2002/09/12
- Re: Help with ODE solver, John W. Eaton, 2002/09/12
- Re: Help with ODE solver, Marco Antoniotti, 2002/09/12
- Re: Help with ODE solver, John W. Eaton, 2002/09/12
- Re: Help with ODE solver, Marco Antoniotti, 2002/09/13
- Re: Help with ODE solver, John W. Eaton, 2002/09/13
- Re: Help with ODE solver, Jon Davis, 2002/09/13
- Re: Help with ODE solver, Marco Antoniotti, 2002/09/13