0

0

Python调用Fortran子程序时ctypes参数类型匹配详解

聖光之護

聖光之護

发布时间:2026-02-15 14:13:26

|

423人浏览过

|

来源于php中文网

原创

Python调用Fortran子程序时ctypes参数类型匹配详解

本文深入解析python通过ctypes调用fortran函数时出现“typeerror: wrong type”错误的根本原因,重点说明fortran默认按引用传递整型参数、字符长度传递约定及ctypes类型声明的正确写法。

本文深入解析python通过ctypes调用fortran函数时出现“typeerror: wrong type”错误的根本原因,重点说明fortran默认按引用传递整型参数、字符长度传递约定及ctypes类型声明的正确写法。

在Python中通过ctypes调用Fortran编译生成的动态链接库(如 .dll 或 .so)是一项常见但易出错的任务。上述报错 ctypes.ArgumentError: argument 3: : wrong type 并非源于数据内容异常,而是函数签名(argtypes)与Fortran实际调用约定不匹配所致。

? 核心问题:Fortran的参数传递机制与ctypes类型映射失配

Fortran子程序默认所有参数均按地址(reference)传递,包括整数(integer)。而Python中ct.c_int(0)创建的是一个值对象(value object),其本身不是指针;若在argtypes中错误声明为ct.c_int,ctypes会在调用时尝试将该值“按值传入”,导致类型校验失败——这正是argument 3: wrong type的直接原因。

✅ 正确做法是:对Fortran中intent(out)或intent(inout)的整型参数,argtypes必须声明为指针类型:

ct.POINTER(ct.c_int)

而非 ct.c_int。

FlowMuse AI
FlowMuse AI

节点式AI视觉创作引擎

下载

立即学习Python免费学习笔记(深入)”;

此外,Fortran中character(*)参数虽以c_char_p(即char*)形式传入,但其长度信息不嵌入字符串本身,也不作为独立参数 interleaved 传入,而是作为额外的、位于参数列表末尾的整型参数(由编译器自动附加)。因此,len(x)应作为最后一个参数传递,且其类型推荐使用平台无关的 ct.c_size_t(对应C标准中的size_t),而非硬编码的ct.c_int。

✅ 修正后的完整调用代码

import ctypes as ct
import sys

lib = ct.CDLL('FTN_lib.dll')
fs = getattr(lib, 'PY_OBJFUN_MOD_mp_GETSIZE')

# ✅ 正确声明 argtypes:
# - c_char_p 对应 character(*)
# - POINTER(c_int) 对应每个 integer intent(out) 参数
# - c_size_t 对应隐式追加的字符串长度(置于末尾)
fs.argtypes = [
    ct.c_char_p,                    # file: character(*), intent(in)
    ct.POINTER(ct.c_int),           # n:   integer, intent(out)
    ct.POINTER(ct.c_int),           # nl:  integer, intent(out)
    ct.POINTER(ct.c_int),           # nnl: integer, intent(out)
    ct.c_size_t                     # len(file): 长度参数,放在最后
]
fs.restype = None

# 构造字节字符串(注意:Fortran通常期望ASCII/ISO-8859-1,避免UTF-8多字节)
x = b'x ' + sys.argv[1].encode('latin-1') + b' debug=900'

# 创建可被修改的整数变量(byref需基于c_int实例)
n   = ct.c_int(0)
nl  = ct.c_int(0)
nnl = ct.c_int(0)

# ✅ 调用顺序必须严格匹配 argtypes 声明:
# 字符串 → n指针 → nl指针 → nnl指针 → 字符串长度
fs(x, ct.byref(n), ct.byref(nl), ct.byref(nnl), len(x))

print(f"n={n.value}, nl={nl.value}, nnl={nnl.value}")

⚠️ 关键注意事项

  • 字符串编码:Fortran运行时通常不支持UTF-8多字节序列。建议统一使用 'latin-1' 编码(1:1字节映射),避免解码异常。
  • 长度参数位置:len(x) 必须作为最后一个参数传入,不可插入在中间。这是gfortran/ifort等主流编译器的ABI约定。
  • ct.byref() vs ct.pointer():ct.byref(x) 更高效且专用于函数调用;ct.pointer(x) 返回新指针对象,一般无需使用。
  • argtypes 不是可选的:虽然省略argtypes可能“偶然成功”,但会禁用类型检查,极易引发内存崩溃(segmentation fault)或静默错误。始终显式声明 argtypes 是生产环境最佳实践。
  • 调试技巧:若仍报错,可用 ctypes.util.find_library() 验证DLL路径,并用 objdump -t FTN_lib.dll | grep GETSIZE 检查符号是否导出(注意名称修饰,如gfortran常为小写+下划线)。

掌握Fortran与C混合编程的ABI细节,是打通科学计算Python化落地的关键一环。精准匹配argtypes不仅是语法要求,更是内存安全与结果可靠性的基石。

热门AI工具

更多
DeepSeek
DeepSeek

幻方量化公司旗下的开源大模型平台

豆包大模型
豆包大模型

字节跳动自主研发的一系列大型语言模型

通义千问
通义千问

阿里巴巴推出的全能AI助手

腾讯元宝
腾讯元宝

腾讯混元平台推出的AI助手

文心一言
文心一言

文心一言是百度开发的AI聊天机器人,通过对话可以生成各种形式的内容。

讯飞写作
讯飞写作

基于讯飞星火大模型的AI写作工具,可以快速生成新闻稿件、品宣文案、工作总结、心得体会等各种文文稿

即梦AI
即梦AI

一站式AI创作平台,免费AI图片和视频生成。

ChatGPT
ChatGPT

最最强大的AI聊天机器人程序,ChatGPT不单是聊天机器人,还能进行撰写邮件、视频脚本、文案、翻译、代码等任务。

相关专题

更多
js 字符串转数组
js 字符串转数组

js字符串转数组的方法:1、使用“split()”方法;2、使用“Array.from()”方法;3、使用for循环遍历;4、使用“Array.split()”方法。本专题为大家提供js字符串转数组的相关的文章、下载、课程内容,供大家免费下载体验。

551

2023.08.03

js截取字符串的方法
js截取字符串的方法

js截取字符串的方法有substring()方法、substr()方法、slice()方法、split()方法和slice()方法。本专题为大家提供字符串相关的文章、下载、课程内容,供大家免费下载体验。

214

2023.09.04

java基础知识汇总
java基础知识汇总

java基础知识有Java的历史和特点、Java的开发环境、Java的基本数据类型、变量和常量、运算符和表达式、控制语句、数组和字符串等等知识点。想要知道更多关于java基础知识的朋友,请阅读本专题下面的的有关文章,欢迎大家来php中文网学习。

1552

2023.10.24

字符串介绍
字符串介绍

字符串是一种数据类型,它可以是任何文本,包括字母、数字、符号等。字符串可以由不同的字符组成,例如空格、标点符号、数字等。在编程中,字符串通常用引号括起来,如单引号、双引号或反引号。想了解更多字符串的相关内容,可以阅读本专题下面的文章。

640

2023.11.24

java读取文件转成字符串的方法
java读取文件转成字符串的方法

Java8引入了新的文件I/O API,使用java.nio.file.Files类读取文件内容更加方便。对于较旧版本的Java,可以使用java.io.FileReader和java.io.BufferedReader来读取文件。在这些方法中,你需要将文件路径替换为你的实际文件路径,并且可能需要处理可能的IOException异常。想了解更多java的相关内容,可以阅读本专题下面的文章。

905

2024.03.22

php中定义字符串的方式
php中定义字符串的方式

php中定义字符串的方式:单引号;双引号;heredoc语法等等。想了解更多字符串的相关内容,可以阅读本专题下面的文章。

856

2024.04.29

go语言字符串相关教程
go语言字符串相关教程

本专题整合了go语言字符串相关教程,阅读专题下面的文章了解更多详细内容。

185

2025.07.29

c++字符串相关教程
c++字符串相关教程

本专题整合了c++字符串相关教程,阅读专题下面的文章了解更多详细内容。

89

2025.08.07

pixiv网页版官网登录与阅读指南_pixiv官网直达入口与在线访问方法
pixiv网页版官网登录与阅读指南_pixiv官网直达入口与在线访问方法

本专题系统整理pixiv网页版官网入口及登录访问方式,涵盖官网登录页面直达路径、在线阅读入口及快速进入方法说明,帮助用户高效找到pixiv官方网站,实现便捷、安全的网页端浏览与账号登录体验。

76

2026.02.13

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
最新Python教程 从入门到精通
最新Python教程 从入门到精通

共4课时 | 22.4万人学习

Django 教程
Django 教程

共28课时 | 4.3万人学习

SciPy 教程
SciPy 教程

共10课时 | 1.6万人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号