Changeset 1412

Timestamp:

07/29/09 14:23:02 (16 years ago)

Author:

Mathieu Morlighem

Message:

Added residue convergence criterion

File:

• issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m (modified) (2 diffs)

issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m

@@ -50 +50 @@
 
                 %Get residue
-                residu=norm(K_ff*soln(count).u_f-p_f,2);
+                residu=norm(K_ff*soln(count).u_f-p_f,2)/norm(p_f,2);
 
                 %Merge back to g set
@@ -72 +72 @@
 
                 %Residue criterion
-                displaystring(m.parameters.debug,'%s%g','      convergence criterion: norm(residue)=',residu);
+                displaystring(m.parameters.debug,'%s%g','      convergence criterion: norm(KU-F)/norm(F)=',residu);
 
                 %Relative criterion