Re: Error in fcmix_serial for Sundials-2.4.0

Date: Mon, 12 Oct 2009 01:11:42 +0200
Message-ID: <3f94f5490910111611v7ccbcac8had910ba5ed37ceea@...>
Subject: Re: [sundials-users] Error in fcmix_serial for Sundials-2.4.0
From: Roberto Lemoine <robertolemoine@...>
To: sundials-users@...
Content-Type: multipart/alternative; boundary=00504502ca5aac98410475b0ef84

--00504502ca5aac98410475b0ef84
Content-Type: text/plain; charset=ISO-8859-1

Hello everyone:

I would appreciate your help. My problem is, that everytime that I start a
simulation using the interface SUNDIALS_TB for Matlab, with the IDA Solver,
using exactly the same parameters, initial estimates, etc., I get very
different situations (sometimes I get Not-A-Numbers, sometimes it converges,
sometimes it uses step sizes ridiculously small, etc), even with consistent
initial conditions.

I suspect that some garbage might be fed in the solver because an erroneous
utilization / combination of IDAInit, IDAReinit and IDAFree. However, I am
not sure. Can somebody help with this issue? Thanks a lot!

Roberto Lemoine

--00504502ca5aac98410475b0ef84
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable

Hello everyone:<br><br>I would appreciate your help. My problem is, that ev=
erytime that I start a simulation using the interface SUNDIALS_TB for Matla=
b, with the IDA Solver, using exactly the same parameters, initial estimate=
s, etc., I get very different situations (sometimes I get Not-A-Numbers, so=
metimes it converges, sometimes it uses step sizes ridiculously small, etc)=
, even with consistent initial conditions.<br>
<br>I suspect that some garbage might be fed in the solver because an erron=
eous utilization / combination of IDAInit, IDAReinit and IDAFree. However, =
I am not sure. Can somebody help with this issue? Thanks a lot!<br><br>
Roberto Lemoine<br><br>

--00504502ca5aac98410475b0ef84--

[Radu Serban-3]

 > I would appreciate your help. My problem is, that everytime that I
 > start a
 > simulation using the interface SUNDIALS_TB for Matlab, with the IDA
 > Solver,
 > using exactly the same parameters, initial estimates, etc., I get very
 > different situations (sometimes I get Not-A-Numbers, sometimes it
 > converges,
 > sometimes it uses step sizes ridiculously small, etc), even with
 > consistent
 > initial conditions.
 >
 > I suspect that some garbage might be fed in the solver because an
 > erroneous
 > utilization / combination of IDAInit, IDAReinit and IDAFree. However,
 > I am
 > not sure. Can somebody help with this issue? Thanks a lot!
 >
 > Roberto Lemoine

Roberto,

Maybe you could post the relevant part of your code showing how you use
these functions?

--Radu

----------------------------------------------------------------------------------------------------------------
Center for Applied Scientific Computing
Lawrence Livermore National Laboratory
P.O. Box 808, L-561
Livermore, CA  94551

[phone]  (Office)            cswoodward@...
[phone]  (Fax)               http://people.llnl.gov/woodward6 

[Alan Hindmarsh]

Hello Roberto-

Have you tried solving the problem with IDA directly, not via Matlab?
That way we would know if the problem is in IDA or in SundialsTB.

-Alan H



On Mon, 12 Oct 2009, by way of Carol Woodward <cswoodward@...> wrote:

> Date: Mon, 12 Oct 2009 01:11:42 +0200
> Message-ID: <3f94f5490910111611v7ccbcac8had910ba5ed37ceea@...>
> Subject: Re: [sundials-users] Error in fcmix_serial for Sundials-2.4.0
> From: Roberto Lemoine <robertolemoine@...>
> To: sundials-users@...
> Content-Type: multipart/alternative; boundary=00504502ca5aac98410475b0ef84
>
> --00504502ca5aac98410475b0ef84
> Content-Type: text/plain; charset=ISO-8859-1
>
> Hello everyone:
>
> I would appreciate your help. My problem is, that everytime that I start a
> simulation using the interface SUNDIALS_TB for Matlab, with the IDA Solver,
> using exactly the same parameters, initial estimates, etc., I get very
> different situations (sometimes I get Not-A-Numbers, sometimes it converges,
> sometimes it uses step sizes ridiculously small, etc), even with consistent
> initial conditions.
>
> I suspect that some garbage might be fed in the solver because an erroneous
> utilization / combination of IDAInit, IDAReinit and IDAFree. However, I am
> not sure. Can somebody help with this issue? Thanks a lot!
>
> Roberto Lemoine
>
> --00504502ca5aac98410475b0ef84
> Content-Type: text/html; charset=ISO-8859-1
> Content-Transfer-Encoding: quoted-printable

k

[Roberto Lemoine]

Hallo Alan:

I fortunately found out that it was an installation error in SundialsTB.
Everything is working fine now. Thanks for the advice!


Roberto Lemoine

--------------------------------------------------
From: "Alan Hindmarsh" <alanh@...>
Sent: Wednesday, December 02, 2009 11:04 PM
To: <sundials-users@...>
Cc: <alanh@...>
Subject: Re: [sundials-users] Error in fcmix_serial for Sundials-2.4.0

 >
 > Hello Roberto-
 >
 > Have you tried solving the problem with IDA directly, not via Matlab?
 > That way we would know if the problem is in IDA or in SundialsTB.
 >
 > -Alan H
 >
 >
 >
 > On Mon, 12 Oct 2009, by way of Carol Woodward <cswoodward@...> wrote:
 >
 >> Date: Mon, 12 Oct 2009 01:11:42 +0200
 >> Message-ID: <3f94f5490910111611v7ccbcac8had910ba5ed37ceea@...>
 >> Subject: Re: [sundials-users] Error in fcmix_serial for Sundials-2.4.0
 >> From: Roberto Lemoine <robertolemoine@...>
 >> To: sundials-users@...
 >> Content-Type: multipart/alternative;
 >> boundary=00504502ca5aac98410475b0ef84
 >>
 >> --00504502ca5aac98410475b0ef84
 >> Content-Type: text/plain; charset=ISO-8859-1
 >>
 >> Hello everyone:
 >>
 >> I would appreciate your help. My problem is, that everytime that I start
 >> a
 >> simulation using the interface SUNDIALS_TB for Matlab, with the IDA
 >> Solver,
 >> using exactly the same parameters, initial estimates, etc., I get very
 >> different situations (sometimes I get Not-A-Numbers, sometimes it
 >> converges,
 >> sometimes it uses step sizes ridiculously small, etc), even with
 >> consistent
 >> initial conditions.
 >>
 >> I suspect that some garbage might be fed in the solver because an
 >> erroneous
 >> utilization / combination of IDAInit, IDAReinit and IDAFree. However, I
 >> am
 >> not sure. Can somebody help with this issue? Thanks a lot!
 >>
 >> Roberto Lemoine
 >>
 >> --00504502ca5aac98410475b0ef84
 >> Content-Type: text/html; charset=ISO-8859-1
 >> Content-Transfer-Encoding: quoted-printable
 > k