pybind11和numpy进行交互的方法


Posted in Python onJuly 04, 2019

使用一个遵循buffer protocol的对象就可以和numpy交互了.

这个buffer_protocol要有哪些东西呢? 要有如下接口:

struct buffer_info {
  void *ptr;
  ssize_t itemsize;
  std::string format;
  ssize_t ndim;
  std::vector<ssize_t> shape;
  std::vector<ssize_t> strides;
};

其实就是一个指向数组的指针+各个维度的信息就可以了. 然后我们就可以用指针+偏移来访问数字中的任意位置上的数字了.

下面是一个可以跑的例子:

#include <pybind11/pybind11.h>
 #include <pybind11/numpy.h>
 namespace py = pybind11;
 py::array_t<double> add_arrays(py::array_t<double> input1, py::array_t<double> input2) {
   py::buffer_info buf1 = input1.request(), buf2 = input2.request();
   if (buf1.ndim != 1 || buf2.ndim != 1)
     throw std::runtime_error("Number of dimensions must be one");
   if (buf1.size != buf2.size)
     throw std::runtime_error("Input shapes must match");
   /* No pointer is passed, so NumPy will allocate the buffer */
   auto result = py::array_t<double>(buf1.size);
   py::buffer_info buf3 = result.request();
   double *ptr1 = (double *) buf1.ptr,
      *ptr2 = (double *) buf2.ptr,
      *ptr3 = (double *) buf3.ptr;
   for (size_t idx = 0; idx < buf1.shape[0]; idx++)
     ptr3[idx] = ptr1[idx] + ptr2[idx];
   return result;
 }
 
 PYBIND11_MODULE(test, m) {
   m.def("add_arrays", &add_arrays, "Add two NumPy arrays");
 }

array_t里的buf就是一个兼容的接口.

buf中可以得到指针和对应数字的维度信息.

为了方便我们甚至可以使用Eigen当作我们兼容numpy的接口:

#include <pybind11/pybind11.h>
 #include <pybind11/eigen.h> 
 #include <Eigen/LU> 
 // N.B. this would equally work with Eigen-types that are not predefined. For example replacing
 // all occurrences of "Eigen::MatrixXd" with "MatD", with the following definition:
 //
 // typedef Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> MatD;
 
 Eigen::MatrixXd inv(const Eigen::MatrixXd &xs)
 {
  return xs.inverse();
 }
 
 double det(const Eigen::MatrixXd &xs)
 {
  return xs.determinant();
 }
 
 namespace py = pybind11;
 
 PYBIND11_MODULE(example,m)
 {
  m.doc() = "pybind11 example plugin";
 
  m.def("inv", &inv);
 
  m.def("det", &det);
 }

更多参考:

https://pybind11.readthedocs.io/en/stable/advanced/pycpp/numpy.html

https://github.com/tdegeus/pybind11_examples

总结

以上所述是小编给大家介绍的pybind11和numpy进行交互的方法,希望对大家有所帮助,如果大家有任何疑问请给我留言,小编会及时回复大家的。在此也非常感谢大家对三水点靠木网站的支持!
如果你觉得本文对你有帮助,欢迎转载,烦请注明出处,谢谢!

Python 相关文章推荐
Python中条件选择和循环语句使用方法介绍
Mar 13 Python
Python3.0与2.X版本的区别实例分析
Aug 25 Python
总结网络IO模型与select模型的Python实例讲解
Jun 27 Python
Python爬虫爬取美剧网站的实现代码
Sep 03 Python
Python编程实现生成特定范围内不重复多个随机数的2种方法
Apr 14 Python
python requests 测试代理ip是否生效
Jul 25 Python
scrapy-redis源码分析之发送POST请求详解
May 15 Python
Python爬虫实现vip电影下载的示例代码
Apr 20 Python
如何在scrapy中捕获并处理各种异常
Sep 28 Python
Python生成pdf目录书签的实例方法
Oct 29 Python
python 实时调取摄像头的示例代码
Nov 25 Python
Python的代理类实现,控制访问和修改属性的权限你都了解吗
Mar 21 Python
pandas计算最大连续间隔的方法
Jul 04 #Python
python SQLAlchemy 中的Engine详解
Jul 04 #Python
Python Pandas实现数据分组求平均值并填充nan的示例
Jul 04 #Python
pybind11在Windows下的使用教程
Jul 04 #Python
Pandas_cum累积计算和rolling滚动计算的用法详解
Jul 04 #Python
Pycharm简单使用教程(入门小结)
Jul 04 #Python
pandas基于时间序列的固定时间间隔求均值的方法
Jul 04 #Python
You might like
php中一个有意思的日期逻辑处理
2012/03/25 PHP
PHP格式化显示时间date()函数代码
2018/10/03 PHP
phpinfo的知识点总结
2019/10/10 PHP
PHP设计模式之迭代器模式Iterator实例分析【对象行为型】
2020/04/26 PHP
关于实现代码语法标亮 dp.SyntaxHighlighter
2007/02/02 Javascript
可以支持多中格式的JS键盘
2007/05/02 Javascript
DOM下的节点属性和操作小结
2009/05/14 Javascript
新鲜出炉的js tips提示效果
2011/04/03 Javascript
JS实现往下不断流动网页背景的方法
2015/02/27 Javascript
Javascript闭包与函数柯里化浅析
2016/06/22 Javascript
AngularJS 与Bootstrap实现表格分页实例代码
2016/10/14 Javascript
详解Js中的模块化是如何实现的
2017/10/18 Javascript
Vue作用域插槽slot-scope实例代码
2018/09/05 Javascript
实例详解带参数的 npm script
2019/05/28 Javascript
微信小程序实现pdf、word等格式文件上传的方法
2019/09/10 Javascript
js正则匹配多个全部数据问题
2019/12/20 Javascript
安装多版本Vue-CLI的实现方法
2020/03/24 Javascript
vue 组件简介
2020/07/31 Javascript
[03:58]兄弟们,回来开黑了!DOTA2昔日战友招募宣传视频
2016/07/17 DOTA
python调用java模块SmartXLS和jpype修改excel文件的方法
2015/04/28 Python
在Python中操作列表之List.pop()方法的使用
2015/05/21 Python
python之当你发现QTimer不能用时的解决方法
2019/06/21 Python
python selenium爬取斗鱼所有直播房间信息过程详解
2019/08/09 Python
Python Sphinx使用实例及问题解决
2020/01/17 Python
Python使用PyQt5/PySide2编写一个极简的音乐播放器功能
2020/02/07 Python
Selenium 滚动页面至元素可见的方法
2020/03/18 Python
基于python获取本地时间并转换时间戳和日期格式
2020/10/27 Python
在Pycharm中安装Pandas库方法(简单易懂)
2021/02/20 Python
HTML5新增属性data-*和js/jquery之间的交互及注意事项
2017/08/08 HTML / CSS
阿玛尼意大利官网:Armani意大利
2018/10/30 全球购物
毕业生护理专业个人求职信范文
2014/01/04 职场文书
《最佳路径》教学反思
2014/04/13 职场文书
团结演讲稿范文
2014/05/23 职场文书
领导班子“四风问题”“整改方案
2014/10/02 职场文书
2014会计年终工作总结
2014/12/20 职场文书
Android自定义双向滑动控件
2022/04/19 Java/Android