Changeset 10333


Ignore:
Timestamp:
10/27/11 09:00:57 (13 years ago)
Author:
Mathieu Morlighem
Message:

More flexible validation.m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/ad/validation/Validation/validation.m

    r10331 r10333  
    55md.miscellaneous.name='validation';
    66
    7 %launch solution directly:
    8 !source validation.queue
     7%run all scripts
     8for i=1:2,
    99
    10 md=loadresultsfromdisk(md,'validation.outbin');
     10        %run queuing script
     11        if i==1,
     12                disp('Checking src/ad/issm.exe ...');
     13                !source validation.queue
     14        elseif i==2,
     15                disp('Checking src/ad/issmROSE.exe ...');
     16                !source validationROSE.queue
     17        end
    1118
    12 s=load('Archive.mat'); h0=s.thickness;
    13 h=md.results.PrognosticSolution.Thickness;
     19        %load solution and remove outbin
     20        md=loadresultsfromdisk(md,'validation.outbin');
     21        !rm validation.outbin
    1422
    15 error_diff=full(max(abs(h-h0))/(max(abs(h0))+eps));
    16 
    17 if (error_diff>tolerance);
    18         disp(sprintf(['\nERROR   difference: %-7.2g > %7.2g \n'],error_diff,tolerance));
    19 else
    20         disp(sprintf(['\nSUCCESS difference: %-7.2g < %7.2g \n'],error_diff,tolerance));
     23        s=load('Archive.mat'); h0=s.thickness;
     24        h=md.results.PrognosticSolution.Thickness;
     25        error_diff=full(max(abs(h-h0))/(max(abs(h0))+eps));
     26        if (error_diff>tolerance);
     27                disp(sprintf(['\nERROR   difference: %-7.2g > %7.2g \n'],error_diff,tolerance));
     28        else
     29                disp(sprintf(['\nSUCCESS difference: %-7.2g < %7.2g \n'],error_diff,tolerance));
     30        end
    2131end
    22 
Note: See TracChangeset for help on using the changeset viewer.