Changeset 79


Ignore:
Timestamp:
04/28/09 13:33:24 (16 years ago)
Author:
Mathieu Morlighem
Message:

Fixed message if eps_abs=NaN

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m

    r78 r79  
    8484
    8585                        %Absolute criterion
    86                         if (isnan(m.parameters.eps_abs) | nduinf<=m.parameters.eps_abs),
    87                                 if m.parameters.debug, disp(sprintf('%s %g %s %g','      convergence criterion: max(du)=',nduinf,' < ',m.parameters.eps_abs)); end
    88                         else
    89                                 if m.parameters.debug, disp(sprintf('%s %g %s %g','      convergence criterion: max(du)=',nduinf,' > ',m.parameters.eps_abs)); end
    90                                 converged=0;
     86                        if ~isnan(m.parameters.eps_abs),
     87                                if (nduinf<=m.parameters.eps_abs),
     88                                        if m.parameters.debug, disp(sprintf('%s %g %s %g','      convergence criterion: max(du)=',nduinf,' < ',m.parameters.eps_abs)); end
     89                                else
     90                                        if m.parameters.debug, disp(sprintf('%s %g %s %g','      convergence criterion: max(du)=',nduinf,' > ',m.parameters.eps_abs)); end
     91                                                converged=0;
     92                                end
    9193                        end
    92                        
     94
    9395                        %rift convergence criterion
    9496                        if ~constraints_converged,
Note: See TracChangeset for help on using the changeset viewer.