2

如果这个问题不完整或难以理解，我很抱歉，我正在尝试自己解决！我有以下 C 代码，在 python 中，我不太擅长这两种语言。

以下 C 代码复制了我认为的 numpy 数组。当我使用它时它不起作用，它使用 SWITCH 语句并打印“默认”答案。有人可以帮忙吗？

是python的版本吗？它似乎适用于旧版本的python。 

  static PyObject *llpy_acopy(PyObject *unused, PyObject *args)
  {
  PyObject *afrom = PyTuple_GET_ITEM(args, 0);
  PyObject *ato = PyTuple_GET_ITEM(args, 1);
  PyArrayObject *tmp = (PyArrayObject *)afrom;
  int nd = tmp->nd;
  int n = 1;
  int i;
  for (i=0; i<nd; i++)
    n *= tmp->dimensions[i];

  switch(tmp->descr->type_num) {
  case PyArray_FLOAT:
      memcpy(((PyArrayObject *)ato)->data, tmp->data,
              (size_t)(n * sizeof(float)));
      break;
  case PyArray_INT:
      memcpy(IDATA(ato), IDATA(afrom), 
           (size_t )(n  * sizeof(int)));
      break;
   case PyArray_CHAR:
      memcpy(CDATA(ato), CDATA(afrom), 
             (size_t )(n * sizeof(char)));
      break;
   default:
      PyErr_SetString(PyExc_TypeError, "Cannot copy array now");
      return NULL;
    }
   Py_INCREF(Py_None);
  return Py_None;
      }

从一个复制到另一个的两个数组是 - 即第一个数组被复制到第二个数组中，替换现有值。

[[ 1.          0.          0.        ]
[ 1.51999998  1.42799997  0.        ]
[ 0.73699999  2.37800002  0.        ]
  ..., 
[-7.55183792 -9.31523228 -0.19252452]
[-6.7193203  -8.40346909  0.58392692]
[-8.47850037 -9.02481842 -0.46105781]]


[[ 1.          0.          0.        ]
 [ 1.51999998  1.42799997  0.        ]
 [ 0.73699999  2.37800002  0.        ]
  ..., 
 [-5.86299992 -9.36299992 -5.0999999 ]
 [-5.61800003 -7.94799995 -4.84499979]
 [-6.12400007 -9.65600014 -6.02799988]]
4

1 回答 1

0

该表达式tmp->descr->type_num将返回一个整数（例如，双精度结果 type_num=12）。但我不知道PyArray_INT, PyArray_FLOAT,PyArray_CHAR等。从 Google 搜索中，这些看起来像旧 numeric_numarray.h 标头中定义的类型。我认为您想要做的是使用此处列出的数据类型检查。参见PyTypeNum_ISFLOAT(num),PyTypeNum_ISINTEGER(num)等。

于 2013-10-08T00:53:50.573 回答