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