Changeset 10330
- Timestamp:
- 10/27/11 08:36:29 (13 years ago)
- Location:
- issm/trunk/src/ad/validation
- Files:
-
- 2 added
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/ad/validation/Update/update.m
r10328 r10330 1 1 tolerance=1e-13; 2 2 3 md=setmesh(model,' Exp/Square.exp',150000);3 md=setmesh(model,'Square.exp',150000); 4 4 md=setmask(md,'all',''); 5 md=parameterize(md,' Par/SquareShelfConstrained.par');5 md=parameterize(md,'SquareShelfConstrained.par'); 6 6 md=setflowequation(md,'macayeal','all'); 7 md.settings.results_on_vertices=1; 7 8 md.cluster=generic('name',oshostname(),'np',3); 9 md.cluster.executionpath=[issmtier() '/src/ad/validation/Validation/']; 8 10 9 md.autodiff. analysis=true;11 md.autodiff.isautodiff=true; 10 12 md.autodiff.forward=true; 11 13 md.verbose=verbose('solution',true); … … 13 15 md=solve(md,PrognosticSolutionEnum); 14 16 15 s=load('Archive/archive.mat'); h0=s.Archive110_field1; 16 h=PatchToVec(md.results.PrognosticSolution.Thickness); 17 18 error_diff=full(max(abs(h-h0))/(max(abs(h0))+eps)); 19 20 if (error_diff>tolerance); 21 disp(sprintf(['\n\n\nERROR difference: %-7.2g > %7.2g \n\n\n'],error_diff,tolerance)); 22 else 23 disp(sprintf(['\n\n\nSUCCESS difference: %-7.2g < %7.2g \n\n\n'],error_diff,tolerance)); 24 end 25 17 thickness=md.results.PrognosticSolution.Thickness; 18 save ../Validation/Archive.mat thickness; -
issm/trunk/src/ad/validation/Validation/validation.m
r10329 r10330 4 4 md.cluster=generic('name',oshostname(),'np',3); 5 5 md.miscellaneous.name='validation'; 6 7 6 8 7 %launch solution directly: … … 12 11 13 12 s=load('Archive/archive.mat'); h0=s.Archive110_field1; 14 h= PatchToVec(md.results.PrognosticSolution.Thickness);13 h=md.results.PrognosticSolution.Thickness; 15 14 16 15 error_diff=full(max(abs(h-h0))/(max(abs(h0))+eps)); -
issm/trunk/src/ad/validation/Validation/validation.queue
r10329 r10330 1 1 #!/bin/sh 2 mpiexec -np 3 $ISSM_TIER/bin/issm.exe PrognosticSolution ./ validation> validation.outlog2 mpiexec -np 3 $ISSM_TIER/bin/issm.exe PrognosticSolution ./ validation > validation.outlog
Note:
See TracChangeset
for help on using the changeset viewer.