Function Math_inverse_fft

Code/Source/PythonAPI/Math_PyModule.cxx:189–240

");

static PyObject *
Math_inverse_fft(PyObject *self, PyObject *args)
{
  auto api = PyUtilApiFunction("Odddi", PyRunTimeErr, __func__);
  PyObject *termsArg;
  double t0 = 0;
  double dt = 0;
  double omega = 0;
  int numPts = 0;

  if (!PyArg_ParseTuple(args, api.format, &termsArg, &t0, &dt, &omega, &numPts)) { return api.argsError();
      return api.argsError();
  }

  // Get an array of points from the pointsArgs list.
  int dim = 2;
  auto terms = GetPointsFromList(api, termsArg, dim, "terms");
  if (terms == nullptr) {
      return nullptr;
  }
  int nlistterms = PyList_Size(termsArg);

  // Perform inverse fft operation.
  //
  auto mathObj = cvMath();
  double **pts = NULL;
  if (mathObj.inverseFFT(terms, nlistterms, t0, dt, omega, numPts, &pts) == SV_ERROR) {
       mathObj.deleteArray(terms,nlistterms,dim);
       api.error("Error calculating the inverse fft.");
       return nullptr;
  }

  // Create result string
  //
  // [TODO:DaveP] This converts results to a string, not
  //    a list list the fft function above.
  //
  char r[2048];
  PyObject *pylist = PyList_New(numPts);
  for (int i = 0; i < numPts; i++) {
      r[0] = '\0';
      sprintf(r,"%.6le %.6le",pts[i][0],pts[i][1]);
      PyObject *rr = PyBytes_FromString(r);
      PyList_SET_ITEM(pylist, i, rr);
  }

  // [TODO:DaveP] this is hideous!
  mathObj.deleteArray(terms,nlistterms,dim);
  mathObj.deleteArray(pts,numPts,dim);

  return pylist;
}

//--------------------------
// Math_compute_womersley