Ignore:
Timestamp:
08/02/12 09:52:47 (13 years ago)
Author:
cborstad
Message:

stripped inversion threshold out of code, not used

Location:
issm/branches/trunk-jpl-damage
Files:
45 edited
1 copied

Legend:

Unmodified
Added
Removed
  • issm/branches/trunk-jpl-damage

  • issm/branches/trunk-jpl-damage/src/c/EnumDefinitions/EnumDefinitions.h

    r12280 r12872  
    8282        InversionStepThresholdEnum,
    8383        InversionThicknessObsEnum,
    84         InversionThreshParametersEnum,
    8584        InversionVelObsEnum,
    8685        InversionVxObsEnum,
  • issm/branches/trunk-jpl-damage/src/c/modules/EnumToStringx/EnumToStringx.cpp

    r12280 r12872  
    8686                case InversionStepThresholdEnum : return "InversionStepThreshold";
    8787                case InversionThicknessObsEnum : return "InversionThicknessObs";
    88                 case InversionThreshParametersEnum : return "InversionThreshParameters";
    8988                case InversionVelObsEnum : return "InversionVelObs";
    9089                case InversionVxObsEnum : return "InversionVxObs";
  • issm/branches/trunk-jpl-damage/src/c/modules/ModelProcessorx/Control/UpdateElementsAndMaterialsControl.cpp

    r11984 r12872  
    4040        iomodel->FetchDataToInput(elements,InversionThicknessObsEnum);
    4141
    42         iomodel->FetchData(5,InversionControlParametersEnum,InversionCostFunctionsCoefficientsEnum,InversionMinParametersEnum,InversionMaxParametersEnum,InversionThreshParametersEnum);
     42        iomodel->FetchData(4,InversionControlParametersEnum,InversionCostFunctionsCoefficientsEnum,InversionMinParametersEnum,InversionMaxParametersEnum);
    4343        for(i=0;i<num_control_type;i++){
    4444                switch((int)iomodel->Data(InversionControlParametersEnum)[i]){
     
    6767       
    6868        /*Free data: */
    69         iomodel->DeleteData(1+5+6,MeshElementsEnum,InversionControlParametersEnum,InversionCostFunctionsCoefficientsEnum,InversionMinParametersEnum,InversionMaxParametersEnum,InversionThreshParametersEnum,BalancethicknessThickeningRateEnum,VxEnum,VyEnum,FrictionCoefficientEnum,MaterialsRheologyBEnum,MaterialsRheologyZEnum);
     69        iomodel->DeleteData(1+4+6,MeshElementsEnum,InversionControlParametersEnum,InversionCostFunctionsCoefficientsEnum,InversionMinParametersEnum,InversionMaxParametersEnum,BalancethicknessThickeningRateEnum,VxEnum,VyEnum,FrictionCoefficientEnum,MaterialsRheologyBEnum,MaterialsRheologyZEnum);
    7070}
  • issm/branches/trunk-jpl-damage/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp

    r12274 r12872  
    4949        if(dim==3)iomodel->FetchData(2,MeshUpperelementsEnum,MeshLowerelementsEnum);
    5050        #endif
    51         if(control_analysis)iomodel->FetchData(4,InversionControlParametersEnum,InversionMinParametersEnum,InversionMaxParametersEnum,InversionThreshParametersEnum);
     51        if(control_analysis)iomodel->FetchData(3,InversionControlParametersEnum,InversionMinParametersEnum,InversionMaxParametersEnum);
    5252       
    5353        /*Create elements and materials: */
     
    6767       
    6868        /*Free data: */
    69         iomodel->DeleteData(11,MeshElementsEnum,MeshElementconnectivityEnum,MeshUpperelementsEnum,MeshLowerelementsEnum,
     69        iomodel->DeleteData(10,MeshElementsEnum,MeshElementconnectivityEnum,MeshUpperelementsEnum,MeshLowerelementsEnum,
    7070                                MaterialsRheologyBEnum,MaterialsRheologyNEnum,MaterialsRheologyZEnum,InversionControlParametersEnum,InversionMinParametersEnum,
    71                                 InversionMaxParametersEnum,InversionThreshParametersEnum);
     71                                InversionMaxParametersEnum);
    7272
    7373        /*Add new constant material property to materials, at the end: */
  • issm/branches/trunk-jpl-damage/src/c/modules/StringToEnumx/StringToEnumx.cpp

    r12280 r12872  
    8787              else if (strcmp(name,"InversionStepThreshold")==0) return InversionStepThresholdEnum;
    8888              else if (strcmp(name,"InversionThicknessObs")==0) return InversionThicknessObsEnum;
    89               else if (strcmp(name,"InversionThreshParameters")==0) return InversionThreshParametersEnum;
    9089              else if (strcmp(name,"InversionVelObs")==0) return InversionVelObsEnum;
    9190              else if (strcmp(name,"InversionVxObs")==0) return InversionVxObsEnum;
     
    137136              else if (strcmp(name,"MiscellaneousName")==0) return MiscellaneousNameEnum;
    138137              else if (strcmp(name,"PrognosticHydrostaticAdjustment")==0) return PrognosticHydrostaticAdjustmentEnum;
     138              else if (strcmp(name,"PrognosticMinThickness")==0) return PrognosticMinThicknessEnum;
    139139         else stage=2;
    140140   }
    141141   if(stage==2){
    142               if (strcmp(name,"PrognosticMinThickness")==0) return PrognosticMinThicknessEnum;
    143               else if (strcmp(name,"PrognosticPenaltyFactor")==0) return PrognosticPenaltyFactorEnum;
     142              if (strcmp(name,"PrognosticPenaltyFactor")==0) return PrognosticPenaltyFactorEnum;
    144143              else if (strcmp(name,"PrognosticSpcthickness")==0) return PrognosticSpcthicknessEnum;
    145144              else if (strcmp(name,"PrognosticStabilization")==0) return PrognosticStabilizationEnum;
     
    260259              else if (strcmp(name,"Input")==0) return InputEnum;
    261260              else if (strcmp(name,"IntInput")==0) return IntInputEnum;
     261              else if (strcmp(name,"IntParam")==0) return IntParamEnum;
    262262         else stage=3;
    263263   }
    264264   if(stage==3){
    265               if (strcmp(name,"IntParam")==0) return IntParamEnum;
    266               else if (strcmp(name,"IntVecParam")==0) return IntVecParamEnum;
     265              if (strcmp(name,"IntVecParam")==0) return IntVecParamEnum;
    267266              else if (strcmp(name,"MacAyeal2dIceFront")==0) return MacAyeal2dIceFrontEnum;
    268267              else if (strcmp(name,"MacAyeal3dIceFront")==0) return MacAyeal3dIceFrontEnum;
     
    383382              else if (strcmp(name,"StressTensoryz")==0) return StressTensoryzEnum;
    384383              else if (strcmp(name,"StressTensorzz")==0) return StressTensorzzEnum;
     384              else if (strcmp(name,"IceVolume")==0) return IceVolumeEnum;
    385385         else stage=4;
    386386   }
    387387   if(stage==4){
    388               if (strcmp(name,"IceVolume")==0) return IceVolumeEnum;
    389               else if (strcmp(name,"P0")==0) return P0Enum;
     388              if (strcmp(name,"P0")==0) return P0Enum;
    390389              else if (strcmp(name,"P1")==0) return P1Enum;
    391390              else if (strcmp(name,"P1DG")==0) return P1DGEnum;
  • issm/branches/trunk-jpl-damage/src/c/objects/Elements/Penta.cpp

    r12278 r12872  
    15331533        double  cmmininputs[6];
    15341534        double  cmmaxinputs[6];
    1535         double  cmthreshinputs[6];
    15361535
    15371536        double  yts;
     
    15661565                                                for(j=0;j<6;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    15671566                                                for(j=0;j<6;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    1568                                                 for(j=0;j<6;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    1569                                                 this->inputs->AddInput(new ControlInput(BalancethicknessThickeningRateEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     1567                                                this->inputs->AddInput(new ControlInput(BalancethicknessThickeningRateEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    15701568                                        }
    15711569                                        break;
     
    15751573                                                for(j=0;j<6;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    15761574                                                for(j=0;j<6;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    1577                                                 for(j=0;j<6;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    1578                                                 this->inputs->AddInput(new ControlInput(VxEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     1575                                                this->inputs->AddInput(new ControlInput(VxEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    15791576                                        }
    15801577                                        break;
     
    15841581                                                for(j=0;j<6;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    15851582                                                for(j=0;j<6;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    1586                                                 for(j=0;j<6;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    1587                                                 this->inputs->AddInput(new ControlInput(VyEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     1583                                                this->inputs->AddInput(new ControlInput(VyEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    15881584                                        }
    15891585                                        break;
     
    15931589                                                for(j=0;j<6;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i];
    15941590                                                for(j=0;j<6;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i];
    1595                                                 for(j=0;j<6;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[(penta_vertex_ids[j]-1)*num_control_type+i]/yts;
    1596                                                 this->inputs->AddInput(new ControlInput(FrictionCoefficientEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     1591                                                this->inputs->AddInput(new ControlInput(FrictionCoefficientEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    15971592                                        }
    15981593                                        break;
  • issm/branches/trunk-jpl-damage/src/c/objects/Elements/Tria.cpp

    r12280 r12872  
    13771377        double cmmininputs[3];
    13781378        double cmmaxinputs[3];
    1379         double cmthreshinputs[3];
    13801379        bool   control_analysis=false;
    13811380        int    num_control_type;
     
    14041403                                                for(j=0;j<3;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    14051404                                                for(j=0;j<3;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    1406                                                 for(j=0;j<3;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    1407                                                 this->inputs->AddInput(new ControlInput(BalancethicknessThickeningRateEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     1405                                                this->inputs->AddInput(new ControlInput(BalancethicknessThickeningRateEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    14081406                                        }
    14091407                                        break;
     
    14131411                                                for(j=0;j<3;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    14141412                                                for(j=0;j<3;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    1415                                                 for(j=0;j<3;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    1416                                                 this->inputs->AddInput(new ControlInput(VxEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     1413                                                this->inputs->AddInput(new ControlInput(VxEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    14171414                                        }
    14181415                                        break;
     
    14221419                                                for(j=0;j<3;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    14231420                                                for(j=0;j<3;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    1424                                                 for(j=0;j<3;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    1425                                                 this->inputs->AddInput(new ControlInput(VyEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     1421                                                this->inputs->AddInput(new ControlInput(VyEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    14261422                                        }
    14271423                                        break;
     
    14311427                                                for(j=0;j<3;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i];
    14321428                                                for(j=0;j<3;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i];
    1433                                                 for(j=0;j<3;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[(tria_vertex_ids[j]-1)*num_control_type+i]/yts;
    1434                                                 this->inputs->AddInput(new ControlInput(FrictionCoefficientEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     1429                                                this->inputs->AddInput(new ControlInput(FrictionCoefficientEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    14351430                                        }
    14361431                                        break;
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/BoolInput.cpp

    r12168 r12872  
    167167/*}}}*/
    168168/*FUNCTION BoolInput::Constrain{{{1*/
    169 void BoolInput::Constrain(double cm_min, double cm_max, double cm_thresh){
     169void BoolInput::Constrain(double cm_min, double cm_max){
    170170
    171171        if(!isnan(cm_min)) if (this->value<cm_min)this->value=cm_min;
    172172        if(!isnan(cm_max)) if (this->value>cm_max)this->value=cm_max;
    173         if(!isnan(cm_thresh)) if (this->value<cm_thresh)this->value=cm_min;
    174173
    175174}
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/BoolInput.h

    r12285 r12872  
    4040                Input* PointwiseMin(Input* inputB){_error_("not implemented yet");};
    4141                Input* PointwiseMax(Input* inputB){_error_("not implemented yet");};
    42                 Input* PointwiseThresh(Input* inputB,Input* inputC){_error_("not implemented yet");};
    4342                ElementResult* SpawnResult(int step, double time);
    4443                void Configure(Parameters* parameters);
     
    7473                void ArtificialNoise(double min,double max){_error_("not implemented yet");};
    7574                void AXPY(Input* xinput,double scalar);
    76                 void Constrain(double cm_min, double cm_max, double cm_thresh);
     75                void Constrain(double cm_min, double cm_max);
    7776                void Extrude(void);
    7877                void VerticallyIntegrate(Input* thickness_input){_error_("not supported yet");};
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/ControlInput.cpp

    r12285 r12872  
    2525        minvalues   = NULL;
    2626        maxvalues   = NULL;
    27         threshvalues   = NULL;
    2827        gradient    = NULL;
    2928}
    3029/*}}}*/
    31 /*FUNCTION ControlInput::ControlInput(int enum_type,int enum_input,double* pvalues,double* pmin,double* pmax,double* pthresh,int id){{{1*/
    32 ControlInput::ControlInput(int in_enum_type,int enum_input,double* pvalues,double* pmin,double* pmax,double* pthresh,int id){
     30/*FUNCTION ControlInput::ControlInput(int enum_type,int enum_input,double* pvalues,double* pmin,double* pmax,int id){{{1*/
     31ControlInput::ControlInput(int in_enum_type,int enum_input,double* pvalues,double* pmin,double* pmax,int id){
    3332
    3433        control_id=id;
     
    4140                        minvalues  =new TriaP1Input(enum_type,pmin);
    4241                        maxvalues  =new TriaP1Input(enum_type,pmax);
    43                         threshvalues  =new TriaP1Input(enum_type,pthresh);
    4442                        break;
    4543                case PentaP1InputEnum:
     
    4846                        minvalues  =new PentaP1Input(enum_type,pmin);
    4947                        maxvalues  =new PentaP1Input(enum_type,pmax);
    50                         threshvalues  =new TriaP1Input(enum_type,pthresh);
    5148                        break;
    5249                default:
     
    6259        delete minvalues;
    6360        delete maxvalues;
    64         delete threshvalues;
    6561        delete gradient;
    6662}
     
    8278        printf("---minvalues: \n");  if (minvalues)   minvalues->Echo();
    8379        printf("---maxvalues: \n");  if (maxvalues)   maxvalues->Echo();
    84         printf("---threshvalues: \n");  if (threshvalues)   threshvalues->Echo();
    8580        printf("---gradient: \n");   if (gradient)    gradient->Echo();
    8681}
     
    160155        }
    161156
    162         /*marshal threshvalues*/
    163         if(!threshvalues){
    164                 flag=0;
    165                 memcpy(marshalled_dataset,&flag,sizeof(flag));marshalled_dataset+=sizeof(flag);
    166         }
    167         else{
    168                 flag=1;
    169                 memcpy(marshalled_dataset,&flag,sizeof(flag));marshalled_dataset+=sizeof(flag);
    170                 this->threshvalues->Marshall(&marshalled_dataset);
    171         }
    172 
    173157        /*marshal gradient*/
    174158        if(!gradient){
     
    200184        if(minvalues)size+=minvalues->MarshallSize();
    201185        if(maxvalues)size+=maxvalues->MarshallSize();
    202         if(threshvalues)size+=threshvalues->MarshallSize();
    203186        if(gradient)   size+=gradient->MarshallSize();
    204187        return size;
     
    291274        }
    292275
    293         /*Demarshal threshvalues*/
    294         memcpy(&flag,marshalled_dataset,sizeof(int));marshalled_dataset+=sizeof(int);
    295         if(flag){
    296                 memcpy(&input_enum_type,marshalled_dataset,sizeof(int)); marshalled_dataset+=sizeof(int);
    297                 if(input_enum_type==PentaP1InputEnum){
    298                         threshvalues=new PentaP1Input();
    299                         threshvalues->Demarshall(&marshalled_dataset);
    300                 }
    301                 else if(input_enum_type==TriaP1InputEnum){
    302                         threshvalues=new TriaP1Input();
    303                         threshvalues->Demarshall(&marshalled_dataset);
    304                 }
    305                 else _error_("Not supported yet");
    306         }
    307         else{
    308                 threshvalues=NULL;
    309         }
    310 
    311276        /*Demarshal gradient*/
    312277        memcpy(&flag,marshalled_dataset,sizeof(int));marshalled_dataset+=sizeof(int);
     
    353318        if(minvalues)   output->minvalues=(Input*)this->minvalues->copy();
    354319        if(maxvalues)   output->maxvalues=(Input*)this->maxvalues->copy();
    355         if(threshvalues)   output->threshvalues=(Input*)this->threshvalues->copy();
    356320        if(gradient)    output->gradient=(Input*)this->gradient->copy();
    357321
     
    379343        newvalues=this->values->PointwiseMax(minvalues);
    380344        delete values; this->values=newvalues;
    381         newvalues=this->values->PointwiseThresh(threshvalues,minvalues);
    382         delete values; this->values=newvalues;
    383 }/*}}}*/
    384 /*FUNCTION ControlInput::Constrain(double min, double max, double thresh){{{1*/
    385 void ControlInput::Constrain(double min, double max, double thresh){
    386            values->Constrain(min,max,thresh);
     345}/*}}}*/
     346/*FUNCTION ControlInput::Constrain(double min, double max){{{1*/
     347void ControlInput::Constrain(double min, double max){
     348           values->Constrain(min,max);
    387349}/*}}}*/
    388350/*FUNCTION ControlInput::Extrude{{{1*/
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/ControlInput.h

    r12285 r12872  
    2323                Input* minvalues;
    2424                Input* maxvalues;
    25                 Input* threshvalues;
    2625                Input* gradient;
    2726
    2827                /*ControlInput constructors, destructors: {{{1*/
    2928                ControlInput();
    30                 ControlInput(int enum_type,int enum_input,double* pvalues,double* pmin,double* pmax,double* pthresh,int id);
     29                ControlInput(int enum_type,int enum_input,double* pvalues,double* pmin,double* pmax,int id);
    3130                ~ControlInput();
    3231                /*}}}*/
     
    4544                Input* PointwiseMin(Input* inputB){_error_("not implemented yet");};
    4645                Input* PointwiseMax(Input* inputB){_error_("not implemented yet");};
    47                 Input* PointwiseThresh(Input* inputB,Input* inputC){_error_("not implemented yet");};
    4846                ElementResult* SpawnResult(int step, double time);
    4947                void AddTimeValues(double* values,int step,double time){_error_("not supported yet");};
     
    7674                void AXPY(Input* xinput,double scalar){_error_("not implemented yet");};
    7775                void Constrain(void);
    78                 void Constrain(double min,double max, double thresh);
     76                void Constrain(double min,double max);
    7977                double InfinityNorm(void){_error_("not implemented yet");};
    8078                double Max(void){_error_("not implemented yet");};
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/DatasetInput.h

    r12285 r12872  
    4040                Input* PointwiseMin(Input* inputB){_error_("not implemented yet");};
    4141                Input* PointwiseMax(Input* inputB){_error_("not implemented yet");};
    42                 Input* PointwiseThresh(Input* inputB,Input* inputC){_error_("not implemented yet");};
    4342                ElementResult* SpawnResult(int step, double time){_error_("not implemented yet");};
    4443                void AddTimeValues(double* values,int step,double time){_error_("not supported yet");};
     
    7069                void AXPY(Input* xinput,double scalar){_error_("not implemented yet");};
    7170                void Constrain(void){_error_("not implemented yet");};
    72                 void Constrain(double min,double max, double thresh){_error_("not implemented yet");};
     71                void Constrain(double min,double max){_error_("not implemented yet");};
    7372                double InfinityNorm(void){_error_("not implemented yet");};
    7473                double Max(void){_error_("not implemented yet");};
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/DoubleInput.cpp

    r12168 r12872  
    223223/*}}}*/
    224224/*FUNCTION DoubleInput::Constrain{{{1*/
    225 void DoubleInput::Constrain(double cm_min, double cm_max, double cm_thresh){
     225void DoubleInput::Constrain(double cm_min, double cm_max){
    226226
    227227        if(!isnan(cm_min)) if (this->value<cm_min)this->value=cm_min;
    228228        if(!isnan(cm_max)) if (this->value>cm_max)this->value=cm_max;
    229         if(!isnan(cm_thresh)) if (this->value<cm_thresh)this->value=cm_min;
    230229
    231230}
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/DoubleInput.h

    r12285 r12872  
    3939                Input* PointwiseMin(Input* inputB);
    4040                Input* PointwiseMax(Input* inputB);
    41                 Input* PointwiseThresh(Input* inputB,Input* inputC){_error_("not supported yet");};
    4241                ElementResult* SpawnResult(int step, double time);
    4342                void AddTimeValues(double* values,int step,double time){_error_("not supported yet");};
     
    6867                void ArtificialNoise(double min,double max){_error_("not implemented yet");};
    6968                void AXPY(Input* xinput,double scalar);
    70                 void Constrain(double cm_min, double cm_max, double cm_thresh);
     69                void Constrain(double cm_min, double cm_max);
    7170                double InfinityNorm(void){_error_("not implemented yet");};
    7271                double Max(void);
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/Input.h

    r12285 r12872  
    5353                virtual void   ArtificialNoise(double min,double max)=0;
    5454                virtual void   AXPY(Input* xinput,double scalar)=0;
    55                 virtual void   Constrain(double cm_min, double cm_max, double cm_thresh)=0;
     55                virtual void   Constrain(double cm_min, double cm_max)=0;
    5656                virtual void   VerticallyIntegrate(Input* thickness_input)=0;
    5757                virtual void   Extrude()=0;
     
    6363                virtual Input* PointwiseMax(Input* inputmax)=0;
    6464                virtual Input* PointwiseMin(Input* inputmin)=0;
    65                 virtual Input* PointwiseThresh(Input* inputthresh,Input* inputmin)=0;
    6665                virtual ElementResult* SpawnResult(int step, double time)=0;
    6766
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/IntInput.cpp

    r12168 r12872  
    172172/*}}}*/
    173173/*FUNCTION IntInput::Constrain{{{1*/
    174 void IntInput::Constrain(double cm_min, double cm_max, double cm_thresh){
     174void IntInput::Constrain(double cm_min, double cm_max){
    175175
    176176        if(!isnan(cm_min)) if (this->value<cm_min)this->value=(int)cm_min;
    177177        if(!isnan(cm_max)) if (this->value>cm_max)this->value=(int)cm_max;
    178         if(!isnan(cm_thresh)) if (this->value<cm_thresh)this->value=(int)cm_min;
    179178
    180179}
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/IntInput.h

    r12285 r12872  
    4040                Input* PointwiseMin(Input* inputB){_error_("not implemented yet");};
    4141                Input* PointwiseMax(Input* inputB){_error_("not implemented yet");};
    42                 Input* PointwiseThresh(Input* inputB,Input* inputC){_error_("not implemented yet");};
    4342                ElementResult* SpawnResult(int step, double time);
    4443                void AddTimeValues(double* values,int step,double time){_error_("not supported yet");};
     
    6968                void ArtificialNoise(double min,double max){_error_("not implemented yet");};
    7069                void AXPY(Input* xinput,double scalar);
    71                 void Constrain(double cm_min, double cm_max, double cm_thresh);
     70                void Constrain(double cm_min, double cm_max);
    7271                double InfinityNorm(void){_error_("InfinityNorm not implemented for integers");};
    7372                double Max(void){_error_("Max not implemented for integers");};
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/PentaP1Input.cpp

    r12285 r12872  
    466466/*}}}*/
    467467/*FUNCTION PentaP1Input::Constrain{{{1*/
    468 void PentaP1Input::Constrain(double cm_min, double cm_max, double cm_thresh){
     468void PentaP1Input::Constrain(double cm_min, double cm_max){
    469469
    470470        int i;
     
    473473        if(!isnan(cm_min)) for(i=0;i<numnodes;i++)if (this->values[i]<cm_min)this->values[i]=cm_min;
    474474        if(!isnan(cm_max)) for(i=0;i<numnodes;i++)if (this->values[i]>cm_max)this->values[i]=cm_max;
    475         if(!isnan(cm_thresh)) for(i=0;i<numnodes;i++)if (this->values[i]<cm_thresh)this->values[i]=cm_min;
    476475
    477476}
     
    609608}
    610609/*}}}*/
    611 /*FUNCTION PentaP1Input::PointwiseThresh{{{1*/
    612 Input* PentaP1Input::PointwiseThresh(Input* inputB,Input* inputC){
    613 
    614         /*Ouput*/
    615         PentaP1Input* outinput=NULL;
    616 
    617         /*Intermediaries*/
    618         int               i;
    619         PentaP1Input *xinputB     = NULL;
    620         PentaP1Input *xinputC     = NULL;
    621         int               B_numvalues;
    622         const int         numnodes    = 6;
    623         double            threshvalues[numnodes];
    624 
    625         /*Check that inputB is of the same type*/
    626         if (inputB->ObjectEnum()!=PentaP1InputEnum) _error_("Operation not permitted because inputB is of type %s",EnumToStringx(inputB->ObjectEnum()));
    627         xinputB=(PentaP1Input*)inputB;
    628 
    629         /*Check that inputB is of the same type*/
    630         if (inputC->ObjectEnum()!=PentaP1InputEnum) _error_("Operation not permitted because inputC is of type %s",EnumToStringx(inputC->ObjectEnum()));
    631         xinputC=(PentaP1Input*)inputC;
    632 
    633         /*Create point wise thresh*/
    634         for(i=0;i<numnodes;i++){
    635                 if(this->values[i] < xinputB->values[i]) threshvalues[i]=xinputC->values[i];
    636                 else threshvalues[i]=this->values[i];
    637         }
    638 
    639         /*Create new Penta vertex input (copy of current input)*/
    640         outinput=new PentaP1Input(this->enum_type,&threshvalues[0]);
    641 
    642         /*Return output pointer*/
    643         return outinput;
    644 
    645 }
    646 /*}}}*/
    647610/*FUNCTION PentaP1Input::GetVectorFromInputs{{{1*/
    648611void PentaP1Input::GetVectorFromInputs(Vector* vector,int* doflist){
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/PentaP1Input.h

    r12285 r12872  
    4040                Input* PointwiseMin(Input* inputB);
    4141                Input* PointwiseMax(Input* inputB);
    42                 Input* PointwiseThresh(Input* inputB,Input* inputC);
    4342                ElementResult* SpawnResult(int step, double time);
    4443                void AddTimeValues(double* values,int step,double time){_error_("not supported yet");};
     
    7069                void ArtificialNoise(double min,double max){_error_("not implemented yet");};
    7170                void AXPY(Input* xinput,double scalar);
    72                 void Constrain(double cm_min, double cm_max, double cm_thresh);
     71                void Constrain(double cm_min, double cm_max);
    7372                double InfinityNorm(void);
    7473                double Max(void);
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/TransientInput.h

    r12285 r12872  
    4343                Input* PointwiseMin(Input* forcingB){_error_("not implemented yet");};
    4444                Input* PointwiseMax(Input* forcingB){_error_("not implemented yet");};
    45                 Input* PointwiseThresh(Input* forcingB,Input* forcingC){_error_("not implemented yet");};
    4645                ElementResult* SpawnResult(int step, double time);
    4746                void Configure(Parameters* parameters);
     
    7271                void ArtificialNoise(double min,double max){_error_("not implemented yet");};
    7372                void AXPY(Input* xforcing,double scalar){_error_("not implemented yet");};
    74                 void Constrain(double cm_min, double cm_max, double cm_thresh){_error_("not implemented yet");};
     73                void Constrain(double cm_min, double cm_max){_error_("not implemented yet");};
    7574                double InfinityNorm(void);
    7675                double Max(void);
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/TriaP1Input.cpp

    r12285 r12872  
    329329/*}}}*/
    330330/*FUNCTION TriaP1Input::Constrain{{{1*/
    331 void TriaP1Input::Constrain(double cm_min, double cm_max, double cm_thresh){
     331void TriaP1Input::Constrain(double cm_min, double cm_max){
    332332
    333333        int i;
     
    336336        if(!isnan(cm_min)) for(i=0;i<numnodes;i++)if (this->values[i]<cm_min)this->values[i]=cm_min;
    337337        if(!isnan(cm_max)) for(i=0;i<numnodes;i++)if (this->values[i]>cm_max)this->values[i]=cm_max;
    338         if(!isnan(cm_thresh)) for(i=0;i<numnodes;i++)if (this->values[i]<cm_thresh)this->values[i]=cm_min;
    339338
    340339}
     
    417416}
    418417/*}}}*/
    419 /*FUNCTION TriaP1Input::PointwiseThresh{{{1*/
    420 Input* TriaP1Input::PointwiseThresh(Input* inputB, Input* inputC){
    421         /*pointwise comparison to inputB.  If less than inputB, set nodal output to input C*/
    422 
    423         /*Ouput*/
    424         TriaP1Input* outinput=NULL;
    425 
    426         /*Intermediaries*/
    427         int               i;
    428         TriaP1Input *xinputB     = NULL;
    429         TriaP1Input *xinputC     = NULL;
    430         int               B_numvalues; /* doesn't look like this is used */
    431         const int         numnodes    = 3;
    432         double            threshvalues[numnodes];
    433 
    434         /*Check that inputB is of the same type*/
    435         if (inputB->ObjectEnum()!=TriaP1InputEnum) _error_("Operation not permitted because inputB is of type %s",EnumToStringx(inputB->ObjectEnum()));
    436         xinputB=(TriaP1Input*)inputB;
    437 
    438         /*Check that inputC is of the same type*/
    439         if (inputC->ObjectEnum()!=TriaP1InputEnum) _error_("Operation not permitted because inputC is of type %s",EnumToStringx(inputC->ObjectEnum()));
    440         xinputC=(TriaP1Input*)inputC;
    441 
    442         /*Create point wise thresh*/
    443         for(i=0;i<numnodes;i++){
    444                 if(this->values[i] < xinputB->values[i]) threshvalues[i]=xinputC->values[i];
    445                 else threshvalues[i]=this->values[i];
    446         }
    447 
    448         /*Create new Tria vertex input (copy of current input)*/
    449         outinput=new TriaP1Input(this->enum_type,&threshvalues[0]);
    450 
    451         /*Return output pointer*/
    452         return outinput;
    453 
    454 }
    455 /*}}}*/
    456418/*FUNCTION TriaP1Input::Configure{{{1*/
    457419void TriaP1Input::Configure(Parameters* parameters){
  • issm/branches/trunk-jpl-damage/src/c/objects/Inputs/TriaP1Input.h

    r12285 r12872  
    4040                Input* PointwiseMin(Input* inputB);
    4141                Input* PointwiseMax(Input* inputB);
    42                 Input* PointwiseThresh(Input* inputB, Input* inputC);
    4342                ElementResult* SpawnResult(int step, double time);
    4443                void AddTimeValues(double* values,int step,double time){_error_("not supported yet");};
     
    7069                void ArtificialNoise(double min,double max);
    7170                void AXPY(Input* xinput,double scalar);
    72                 void Constrain(double cm_min, double cm_max, double cm_thresh);
     71                void Constrain(double cm_min, double cm_max);
    7372                double InfinityNorm(void);
    7473                double Max(void);
  • issm/branches/trunk-jpl-damage/src/c/objects/Materials/Matice.cpp

    r12275 r12872  
    758758                double    cmmininputs[num_vertices];
    759759                double    cmmaxinputs[num_vertices];
    760                 double    cmthreshinputs[num_vertices];
    761760
    762761                /*Get B*/
     
    789788                                                        for(j=0;j<num_vertices;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    790789                                                        for(j=0;j<num_vertices;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    791                                                         for(j=0;j<num_vertices;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    792                                                         this->inputs->AddInput(new ControlInput(MaterialsRheologyBbarEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     790                                                        this->inputs->AddInput(new ControlInput(MaterialsRheologyBbarEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    793791                                                }
    794792                                                break;
     
    799797                                                        for(j=0;j<num_vertices;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    800798                                                        for(j=0;j<num_vertices;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    801                                                         for(j=0;j<num_vertices;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    802                                                         this->inputs->AddInput(new ControlInput(MaterialsRheologyZbarEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     799                                                        this->inputs->AddInput(new ControlInput(MaterialsRheologyZbarEnum,TriaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    803800                                                }
    804801                                                break;
     
    818815                double    cmmininputs[num_vertices];
    819816                double    cmmaxinputs[num_vertices];
    820                 double    cmthreshinputs[num_vertices];
    821817
    822818                /*Get B*/
     
    849845                                                        for(j=0;j<num_vertices;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    850846                                                        for(j=0;j<num_vertices;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    851                                                         for(j=0;j<num_vertices;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    852                                                         this->inputs->AddInput(new ControlInput(MaterialsRheologyBEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     847                                                        this->inputs->AddInput(new ControlInput(MaterialsRheologyBEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    853848                                                }
    854849                                                break;
     
    859854                                                        for(j=0;j<num_vertices;j++)cmmininputs[j]=iomodel->Data(InversionMinParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    860855                                                        for(j=0;j<num_vertices;j++)cmmaxinputs[j]=iomodel->Data(InversionMaxParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    861                                                         for(j=0;j<num_vertices;j++)cmthreshinputs[j]=iomodel->Data(InversionThreshParametersEnum)[int(iomodel->Data(MeshElementsEnum)[num_vertices*index+j]-1)*num_control_type+i];
    862                                                         this->inputs->AddInput(new ControlInput(MaterialsRheologyZEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,cmthreshinputs,i+1));
     856                                                        this->inputs->AddInput(new ControlInput(MaterialsRheologyZEnum,PentaP1InputEnum,nodeinputs,cmmininputs,cmmaxinputs,i+1));
    863857                                                }
    864858                                                break;
  • issm/branches/trunk-jpl-damage/src/m/classes/inversion.m

    r12286 r12872  
    1818                min_parameters              = NaN
    1919                max_parameters              = NaN
    20                 thresh_parameters           = NaN
    2120                step_threshold              = NaN
    2221                gradient_only               = 0
     
    9796                        checkfield(md,'inversion.min_parameters','size',[md.mesh.numberofvertices num_controls]);
    9897                        checkfield(md,'inversion.max_parameters','size',[md.mesh.numberofvertices num_controls]);
    99                         checkfield(md,'inversion.thresh_parameters','size',[md.mesh.numberofvertices num_controls]);
    10098
    10199                        if solution==BalancethicknessSolutionEnum
     
    119117                        fielddisplay(obj,'min_parameters','absolute minimum acceptable value of the inversed parameter on each vertex');
    120118                        fielddisplay(obj,'max_parameters','absolute maximum acceptable value of the inversed parameter on each vertex');
    121                         fielddisplay(obj,'thresh_parameters','threshold value of the inversed parameter, below which it jumps to min_parameters. Set to min_parameters to disable');
    122119                        fielddisplay(obj,'gradient_only','stop control method solution at gradient');
    123120                        fielddisplay(obj,'vx_obs','observed velocity x component [m/a]');
     
    149146                        WriteData(fid,'object',obj,'fieldname','min_parameters','format','DoubleMat','mattype',3);
    150147                        WriteData(fid,'object',obj,'fieldname','max_parameters','format','DoubleMat','mattype',3);
    151                         WriteData(fid,'object',obj,'fieldname','thresh_parameters','format','DoubleMat','mattype',3);
    152148                        WriteData(fid,'object',obj,'fieldname','step_threshold','format','DoubleMat','mattype',3);
    153149                        WriteData(fid,'object',obj,'fieldname','gradient_only','format','Boolean');
  • issm/branches/trunk-jpl-damage/src/m/classes/model/model.m

    r12168 r12872  
    189189                         if isfield(structmd,'cm_min'), md.inversion.min_parameters=structmd.cm_min; end
    190190                         if isfield(structmd,'cm_max'), md.inversion.max_parameters=structmd.cm_max; end
    191                          if isfield(structmd,'cm_thresh'), md.inversion.thresh_parameters=structmd.cm_thresh; end
    192191                         if isfield(structmd,'vx_obs'), md.inversion.vx_obs=structmd.vx_obs; end
    193192                         if isfield(structmd,'vy_obs'), md.inversion.vy_obs=structmd.vy_obs; end
  • issm/branches/trunk-jpl-damage/src/m/enum/MaximumNumberOfEnums.m

    r12280 r12872  
    99%      macro=MaximumNumberOfEnums()
    1010
    11 macro=440;
     11macro=439;
  • issm/branches/trunk-jpl-damage/src/m/model/parameterization/parametercontrolZ.m

    r11990 r12872  
    5757else
    5858        md.inversion.max_parameters=cm_max;
    59 end
    60 
    61 %cm_thresh
    62 cm_thresh=getfieldvalue(options,'cm_thresh',1e-4*ones(md.mesh.numberofvertices,1));
    63 if (length(cm_thresh)==1)
    64         md.inversion.thresh_parameters=cm_thresh*ones(md.mesh.numberofvertices,1);
    65 elseif (length(cm_thresh)==md.mesh.numberofvertices)
    66         md.inversion.thresh_parameters=cm_thresh;
    67 else
    68         md.inversion.thresh_parameters=cm_thresh;
    6959end
    7060
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test224.m

    r12276 r12872  
    99md.inversion.min_parameters=10^6*ones(md.mesh.numberofvertices,1);
    1010md.inversion.max_parameters=2*10^9*ones(md.mesh.numberofvertices,1);
    11 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1211md.inversion.nsteps=2;
    1312md.inversion.cost_functions=101*ones(md.inversion.nsteps,1);
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test226.m

    r12276 r12872  
    1010md.inversion.min_parameters=10^6*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=2*10^9*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=101*ones(md.inversion.nsteps,1);
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test228.m

    r12276 r12872  
    1010md.inversion.min_parameters=10^6*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=2*10^9*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=101*ones(md.inversion.nsteps,1);
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test230.m

    r12276 r12872  
    1010md.inversion.min_parameters=10^6*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=2*10^9*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=101*ones(md.inversion.nsteps,1);
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test272.m

    r12259 r12872  
    1010md.inversion.min_parameters=10^-13*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=101*ones(md.inversion.nsteps,1);
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test338.m

    r12276 r12872  
    99md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
    1010md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
    11 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1211md.inversion.nsteps=2;
    1312md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1)  501*ones(md.inversion.nsteps,1)];
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test340.m

    r12276 r12872  
    1010md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1)  501*ones(md.inversion.nsteps,1)];
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test342.m

    r12276 r12872  
    1010md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=[102*ones(md.inversion.nsteps,1)  501*ones(md.inversion.nsteps,1)];
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test344.m

    r12276 r12872  
    1010md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=[104*ones(md.inversion.nsteps,1)  501*ones(md.inversion.nsteps,1)];
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test430.m

    r12276 r12872  
    1010md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1)  501*ones(md.inversion.nsteps,1)];
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test432.m

    r12276 r12872  
    1010md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=[102*ones(md.inversion.nsteps,1)  501*ones(md.inversion.nsteps,1)];
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test522.m

    r12276 r12872  
    1515md.inversion.min_parameters=10^6*ones(md.mesh.numberofvertices,1);
    1616md.inversion.max_parameters=2*10^9*ones(md.mesh.numberofvertices,1);
    17 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1817md.inversion.nsteps=2;
    1918md.inversion.cost_functions=101*ones(md.inversion.nsteps,1);
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test524.m

    r12276 r12872  
    1010md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1)  501*ones(md.inversion.nsteps,1)];
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test526.m

    r12276 r12872  
    1010md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
    1111md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
    12 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1312md.inversion.nsteps=2;
    1413md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1)  501*ones(md.inversion.nsteps,1)];
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test622.m

    r12276 r12872  
    1212md.inversion.min_parameters=-50*ones(md.mesh.numberofvertices,1);
    1313md.inversion.max_parameters=50*ones(md.mesh.numberofvertices,1);
    14 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1514md.inversion.cost_functions=201*ones(md.inversion.nsteps,1);
    1615md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1);
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test624.m

    r12276 r12872  
    1313md.inversion.min_parameters=-50*ones(md.mesh.numberofvertices,1);
    1414md.inversion.max_parameters=50*ones(md.mesh.numberofvertices,1);
    15 md.inversion.thresh_parameters=md.inversion.min_parameters;
    1615md.inversion.cost_functions=201*ones(md.inversion.nsteps,1);
    1716md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1);
  • issm/branches/trunk-jpl-damage/test/NightlyRun/test626.m

    r12276 r12872  
    2020md.inversion.min_parameters=[-2000*ones(md.mesh.numberofvertices,1) -2000*ones(md.mesh.numberofvertices,1)];
    2121md.inversion.max_parameters=[+2000*ones(md.mesh.numberofvertices,1) +2000*ones(md.mesh.numberofvertices,1)];
    22 md.inversion.thresh_parameters=md.inversion.min_parameters;
    2322md.inversion.cost_functions=201*ones(md.inversion.nsteps,1);
    2423md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1);
Note: See TracChangeset for help on using the changeset viewer.