Changeset 26529


Ignore:
Timestamp:
11/02/21 10:02:15 (3 years ago)
Author:
Mathieu Morlighem
Message:

BUG: should not divide by 0 (detected with CoDiPack)

Location:
issm/trunk-jpl/src/c/classes/Elements
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

    r26503 r26529  
    39223922
    39233923        /*CFL criterion: */
    3924         IssmDouble dt = C/(maxabsvx/dx+maxabsvy/dy+maxabsvz/dz);
     3924        IssmDouble dt = C/(maxabsvx/dx+maxabsvy/dy+maxabsvz/dz + 1.e-18);
    39253925
    39263926        return dt;

  • issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r26487 r26529  
    50725072
    50735073        /*CFL criterion: */
    5074         IssmDouble dt = C/(maxabsvx/dx+maxabsvy/dy);
     5074        IssmDouble dt = C/(maxabsvx/dx+maxabsvy/dy + 1.e-18);
    50755075
    50765076        return dt;
Note: See TracChangeset for help on using the changeset viewer.