Ignore:
Timestamp:
Jun 19, 2011, 3:54:42 PM (14 years ago)
Author:
paul
Message:

Added state evolved quantity to pipe domain

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/anuga_work/development/2010-projects/anuga_1d/pipe/pipe_domain_ext.c

    r8187 r8188  
    230230                                double* width_edge_values,
    231231                                double* top_edge_values,
     232                                double* state_edge_values,
    232233                                double* area_boundary_values,
    233234                                double* discharge_boundary_values,
    … …  
    237238                                double* width_boundary_values,
    238239                                double* top_boundary_values,
     240                                double* state_boundary_values,
    239241                                double* area_explicit_update,
    240242                                double* discharge_explicit_update,
    … …  
    242244                                double* max_speed_array) {
    243245
    244   double flux[2], qlm[7], qlp[7], qrm[7], qrp[7], edgeflux[2];
     246  double flux[2], qlm[8], qlp[8], qrm[8], qrp[8], edgeflux[2];
    245247  double max_speed;
    246248  int k, ki, n, m, nm=0;
    … …  
    264266      qlm[5] = width_boundary_values[m];
    265267      qlm[6] = top_boundary_values[m];
     268      qlm[7] = state_boundary_values[m];
    266269
    267270    } else {
    … …  
    277280      qlm[5] = width_edge_values[nm];
    278281      qlm[6] = top_edge_values[nm];
     282      qlm[7] = state_edge_values[nm];
    279283    }
    280284
    … …  
    286290    qlp[5] = width_edge_values[ki];
    287291    qlp[6] = top_edge_values[ki];
     292    qlp[7] = state_edge_values[ki];
    288293
    289294    ki = k*2+1;
    … …  
    299304      qrp[5] = width_boundary_values[m];
    300305      qrp[6] = top_boundary_values[m];
     306      qrp[7] = state_boundary_values[m];
    301307    } else {
    302308      m = neighbour_vertices[ki];
    … …  
    311317      qrp[5] = width_edge_values[nm];
    312318      qrp[6] = top_edge_values[nm];
     319      qrp[7] = state_edge_values[nm];
    313320    }
    314321
    … …  
    320327    qrm[5] = width_edge_values[ki];
    321328    qrm[6] = top_edge_values[ki];
     329    qrm[7] = state_edge_values[ki];
    322330
    323331    _flux_function_pipe(qlm,qlp,qrm,qrp,g,epsilon,h0,edgeflux,&max_speed);
    … …  
    367375    *velocity,
    368376    *width,
    369     *top;
     377    *top,
     378    *state;
    370379
    371380  PyArrayObject
    … …  
    381390    *width_vertex_values,
    382391    *top_vertex_values,
     392    *state_vertex_values,
    383393    *area_boundary_values,
    384394    *discharge_boundary_values,
    … …  
    388398    *width_boundary_values,
    389399    *top_boundary_values,
     400    *state_boundary_values,
    390401    *area_explicit_update,
    391402    *discharge_explicit_update,
    … …  
    397408
    398409  // Convert Python arguments to C
    399   if (!PyArg_ParseTuple(args, "dOOOOOOOO",
     410  if (!PyArg_ParseTuple(args, "dOOOOOOOOO",
    400411                        &timestep,
    401412                        &domain,
    … …  
    406417                        &velocity,
    407418                        &width,
    408                         &top)) {
     419                        &top,
     420                        &state)) {
    409421    PyErr_SetString(PyExc_RuntimeError, "comp_flux_pipe_ext.c: compute_fluxes_pipe_ext could not parse input");
    410422    return NULL;
    … …  
    431443  width_vertex_values      = get_consecutive_array(width, "vertex_values");
    432444  top_vertex_values      = get_consecutive_array(top, "vertex_values");
     445  state_vertex_values      = get_consecutive_array(state, "vertex_values");
    433446
    434447  area_boundary_values     = get_consecutive_array(area,     "boundary_values");
    … …  
    439452  width_boundary_values       = get_consecutive_array(width,     "boundary_values");
    440453  top_boundary_values       = get_consecutive_array(top,     "boundary_values");
     454  state_boundary_values       = get_consecutive_array(state,     "boundary_values");
    441455
    442456
    … …  
    464478                                      (double*) width_vertex_values -> data,
    465479                                      (double*) top_vertex_values -> data,
     480                                      (double*) state_vertex_values -> data,
    466481                                      (double*) area_boundary_values -> data,
    467482                                      (double*) discharge_boundary_values -> data,
    … …  
    471486                                      (double*) width_boundary_values -> data,
    472487                                      (double*) top_boundary_values -> data,
     488                                      (double*) state_boundary_values -> data,
    473489                                      (double*) area_explicit_update -> data,
    474490                                      (double*) discharge_explicit_update -> data,
    … …  
    488504  Py_DECREF(width_vertex_values);
    489505  Py_DECREF(top_vertex_values);
     506  Py_DECREF(state_vertex_values);
    490507  Py_DECREF(area_boundary_values);
    491508  Py_DECREF(discharge_boundary_values);
    … …  
    495512  Py_DECREF(width_boundary_values);
    496513  Py_DECREF(top_boundary_values);
     514  Py_DECREF(state_boundary_values);
    497515  Py_DECREF(area_explicit_update);
    498516  Py_DECREF(discharge_explicit_update);
Note: See TracChangeset for help on using the changeset viewer.