Changeset 1919


Ignore:
Timestamp:
08/26/09 10:11:54 (16 years ago)
Author:
Mathieu Morlighem
Message:

updated updatearchive of CM

Location:
issm/trunk/test/Verification
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/test/Verification/PigControlMethodBM2d_21/runme.m

    r1912 r1919  
    7272
    7373                %compute fields to be checked
    74                 fields=testsgetfields(md.type,'control');
     74                fields=testsgetfields(md.type,solution);
    7575
    7676                %load archive

  • issm/trunk/test/Verification/PigControlMethodBM2d_21/updatearchive.m

    r1908 r1919  
    2626end
    2727
    28 %packages and misfits to be tested
     28%packages and solutions to be tested
    2929if nargin==2,
    30         misfits=varargin{2};
     30        solutions=varargin{2};
    3131else
    32         misfits={'absolute','relative','logarithmic'};
     32        solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
    3333end
    3434if nargin,
    … …  
    4545        package=packages{i};
    4646
    47         for j=1:length(misfits),
    48                 misfit=misfits{j};
     47        for j=1:length(solutions),
     48                solution=solutions{j};
    4949
    5050                %check solution requested
    51                 if ((strcmpi(package,'macayeal') & ~strcmpi(misfit,'absolute')) | strcmpi(package,'ice')),
    52                         disp(['package: ' package  ' and misfit: ' misfit  ' in test: ' testname  ' not supported yet']);
     51                if ((strcmpi(package,'macayeal') & ~strcmpi(solution,'absolute')) | strcmpi(package,'ice')),
     52                        disp(['package: ' package  ' and solution: ' solution  ' in test: ' testname  ' not supported yet']);
    5353                        continue
    5454                end
    … …  
    6060                md=parameterize(md,'Pig.par');
    6161                md=setelementstype(md,'macayeal','all');
    62                 if strcmpi(misfit,'absolute'),
    63                         md.fit(:)=0;
    64                 elseif strcmpi(misfit,'relative'),
    65                         md.fit(:)=1;
    66                 elseif strcmpi(misfit,'logarithmic'),
    67                         md.fit(:)=2;
    68                 end
    6962                md=modelextract(md,md.elementoniceshelf);
    7063
    71                 %launch control method
    72                 if ~testscheckmisfit(misfit), continue; end
    73                 md=solve(md,'analysis_type','control','package',packagefinal);
     64                %compute solution
     65                [md analysis_type sub_analysis_type packagefinal]=testsparameterization(md,solution,package);
     66                if ~testschecksolution(md.type,solution), continue; end
     67                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',packagefinal);
    7468
    7569                %compute fields to be checked
    76                 fields=testsgetfields(md.type,'control');
     70                fields=testsgetfields(md.type,solution);
    7771
    7872                %save new archive
    7973                for k=1:length(fields),
    8074                        field=fields{k};
    81                         eval(['Archive' package misfit '_field' num2str(k) '=md.results.' field  ';']);
     75                        eval(['Archive' package solution '_field' num2str(k) '=md.results.' field  ';']);
    8276                end
    83                 eval(['save Archive' package misfit ' Archive' package misfit '_field*']);
     77                eval(['save Archive' package solution ' Archive' package solution '_field*']);
    8478        end
    8579end

  • issm/trunk/test/Verification/PigControlMethodDragM2d_20/runme.m

    r1912 r1919  
    7171
    7272                %compute fields to be checked
    73                 fields=testsgetfields(md.type,'control');
     73                fields=testsgetfields(md.type,solution);
    7474
    7575                %load archive

  • issm/trunk/test/Verification/PigControlMethodDragM2d_20/updatearchive.m

    r1908 r1919  
    2626end
    2727
    28 %packages and misfits to be tested
     28%packages and solutions to be tested
    2929if nargin==2,
    30         misfits=varargin{2};
     30        solutions=varargin{2};
    3131else
    32         misfits={'absolute','relative','logarithmic'};
     32        solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
    3333end
    3434if nargin,
    … …  
    4545        package=packages{i};
    4646
    47         for j=1:length(misfits),
    48                 misfit=misfits{j};
     47        for j=1:length(solutions),
     48                solution=solutions{j};
    4949
    5050                %check package
    … …  
    6060                md=parameterize(md,'Pig.par');
    6161                md=setelementstype(md,'macayeal','all');
    62                 if strcmpi(misfit,'absolute'),
    63                         md.fit(:)=0;
    64                 elseif strcmpi(misfit,'relative'),
    65                         md.fit(:)=1;
    66                 elseif strcmpi(misfit,'logarithmic'),
    67                         md.fit(:)=2;
    68                 end
    6962
    70                 %launch control method
    71                 if ~testscheckmisfit(misfit), continue; end
    72                 md=solve(md,'analysis_type','control','package',packagefinal);
     63                %compute solution
     64                [md analysis_type sub_analysis_type packagefinal]=testsparameterization(md,solution,package);
     65                if ~testschecksolution(md.type,solution), continue; end
     66                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',packagefinal);
    7367
    7468                %compute fields to be checked
    75                 fields=testsgetfields(md.type,'control');
     69                fields=testsgetfields(md.type,solution);
    7670
    7771                %save new archive
    7872                for k=1:length(fields),
    7973                        field=fields{k};
    80                         eval(['Archive' package misfit '_field' num2str(k) '=md.results.' field  ';']);
     74                        eval(['Archive' package solution '_field' num2str(k) '=md.results.' field  ';']);
    8175                end
    82                 eval(['save Archive' package misfit ' Archive' package misfit '_field*']);
     76                eval(['save Archive' package solution ' Archive' package solution '_field*']);
    8377        end
    8478end

  • issm/trunk/test/Verification/PigControlMethodDragP3d_22/runme.m

    r1912 r1919  
    7272
    7373                %compute fields to be checked
    74                 fields=testsgetfields(md.type,'control');
     74                fields=testsgetfields(md.type,solution);
    7575
    7676                %load archive

  • issm/trunk/test/Verification/PigControlMethodDragP3d_22/updatearchive.m

    r1908 r1919  
    2626end
    2727
    28 %packages and misfits to be tested
     28%packages and solutions to be tested
    2929if nargin==2,
    30         misfits=varargin{2};
     30        solutions=varargin{2};
    3131else
    32         misfits={'absolute','relative','logarithmic'};
     32        solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
    3333end
    3434if nargin,
    … …  
    4545        package=packages{i};
    4646
    47         for j=1:length(misfits),
    48                 misfit=misfits{j};
     47        for j=1:length(solutions),
     48                solution=solutions{j};
    4949
    5050                %check package
    … …  
    6161                md=extrude(md,4,2);
    6262                md=setelementstype(md,'pattyn','all');
    63                 if strcmpi(misfit,'absolute'),
    64                         md.fit(:)=0;
    65                 elseif strcmpi(misfit,'relative'),
    66                         md.fit(:)=1;
    67                 elseif strcmpi(misfit,'logarithmic'),
    68                         md.fit(:)=2;
    69                 end
    7063
    71                 %launch control method
    72                 if ~testscheckmisfit(misfit), continue; end
    73                 md=solve(md,'analysis_type','control','package',packagefinal);
     64                %compute solution
     65                [md analysis_type sub_analysis_type packagefinal]=testsparameterization(md,solution,package);
     66                if ~testschecksolution(md.type,solution), continue; end
     67                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',packagefinal);
    7468
    7569                %compute fields to be checked
    76                 fields=testsgetfields(md.type,'control');
     70                fields=testsgetfields(md.type,solution);
    7771
    7872                %save new archive
    7973                for k=1:length(fields),
    8074                        field=fields{k};
    81                         eval(['Archive' package misfit '_field' num2str(k) '=md.results.' field  ';']);
     75                        eval(['Archive' package solution '_field' num2str(k) '=md.results.' field  ';']);
    8276                end
    83                 eval(['save Archive' package misfit ' Archive' package misfit '_field*']);
     77                eval(['save Archive' package solution ' Archive' package solution '_field*']);
    8478        end
    8579end

  • issm/trunk/test/Verification/PigControlMethodDragS3d_23/runme.m

    r1912 r1919  
    7272
    7373                %compute fields to be checked
    74                 fields=testsgetfields(md.type,'control');
     74                fields=testsgetfields(md.type,solution);
    7575
    7676                %load archive

  • issm/trunk/test/Verification/PigControlMethodDragS3d_23/updatearchive.m

    r1908 r1919  
    2626end
    2727
    28 %packages and misfits to be tested
     28%packages and solutions to be tested
    2929if nargin==2,
    30         misfits=varargin{2};
     30        solutions=varargin{2};
    3131else
    32         misfits={'absolute','relative','logarithmic'};
     32        solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
    3333end
    3434if nargin,
    … …  
    4545        package=packages{i};
    4646
    47         for j=1:length(misfits),
    48                 misfit=misfits{j};
     47        for j=1:length(solutions),
     48                solution=solutions{j};
    4949
    5050                %check package
    … …  
    6161                md=extrude(md,4,3);
    6262                md=setelementstype(md,'pattyn','all','stokes','all');
    63                 if strcmpi(misfit,'absolute'),
    64                         md.fit(:)=0;
    65                 elseif strcmpi(misfit,'relative'),
    66                         md.fit(:)=1;
    67                 elseif strcmpi(misfit,'logarithmic'),
    68                         md.fit(:)=2;
    69                 end
    7063
    71                 %launch control method
    72                 if ~testscheckmisfit(misfit), continue; end
    73                 md=solve(md,'analysis_type','control','package',packagefinal);
     64                %compute solution
     65                [md analysis_type sub_analysis_type packagefinal]=testsparameterization(md,solution,package);
     66                if ~testschecksolution(md.type,solution), continue; end
     67                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',packagefinal);
    7468
    7569                %compute fields to be checked
    76                 fields=testsgetfields(md.type,'control');
     70                fields=testsgetfields(md.type,solution);
    7771
    7872                %save new archive
    7973                for k=1:length(fields),
    8074                        field=fields{k};
    81                         eval(['Archive' package misfit '_field' num2str(k) '=md.results.' field  ';']);
     75                        eval(['Archive' package solution '_field' num2str(k) '=md.results.' field  ';']);
    8276                end
    83                 eval(['save Archive' package misfit ' Archive' package misfit '_field*']);
     77                eval(['save Archive' package solution ' Archive' package solution '_field*']);
    8478        end
    8579end
Note: See TracChangeset for help on using the changeset viewer.