Re: F2CL under LWW a few warnings
Raymond Toy <[email protected]> Sat, 14 Aug 2004 11:38:02 -0400
| Newsgroups | gmane.lisp.clocc.general |
|---|---|
| Message-ID | <[email protected]> |
Marco Antoniotti wrote:
> ? (sorry if the above is not all that correct; you get the idea)
>
> It is tedious, but it can be done.
With the latest changes to f2cl, f2cl seems able to compile odepack and
can run some of the demo programs.
I didn't have to modify any of the ODEPACK routines, but did have to
modify slightly the demo driver program. This was because ODEPACK
assumes some variables are actually arrays but calls the routine with a
single variable. Not a problem in Fortran, but Lisp doesn't like that.
Anyway, here is sample output from the first demo program. It doesn't
finish the entire first demo because of some errors, but at least
there's hope.
Ray
CL-USER> (*main*)
Demonstration program for DLSODE package
Problem 1: Van der Pol oscillator:
xdotdot - 3*(1 - x**2)*xdot + x = 0, x(0) = 2, xdot(0) = 0
neq = 2
itol = 1 rtol = 0.0D+01 atol = 1.0D-06
Solution with mf = 10
t x xdot nq h
0.13928D+01 0.16801D+01 -0.291D+00 3 0.123D+00
0.36076D+01 -0.77986D-04 -0.317D+01 5 0.217D-01
0.58224D+01 -0.16801D+01 0.291D+00 3 0.475D-01
0.80372D+01 0.11669D-03 0.317D+01 5 0.234D-01
Final statistics for this run:
rwork size = 52 iwork size = 20
number of steps = 297
number of f-s = 352
(excluding J-s) = 352
number of J-s = 0
error overrun = 0.12D+03
-------------------------------------------------------
SF.Net email is sponsored by Shop4tech.com-Lowest price on Blank Media
100pk Sonic DVD-R 4x for only $29 -100pk Sonic DVD+R for only $33
Save 50% off Retail on Ink & Toner - Free Shipping and Free Gift.
http://www.shop4tech.com/z/Inkjet_Cartridges/9_108_r285