Changeset 18357

Timestamp:

08/08/14 20:35:30 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added check on Qminus and Qplus

File:

• issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp ¶

@@ -236 +236 @@
                 //double Qi_plus  = ml_serial[row]/deltat*(ulmax[row] - u[row]);
                 //double Qi_minus = ml_serial[row]/deltat*(ulmin[row] - u[row]);
+                _assert_(Qi_plus  >= 0.);
+                _assert_(Qi_minus <= 0.);
                 d = 1.;
                 if(Pi_plus!=0.) d = min(1.,Qi_plus/Pi_plus);