f90mod_rules.py 9.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270
  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 isvariable(var)) 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 len(onlyvars) == 0 and len(notvars) == 1 and m['name'] in notvars:
  99. outmess(f"\t\t\tSkipping {m['name']} since there are no public vars/func in this module...\n")
  100. continue
  101. # gh-25186
  102. if m['name'] in usenames and containscommon(m):
  103. outmess(f"\t\t\tSkipping {m['name']} since it is in 'use' and contains a common block...\n")
  104. continue
  105. if onlyvars:
  106. outmess('\t\t Variables: %s\n' % (' '.join(onlyvars)))
  107. chooks = ['']
  108. def cadd(line, s=chooks):
  109. s[0] = '%s\n%s' % (s[0], line)
  110. ihooks = ['']
  111. def iadd(line, s=ihooks):
  112. s[0] = '%s\n%s' % (s[0], line)
  113. vrd = capi_maps.modsign2map(m)
  114. cadd('static FortranDataDef f2py_%s_def[] = {' % (m['name']))
  115. dadd('\\subsection{Fortran 90/95 module \\texttt{%s}}\n' % (m['name']))
  116. if hasnote(m):
  117. note = m['note']
  118. if isinstance(note, list):
  119. note = '\n'.join(note)
  120. dadd(note)
  121. if onlyvars:
  122. dadd('\\begin{description}')
  123. for n in onlyvars:
  124. var = m['vars'][n]
  125. modobjs.append(n)
  126. ct = capi_maps.getctype(var)
  127. at = capi_maps.c2capi_map[ct]
  128. dm = capi_maps.getarrdims(n, var)
  129. dms = dm['dims'].replace('*', '-1').strip()
  130. dms = dms.replace(':', '-1').strip()
  131. if not dms:
  132. dms = '-1'
  133. use_fgetdims2 = fgetdims2
  134. cadd('\t{"%s",%s,{{%s}},%s, %s},' %
  135. (undo_rmbadname1(n), dm['rank'], dms, at,
  136. capi_maps.get_elsize(var)))
  137. dadd('\\item[]{{}\\verb@%s@{}}' %
  138. (capi_maps.getarrdocsign(n, var)))
  139. if hasnote(var):
  140. note = var['note']
  141. if isinstance(note, list):
  142. note = '\n'.join(note)
  143. dadd('--- %s' % (note))
  144. if isallocatable(var):
  145. fargs.append('f2py_%s_getdims_%s' % (m['name'], n))
  146. efargs.append(fargs[-1])
  147. sargs.append(
  148. 'void (*%s)(int*,npy_intp*,void(*)(char*,npy_intp*),int*)' % (n))
  149. sargsp.append('void (*)(int*,npy_intp*,void(*)(char*,npy_intp*),int*)')
  150. iadd('\tf2py_%s_def[i_f2py++].func = %s;' % (m['name'], n))
  151. fadd('subroutine %s(r,s,f2pysetdata,flag)' % (fargs[-1]))
  152. fadd('use %s, only: d => %s\n' %
  153. (m['name'], undo_rmbadname1(n)))
  154. fadd('integer flag\n')
  155. fhooks[0] = fhooks[0] + fgetdims1
  156. dms = range(1, int(dm['rank']) + 1)
  157. fadd(' allocate(d(%s))\n' %
  158. (','.join(['s(%s)' % i for i in dms])))
  159. fhooks[0] = fhooks[0] + use_fgetdims2
  160. fadd('end subroutine %s' % (fargs[-1]))
  161. else:
  162. fargs.append(n)
  163. sargs.append('char *%s' % (n))
  164. sargsp.append('char*')
  165. iadd('\tf2py_%s_def[i_f2py++].data = %s;' % (m['name'], n))
  166. if onlyvars:
  167. dadd('\\end{description}')
  168. if hasbody(m):
  169. for b in m['body']:
  170. if not isroutine(b):
  171. outmess("f90mod_rules.buildhooks:"
  172. f" skipping {b['block']} {b['name']}\n")
  173. continue
  174. modobjs.append('%s()' % (b['name']))
  175. b['modulename'] = m['name']
  176. api, wrap = rules.buildapi(b)
  177. if isfunction(b):
  178. fhooks[0] = fhooks[0] + wrap
  179. fargs.append('f2pywrap_%s_%s' % (m['name'], b['name']))
  180. ifargs.append(func2subr.createfuncwrapper(b, signature=1))
  181. else:
  182. if wrap:
  183. fhooks[0] = fhooks[0] + wrap
  184. fargs.append('f2pywrap_%s_%s' % (m['name'], b['name']))
  185. ifargs.append(
  186. func2subr.createsubrwrapper(b, signature=1))
  187. else:
  188. fargs.append(b['name'])
  189. mfargs.append(fargs[-1])
  190. api['externroutines'] = []
  191. ar = applyrules(api, vrd)
  192. ar['docs'] = []
  193. ar['docshort'] = []
  194. ret = dictappend(ret, ar)
  195. cadd(('\t{"%s",-1,{{-1}},0,0,NULL,(void *)'
  196. 'f2py_rout_#modulename#_%s_%s,'
  197. 'doc_f2py_rout_#modulename#_%s_%s},')
  198. % (b['name'], m['name'], b['name'], m['name'], b['name']))
  199. sargs.append('char *%s' % (b['name']))
  200. sargsp.append('char *')
  201. iadd('\tf2py_%s_def[i_f2py++].data = %s;' %
  202. (m['name'], b['name']))
  203. cadd('\t{NULL}\n};\n')
  204. iadd('}')
  205. ihooks[0] = 'static void f2py_setup_%s(%s) {\n\tint i_f2py=0;%s' % (
  206. m['name'], ','.join(sargs), ihooks[0])
  207. if '_' in m['name']:
  208. F_FUNC = 'F_FUNC_US'
  209. else:
  210. F_FUNC = 'F_FUNC'
  211. iadd('extern void %s(f2pyinit%s,F2PYINIT%s)(void (*)(%s));'
  212. % (F_FUNC, m['name'], m['name'].upper(), ','.join(sargsp)))
  213. iadd('static void f2py_init_%s(void) {' % (m['name']))
  214. iadd('\t%s(f2pyinit%s,F2PYINIT%s)(f2py_setup_%s);'
  215. % (F_FUNC, m['name'], m['name'].upper(), m['name']))
  216. iadd('}\n')
  217. ret['f90modhooks'] = ret['f90modhooks'] + chooks + ihooks
  218. ret['initf90modhooks'] = ['\tPyDict_SetItemString(d, "%s", PyFortranObject_New(f2py_%s_def,f2py_init_%s));' % (
  219. m['name'], m['name'], m['name'])] + ret['initf90modhooks']
  220. fadd('')
  221. fadd('subroutine f2pyinit%s(f2pysetupfunc)' % (m['name']))
  222. if mfargs:
  223. for a in undo_rmbadname(mfargs):
  224. fadd('use %s, only : %s' % (m['name'], a))
  225. if ifargs:
  226. fadd(' '.join(['interface'] + ifargs))
  227. fadd('end interface')
  228. fadd('external f2pysetupfunc')
  229. if efargs:
  230. for a in undo_rmbadname(efargs):
  231. fadd('external %s' % (a))
  232. fadd('call f2pysetupfunc(%s)' % (','.join(undo_rmbadname(fargs))))
  233. fadd('end subroutine f2pyinit%s\n' % (m['name']))
  234. dadd('\n'.join(ret['latexdoc']).replace(
  235. r'\subsection{', r'\subsubsection{'))
  236. ret['latexdoc'] = []
  237. ret['docs'].append('"\t%s --- %s"' % (m['name'],
  238. ','.join(undo_rmbadname(modobjs))))
  239. ret['routine_defs'] = ''
  240. ret['doc'] = []
  241. ret['docshort'] = []
  242. ret['latexdoc'] = doc[0]
  243. if len(ret['docs']) <= 1:
  244. ret['docs'] = ''
  245. return ret, fhooks[0]