From 5d2ef4ffd05593f91f9b00546be76de4b8075ae3 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Thu, 21 Mar 2019 16:49:34 +0100 Subject: ndarray: add new API for n-dimensional arrays --- ndarray.h | 55 +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 55 insertions(+) create mode 100644 ndarray.h (limited to 'ndarray.h') diff --git a/ndarray.h b/ndarray.h new file mode 100644 index 0000000..e264fba --- /dev/null +++ b/ndarray.h @@ -0,0 +1,55 @@ +/* + * N-dimensional strided array + * Copyright 2018 Anton Khirnov + * + * This program is free software: you can redistribute it and/or modify + * it under the terms of the GNU General Public License as published by + * the Free Software Foundation, either version 3 of the License, or + * (at your option) any later version. + * + * This program is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + * GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with this program. If not, see . + */ + +#ifndef MG2D_NDARRAY_H +#define MG2D_NDARRAY_H + +#include +#include + +typedef struct Slice { + ptrdiff_t start; + ptrdiff_t end; + size_t step; +} Slice; + +#define SLICE(start, end, step) ((Slice){ start, end, step }) +#define SLICE_NULL ((Slice){ PTRDIFF_MAX }) + +typedef struct NDArrayInternal NDArrayInternal; + +typedef struct NDArray { + NDArrayInternal *priv; + + unsigned int dims; + + double *data; + const size_t *shape; + const ptrdiff_t *stride; +} NDArray; + +#define NDARRAY_ALLOC_ZERO (1 << 0) + +int mg2di_ndarray_alloc(NDArray **result, unsigned int dims, + const size_t * const size, unsigned int flags); +void mg2di_ndarray_free(NDArray **a); + +int mg2di_ndarray_slice(NDArray **result, NDArray *src, + const Slice * const slice); + +#endif // MG2D_ARRAY_H -- cgit v1.2.3