I decided to give wrapping this code a try:
http://morrislab.med.utoronto.ca/~dwf/GLMnet.f90
I'm afraid my Fortran skills are fairly limited, but I do know that
gfortran compiles it fine. f2py run on this file produces lots of
errors of the form,
Reading fortran codes...
Reading file 'GLMnet.f90' (format:fix)
Line #263 in GLMnet.f90:" real
x(no,ni),y(no),w(no),vp(ni),ca(nx,nlam) 353"
updatevars: could not crack entity declaration "ca(nx,nlam)353".
Ignoring.
Line #264 in GLMnet.f90:" real
ulam(nlam),a0(nlam),rsq(nlam),alm(nlam) 354"
updatevars: could not crack entity declaration "alm(nlam)354".
Ignoring.
Line #265 in GLMnet.f90:" integer
jd(*),ia(nx),nin(nlam) 355"
updatevars: could not crack entity declaration "nin(nlam)355".
Ignoring.
Line #289 in GLMnet.f90:" real
x(no,ni),y(no),w(no),vp(ni),ulam(nlam) 378"
updatevars: could not crack entity declaration "ulam(nlam)378".
Ignoring.
Line #290 in GLMnet.f90:" real
ca(nx,nlam),a0(nlam),rsq(nlam),alm(nlam) 379"
updatevars: could not crack entity declaration "alm(nlam)379".
Ignoring.
Line #291 in GLMnet.f90:" integer
jd(*),ia(nx),nin(nlam) 380"
updatevars: could not crack entity declaration "nin(nlam)380".
Ignoring.
Line #306 in GLMnet.f90:" call
chkvars(no,ni,x,ju) 392"
analyzeline: No name/args pattern found for li
Is it the numbers that it is objecting to (I'm assuming these are some
sort of punchcard thing)? Do I need to modify the code in some way to
make it f2py-friendly?
Thanks,
David
_______________________________________________
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion
