ceps.cpp 20.9 KB
Newer Older
André Anjos's avatar
André Anjos committed
1 2
/**
 * @author Andre Anjos <andre.anjos@idiap.ch>
3
 * @author Pavel Korshunov <Pavel.Korshunov@idiap.ch>
André Anjos's avatar
André Anjos committed
4 5 6 7 8
 * @date Thu  6 Feb 09:00:05 2014
 *
 * @brief Bindings to the base class bob::ap::Ceps
 */

André Anjos's avatar
André Anjos committed
9 10
#include <bob.blitz/cppapi.h>
#include <bob.blitz/cleanup.h>
11
#include <bob.extension/defines.h>
André Anjos's avatar
André Anjos committed
12 13
#include "types.h"

André Anjos's avatar
André Anjos committed
14
PyDoc_STRVAR(s_ceps_str, BOB_EXT_MODULE_PREFIX ".Ceps");
André Anjos's avatar
André Anjos committed
15 16

PyDoc_STRVAR(s_ceps_doc,
17
"Ceps(sampling_frequency, [win_length_ms=20., [win_shift_ms=10., [n_filters=24, [n_ceps=19, [f_min=0., [f_max=4000., [delta_win=2, [pre_emphasis_coeff=0.95, [mel_scale=True, [dct_norm=True, [normalize_mean=True, [rect_filter=False, [inverse_filter=False, [normalize_spectrum=False, [ssfc_features=False, [scfc_features=False, [scmc_features=False]]]]]]]]]]]]]]]]]) -> new Ceps\n\
André Anjos's avatar
André Anjos committed
18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Ceps(other) -> new Ceps\n\
\n\
Objects of this class, after configuration, can extract the\n\
cepstral coefficients from 1D audio array/signals.\n\
\n\
Parameters:\n\
\n\
sampling_frequency\n\
  [float] the sampling frequency/frequency rate\n\
\n\
win_length_ms\n\
  [float] the window length in miliseconds\n\
\n\
win_shift_ms\n\
  [float] the window shift in miliseconds\n\
\n\
n_filters\n\
  [int] the number of filter bands\n\
\n\
n_ceps\n\
  [int] the number of cepstral coefficients\n\
\n\
f_min\n\
  [double] the minimum frequency of the filter bank\n\
\n\
f_max\n\
  [double] the maximum frequency of the filter bank\n\
\n\
delta_win\n\
  [int] The integer delta value used for computing the\n\
  first and second order derivatives\n\
\n\
pre_emphasis_coeff\n\
  [double] the coefficient used for the pre-emphasis\n\
\n\
mel_scale\n\
  [bool] tells whether cepstral features are extracted\n\
  on a linear (LFCC, set it to ``False``) or Mel (MFCC,\n\
  set it to ``True`` - the default)\n\
\n\
dct_norm\n\
  [bool] A factor by which the cepstral coefficients are\n\
  multiplied\n\
61 62 63 64
normalize_mean\n\
  [bool] Tells whether frame should be normalized \n\
  by subtracting mean (True) or dividing by max_range (False)\n\
  ``True`` is the default value.\n\
65 66 67 68 69 70 71 72 73 74 75 76 77 78
\n\
rect_filter\n\
  [bool] tells whether to apply the filter in the\n\
  inversed order, i.e., from high frequencies to low\n\
  (set it to ``True''). ``False`` is the default value.\n\
\n\
inverse_filter\n\
  [bool] tells whether cepstral features are extracted\n\
  using a rectungular filter (set it to ``True``), i.e., RFCC features,\n\
  instead of the default filter (the default value is ``False``)\n\
\n\
normalize_spectrum\n\
  [bool] Tells whether to normalize the power spectrum of the signal.\n\
  The default value is ``False``.\n\
André Anjos's avatar
André Anjos committed
79
\n\
80 81 82 83 84 85 86 87 88 89 90 91 92 93 94
ssfc_features\n\
  [bool] Set to true if you want to compute\n\
  Subband Spectral Flux Coefficients (SSFC), which measures\n\
  the frame-by-frame change in the power spectrum\n\
\n\
scfc_features\n\
  [bool] Set to true if you want to compute\n\
  Spectral Centroid Frequency Coefficients (SCFC), which\n\
  capture detailed information about subbands similar to formant frequencies\n\
\n\
scmc_features\n\
  [bool] Set to true if you want to compute\n\
  Spectral Centroid Magnitude  Coefficients (SCMC), which\n\
  capture detailed information about subbands similar to SCFC features\n\
\n\
André Anjos's avatar
André Anjos committed
95 96 97 98 99 100 101 102 103 104 105 106 107
other\n\
  [Ceps] an object of which is or inherits from ``Ceps``\n\
  that will be deep-copied into a new instance.\n\
\n\
"
);

int PyBobApCeps_Check(PyObject* o) {
  return PyObject_IsInstance(o, reinterpret_cast<PyObject*>(&PyBobApCeps_Type));
}

static void PyBobApCeps_Delete (PyBobApCepsObject* o) {

108 109 110
  o->parent.parent.parent.cxx = 0; // FrameExtractor
  o->parent.parent.cxx = 0; // Energy
  o->parent.cxx = 0; // Spectrogram
André Anjos's avatar
André Anjos committed
111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161
  delete o->cxx;
  Py_TYPE(o)->tp_free((PyObject*)o);

}

static int PyBobApCeps_InitCopy
(PyBobApCepsObject* self, PyObject* args, PyObject* kwds) {

  /* Parses input arguments in a single shot */
  static const char* const_kwlist[] = {"other", 0};
  static char** kwlist = const_cast<char**>(const_kwlist);

  PyObject* other = 0;

  if (!PyArg_ParseTupleAndKeywords(args, kwds, "O!", kwlist,
        &PyBobApCeps_Type, &other)) return -1;

  auto copy = reinterpret_cast<PyBobApCepsObject*>(other);

  try {
    self->cxx = new bob::ap::Ceps(*(copy->cxx));
    if (!self->cxx) {
      PyErr_Format(PyExc_MemoryError, "cannot create new object of type `%s' - no more memory", Py_TYPE(self)->tp_name);
      return -1;
    }
    self->parent.parent.parent.cxx = self->cxx;
    self->parent.parent.cxx = self->cxx;
    self->parent.cxx = self->cxx;
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return -1;
  }
  catch (...) {
    PyErr_Format(PyExc_RuntimeError, "cannot create new object of type `%s' - unknown exception thrown", Py_TYPE(self)->tp_name);
    return -1;
  }

  return 0;

}

static int PyBobApCeps_InitParameters
(PyBobApCepsObject* self, PyObject *args, PyObject* kwds) {

  /* Parses input arguments in a single shot */
  static const char* const_kwlist[] = {
    "sampling_frequency",
    "win_length_ms",
    "win_shift_ms",
    "n_filters",
162
    "n_ceps",
André Anjos's avatar
André Anjos committed
163 164
    "f_min",
    "f_max",
165
    "delta_win",
André Anjos's avatar
André Anjos committed
166 167
    "pre_emphasis_coeff",
    "mel_scale",
168 169
    "dct_norm",
    "normalize_mean",
170 171 172 173 174 175
    "rect_filter",
    "inverse_filter",
    "normalize_spectrum",
    "ssfc_features",
    "scfc_features",
    "scmc_features",
André Anjos's avatar
André Anjos committed
176 177 178 179 180 181 182
    0};
  static char** kwlist = const_cast<char**>(const_kwlist);

  double sampling_frequency = 0.;
  double win_length_ms = 20.;
  double win_shift_ms = 10.;
  Py_ssize_t n_filters = 24;
183
  Py_ssize_t n_ceps = 19;
André Anjos's avatar
André Anjos committed
184
  double f_min = 0.;
185
  double f_max = 8000.;
186 187
  Py_ssize_t delta_win = 2;
  double pre_emphasis_coeff = 0.95;
André Anjos's avatar
André Anjos committed
188
  PyObject* mel_scale = Py_True;
189 190
  PyObject* dct_norm = Py_False;
  PyObject* normalize_mean = Py_True;
191 192 193 194 195 196
  PyObject* rect_filter = Py_False;
  PyObject* inverse_filter = Py_False;
  PyObject* normalize_spectrum = Py_False;
  PyObject* ssfc_features = Py_False;
  PyObject* scfc_features = Py_False;
  PyObject* scmc_features = Py_False;
197 198
  if (!PyArg_ParseTupleAndKeywords(args, kwds, "d|ddnnddndOOOOOOOOO", kwlist,
        &sampling_frequency, &win_length_ms, &win_shift_ms, &n_filters,
199
        &n_ceps, &f_min, &f_max, &delta_win, &pre_emphasis_coeff,
200 201
        &mel_scale, &dct_norm, &normalize_mean, &rect_filter,
        &inverse_filter, &normalize_spectrum,
202
        &ssfc_features, &scfc_features, &scmc_features))
André Anjos's avatar
André Anjos committed
203 204 205
    return -1;

  bool mel_scale_ = PyObject_IsTrue(mel_scale);
206 207
  bool dct_norm_ = PyObject_IsTrue(dct_norm);
  bool normalize_mean_ = PyObject_IsTrue(normalize_mean);
208 209 210 211 212 213
  bool rect_filter_ = PyObject_IsTrue(rect_filter);
  bool inverse_filter_ = PyObject_IsTrue(inverse_filter);
  bool normalize_spectrum_ = PyObject_IsTrue(normalize_spectrum);
  bool ssfc_features_ = PyObject_IsTrue(ssfc_features);
  bool scfc_features_ = PyObject_IsTrue(scfc_features);
  bool scmc_features_ = PyObject_IsTrue(scmc_features);
André Anjos's avatar
André Anjos committed
214 215 216

  try {
    self->cxx = new bob::ap::Ceps(sampling_frequency,
217 218 219 220
        win_length_ms, win_shift_ms, n_filters, n_ceps, f_min, f_max,
        delta_win, pre_emphasis_coeff, mel_scale_, dct_norm_, normalize_mean_,
        rect_filter_, inverse_filter_, normalize_spectrum_,
        ssfc_features_, scfc_features_, scmc_features_);
André Anjos's avatar
André Anjos committed
221 222 223 224
    if (!self->cxx) {
      PyErr_Format(PyExc_MemoryError, "cannot create new object of type `%s' - no more memory", Py_TYPE(self)->tp_name);
      return -1;
    }
225
    self->parent.parent.parent.cxx = self->cxx;
André Anjos's avatar
André Anjos committed
226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285
    self->parent.parent.cxx = self->cxx;
    self->parent.cxx = self->cxx;
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return -1;
  }
  catch (...) {
    PyErr_Format(PyExc_RuntimeError, "cannot create new object of type `%s' - unknown exception thrown", Py_TYPE(self)->tp_name);
    return -1;
  }

  return 0; ///< SUCCESS

}

static int PyBobApCeps_Init(PyBobApCepsObject* self,
    PyObject* args, PyObject* kwds) {

  Py_ssize_t nargs = (args?PyTuple_Size(args):0) + (kwds?PyDict_Size(kwds):0);

  switch (nargs) {

    case 1:

      {

        PyObject* arg = 0; ///< borrowed (don't delete)
        if (PyTuple_Size(args)) arg = PyTuple_GET_ITEM(args, 0);
        else {
          PyObject* tmp = PyDict_Values(kwds);
          auto tmp_ = make_safe(tmp);
          arg = PyList_GET_ITEM(tmp, 0);
        }

        if (PyBobApCeps_Check(arg)) {
          return PyBobApCeps_InitCopy(self, args, kwds);
        }

        else {
          return PyBobApCeps_InitParameters(self, args, kwds);
        }

        PyErr_Format(PyExc_TypeError, "cannot initialize `%s' with `%s' (see help)", Py_TYPE(self)->tp_name, Py_TYPE(arg)->tp_name);

      }

      break;

    default:

      return PyBobApCeps_InitParameters(self, args, kwds);

  }

  return -1;

}

static PyObject* PyBobApCeps_Repr(PyBobApCepsObject* self) {
286 287
  static const int MAXSIZE = 256;
  char buffer[MAXSIZE];
André Anjos's avatar
André Anjos committed
288
  Py_ssize_t n_filters = self->cxx->getNFilters();
289 290
  Py_ssize_t n_ceps = self->cxx->getNCeps();
  Py_ssize_t delta_win = self->cxx->getDeltaWin();
291
  auto count = std::snprintf(buffer, MAXSIZE, "%s(sampling_frequency=%f, win_length_ms=%f, win_shift_ms=%f, n_filters=%" PY_FORMAT_SIZE_T "d, n_ceps=%" PY_FORMAT_SIZE_T "d, f_min=%f, f_max=%f, delta_win=%" PY_FORMAT_SIZE_T "d, pre_emphasis_coeff=%f, mel_scale=%s, dct_norm=%s, normalize_mean=%s, rect_filter=%s, inverse_filter=%s, normalize_spectrum=%s, ssfc_features=%s, scfc_features=%s, scmc_features=%s)", Py_TYPE(self)->tp_name, self->cxx->getSamplingFrequency(), self->cxx->getWinLengthMs(), self->cxx->getWinShiftMs(), n_filters, n_ceps, self->cxx->getFMin(), self->cxx->getFMax(), delta_win, self->cxx->getPreEmphasisCoeff(), self->cxx->getMelScale()?"True":"False", self->cxx->getDctNorm()?"True":"False", self->cxx->getNormalizeMean()?"True":"False", self->cxx->getRectangularFilter()?"True":"False", self->cxx->getInverseFilter()?"True":"False", self->cxx->getNormalizeSpectrum()?"True":"False", self->cxx->getSSFCFeatures()?"True":"False", self->cxx->getSCFCFeatures()?"True":"False", self->cxx->getSCMCFeatures()?"True":"False");
André Anjos's avatar
André Anjos committed
292 293
  return
# if PY_VERSION_HEX >= 0x03000000
294
  PyUnicode_FromStringAndSize
André Anjos's avatar
André Anjos committed
295
# else
296
  PyString_FromStringAndSize
André Anjos's avatar
André Anjos committed
297
# endif
298
    (buffer, (count<=MAXSIZE)?count:MAXSIZE);
André Anjos's avatar
André Anjos committed
299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327
}

static PyObject* PyBobApCeps_RichCompare (PyBobApCepsObject* self,
    PyObject* other, int op) {

  if (!PyBobApCeps_Check(other)) {
    PyErr_Format(PyExc_TypeError, "cannot compare `%s' with `%s'",
        Py_TYPE(self)->tp_name, Py_TYPE(other)->tp_name);
    return 0;
  }

  auto other_ = reinterpret_cast<PyBobApCepsObject*>(other);

  switch (op) {
    case Py_EQ:
      if (self->cxx->operator==(*other_->cxx)) Py_RETURN_TRUE;
      Py_RETURN_FALSE;
      break;
    case Py_NE:
      if (self->cxx->operator!=(*other_->cxx)) Py_RETURN_TRUE;
      Py_RETURN_FALSE;
      break;
    default:
      Py_INCREF(Py_NotImplemented);
      return Py_NotImplemented;
  }

}

328 329 330
PyDoc_STRVAR(s_n_ceps_str, "n_ceps");
PyDoc_STRVAR(s_n_ceps_doc,
"The number of cepstral coefficients"
André Anjos's avatar
André Anjos committed
331 332
);

333
static PyObject* PyBobApCeps_GetNCeps
André Anjos's avatar
André Anjos committed
334
(PyBobApCepsObject* self, void* /*closure*/) {
335
  return Py_BuildValue("n", self->cxx->getNCeps());
André Anjos's avatar
André Anjos committed
336 337
}

338
static int PyBobApCeps_SetNCeps
André Anjos's avatar
André Anjos committed
339 340
(PyBobApCepsObject* self, PyObject* o, void* /*closure*/) {

341
  if (!PyBob_NumberCheck(o)) {
342
    PyErr_Format(PyExc_TypeError, "`%s' n_ceps can only be set using a number, not `%s'", Py_TYPE(self)->tp_name, Py_TYPE(o)->tp_name);
André Anjos's avatar
André Anjos committed
343 344 345 346 347 348 349
    return -1;
  }

  Py_ssize_t n = PyNumber_AsSsize_t(o, PyExc_OverflowError);
  if (PyErr_Occurred()) return -1;

  try {
350
    self->cxx->setNCeps(n);
André Anjos's avatar
André Anjos committed
351 352 353 354 355 356
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return -1;
  }
  catch (...) {
357
    PyErr_Format(PyExc_RuntimeError, "cannot reset `n_ceps' of %s: unknown exception caught", Py_TYPE(self)->tp_name);
André Anjos's avatar
André Anjos committed
358 359 360 361 362 363 364
    return -1;
  }

  return 0;

}

365 366 367 368
PyDoc_STRVAR(s_delta_win_str, "delta_win");
PyDoc_STRVAR(s_delta_win_doc,
"The integer delta value used for computing the first and\n\
second order derivatives"
André Anjos's avatar
André Anjos committed
369 370
);

371
static PyObject* PyBobApCeps_GetDeltaWin
André Anjos's avatar
André Anjos committed
372
(PyBobApCepsObject* self, void* /*closure*/) {
373
  return Py_BuildValue("n", self->cxx->getDeltaWin());
André Anjos's avatar
André Anjos committed
374 375
}

376
static int PyBobApCeps_SetDeltaWin
André Anjos's avatar
André Anjos committed
377 378
(PyBobApCepsObject* self, PyObject* o, void* /*closure*/) {

379
  if (!PyBob_NumberCheck(o)) {
380
    PyErr_Format(PyExc_TypeError, "`%s' delta_win can only be set using a number, not `%s'", Py_TYPE(self)->tp_name, Py_TYPE(o)->tp_name);
André Anjos's avatar
André Anjos committed
381 382 383
    return -1;
  }

384
  Py_ssize_t n = PyNumber_AsSsize_t(o, PyExc_OverflowError);
André Anjos's avatar
André Anjos committed
385 386 387
  if (PyErr_Occurred()) return -1;

  try {
388
    self->cxx->setDeltaWin(n);
André Anjos's avatar
André Anjos committed
389 390 391 392 393 394
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return -1;
  }
  catch (...) {
395
    PyErr_Format(PyExc_RuntimeError, "cannot reset `delta_win' of %s: unknown exception caught", Py_TYPE(self)->tp_name);
André Anjos's avatar
André Anjos committed
396 397 398 399 400 401 402
    return -1;
  }

  return 0;

}

403 404 405
PyDoc_STRVAR(s_dct_norm_str, "dct_norm");
PyDoc_STRVAR(s_dct_norm_doc,
"A factor by which the cepstral coefficients are multiplied"
André Anjos's avatar
André Anjos committed
406 407
);

408
static PyObject* PyBobApCeps_GetDctNorm
André Anjos's avatar
André Anjos committed
409
(PyBobApCepsObject* self, void* /*closure*/) {
410
  if (self->cxx->getDctNorm()) Py_RETURN_TRUE;
André Anjos's avatar
André Anjos committed
411 412 413
  else Py_RETURN_FALSE;
}

414
static int PyBobApCeps_SetDctNorm
André Anjos's avatar
André Anjos committed
415 416 417 418 419 420
(PyBobApCepsObject* self, PyObject* o, void* /*closure*/) {

  bool b = PyObject_IsTrue(o);
  if (PyErr_Occurred()) return -1;

  try {
421
    self->cxx->setDctNorm(b);
André Anjos's avatar
André Anjos committed
422 423 424 425 426 427
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return -1;
  }
  catch (...) {
428
    PyErr_Format(PyExc_RuntimeError, "cannot reset `dct_norm' of %s: unknown exception caught", Py_TYPE(self)->tp_name);
André Anjos's avatar
André Anjos committed
429 430 431 432 433 434 435
    return -1;
  }

  return 0;

}

436 437 438
PyDoc_STRVAR(s_with_energy_str, "with_energy");
PyDoc_STRVAR(s_with_energy_doc,
"Tells if we add the energy to the output feature"
André Anjos's avatar
André Anjos committed
439 440
);

441
static PyObject* PyBobApCeps_GetWithEnergy
André Anjos's avatar
André Anjos committed
442
(PyBobApCepsObject* self, void* /*closure*/) {
443
  if (self->cxx->getWithEnergy()) Py_RETURN_TRUE;
André Anjos's avatar
André Anjos committed
444 445 446
  else Py_RETURN_FALSE;
}

447
static int PyBobApCeps_SetWithEnergy
André Anjos's avatar
André Anjos committed
448 449 450 451 452 453
(PyBobApCepsObject* self, PyObject* o, void* /*closure*/) {

  bool b = PyObject_IsTrue(o);
  if (PyErr_Occurred()) return -1;

  try {
454
    self->cxx->setWithEnergy(b);
André Anjos's avatar
André Anjos committed
455 456 457 458 459 460
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return -1;
  }
  catch (...) {
461
    PyErr_Format(PyExc_RuntimeError, "cannot reset `with_energy' of %s: unknown exception caught", Py_TYPE(self)->tp_name);
André Anjos's avatar
André Anjos committed
462 463 464 465 466 467 468
    return -1;
  }

  return 0;

}

469 470 471
PyDoc_STRVAR(s_with_delta_str, "with_delta");
PyDoc_STRVAR(s_with_delta_doc,
"Tells if we add the first derivatives to the output feature"
André Anjos's avatar
André Anjos committed
472 473
);

474
static PyObject* PyBobApCeps_GetWithDelta
André Anjos's avatar
André Anjos committed
475
(PyBobApCepsObject* self, void* /*closure*/) {
476
  if (self->cxx->getWithDelta()) Py_RETURN_TRUE;
André Anjos's avatar
André Anjos committed
477 478 479
  else Py_RETURN_FALSE;
}

480
static int PyBobApCeps_SetWithDelta
André Anjos's avatar
André Anjos committed
481 482 483 484 485 486
(PyBobApCepsObject* self, PyObject* o, void* /*closure*/) {

  bool b = PyObject_IsTrue(o);
  if (PyErr_Occurred()) return -1;

  try {
487
    self->cxx->setWithDelta(b);
André Anjos's avatar
André Anjos committed
488 489 490 491 492 493
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return -1;
  }
  catch (...) {
494
    PyErr_Format(PyExc_RuntimeError, "cannot reset `with_delta' of %s: unknown exception caught", Py_TYPE(self)->tp_name);
André Anjos's avatar
André Anjos committed
495 496 497 498 499 500 501
    return -1;
  }

  return 0;

}

502 503 504
PyDoc_STRVAR(s_with_delta_delta_str, "with_delta_delta");
PyDoc_STRVAR(s_with_delta_delta_doc,
"Tells if we add the second derivatives to the output feature"
André Anjos's avatar
André Anjos committed
505 506
);

507
static PyObject* PyBobApCeps_GetWithDeltaDelta
André Anjos's avatar
André Anjos committed
508
(PyBobApCepsObject* self, void* /*closure*/) {
509
  if (self->cxx->getWithDeltaDelta()) Py_RETURN_TRUE;
André Anjos's avatar
André Anjos committed
510 511 512
  else Py_RETURN_FALSE;
}

513
static int PyBobApCeps_SetWithDeltaDelta
André Anjos's avatar
André Anjos committed
514 515 516 517 518 519
(PyBobApCepsObject* self, PyObject* o, void* /*closure*/) {

  bool b = PyObject_IsTrue(o);
  if (PyErr_Occurred()) return -1;

  try {
520
    self->cxx->setWithDeltaDelta(b);
André Anjos's avatar
André Anjos committed
521 522 523 524 525 526
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return -1;
  }
  catch (...) {
527
    PyErr_Format(PyExc_RuntimeError, "cannot reset `with_delta_delta' of %s: unknown exception caught", Py_TYPE(self)->tp_name);
André Anjos's avatar
André Anjos committed
528 529 530 531 532 533 534 535 536
    return -1;
  }

  return 0;

}

static PyGetSetDef PyBobApCeps_getseters[] = {
    {
537 538 539 540
      s_n_ceps_str,
      (getter)PyBobApCeps_GetNCeps,
      (setter)PyBobApCeps_SetNCeps,
      s_n_ceps_doc,
André Anjos's avatar
André Anjos committed
541 542 543
      0
    },
    {
544 545 546 547
      s_delta_win_str,
      (getter)PyBobApCeps_GetDeltaWin,
      (setter)PyBobApCeps_SetDeltaWin,
      s_delta_win_doc,
André Anjos's avatar
André Anjos committed
548 549 550
      0
    },
    {
551 552 553 554
      s_dct_norm_str,
      (getter)PyBobApCeps_GetDctNorm,
      (setter)PyBobApCeps_SetDctNorm,
      s_dct_norm_doc,
André Anjos's avatar
André Anjos committed
555 556 557
      0
    },
    {
558 559 560 561
      s_with_energy_str,
      (getter)PyBobApCeps_GetWithEnergy,
      (setter)PyBobApCeps_SetWithEnergy,
      s_with_energy_doc,
André Anjos's avatar
André Anjos committed
562 563 564
      0
    },
    {
565 566 567 568
      s_with_delta_str,
      (getter)PyBobApCeps_GetWithDelta,
      (setter)PyBobApCeps_SetWithDelta,
      s_with_delta_doc,
André Anjos's avatar
André Anjos committed
569 570 571
      0
    },
    {
572 573 574 575
      s_with_delta_delta_str,
      (getter)PyBobApCeps_GetWithDeltaDelta,
      (setter)PyBobApCeps_SetWithDeltaDelta,
      s_with_delta_delta_doc,
André Anjos's avatar
André Anjos committed
576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649
      0
    },
    {0}  /* Sentinel */
};

static PyObject* PyBobApCeps_Call
(PyBobApCepsObject* self, PyObject *args, PyObject* kwds) {

  /* Parses input arguments in a single shot */
  static const char* const_kwlist[] = {"input", "output", 0};
  static char** kwlist = const_cast<char**>(const_kwlist);

  PyBlitzArrayObject* input = 0;
  PyBlitzArrayObject* output = 0;
  if (!PyArg_ParseTupleAndKeywords(args, kwds, "O&|O&", kwlist,
        &PyBlitzArray_Converter, &input,
        &PyBlitzArray_OutputConverter, &output
        )) return 0;

  auto input_ = make_safe(input);
  auto output_ = make_xsafe(output);

  if (input->type_num != NPY_FLOAT64) {
    PyErr_Format(PyExc_TypeError, "`%s' only supports 1D 64-bit float arrays for input array `input'", Py_TYPE(self)->tp_name);
    return 0;
  }

  if (input->ndim != 1) {
    PyErr_Format(PyExc_TypeError, "`%s' only supports 1D 64-bit float arrays for input array `input'", Py_TYPE(self)->tp_name);
    return 0;
  }

  auto bz_input = PyBlitzArrayCxx_AsBlitz<double,1>(input);

  if (output) {

    if (output->type_num != NPY_FLOAT64) {
      PyErr_Format(PyExc_TypeError, "`%s' only supports 2D 64-bit float arrays for output array `output'", Py_TYPE(self)->tp_name);
      return 0;
    }

    if (output->ndim != 2) {
      PyErr_Format(PyExc_TypeError, "`%s' only supports 2D 64-bit float arrays for output array `output'", Py_TYPE(self)->tp_name);
      return 0;
    }

  }

  else {

    Py_ssize_t length[2];
    auto s = self->cxx->getShape(*bz_input);
    length[0] = s(0);
    length[1] = s(1);
    output = (PyBlitzArrayObject*)PyBlitzArray_SimpleNew(NPY_FLOAT64, 2, length);
    if (!output) return 0;
    output_ = make_safe(output);

  }

  auto bz_output = PyBlitzArrayCxx_AsBlitz<double,2>(output);

  try {
    self->cxx->operator()(*bz_input, *bz_output);
  }
  catch (std::exception& ex) {
    PyErr_SetString(PyExc_RuntimeError, ex.what());
    return 0;
  }
  catch (...) {
    PyErr_Format(PyExc_RuntimeError, "cannot call object of type `%s' - unknown exception thrown", Py_TYPE(self)->tp_name);
    return 0;
  }

650
  return PyBlitzArray_NUMPY_WRAP(Py_BuildValue("O", output));
André Anjos's avatar
André Anjos committed
651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693

}

PyTypeObject PyBobApCeps_Type = {
    PyVarObject_HEAD_INIT(0, 0)
    s_ceps_str,                               /*tp_name*/
    sizeof(PyBobApCepsObject),                /*tp_basicsize*/
    0,                                        /*tp_itemsize*/
    (destructor)PyBobApCeps_Delete,           /*tp_dealloc*/
    0,                                        /*tp_print*/
    0,                                        /*tp_getattr*/
    0,                                        /*tp_setattr*/
    0,                                        /*tp_compare*/
    (reprfunc)PyBobApCeps_Repr,               /*tp_repr*/
    0,                                        /*tp_as_number*/
    0,                                        /*tp_as_sequence*/
    0,                                        /*tp_as_mapping*/
    0,                                        /*tp_hash */
    (ternaryfunc)PyBobApCeps_Call,            /* tp_call */
    (reprfunc)PyBobApCeps_Repr,               /*tp_str*/
    0,                                        /*tp_getattro*/
    0,                                        /*tp_setattro*/
    0,                                        /*tp_as_buffer*/
    Py_TPFLAGS_DEFAULT | Py_TPFLAGS_BASETYPE, /*tp_flags*/
    s_ceps_doc,                               /* tp_doc */
    0,		                                    /* tp_traverse */
    0,		                                    /* tp_clear */
    (richcmpfunc)PyBobApCeps_RichCompare,     /* tp_richcompare */
    0,		                                    /* tp_weaklistoffset */
    0,		                                    /* tp_iter */
    0,		                                    /* tp_iternext */
    0,                                        /* tp_methods */
    0,                                        /* tp_members */
    PyBobApCeps_getseters,                    /* tp_getset */
    0,                                        /* tp_base */
    0,                                        /* tp_dict */
    0,                                        /* tp_descr_get */
    0,                                        /* tp_descr_set */
    0,                                        /* tp_dictoffset */
    (initproc)PyBobApCeps_Init,               /* tp_init */
    0,                                        /* tp_alloc */
    0,                                        /* tp_new */
};