0

0

如何在 NumPy 中从 Python 调用 C 函数:源码级调用机制解析

花韻仙語

花韻仙語

发布时间:2026-01-18 14:13:40

|

855人浏览过

|

来源于php中文网

原创

如何在 NumPy 中从 Python 调用 C 函数:源码级调用机制解析

numpy 通过 c 扩展模块(如 `_multiarray_umath`)将 python 接口直接绑定到底层 c 函数,调用 `np.array()` 等函数时,python 解释器会跳过字节码执行,直接动态链接并调用编译后的共享库函数。

NumPy 的高性能核心源于其大量关键操作(如数组创建、内存布局管理、基础 ufunc 计算)均由 C 语言实现。但这些 C 函数并非通过 ctypes 或 cffi 等外部桥接方式调用,而是以 原生 C 扩展模块(CPython Extension Module) 的形式深度集成进 Python 解释器——这是最高效、最标准的 Python/C 互操作机制。

核心机制:C 扩展模块注册与方法绑定

当 NumPy 安装完成,其核心模块 numpy.core._multiarray_umath 实际对应一个已编译的共享对象文件(Linux/macOS 下为 .so,Windows 下为 .pyd),例如:

site-packages/numpy/core/_multiarray_umath.cpython-311-darwin.so

该二进制文件由 C 源码(主要位于 numpy/_core/src/multiarray/multiarraymodule.c)经 CPython C API 编译生成。其中最关键的结构是 PyMethodDef 数组,它显式声明了哪些 C 函数将暴露为 Python 可调用的方法:

// numpy/_core/src/multiarray/multiarraymodule.c
static PyMethodDef array_module_methods[] = {
    // ...
    {"array",
        (PyCFunction)array_array,           // 绑定到 C 函数 array_array()
        METH_FASTCALL | METH_KEYWORDS,      // 支持快速调用 + 关键字参数
        NULL},
    // ...
    {NULL, NULL, 0, NULL}  // 结束标记
};

随后,该方法表被传入模块定义结构体,并在模块初始化时注册:

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

人声去除
人声去除

用强大的AI算法将声音从音乐中分离出来

下载
static struct PyModuleDef moduledef = {
    PyModuleDef_HEAD_INIT,
    "_multiarray_umath",     // 模块名(Python 中 import 的名称)
    NULL,
    -1,
    array_module_methods,    // 关键:绑定所有 C 函数
    NULL, NULL, NULL, NULL
};

PyMODINIT_FUNC PyInit__multiarray_umath(void) {
    PyObject *m;
    m = PyModule_Create(&moduledef);
    if (m == NULL) return NULL;
    // ... 其他初始化逻辑(如添加常量、类型等)
    return m;
}

当 Python 执行 import numpy.core._multiarray_umath 时,解释器自动加载对应 .so 文件,调用 PyInit__multiarray_umath 初始化模块,并将 array_module_methods 中声明的每个条目(如 "array")作为模块属性注入——此时 numpy.core._multiarray_umath.array 即是一个指向 array_array C 函数的 builtin_function_or_method 对象。

验证:Python 层面的直接映射

你可以通过以下代码验证该绑定关系的真实性:

import numpy as np
import numpy.core._multiarray_umath as mu

# 二者完全等价,是同一对象
print(mu.array is np.array)  # True

# 尝试反汇编:无法获取字节码 → 确认为内置 C 函数
import dis
try:
    dis.dis(np.array)
except TypeError as e:
    print(e)  # "don't know how to disassemble builtin_function_or_method objects"
✅ 注意:np.array 并非 Python 函数包装器(如 @array_function_from_c_func_and_dispatcher 装饰的 empty_like),而是直连 C 函数的原始入口。装饰器方案用于支持 __array_function__ 协议的高级封装函数,而 array、zeros、copy 等基础构造函数走的是更底层、零开销的 C 扩展路径。

补充说明:从 _multiarray_umath.array 到 np.array 的导入链

虽然调用发生在 _multiarray_umath 模块,但用户最终使用的是 np.array。这一映射通过多层 __init__.py 导入和 sys.modules 动态挂载实现(例如 numpy/__init__.py 中 from numpy.core import multiarray,再经 multiarray.py 中 from numpy.core._multiarray_umath import array),但该过程仅影响命名空间组织,不引入任何额外调用开销——np.array() 调用仍直接路由至 C 函数。

总结与建议

  • 本质:NumPy 的核心函数是标准 CPython C 扩展,遵循 Extending Python with C 规范;
  • 性能关键:METH_FASTCALL | METH_KEYWORDS 标志启用 Python 3.6+ 最快的参数传递协议,避免传统 *args, **kwargs 解包开销;
  • ⚠️ 调试提示:若需定位某函数的 C 实现,可先查 multiarraymodule.c 中的 PyMethodDef 条目,再按函数名(如 array_array)搜索具体实现;
  • ? 进阶参考:PyCFunction 类型、PyModule_Create 流程、CPython 的 Objects/methodobject.c 是理解底层调用的关键源码节点。

掌握这一机制,不仅能深入理解 NumPy 的性能根基,也为开发高性能科学计算扩展提供了标准范式。

热门AI工具

更多
DeepSeek
DeepSeek

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

豆包大模型
豆包大模型

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

通义千问
通义千问

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

腾讯元宝
腾讯元宝

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

文心一言
文心一言

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

讯飞写作
讯飞写作

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

即梦AI
即梦AI

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

ChatGPT
ChatGPT

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

相关专题

更多
golang结构体相关大全
golang结构体相关大全

本专题整合了golang结构体相关大全,想了解更多内容,请阅读专题下面的文章。

470

2025.06.09

golang结构体方法
golang结构体方法

本专题整合了golang结构体相关内容,请阅读专题下面的文章了解更多。

201

2025.07.04

硬盘接口类型介绍
硬盘接口类型介绍

硬盘接口类型有IDE、SATA、SCSI、Fibre Channel、USB、eSATA、mSATA、PCIe等等。详细介绍:1、IDE接口是一种并行接口,主要用于连接硬盘和光驱等设备,它主要有两种类型:ATA和ATAPI,IDE接口已经逐渐被SATA接口;2、SATA接口是一种串行接口,相较于IDE接口,它具有更高的传输速度、更低的功耗和更小的体积;3、SCSI接口等等。

1846

2023.10.19

PHP接口编写教程
PHP接口编写教程

本专题整合了PHP接口编写教程,阅读专题下面的文章了解更多详细内容。

614

2025.10.17

php8.4实现接口限流的教程
php8.4实现接口限流的教程

PHP8.4本身不内置限流功能,需借助Redis(令牌桶)或Swoole(漏桶)实现;文件锁因I/O瓶颈、无跨机共享、秒级精度等缺陷不适用高并发场景。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

2354

2025.12.29

java接口相关教程
java接口相关教程

本专题整合了java接口相关内容,阅读专题下面的文章了解更多详细内容。

47

2026.01.19

堆和栈的区别
堆和栈的区别

堆和栈的区别:1、内存分配方式不同;2、大小不同;3、数据访问方式不同;4、数据的生命周期。本专题为大家提供堆和栈的区别的相关的文章、下载、课程内容,供大家免费下载体验。

435

2023.07.18

堆和栈区别
堆和栈区别

堆(Heap)和栈(Stack)是计算机中两种常见的内存分配机制。它们在内存管理的方式、分配方式以及使用场景上有很大的区别。本文将详细介绍堆和栈的特点、区别以及各自的使用场景。php中文网给大家带来了相关的教程以及文章欢迎大家前来学习阅读。

601

2023.08.10

JavaScript浏览器渲染机制与前端性能优化实践
JavaScript浏览器渲染机制与前端性能优化实践

本专题围绕 JavaScript 在浏览器中的执行与渲染机制展开,系统讲解 DOM 构建、CSSOM 解析、重排与重绘原理,以及关键渲染路径优化方法。内容涵盖事件循环机制、异步任务调度、资源加载优化、代码拆分与懒加载等性能优化策略。通过真实前端项目案例,帮助开发者理解浏览器底层工作原理,并掌握提升网页加载速度与交互体验的实用技巧。

23

2026.03.06

热门下载

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

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
PostgreSQL 教程
PostgreSQL 教程

共48课时 | 10.3万人学习

Git 教程
Git 教程

共21课时 | 4.1万人学习

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

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