Changeset 7782


Ignore:
Timestamp:
Jun 5, 2010, 7:03:20 PM (12 years ago)
Author:
steve
Message:

cleaning up shallow water domain

File:
1 copied

Legend:

Unmodified
Added
Removed
  • anuga_work/development/pipeflow/swb_comp_flux_ext.c

    r7780 r7782  
    77
    88//Innermost flux function (using w=z+h)
    9 int _flux_function_wellbalanced(double *q_left,
     9int _flux_function(double *q_left,
    1010                                                                double *q_right,
    1111                                                                double normals,
     
    1515                                                                double *max_speed) {           
    1616                int i;
    17         double z_left, z_right;
     17                double z_left, z_right;
    1818                double flux_left[2], flux_right[2];
    1919                double z_star, w_left, h_left, h_left_star, uh_left, soundspeed_left, u_left;
     
    2424                uh_left = q_left[1];
    2525                uh_left = uh_left*normals;
    26         z_left  = q_left[2];
     26                z_left  = q_left[2];
    2727                h_left  = q_left[3];
    2828                u_left  = q_left[4];
     
    3131                uh_right = q_right[1];
    3232                uh_right = uh_right*normals;
    33         z_right  = q_right[2];
     33                z_right  = q_right[2];
    3434                h_right  = q_right[3];
    3535                u_right  = q_right[4];
     
    108108       
    109109// Computational function for flux computation
    110 double _compute_fluxes_ext_wellbalanced(double timestep,
     110double _compute_fluxes_ext(double timestep,
    111111                double epsilon,
    112112                double g,
     
    152152                                ql[1] = xmom_edge_values[ki];
    153153                                ql[2] = bed_edge_values[ki];
    154                 ql[3] = height_edge_values[ki];
     154                                ql[3] = height_edge_values[ki];
    155155                                ql[4] = vel_edge_values[ki];
    156156               
     
    176176                               
    177177                                normal = normals[ki];
    178                                 _flux_function_wellbalanced(ql, qr, normal, g, epsilon, edgeflux, &max_speed);
     178                                _flux_function(ql, qr, normal, g, epsilon, edgeflux, &max_speed);
    179179                                flux[0] -= edgeflux[0];
    180180                                flux[1] -= edgeflux[1];
     
    209209// Python Glue
    210210//=========================================================================
    211 PyObject *compute_fluxes_ext_wellbalanced(PyObject *self, PyObject *args) {
     211PyObject *compute_fluxes_ext(PyObject *self, PyObject *args) {
    212212 
    213213    PyArrayObject
     
    271271                        &number_of_elements,
    272272                        &max_speed_array)) {
    273     PyErr_SetString(PyExc_RuntimeError, "comp_flux_ext_wellbalanced.c: compute_fluxes_ext_wellbalanced could not parse input");
     273    PyErr_SetString(PyExc_RuntimeError, "swb_comp_flux_ext.c: compute_fluxes_ext could not parse input");
    274274    return NULL;
    275275  }
     
    278278  // Call underlying flux computation routine and update
    279279  // the explicit update arrays
    280   timestep = _compute_fluxes_ext_wellbalanced(timestep,
     280  timestep = _compute_fluxes_ext(timestep,
    281281                                 epsilon,
    282282                                 g,
     
    319319
    320320static struct PyMethodDef MethodTable[] = {
    321   {"compute_fluxes_ext_wellbalanced", compute_fluxes_ext_wellbalanced, METH_VARARGS, "Print out"},
     321  {"compute_fluxes_ext", compute_fluxes_ext, METH_VARARGS, "Print out"},
    322322  {NULL, NULL}
    323323};
    324324
    325325// Module initialisation
    326 void initcomp_flux_ext_wellbalanced(void){
    327   Py_InitModule("comp_flux_ext_wellbalanced", MethodTable);
     326void initswb_comp_flux_ext(void){
     327  Py_InitModule("swb_comp_flux_ext", MethodTable);
    328328  import_array();
    329329}
Note: See TracChangeset for help on using the changeset viewer.