Changeset 18824


Ignore:
Timestamp:
11/21/14 15:15:37 (10 years ago)
Author:
hongjuy
Message:

NEW: (1) added one new cost function (507) for rheology B regularization
(2) added one control parameter (RheologyB not depth averaged)
(3) cosmetics

Location:
issm/trunk-jpl/src/m/classes
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/classes/clusters/greenplanet.m

    r18305 r18824  
    99    properties (SetAccess=public) 
    1010                 % {{{
    11                  name='greenplanet'
    12                  login='';
    13                  numnodes=20;
    14                  cpuspernode=8;
    15                  port=8000;
    16                  queue='c6145';
    17                  codepath='';
    18                  executionpath='';
    19                  interactive=0;
    20                  time=24*60;
    21                  memory=2;
     11                 name          = 'greenplanet'
     12                 login         = '';
     13                 numnodes      = 20;
     14                 cpuspernode   = 8;
     15                 port          = 8000;
     16                 queue         = 'c6145';
     17                 codepath      = '';
     18                 executionpath = '';
     19                 interactive   = 0;
     20                 time          = 24*60;
     21                 memory        = 2;
    2222         end
    2323         properties (SetAccess=private)

  • issm/trunk-jpl/src/m/classes/inversion.m

    r18771 r18824  
    134134                        md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',...
    135135                                {'BalancethicknessThickeningRate' 'FrictionCoefficient' 'MaterialsRheologyBbar' 'DamageDbar' 'Vx' 'Vy' 'Thickness',...
    136                                 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance'});
     136                                'BalancethicknessOmega' 'BalancethicknessApparentMassbalance','MaterialsRheologyB'});
    137137                        md = checkfield(md,'fieldname','inversion.nsteps','numel',1,'>=',0);
    138138                        md = checkfield(md,'fieldname','inversion.maxiter_per_step','size',[md.inversion.nsteps 1],'>=',0);
    139139                        md = checkfield(md,'fieldname','inversion.step_threshold','size',[md.inversion.nsteps 1]);
    140                         md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:506 601:604]);
     140                        md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:507 601:604]);
    141141                        md = checkfield(md,'fieldname','inversion.cost_functions_coefficients','size',[md.mesh.numberofvertices num_costfunc],'>=',0);
    142142                        md = checkfield(md,'fieldname','inversion.gradient_scaling','size',[md.inversion.nsteps num_controls]);
    … …  
    242242                        pos=find(obj.cost_functions==505); data(pos)=ThicknessAcrossGradientEnum();
    243243                        pos=find(obj.cost_functions==506); data(pos)=BalancethicknessMisfitEnum();
     244                        pos=find(obj.cost_functions==507); data(pos)=RheologyBAbsGradientEnum();
    244245                        pos=find(obj.cost_functions==601); data(pos)=SurfaceAbsMisfitEnum();
    245246                        WriteData(fid,'data',data,'enum',InversionCostFunctionsEnum(),'format','DoubleMat','mattype',3);

  • issm/trunk-jpl/src/m/classes/inversionvalidation.m

    r18617 r18824  
    5959                        md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',...
    6060                                {'BalancethicknessThickeningRate' 'FrictionCoefficient' 'MaterialsRheologyBbar' 'DamageDbar' 'Vx' 'Vy' 'Thickness',...
    61                                 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance'});
     61                                'BalancethicknessOmega' 'BalancethicknessApparentMassbalance','MaterialsRheologyB'});
    6262                        md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1);
    63                         md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:506 601:604]);
     63                        md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:507 601:604]);
    6464                        md = checkfield(md,'fieldname','inversion.cost_functions_coefficients','size',[md.mesh.numberofvertices num_costfunc],'>=',0);
    6565                        md = checkfield(md,'fieldname','inversion.min_parameters','size',[md.mesh.numberofvertices num_controls]);
    … …  
    150150                        pos=find(obj.cost_functions==505); data(pos)=ThicknessAcrossGradientEnum();
    151151                        pos=find(obj.cost_functions==506); data(pos)=BalancethicknessMisfitEnum();
     152                        pos=find(obj.cost_functions==507); data(pos)=RheologyBAbsGradientEnum();
    152153                        pos=find(obj.cost_functions==601); data(pos)=SurfaceAbsMisfitEnum();
    153154                        WriteData(fid,'data',data,'enum',InversionCostFunctionsEnum(),'format','DoubleMat','mattype',3);

  • issm/trunk-jpl/src/m/classes/m1qn3inversion.m

    r18615 r18824  
    6969                                md = checkmessage(md,['M1QN3 has not been installed, ISSM needs to be reconfigured and recompiled with M1QN3']);
    7070                        end
    71 
    7271                        num_controls=numel(md.inversion.control_parameters);
    7372                        num_costfunc=size(md.inversion.cost_functions,2);
    … …  
    7776                        md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',...
    7877                                {'BalancethicknessThickeningRate' 'FrictionCoefficient' 'MaterialsRheologyBbar' 'DamageDbar',...
    79                                 'Vx' 'Vy' 'Thickness' 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance'});
     78                                'Vx' 'Vy' 'Thickness' 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance','MaterialsRheologyB'});
    8079                        md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1);
    8180                        md = checkfield(md,'fieldname','inversion.maxsteps','numel',1,'>=',0);
    … …  
    8382                        md = checkfield(md,'fieldname','inversion.dxmin','numel',1,'>',0);
    8483                        md = checkfield(md,'fieldname','inversion.gttol','numel',1,'>',0);
    85                         md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:506 601:604]);
     84                        md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:507 601:604]);
    8685                        md = checkfield(md,'fieldname','inversion.cost_functions_coefficients','size',[md.mesh.numberofvertices num_costfunc],'>=',0);
    8786                        md = checkfield(md,'fieldname','inversion.min_parameters','size',[md.mesh.numberofvertices num_controls]);
    … …  
    181180                        pos=find(obj.cost_functions==505); data(pos)=ThicknessAcrossGradientEnum();
    182181                        pos=find(obj.cost_functions==506); data(pos)=BalancethicknessMisfitEnum();
     182                        pos=find(obj.cost_functions==507); data(pos)=RheologyBAbsGradientEnum();
    183183                        pos=find(obj.cost_functions==601); data(pos)=SurfaceAbsMisfitEnum();
    184184                        WriteData(fid,'data',data,'enum',InversionCostFunctionsEnum(),'format','DoubleMat','mattype',3);
Note: See TracChangeset for help on using the changeset viewer.