npy_intp sz = npy_intp(n), sz0 = 0, stride1 = sizeof(double);
PyObject *xpy = PyArray_New(&PyArray_Type, 1, &sz, NPY_DOUBLE, &stride1,
const_cast<double*>(x), // not NPY_WRITEABLE
- 0, NPY_ARRAY_IN_ARRAY, NULL);
+ 0, NPY_ARRAY_C_CONTIGUOUS | NPY_ARRAY_ALIGNED, NULL);
PyObject *gradpy = grad
? PyArray_SimpleNewFromData(1, &sz, NPY_DOUBLE, grad)
: PyArray_SimpleNew(1, &sz0, NPY_DOUBLE);
npy_intp sz0 = 0, stride1 = sizeof(double);
PyObject *xpy = PyArray_New(&PyArray_Type, 1, &nsz, NPY_DOUBLE, &stride1,
const_cast<double*>(x), // not NPY_WRITEABLE
- 0, NPY_ARRAY_IN_ARRAY, NULL);
+ 0, NPY_ARRAY_C_CONTIGUOUS | NPY_ARRAY_ALIGNED, NULL);
PyObject *rpy = PyArray_SimpleNewFromData(1, &msz, NPY_DOUBLE, result);
PyObject *gradpy = grad
? PyArray_SimpleNewFromData(2, mnsz, NPY_DOUBLE, grad)