Ignore:
Timestamp:
Feb 5, 2010, 5:24:39 PM (13 years ago)
Author:
steve
Message:

Updating to relocated repository

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/abstract_2d_finite_volumes/quantity_ext.c

    r7573 r7616  
    362362   
    363363    sign = 0.0;
    364     if (qmin > 0) {
     364    if (qmin > 0.0) {
    365365      sign = 1.0;
    366366    } else if (qmax < 0) {
     
    373373      dqa[i] = dq;                      //Save dq for use in updating vertex values 
    374374     
    375       // FIXME SR 20091125 This caused problems in shallow_water_balanced
    376       // commenting out problem
     375
    377376      // Just limit non boundary edges so that we can reconstruct a linear function
    378       if (neighbours[k3+i] >= 0) {
     377      // FIXME Problem with stability on edges
     378      //if (neighbours[k3+i] >= 0) {
    379379        r = 1.0;
    380380     
     
    383383           
    384384        phi = min( min(r*beta, 1.0), phi);
    385         }
    386 
    387       if (neighbours[k3+i] < 0) {
    388         r = 1.0;
    389      
    390         if (dq > 0.0 && sign == -1.0 ) r = (0.0 - qc)/dq;
    391         if (dq < 0.0 && sign ==  1.0 ) r = (0.0 - qc)/dq;     
     385        //      }
     386
     387      //
     388      /* if (neighbours[k3+i] < 0) { */
     389      /*        r = 1.0; */
     390     
     391      /*        if (dq > 0.0 && (sign == -1.0 || sign == 0.0 )) r = (0.0 - qc)/dq; */
     392      /*        if (dq < 0.0 && (sign ==  1.0 || sign == 0.0 )) r = (0.0 - qc)/dq; */
    392393           
    393         phi = min( min(r*beta, 1.0), phi);
    394         }
     394      /*        phi = min( min(r*beta, 1.0), phi); */
     395      /*        } */
    395396   
    396397    }
     
    11941195  if (!PyArg_ParseTuple(args, "O",&quantity)) {
    11951196      PyErr_SetString(PyExc_RuntimeError,
    1196                       "quantity_ext.c: extrapolate_second_order_and_limit could not parse input");     
     1197                      "quantity_ext.c: extrapolate_second_order_and_limit_by_edge could not parse input");     
    11971198      return NULL;
    11981199  }
     
    12011202  if (!domain) {
    12021203      PyErr_SetString(PyExc_RuntimeError,
    1203                       "quantity_ext.c: extrapolate_second_order_and_limit could not obtain domain object from quantity");       
     1204                      "quantity_ext.c: extrapolate_second_order_and_limit_by_edge could not obtain domain object from quantity");       
    12041205      return NULL;
    12051206  }
Note: See TracChangeset for help on using the changeset viewer.