1 /****************************************************************************** 2 * 3 * Copyright 2022 Google LLC 4 * 5 * Licensed under the Apache License, Version 2.0 (the "License"); 6 * you may not use this file except in compliance with the License. 7 * You may obtain a copy of the License at: 8 * 9 * http://www.apache.org/licenses/LICENSE-2.0 10 * 11 * Unless required by applicable law or agreed to in writing, software 12 * distributed under the License is distributed on an "AS IS" BASIS, 13 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 14 * See the License for the specific language governing permissions and 15 * limitations under the License. 16 * 17 ******************************************************************************/ 18 19 #include <Python.h> 20 #include <numpy/ndarrayobject.h> 21 22 #include <mdct.c> 23 #include "ctypes.h" 24 25 26 static PyObject *mdct_forward_py(PyObject *m, PyObject *args) 27 { 28 PyObject *x_obj, *xd_obj, *y_obj, *d_obj; 29 enum lc3_dt dt; 30 enum lc3_srate sr; 31 float *x, *xd, *y, *d; 32 33 if (!PyArg_ParseTuple(args, "iiOO", &dt, &sr, &x_obj, &xd_obj)) 34 return NULL; 35 36 CTYPES_CHECK("dt", (unsigned)dt < LC3_NUM_DT); 37 CTYPES_CHECK("sr", (unsigned)sr < LC3_NUM_SRATE); 38 39 int ns = LC3_NS(dt, sr), nd = LC3_ND(dt, sr); 40 41 CTYPES_CHECK("x", to_1d_ptr(x_obj, NPY_FLOAT, ns, &x)); 42 CTYPES_CHECK("xd", to_1d_ptr(xd_obj, NPY_FLOAT, nd, &xd)); 43 d_obj = new_1d_ptr(NPY_FLOAT, nd, &d); 44 y_obj = new_1d_ptr(NPY_FLOAT, ns, &y); 45 46 memcpy(d, xd, nd * sizeof(float)); 47 48 lc3_mdct_forward(dt, sr, sr, x, d, y); 49 50 return Py_BuildValue("NN", y_obj, d_obj); 51 } 52 53 static PyObject *mdct_inverse_py(PyObject *m, PyObject *args) 54 { 55 PyObject *x_obj, *xd_obj, *d_obj, *y_obj; 56 enum lc3_dt dt; 57 enum lc3_srate sr; 58 float *x, *xd, *d, *y; 59 60 if (!PyArg_ParseTuple(args, "iiOO", &dt, &sr, &x_obj, &xd_obj)) 61 return NULL; 62 63 CTYPES_CHECK("dt", (unsigned)dt < LC3_NUM_DT); 64 CTYPES_CHECK("sr", (unsigned)sr < LC3_NUM_SRATE); 65 66 int ns = LC3_NS(dt, sr), nd = LC3_ND(dt, sr); 67 68 CTYPES_CHECK("x", to_1d_ptr(x_obj, NPY_FLOAT, ns, &x)); 69 CTYPES_CHECK("xd", to_1d_ptr(xd_obj, NPY_FLOAT, nd, &xd)); 70 d_obj = new_1d_ptr(NPY_FLOAT, nd, &d); 71 y_obj = new_1d_ptr(NPY_FLOAT, ns, &y); 72 73 memcpy(d, xd, nd * sizeof(float)); 74 75 lc3_mdct_inverse(dt, sr, sr, x, d, y); 76 77 return Py_BuildValue("NN", y_obj, d_obj); 78 } 79 80 static PyMethodDef methods[] = { 81 { "mdct_forward", mdct_forward_py, METH_VARARGS }, 82 { "mdct_inverse", mdct_inverse_py, METH_VARARGS }, 83 { NULL }, 84 }; 85 86 PyMODINIT_FUNC lc3_mdct_py_init(PyObject *m) 87 { 88 import_array(); 89 90 PyModule_AddFunctions(m, methods); 91 92 return m; 93 } 94