func2subr.py 10 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323
  1. """
  2. Rules for building C/API module with 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. import copy
  10. from .auxfuncs import (
  11. getfortranname, isexternal, isfunction, isfunction_wrap, isintent_in,
  12. isintent_out, islogicalfunction, ismoduleroutine, isscalar,
  13. issubroutine, issubroutine_wrap, outmess, show
  14. )
  15. from ._isocbind import isoc_kindmap
  16. def var2fixfortran(vars, a, fa=None, f90mode=None):
  17. if fa is None:
  18. fa = a
  19. if a not in vars:
  20. show(vars)
  21. outmess('var2fixfortran: No definition for argument "%s".\n' % a)
  22. return ''
  23. if 'typespec' not in vars[a]:
  24. show(vars[a])
  25. outmess('var2fixfortran: No typespec for argument "%s".\n' % a)
  26. return ''
  27. vardef = vars[a]['typespec']
  28. if vardef == 'type' and 'typename' in vars[a]:
  29. vardef = '%s(%s)' % (vardef, vars[a]['typename'])
  30. selector = {}
  31. lk = ''
  32. if 'kindselector' in vars[a]:
  33. selector = vars[a]['kindselector']
  34. lk = 'kind'
  35. elif 'charselector' in vars[a]:
  36. selector = vars[a]['charselector']
  37. lk = 'len'
  38. if '*' in selector:
  39. if f90mode:
  40. if selector['*'] in ['*', ':', '(*)']:
  41. vardef = '%s(len=*)' % (vardef)
  42. else:
  43. vardef = '%s(%s=%s)' % (vardef, lk, selector['*'])
  44. else:
  45. if selector['*'] in ['*', ':']:
  46. vardef = '%s*(%s)' % (vardef, selector['*'])
  47. else:
  48. vardef = '%s*%s' % (vardef, selector['*'])
  49. else:
  50. if 'len' in selector:
  51. vardef = '%s(len=%s' % (vardef, selector['len'])
  52. if 'kind' in selector:
  53. vardef = '%s,kind=%s)' % (vardef, selector['kind'])
  54. else:
  55. vardef = '%s)' % (vardef)
  56. elif 'kind' in selector:
  57. vardef = '%s(kind=%s)' % (vardef, selector['kind'])
  58. vardef = '%s %s' % (vardef, fa)
  59. if 'dimension' in vars[a]:
  60. vardef = '%s(%s)' % (vardef, ','.join(vars[a]['dimension']))
  61. return vardef
  62. def useiso_c_binding(rout):
  63. useisoc = False
  64. for key, value in rout['vars'].items():
  65. kind_value = value.get('kindselector', {}).get('kind')
  66. if kind_value in isoc_kindmap:
  67. return True
  68. return useisoc
  69. def createfuncwrapper(rout, signature=0):
  70. assert isfunction(rout)
  71. extra_args = []
  72. vars = rout['vars']
  73. for a in rout['args']:
  74. v = rout['vars'][a]
  75. for i, d in enumerate(v.get('dimension', [])):
  76. if d == ':':
  77. dn = 'f2py_%s_d%s' % (a, i)
  78. dv = dict(typespec='integer', intent=['hide'])
  79. dv['='] = 'shape(%s, %s)' % (a, i)
  80. extra_args.append(dn)
  81. vars[dn] = dv
  82. v['dimension'][i] = dn
  83. rout['args'].extend(extra_args)
  84. need_interface = bool(extra_args)
  85. ret = ['']
  86. def add(line, ret=ret):
  87. ret[0] = '%s\n %s' % (ret[0], line)
  88. name = rout['name']
  89. fortranname = getfortranname(rout)
  90. f90mode = ismoduleroutine(rout)
  91. newname = '%sf2pywrap' % (name)
  92. if newname not in vars:
  93. vars[newname] = vars[name]
  94. args = [newname] + rout['args'][1:]
  95. else:
  96. args = [newname] + rout['args']
  97. l_tmpl = var2fixfortran(vars, name, '@@@NAME@@@', f90mode)
  98. if l_tmpl[:13] == 'character*(*)':
  99. if f90mode:
  100. l_tmpl = 'character(len=10)' + l_tmpl[13:]
  101. else:
  102. l_tmpl = 'character*10' + l_tmpl[13:]
  103. charselect = vars[name]['charselector']
  104. if charselect.get('*', '') == '(*)':
  105. charselect['*'] = '10'
  106. l1 = l_tmpl.replace('@@@NAME@@@', newname)
  107. rl = None
  108. useisoc = useiso_c_binding(rout)
  109. sargs = ', '.join(args)
  110. if f90mode:
  111. # gh-23598 fix warning
  112. # Essentially, this gets called again with modules where the name of the
  113. # function is added to the arguments, which is not required, and removed
  114. sargs = sargs.replace(f"{name}, ", '')
  115. args = [arg for arg in args if arg != name]
  116. rout['args'] = args
  117. add('subroutine f2pywrap_%s_%s (%s)' %
  118. (rout['modulename'], name, sargs))
  119. if not signature:
  120. add('use %s, only : %s' % (rout['modulename'], fortranname))
  121. if useisoc:
  122. add('use iso_c_binding')
  123. else:
  124. add('subroutine f2pywrap%s (%s)' % (name, sargs))
  125. if useisoc:
  126. add('use iso_c_binding')
  127. if not need_interface:
  128. add('external %s' % (fortranname))
  129. rl = l_tmpl.replace('@@@NAME@@@', '') + ' ' + fortranname
  130. if need_interface:
  131. for line in rout['saved_interface'].split('\n'):
  132. if line.lstrip().startswith('use ') and '__user__' not in line:
  133. add(line)
  134. args = args[1:]
  135. dumped_args = []
  136. for a in args:
  137. if isexternal(vars[a]):
  138. add('external %s' % (a))
  139. dumped_args.append(a)
  140. for a in args:
  141. if a in dumped_args:
  142. continue
  143. if isscalar(vars[a]):
  144. add(var2fixfortran(vars, a, f90mode=f90mode))
  145. dumped_args.append(a)
  146. for a in args:
  147. if a in dumped_args:
  148. continue
  149. if isintent_in(vars[a]):
  150. add(var2fixfortran(vars, a, f90mode=f90mode))
  151. dumped_args.append(a)
  152. for a in args:
  153. if a in dumped_args:
  154. continue
  155. add(var2fixfortran(vars, a, f90mode=f90mode))
  156. add(l1)
  157. if rl is not None:
  158. add(rl)
  159. if need_interface:
  160. if f90mode:
  161. # f90 module already defines needed interface
  162. pass
  163. else:
  164. add('interface')
  165. add(rout['saved_interface'].lstrip())
  166. add('end interface')
  167. sargs = ', '.join([a for a in args if a not in extra_args])
  168. if not signature:
  169. if islogicalfunction(rout):
  170. add('%s = .not.(.not.%s(%s))' % (newname, fortranname, sargs))
  171. else:
  172. add('%s = %s(%s)' % (newname, fortranname, sargs))
  173. if f90mode:
  174. add('end subroutine f2pywrap_%s_%s' % (rout['modulename'], name))
  175. else:
  176. add('end')
  177. return ret[0]
  178. def createsubrwrapper(rout, signature=0):
  179. assert issubroutine(rout)
  180. extra_args = []
  181. vars = rout['vars']
  182. for a in rout['args']:
  183. v = rout['vars'][a]
  184. for i, d in enumerate(v.get('dimension', [])):
  185. if d == ':':
  186. dn = 'f2py_%s_d%s' % (a, i)
  187. dv = dict(typespec='integer', intent=['hide'])
  188. dv['='] = 'shape(%s, %s)' % (a, i)
  189. extra_args.append(dn)
  190. vars[dn] = dv
  191. v['dimension'][i] = dn
  192. rout['args'].extend(extra_args)
  193. need_interface = bool(extra_args)
  194. ret = ['']
  195. def add(line, ret=ret):
  196. ret[0] = '%s\n %s' % (ret[0], line)
  197. name = rout['name']
  198. fortranname = getfortranname(rout)
  199. f90mode = ismoduleroutine(rout)
  200. args = rout['args']
  201. useisoc = useiso_c_binding(rout)
  202. sargs = ', '.join(args)
  203. if f90mode:
  204. add('subroutine f2pywrap_%s_%s (%s)' %
  205. (rout['modulename'], name, sargs))
  206. if useisoc:
  207. add('use iso_c_binding')
  208. if not signature:
  209. add('use %s, only : %s' % (rout['modulename'], fortranname))
  210. else:
  211. add('subroutine f2pywrap%s (%s)' % (name, sargs))
  212. if useisoc:
  213. add('use iso_c_binding')
  214. if not need_interface:
  215. add('external %s' % (fortranname))
  216. if need_interface:
  217. for line in rout['saved_interface'].split('\n'):
  218. if line.lstrip().startswith('use ') and '__user__' not in line:
  219. add(line)
  220. dumped_args = []
  221. for a in args:
  222. if isexternal(vars[a]):
  223. add('external %s' % (a))
  224. dumped_args.append(a)
  225. for a in args:
  226. if a in dumped_args:
  227. continue
  228. if isscalar(vars[a]):
  229. add(var2fixfortran(vars, a, f90mode=f90mode))
  230. dumped_args.append(a)
  231. for a in args:
  232. if a in dumped_args:
  233. continue
  234. add(var2fixfortran(vars, a, f90mode=f90mode))
  235. if need_interface:
  236. if f90mode:
  237. # f90 module already defines needed interface
  238. pass
  239. else:
  240. add('interface')
  241. for line in rout['saved_interface'].split('\n'):
  242. if line.lstrip().startswith('use ') and '__user__' in line:
  243. continue
  244. add(line)
  245. add('end interface')
  246. sargs = ', '.join([a for a in args if a not in extra_args])
  247. if not signature:
  248. add('call %s(%s)' % (fortranname, sargs))
  249. if f90mode:
  250. add('end subroutine f2pywrap_%s_%s' % (rout['modulename'], name))
  251. else:
  252. add('end')
  253. return ret[0]
  254. def assubr(rout):
  255. if isfunction_wrap(rout):
  256. fortranname = getfortranname(rout)
  257. name = rout['name']
  258. outmess('\t\tCreating wrapper for Fortran function "%s"("%s")...\n' % (
  259. name, fortranname))
  260. rout = copy.copy(rout)
  261. fname = name
  262. rname = fname
  263. if 'result' in rout:
  264. rname = rout['result']
  265. rout['vars'][fname] = rout['vars'][rname]
  266. fvar = rout['vars'][fname]
  267. if not isintent_out(fvar):
  268. if 'intent' not in fvar:
  269. fvar['intent'] = []
  270. fvar['intent'].append('out')
  271. flag = 1
  272. for i in fvar['intent']:
  273. if i.startswith('out='):
  274. flag = 0
  275. break
  276. if flag:
  277. fvar['intent'].append('out=%s' % (rname))
  278. rout['args'][:] = [fname] + rout['args']
  279. return rout, createfuncwrapper(rout)
  280. if issubroutine_wrap(rout):
  281. fortranname = getfortranname(rout)
  282. name = rout['name']
  283. outmess('\t\tCreating wrapper for Fortran subroutine "%s"("%s")...\n'
  284. % (name, fortranname))
  285. rout = copy.copy(rout)
  286. return rout, createsubrwrapper(rout)
  287. return rout, ''