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