Changeset 10333
- Timestamp:
- 10/27/11 09:00:57 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/ad/validation/Validation/validation.m
r10331 r10333 5 5 md.miscellaneous.name='validation'; 6 6 7 % launch solution directly:8 !source validation.queue 7 %run all scripts 8 for i=1:2, 9 9 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 11 18 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 14 22 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 21 31 end 22
Note:
See TracChangeset
for help on using the changeset viewer.