source: trunk/anuga_work/development/2010-projects/anuga_1d/sqpipe/sqpipe_pressurised_fixed_width_a_d_comp_flux.c @ 8239

#include "Python.h"
#include "numpy/arrayobject.h"
#include "math.h"
#include <stdio.h>

const double pi = 3.14159265358979;

// Shared code snippets
#include "util_ext.h"
#include "solver.h"
#include "python_solver.h"
#include "sqpipe_pressurised_fixed_width_a_d.h"
#include "sqpipe_python.h"

//=========================================================================
// Python Glue
//=========================================================================

PyObject *compute_fluxes_ext(PyObject *self, PyObject *args) {
  struct domain *D = (struct domain *) sqpipe_pressurised_fixed_width_a_d_domain_new();
  double timestep;
  PyObject *domain;

  // Convert Python arguments to C
  if (!PyArg_ParseTuple(args, "dO", &timestep, &domain)) {
    PyErr_SetString(PyExc_RuntimeError, "comp_flux_ext.c: compute_fluxes_ext could not parse input");
    return NULL;
  }

  D = (struct domain *) sqpipe_domain_python_get((struct sqpipe_domain *)D, domain, timestep);

  timestep = domain_compute_fluxes(D);
  domain_destroy(D);
  
  return Py_BuildValue("d", timestep);
}

//-------------------------------
// Method table for python module
//-------------------------------

static struct PyMethodDef MethodTable[] = {
  {"compute_fluxes_ext", compute_fluxes_ext, METH_VARARGS, "Print out"},
  {NULL, NULL}
};

// Module initialisation
void initsqpipe_pressurised_fixed_width_a_d_comp_flux(void){
  Py_InitModule("sqpipe_pressurised_fixed_width_a_d_comp_flux", MethodTable);
  import_array();
}