change function to create array

fftma_module/gen/lib_src/clean_real.c

@@ -15,13 +15,13 @@ void clean_real(struct realization_mod* realin, int n[3], struct grid_mod grid,
     /*is the output realization already allocated?*/
     /*if not, memory allocation*/
 
-    if (realout->vector == NULL || realout->n != realin->n) {
+    /*if (realout->vector == NULL || realout->n != realin->n) {
         realout->vector = (double*)malloc(realin->n * sizeof(double));
         if (realout->vector == NULL) {
             printf("Clean_real.c:  No memory available\n");
             exit;
         }
-    }
+    }*/
 
     realout->n = realin->n;
     realout->code = 1;

fftma_module/gen/moduleFFTMA.c

@@ -28,9 +28,12 @@ static PyObject* genFunc(PyObject* self, PyObject* args)
     struct vario_mod variogram;
     long seed;
     struct statistic_mod stat;
-    PyObject* out_array;
+    PyArrayObject* out_array;
     npy_intp out_dims[NPY_MAXDIMS];
 
+    out_array = (PyArrayObject *) PyArray_FromDims(NDIMENSIONS,out_dims,NPY_DOUBLE);
+    Y.vector = (double *) out_array->data;
+
     if (!Py_getvalues(args, &seed, &grid, &variogram, &stat))
         return NULL;
     Py_kgeneration(seed, grid, stat, variogram, &Z, &Y, n);
@@ -39,11 +42,12 @@ static PyObject* genFunc(PyObject* self, PyObject* args)
     out_dims[1] = grid.NY;
     out_dims[2] = grid.NX;
 
-    out_array = PyArray_SimpleNewFromData(NDIMENSIONS, out_dims, NPY_DOUBLE, Y.vector);
+    out_array = (PyArrayObject*)PyArray_SimpleNewFromData(NDIMENSIONS, out_dims, NPY_DOUBLE, Y.vector);
     if (out_array == NULL)
         return NULL;
+    // PyArray_ENABLEFLAGS(out_array, NPY_ARRAY_OWNDATA);
 
-    PyArray_ENABLEFLAGS(out_array, NPY_ARRAY_OWNDATA);
+    // out_array->flags |= NPY_ARRAY_OWNDATA;
 
     free(stat.mean);
     free(stat.variance);