Re: F2CL under LWW a few warnings

Marco Antoniotti <[email protected]> Tue, 10 Aug 2004 17:11:26 -0400
Newsgroups gmane.lisp.clocc.general
Organization NYU Courant Bioinformatics Group
Message-ID <[email protected]>
More on F2CL

I am trying to convert the odepack-lsode library from netlib.

The library comes in two flavors: single and double precision.

The double precision library consists of three files.

opkda1.f
opkda2.f
opkdmain.f

The first two get converted with a few warnings, but the thirds 
generates a coercion error.

Here is a report.

Cheers
Marco


CL-USER 8 > (f2cl:f2cl "opkdmain.f" :verbose t)
;; "opkdmain.f" -> #P"opkdmain.lisp"
preprocessing begins ...
       SUBROUTINE DLSODE (F, NEQ, Y, T, TOUT, ITOL, RTOL, ATOL, ITASK,
      1                  ISTATE, IOPT, RWORK, LRW, IWORK, LIW, JAC, MF)
       EXTERNAL F, JAC
       INTEGER NEQ, ITOL, ITASK, ISTATE, IOPT, LRW, IWORK, LIW, MF
       DOUBLE PRECISION Y, T, TOUT, RTOL, ATOL, RWORK
       DIMENSION NEQ(*), Y(*), RTOL(*), ATOL(*), RWORK(LRW), IWORK(LIW)
       EXTERNAL DPREPJ, DSOLSY
       DOUBLE PRECISION DUMACH, DVNORM
       INTEGER INIT, MXSTEP, MXHNIL, NHNIL, NSLAST, NYH, IOWNS,

....

(call dintdy (tout \, 0 \, rwork (lyh) \, nyh \, y \, iflag))
Warning: Generating call to unknown function (DINTDY).  Check generated 
call!(if (iflag >< 0) go to 627)
(go to 627)
(t = tout)
(go to 420)
(tcrit = rwork (1))
(if ((tn - tcrit) * h > 0.0) go to 624)
(go to 624)
(hmx = abs (tn) + abs (h))
(ihit = abs (tn - tcrit) <= 100.0 * uround * hmx)

Error: Cannot coerce 0 to type LOGICAL.
   1 (abort) Return to level 0.
   2 Return to top loop level 0.

Type :b for backtrace, :c <option number> to proceed,  or :? for other 
options

F2CL 9 : 1 >



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