C++扩展python(三)boost.Python实现(附swig比较)

2023-12-27 02:32

本文主要是介绍C++扩展python(三)boost.Python实现(附swig比较),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

本文98%内容转载自博客:https://www.jianshu.com/p/1b4235c90567

除了可以使用 SWIG 包装 C/C++ 之外,还可以使用 C++ 的 boost.Python 库包装 C/C++ MPI 程序以供 mpi4py 调用。

关于各种包装方式的优缺点可以看这一篇博客,或者这一篇。

 

话不多说下面我们以包装 C++ MPI 程序为例对此进行介绍。

Boost.Python 是 C++ 的 Boost 库的一部分,用来在 C++ 代码中调用 Python 代码以及在 Python 代码中调用 C++ 代码,并且避免用户直接操作指针。Boost.Python 库可以让 C++ 和 Python 几乎无缝结合。这里不介绍 Boost.Python 的详细使用方法,感兴趣的读者可以参考其文档。我们只会以简单的例子展示如何使用 Boost.Python 库包装 C/C++ MPI 程序以供 mpi4py 调用。

假设我们有以下使用了 Boost.Python 的 C++ 程序文件 helloworld.cxx,其中定义了函数 sayhello,其接受一个 MPI 通信子作为参数。

// helloworld.cxx#include <mpi.h>
#include <iostream>static void sayhello(MPI_Comm comm)
{if (comm == MPI_COMM_NULL) {std::cout << "You passed MPI_COMM_NULL !!!" << std::endl;return;}int size;MPI_Comm_size(comm, &size);int rank;MPI_Comm_rank(comm, &rank);int plen; char pname[MPI_MAX_PROCESSOR_NAME];MPI_Get_processor_name(pname, &plen);std::cout <<"Hello, World! " <<"I am process "  << rank  <<" of "           << size  <<" on  "          << pname <<"."              << std::endl;
}#include <boost/python.hpp>
#include <mpi4py/mpi4py.h>
using namespace boost::python;static void hw_sayhello(object py_comm)
{PyObject* py_obj = py_comm.ptr();MPI_Comm *comm_p = PyMPIComm_Get(py_obj);if (comm_p == NULL) throw_error_already_set();sayhello(*comm_p);
}BOOST_PYTHON_MODULE(helloworld)
{if (import_mpi4py() < 0) return; /* Python 2.X */def("sayhello", hw_sayhello);
}/** Local Variables:* mode: C++* End:*/

我们要将其编译成扩展模块 helloworld.so 以供 mpi4py 程序调用,可以使用类似于以下的编译命令(注意将其中的头文件路径改成你的系统中实际的路径):

$ mpicxx -I/path/to/python/include/python2.7  -I/path/to/python/lib/python2.7/site-packages/mpi4py/include -o helloworld.so helloworld.cxx -fPIC -shared -lpthread -ldl -lutil -lm -lpython2.7 -lboost_python -lboost_python-mt

这编译命令与前两篇中使用的编译命令非常类似,不同的是,对 C++ 程序,这里使用的是 mpicxx 编译器,另外除了指定 Python.h 和 mpi4py.h 的头文件路径及链接 python2.7 等库文件外,还需额外链接 boost_python 和 boost_python-mt 库文件,因为这里我们是使用 Boost.Python 库来包装的。

编译成功后会生成扩展模块 helloworld.so,然后就可以在我们的 mpi4py 程序中像使用其它 Python 模块一样导入该模块并调用该模块中定义的 sayhello 函数,可以向此函数传递一个 mpi4py 中定义的通信子,如 MPI.COMM_WORLD 或者其它通信子对象。例如使用以下的 test 例程:

# test.pyfrom mpi4py import MPI
import helloworld as hwnull = MPI.COMM_NULL
hw.sayhello(null)comm = MPI.COMM_WORLD
hw.sayhello(comm)try:hw.sayhello(list())
except:pass
else:assert 0, "exception not raised"

执行结果如下:

 

$ mpiexec -n 4 python test.py
You passed MPI_COMM_NULL !!!
Hello, World! I am process 1 of 4 on  node2.
You passed MPI_COMM_NULL !!!
Hello, World! I am process 2 of 4 on  node2.
You passed MPI_COMM_NULL !!!
Hello, World! I am process 3 of 4 on  node2.
You passed MPI_COMM_NULL !!!
Hello, World! I am process 0 of 4 on  node2.

为了方便,我们也可以编写如下 Makefile 以简化上述操作(注意其中使用了在上一篇中介绍的 python-config 文件):

 

# Makefile.PHONY: default
default: build test cleanPYTHON = python
PYTHON_CONFIG = ${PYTHON} ./python-config
MPI4PY_INCLUDE = ${shell ${PYTHON} -c 'import mpi4py; print( mpi4py.get_include() )'}
BOOST_INCS =
BOOST_LIBS = -lboost_python -lboost_python-mtMPICXX = mpicxx
CXXFLAGS = -fPIC ${shell ${PYTHON_CONFIG} --includes} ${BOOST_INCS}
LDFLAGS = -shared ${shell ${PYTHON_CONFIG} --libs} ${BOOST_LIBS}
SO = ${shell ${PYTHON_CONFIG} --extension-suffix}
.PHONY: build
build: helloworld${SO}
helloworld${SO}: helloworld.cxx${MPICXX} ${CXXFLAGS} -I${MPI4PY_INCLUDE} -o $@ $< ${LDFLAGS}MPIEXEC = mpiexec
NP_FLAG = -n
NP = 5
.PHONY: test
test: build${MPIEXEC} ${NP_FLAG} ${NP} ${PYTHON} test.py.PHONY: clean
clean:${RM} helloworld${SO}

编译扩展库,执行程序及清理可以分别使用如下命令:

 

$ make build
$ make test
$ make clean

 

ps:虽然boost.Python据说扩展性不错,但是鉴于我目前的水平,还是最终选择swig来实现吧,我看openmm其实也使用了swig,毕竟以后跟它一起走,那就和(tou)谐(lan)一点哈哈哈。

 

感谢原文博主!

这篇关于C++扩展python(三)boost.Python实现(附swig比较)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/541618

相关文章

C++中unordered_set哈希集合的实现

《C++中unordered_set哈希集合的实现》std::unordered_set是C++标准库中的无序关联容器,基于哈希表实现,具有元素唯一性和无序性特点,本文就来详细的介绍一下unorder... 目录一、概述二、头文件与命名空间三、常用方法与示例1. 构造与析构2. 迭代器与遍历3. 容量相关4

C++中悬垂引用(Dangling Reference) 的实现

《C++中悬垂引用(DanglingReference)的实现》C++中的悬垂引用指引用绑定的对象被销毁后引用仍存在的情况,会导致访问无效内存,下面就来详细的介绍一下产生的原因以及如何避免,感兴趣... 目录悬垂引用的产生原因1. 引用绑定到局部变量,变量超出作用域后销毁2. 引用绑定到动态分配的对象,对象

SpringBoot基于注解实现数据库字段回填的完整方案

《SpringBoot基于注解实现数据库字段回填的完整方案》这篇文章主要为大家详细介绍了SpringBoot如何基于注解实现数据库字段回填的相关方法,文中的示例代码讲解详细,感兴趣的小伙伴可以了解... 目录数据库表pom.XMLRelationFieldRelationFieldMapping基础的一些代

Java HashMap的底层实现原理深度解析

《JavaHashMap的底层实现原理深度解析》HashMap基于数组+链表+红黑树结构,通过哈希算法和扩容机制优化性能,负载因子与树化阈值平衡效率,是Java开发必备的高效数据结构,本文给大家介绍... 目录一、概述:HashMap的宏观结构二、核心数据结构解析1. 数组(桶数组)2. 链表节点(Node

Java AOP面向切面编程的概念和实现方式

《JavaAOP面向切面编程的概念和实现方式》AOP是面向切面编程,通过动态代理将横切关注点(如日志、事务)与核心业务逻辑分离,提升代码复用性和可维护性,本文给大家介绍JavaAOP面向切面编程的概... 目录一、AOP 是什么?二、AOP 的核心概念与实现方式核心概念实现方式三、Spring AOP 的关

Python版本信息获取方法详解与实战

《Python版本信息获取方法详解与实战》在Python开发中,获取Python版本号是调试、兼容性检查和版本控制的重要基础操作,本文详细介绍了如何使用sys和platform模块获取Python的主... 目录1. python版本号获取基础2. 使用sys模块获取版本信息2.1 sys模块概述2.1.1

一文详解Python如何开发游戏

《一文详解Python如何开发游戏》Python是一种非常流行的编程语言,也可以用来开发游戏模组,:本文主要介绍Python如何开发游戏的相关资料,文中通过代码介绍的非常详细,需要的朋友可以参考下... 目录一、python简介二、Python 开发 2D 游戏的优劣势优势缺点三、Python 开发 3D

Python函数作用域与闭包举例深度解析

《Python函数作用域与闭包举例深度解析》Python函数的作用域规则和闭包是编程中的关键概念,它们决定了变量的访问和生命周期,:本文主要介绍Python函数作用域与闭包的相关资料,文中通过代码... 目录1. 基础作用域访问示例1:访问全局变量示例2:访问外层函数变量2. 闭包基础示例3:简单闭包示例4

Python实现字典转字符串的五种方法

《Python实现字典转字符串的五种方法》本文介绍了在Python中如何将字典数据结构转换为字符串格式的多种方法,首先可以通过内置的str()函数进行简单转换;其次利用ison.dumps()函数能够... 目录1、使用json模块的dumps方法:2、使用str方法:3、使用循环和字符串拼接:4、使用字符

Python版本与package版本兼容性检查方法总结

《Python版本与package版本兼容性检查方法总结》:本文主要介绍Python版本与package版本兼容性检查方法的相关资料,文中提供四种检查方法,分别是pip查询、conda管理、PyP... 目录引言为什么会出现兼容性问题方法一:用 pip 官方命令查询可用版本方法二:conda 管理包环境方法