Function Math_compute_womersley

Code/Source/PythonAPI/Math_PyModule.cxx:257–296

");

static PyObject *
Math_compute_womersley(PyObject *self, PyObject *args)
{
  auto api = PyUtilApiFunction("Odddddd", PyRunTimeErr, __func__);
  PyObject *termsArg;
  double time = 0;
  double viscosity = 0;
  double omega = 0;
  double density = 0;
  double radmax = 0;
  double radius = 0;

  if (!PyArg_ParseTuple(args, api.format, &termsArg, &time, &viscosity, &omega, &density, &radmax, &radius)) {
      return api.argsError();
  }

  // TODO:DaveP] check other args.

  // 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 the womersley operation.
  //
  auto mathObj = cvMath();
  double velocity = 0;
  if (mathObj.compute_v_womersley(terms, nlistterms, viscosity, density, omega, radmax, radius, time, &velocity) == SV_ERROR) {
       mathObj.deleteArray(terms,nlistterms,dim);
       api.error("Error calculating the womersley velocity.");
       return nullptr;
  }

  mathObj.deleteArray(terms,nlistterms,dim);

  return Py_BuildValue("d",velocity);
}

//---------------------------
// Math_linear_interpolate