From bug-request at octave dot org Wed May 12 14:13:29 2004
Subject: Re: DASSL in 2.1.57 and 2.1.35
From: Shuo-Huan Hsu <shsu at ecn dot purdue dot edu>
To: "John W. Eaton" <jwe at bevo dot che dot wisc dot edu>
cc: bug at octave dot org
Date: Wed, 12 May 2004 14:13:04 -0500 (EST)

Hi, John,

If I use octave environment directly, I don't have any problem to solve
it.  My problem is when I try to build my own C++ application by linking
liboctave in ver 2.1.57 the following error shows up. However, it works
fine if I use the one in ver 2.1.35. I guess the problem comes from gcc
2.96 as Dmtiri mentioned. I'll upgrade my gcc and try it again.

Thanks a lot,
Shuo-Huan

On Wed, 12 May 2004, John W. Eaton wrote:

> On 12-May-2004, Shuo-Huan Hsu <shsu at ecn dot purdue dot edu> wrote:
>
> | Hi, All,
> |
> | I don't know if it is an adequate question to ask in this mail list. I try
> | to use DASSL to solve my DAE system using liboctave. When I use 2.1.35, it
> | works fine, but when I upgrade my octave to 2.1.57, my source code is
> | compiled but the system cannot be solved. DASSL gives me some error. I
> | attached my source code and here is the error message from 2.1.57.
> |
> |  ***MESSAGE FROM ROUTINE DDASSL IN LIBRARY SLATEC.
> |  ***POTENTIALLY RECOVERABLE ERROR, PROG ABORTED, TRACEBACK REQUESTED
> |  *  AT T =    0.000000E+00  AND STEPSIZE H =    9.536743E-11  THE INITIAL
> |  *  YPRIME COULD NOT BE COMPUTED
> |  *  ERROR NUMBER = -12
> |  *
> |  ***END OF MESSAGE
> |
> |  ***JOB ABORT DUE TO UNRECOVERED ERROR.
> | 0          ERROR MESSAGE SUMMARY
> |  LIBRARY    SUBROUTINE MESSAGE START             NERR     LEVEL     COUNT
> |  SLATEC     DDASSL     AT T =    0.000000E+       -12         1         1
> |
> | fatal: exception encountered in Fortran subroutine ddassl_
> |
> | P.S. I use RedHead Linux 7.3 and gcc 2.96.
> |
> | Is there any bug or am I doing something wrong?
>
> I am unable to reproduce your problem using the Octave code
>
>   function res = f (x, xdot, t)
>     k1 = 1;
>     k2 = 0.5;
>     res = [-xdot(1)-k1*x(1); k1*x(1)-k2*x(2); -xdot(3)+k2*x(2)];
>   endfunction
>
>   dassl ("f", [0.75; 0; 0], [0; 0; 0], 0:0.1:2)
>
> or when running your example code using
>
>   mkoctfile --link-stand-alone foo.cc
>
>
> jwe
>
>



-------------------------------------------------------------
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
-------------------------------------------------------------