Changeset 2978

Timestamp:

02/08/10 08:53:00 (15 years ago)

Author:

Mathieu Morlighem

Message:

Check that nu!=0 to avoid Matlab's annoying warning

File:

• issm/trunk/src/m/solutions/jpl/convergence.m (modified) (1 diff)

issm/trunk/src/m/solutions/jpl/convergence.m

@@ -45 +45 @@
         %print criterion
         if ~isnan(eps_rel),
-                if (ndu/nu<=eps_rel),
-                        displaystring(verbose,'%-60s%g%s%g%s','      relative convergence criterion: norm(du)/norm(u)',ndu/nu*100,' < ',eps_rel*100,' %');
+                if nu, %avoid "dividing by zero" warning
+                        if (ndu/nu<=eps_rel),
+                                displaystring(verbose,'%-60s%g%s%g%s','      relative convergence criterion: norm(du)/norm(u)',ndu/nu*100,' < ',eps_rel*100,' %');
+                        else
+                                displaystring(verbose,'%-60s%g%s%g%s','      relative convergence criterion: norm(du)/norm(u)',ndu/nu*100,' > ',eps_rel*100,' %');
+                                converged=0;
+                        end
                 else
-                        displaystring(verbose,'%-60s%g%s%g%s','      relative convergence criterion: norm(du)/norm(u)',ndu/nu*100,' > ',eps_rel*100,' %');
                         converged=0;
                 end
         else
-                displaystring(verbose,'%-60s%g%s','      relative convergence criterion: norm(du)/norm(u)',ndu/nu*100,' %');
+                if nu, %avoid "dividing by zero" warning
+                        displaystring(verbose,'%-60s%g%s','      relative convergence criterion: norm(du)/norm(u)',ndu/nu*100,' %');
+                end
         end