f90mod_rules.py 9.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264
  1. """
  2. Build F90 module support for f2py2e.
  3. Copyright 1999 -- 2011 Pearu Peterson all rights reserved.
  4. Copyright 2011 -- present NumPy Developers.
  5. Permission to use, modify, and distribute this software is given under the
  6. terms of the NumPy License.
  7. NO WARRANTY IS EXPRESSED OR IMPLIED. USE AT YOUR OWN RISK.
  8. """
  9. __version__ = "$Revision: 1.27 $"[10:-1]
  10. f2py_version = 'See `f2py -v`'
  11. import numpy as np
  12. from . import capi_maps
  13. from . import func2subr
  14. from .crackfortran import undo_rmbadname, undo_rmbadname1
  15. # The environment provided by auxfuncs.py is needed for some calls to eval.
  16. # As the needed functions cannot be determined by static inspection of the
  17. # code, it is safest to use import * pending a major refactoring of f2py.
  18. from .auxfuncs import *
  19. options = {}
  20. def findf90modules(m):
  21. if ismodule(m):
  22. return [m]
  23. if not hasbody(m):
  24. return []
  25. ret = []
  26. for b in m['body']:
  27. if ismodule(b):
  28. ret.append(b)
  29. else:
  30. ret = ret + findf90modules(b)
  31. return ret
  32. fgetdims1 = """\
  33. external f2pysetdata
  34. logical ns
  35. integer r,i
  36. integer(%d) s(*)
  37. ns = .FALSE.
  38. if (allocated(d)) then
  39. do i=1,r
  40. if ((size(d,i).ne.s(i)).and.(s(i).ge.0)) then
  41. ns = .TRUE.
  42. end if
  43. end do
  44. if (ns) then
  45. deallocate(d)
  46. end if
  47. end if
  48. if ((.not.allocated(d)).and.(s(1).ge.1)) then""" % np.intp().itemsize
  49. fgetdims2 = """\
  50. end if
  51. if (allocated(d)) then
  52. do i=1,r
  53. s(i) = size(d,i)
  54. end do
  55. end if
  56. flag = 1
  57. call f2pysetdata(d,allocated(d))"""
  58. fgetdims2_sa = """\
  59. end if
  60. if (allocated(d)) then
  61. do i=1,r
  62. s(i) = size(d,i)
  63. end do
  64. !s(r) must be equal to len(d(1))
  65. end if
  66. flag = 2
  67. call f2pysetdata(d,allocated(d))"""
  68. def buildhooks(pymod):
  69. from . import rules
  70. ret = {'f90modhooks': [], 'initf90modhooks': [], 'body': [],
  71. 'need': ['F_FUNC', 'arrayobject.h'],
  72. 'separatorsfor': {'includes0': '\n', 'includes': '\n'},
  73. 'docs': ['"Fortran 90/95 modules:\\n"'],
  74. 'latexdoc': []}
  75. fhooks = ['']
  76. def fadd(line, s=fhooks):
  77. s[0] = '%s\n %s' % (s[0], line)
  78. doc = ['']
  79. def dadd(line, s=doc):
  80. s[0] = '%s\n%s' % (s[0], line)
  81. usenames = getuseblocks(pymod)
  82. for m in findf90modules(pymod):
  83. sargs, fargs, efargs, modobjs, notvars, onlyvars = [], [], [], [], [
  84. m['name']], []
  85. sargsp = []
  86. ifargs = []
  87. mfargs = []
  88. if hasbody(m):
  89. for b in m['body']:
  90. notvars.append(b['name'])
  91. for n in m['vars'].keys():
  92. var = m['vars'][n]
  93. if (n not in notvars) and (not l_or(isintent_hide, isprivate)(var)):
  94. onlyvars.append(n)
  95. mfargs.append(n)
  96. outmess('\t\tConstructing F90 module support for "%s"...\n' %
  97. (m['name']))
  98. if m['name'] in usenames and not onlyvars:
  99. outmess(f"\t\t\tSkipping {m['name']} since it is in 'use'...\n")
  100. continue
  101. if onlyvars:
  102. outmess('\t\t Variables: %s\n' % (' '.join(onlyvars)))
  103. chooks = ['']
  104. def cadd(line, s=chooks):
  105. s[0] = '%s\n%s' % (s[0], line)
  106. ihooks = ['']
  107. def iadd(line, s=ihooks):
  108. s[0] = '%s\n%s' % (s[0], line)
  109. vrd = capi_maps.modsign2map(m)
  110. cadd('static FortranDataDef f2py_%s_def[] = {' % (m['name']))
  111. dadd('\\subsection{Fortran 90/95 module \\texttt{%s}}\n' % (m['name']))
  112. if hasnote(m):
  113. note = m['note']
  114. if isinstance(note, list):
  115. note = '\n'.join(note)
  116. dadd(note)
  117. if onlyvars:
  118. dadd('\\begin{description}')
  119. for n in onlyvars:
  120. var = m['vars'][n]
  121. modobjs.append(n)
  122. ct = capi_maps.getctype(var)
  123. at = capi_maps.c2capi_map[ct]
  124. dm = capi_maps.getarrdims(n, var)
  125. dms = dm['dims'].replace('*', '-1').strip()
  126. dms = dms.replace(':', '-1').strip()
  127. if not dms:
  128. dms = '-1'
  129. use_fgetdims2 = fgetdims2
  130. cadd('\t{"%s",%s,{{%s}},%s, %s},' %
  131. (undo_rmbadname1(n), dm['rank'], dms, at,
  132. capi_maps.get_elsize(var)))
  133. dadd('\\item[]{{}\\verb@%s@{}}' %
  134. (capi_maps.getarrdocsign(n, var)))
  135. if hasnote(var):
  136. note = var['note']
  137. if isinstance(note, list):
  138. note = '\n'.join(note)
  139. dadd('--- %s' % (note))
  140. if isallocatable(var):
  141. fargs.append('f2py_%s_getdims_%s' % (m['name'], n))
  142. efargs.append(fargs[-1])
  143. sargs.append(
  144. 'void (*%s)(int*,npy_intp*,void(*)(char*,npy_intp*),int*)' % (n))
  145. sargsp.append('void (*)(int*,npy_intp*,void(*)(char*,npy_intp*),int*)')
  146. iadd('\tf2py_%s_def[i_f2py++].func = %s;' % (m['name'], n))
  147. fadd('subroutine %s(r,s,f2pysetdata,flag)' % (fargs[-1]))
  148. fadd('use %s, only: d => %s\n' %
  149. (m['name'], undo_rmbadname1(n)))
  150. fadd('integer flag\n')
  151. fhooks[0] = fhooks[0] + fgetdims1
  152. dms = range(1, int(dm['rank']) + 1)
  153. fadd(' allocate(d(%s))\n' %
  154. (','.join(['s(%s)' % i for i in dms])))
  155. fhooks[0] = fhooks[0] + use_fgetdims2
  156. fadd('end subroutine %s' % (fargs[-1]))
  157. else:
  158. fargs.append(n)
  159. sargs.append('char *%s' % (n))
  160. sargsp.append('char*')
  161. iadd('\tf2py_%s_def[i_f2py++].data = %s;' % (m['name'], n))
  162. if onlyvars:
  163. dadd('\\end{description}')
  164. if hasbody(m):
  165. for b in m['body']:
  166. if not isroutine(b):
  167. outmess("f90mod_rules.buildhooks:"
  168. f" skipping {b['block']} {b['name']}\n")
  169. continue
  170. modobjs.append('%s()' % (b['name']))
  171. b['modulename'] = m['name']
  172. api, wrap = rules.buildapi(b)
  173. if isfunction(b):
  174. fhooks[0] = fhooks[0] + wrap
  175. fargs.append('f2pywrap_%s_%s' % (m['name'], b['name']))
  176. ifargs.append(func2subr.createfuncwrapper(b, signature=1))
  177. else:
  178. if wrap:
  179. fhooks[0] = fhooks[0] + wrap
  180. fargs.append('f2pywrap_%s_%s' % (m['name'], b['name']))
  181. ifargs.append(
  182. func2subr.createsubrwrapper(b, signature=1))
  183. else:
  184. fargs.append(b['name'])
  185. mfargs.append(fargs[-1])
  186. api['externroutines'] = []
  187. ar = applyrules(api, vrd)
  188. ar['docs'] = []
  189. ar['docshort'] = []
  190. ret = dictappend(ret, ar)
  191. cadd(('\t{"%s",-1,{{-1}},0,0,NULL,(void *)'
  192. 'f2py_rout_#modulename#_%s_%s,'
  193. 'doc_f2py_rout_#modulename#_%s_%s},')
  194. % (b['name'], m['name'], b['name'], m['name'], b['name']))
  195. sargs.append('char *%s' % (b['name']))
  196. sargsp.append('char *')
  197. iadd('\tf2py_%s_def[i_f2py++].data = %s;' %
  198. (m['name'], b['name']))
  199. cadd('\t{NULL}\n};\n')
  200. iadd('}')
  201. ihooks[0] = 'static void f2py_setup_%s(%s) {\n\tint i_f2py=0;%s' % (
  202. m['name'], ','.join(sargs), ihooks[0])
  203. if '_' in m['name']:
  204. F_FUNC = 'F_FUNC_US'
  205. else:
  206. F_FUNC = 'F_FUNC'
  207. iadd('extern void %s(f2pyinit%s,F2PYINIT%s)(void (*)(%s));'
  208. % (F_FUNC, m['name'], m['name'].upper(), ','.join(sargsp)))
  209. iadd('static void f2py_init_%s(void) {' % (m['name']))
  210. iadd('\t%s(f2pyinit%s,F2PYINIT%s)(f2py_setup_%s);'
  211. % (F_FUNC, m['name'], m['name'].upper(), m['name']))
  212. iadd('}\n')
  213. ret['f90modhooks'] = ret['f90modhooks'] + chooks + ihooks
  214. ret['initf90modhooks'] = ['\tPyDict_SetItemString(d, "%s", PyFortranObject_New(f2py_%s_def,f2py_init_%s));' % (
  215. m['name'], m['name'], m['name'])] + ret['initf90modhooks']
  216. fadd('')
  217. fadd('subroutine f2pyinit%s(f2pysetupfunc)' % (m['name']))
  218. if mfargs:
  219. for a in undo_rmbadname(mfargs):
  220. fadd('use %s, only : %s' % (m['name'], a))
  221. if ifargs:
  222. fadd(' '.join(['interface'] + ifargs))
  223. fadd('end interface')
  224. fadd('external f2pysetupfunc')
  225. if efargs:
  226. for a in undo_rmbadname(efargs):
  227. fadd('external %s' % (a))
  228. fadd('call f2pysetupfunc(%s)' % (','.join(undo_rmbadname(fargs))))
  229. fadd('end subroutine f2pyinit%s\n' % (m['name']))
  230. dadd('\n'.join(ret['latexdoc']).replace(
  231. r'\subsection{', r'\subsubsection{'))
  232. ret['latexdoc'] = []
  233. ret['docs'].append('"\t%s --- %s"' % (m['name'],
  234. ','.join(undo_rmbadname(modobjs))))
  235. ret['routine_defs'] = ''
  236. ret['doc'] = []
  237. ret['docshort'] = []
  238. ret['latexdoc'] = doc[0]
  239. if len(ret['docs']) <= 1:
  240. ret['docs'] = ''
  241. return ret, fhooks[0]