Changeset 23588


Ignore:
Timestamp:
01/02/19 20:42:58 (6 years ago)
Author:
Mathieu Morlighem
Message:

CHG: removing IsAnalysis now that each constraint, node and load is in an analysis specific dataset

Location:
issm/trunk-jpl/src/c
Files:
31 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/classes/Constraints/Constraint.h

    r22200 r23588  
    2222                virtual void ActivatePenaltyMethod(void)=0;
    2323                virtual void ConstrainNode(Nodes* nodes,Parameters* parameters)=0;
    24                 virtual bool InAnalysis(int analysis_type)=0;
    2524                virtual void PenaltyDofAndValue(int* dof,IssmDouble* value,Nodes* nodes,Parameters* parameters)=0;
    2625                virtual void InputUpdateFromVectorDakota(IssmDouble* vector,Nodes* nodes,int name,int type) = 0;

  • issm/trunk-jpl/src/c/classes/Constraints/Constraints.cpp

    r23533 r23588  
    2424        for(int i=0;i<this->Size();i++){
    2525                Constraint* constraint=(Constraint*)this->GetObjectByOffset(i);
    26                 if(constraint->InAnalysis(in_analysis)){
    27                         constraint->ActivatePenaltyMethod();
    28                 }
     26                constraint->ActivatePenaltyMethod();
    2927        }
    3028

  • issm/trunk-jpl/src/c/classes/Constraints/SpcDynamic.cpp

    r22361 r23588  
    112112}
    113113/*}}}*/
    114 bool SpcDynamic::InAnalysis(int in_analysis_type){/*{{{*/
    115         if (in_analysis_type==this->analysis_type) return true;
    116         else return false;
    117 }
    118 /*}}}*/
    119114
    120115/*SpcDynamic functions*/

  • issm/trunk-jpl/src/c/classes/Constraints/SpcDynamic.h

    r22361 r23588  
    4040                void ActivatePenaltyMethod(void);
    4141                void ConstrainNode(Nodes* nodes,Parameters* parameters);
    42                 bool InAnalysis(int analysis_type);
    4342                void PenaltyDofAndValue(int* dof,IssmDouble* value,Nodes* nodes,Parameters* parameters){_error_("not implemented yet");};
    4443                void InputUpdateFromVectorDakota(IssmDouble* vector,Nodes* nodes,int name,int type){_error_("not implemented yet");};

  • issm/trunk-jpl/src/c/classes/Constraints/SpcStatic.cpp

    r23066 r23588  
    109109}
    110110/*}}}*/
    111 bool SpcStatic::InAnalysis(int in_analysis_type){/*{{{*/
    112         if (in_analysis_type==this->analysis_type) return true;
    113         else return false;
    114 }
    115 /*}}}*/
    116111void SpcStatic::InputUpdateFromVectorDakota(IssmDouble* vector,Nodes* nodes,int name,int type){/*{{{*/
    117112

  • issm/trunk-jpl/src/c/classes/Constraints/SpcStatic.h

    r22361 r23588  
    3939                void ActivatePenaltyMethod(void);
    4040                void ConstrainNode(Nodes* nodes,Parameters* parameters);
    41                 bool InAnalysis(int analysis_type);
    4241                void PenaltyDofAndValue(int* dof,IssmDouble* value,Nodes* nodes,Parameters* parameters){_error_("not implemented yet");};
    4342                void InputUpdateFromVectorDakota(IssmDouble* vector,Nodes* nodes,int name,int type);

  • issm/trunk-jpl/src/c/classes/Constraints/SpcTransient.cpp

    r22361 r23588  
    159159}
    160160/*}}}*/
    161 bool SpcTransient::InAnalysis(int in_analysis_type){/*{{{*/
    162 
    163         if (in_analysis_type==this->analysis_type) return true;
    164         else return false;
    165 }
    166 /*}}}*/
    167161void SpcTransient::PenaltyDofAndValue(int* pdof,IssmDouble* pvalue,Nodes* nodes,Parameters* parameters){/*{{{*/
    168162

  • issm/trunk-jpl/src/c/classes/Constraints/SpcTransient.h

    r22361 r23588  
    4141                void   ActivatePenaltyMethod(void);
    4242                void   ConstrainNode(Nodes* nodes,Parameters* parameters);
    43                 bool   InAnalysis(int analysis_type);
    4443                void   PenaltyDofAndValue(int* dof,IssmDouble* value,Nodes* nodes,Parameters* parameters);
    4544                void InputUpdateFromVectorDakota(IssmDouble* vector,Nodes* nodes,int name,int type){_error_("not implemented yet");};

  • issm/trunk-jpl/src/c/classes/FemModel.cpp

    r23587 r23588  
    34433443
    34443444                Constraint* constraint=(Constraint*)this->constraints_list[analysis_index]->GetObjectByOffset(i);
    3445                 if(!constraint->InAnalysis(analysis_enum)) _error_("AMR create constraints for "<<EnumToStringx(analysis_enum)<<" not supported yet!\n");
    34463445
    34473446                SpcStatic* spcstatic = xDynamicCast<SpcStatic*>(constraint);
    … …  
    45434542        for (int i=0;i<nodes->Size();i++){
    45444543                Node* node=xDynamicCast<Node*>(nodes->GetObjectByOffset(i));
    4545                 _assert_(node->InAnalysis(HydrologyDCEfficientAnalysisEnum));
    45464544                if(serial_active[node->Sid()]==1.){
    45474545                        node->Activate();
    … …  
    46164614        for (int i=0;i<nodes->Size();i++){
    46174615                Node* node=xDynamicCast<Node*>(nodes->GetObjectByOffset(i));
    4618                 _assert_(node->InAnalysis(HydrologyDCInefficientAnalysisEnum));
    46194616                if(serial_active[node->Sid()]==1.){
    46204617                        node->Activate();
    … …  
    46654662        for (int i=0;i<nodes->Size();i++){
    46664663                Node* node=xDynamicCast<Node*>(nodes->GetObjectByOffset(i));
    4667                 _assert_(node->InAnalysis(L2ProjectionEPLAnalysisEnum));
    46684664                if(serial_active[node->Sid()]==1.){
    46694665                        node->Activate();

  • issm/trunk-jpl/src/c/classes/Loads/Load.h

    r23524 r23588  
    3131                virtual void  GetNodesSidList(int* sidlist)=0;
    3232                virtual int   GetNumberOfNodes(void)=0;
    33                 virtual bool  InAnalysis(int analysis_type)=0;
    3433                virtual bool  IsPenalty(void)=0;
    3534                virtual void  PenaltyCreateJacobianMatrix(Matrix<IssmDouble>* Jff,IssmDouble kmax)=0;

  • issm/trunk-jpl/src/c/classes/Loads/Moulin.cpp

    r23532 r23588  
    219219
    220220        return NUMVERTICES;
    221 }
    222 /*}}}*/
    223 bool  Moulin::InAnalysis(int in_analysis_type){/*{{{*/
    224         if (in_analysis_type==this->analysis_type)return true;
    225         else return false;
    226221}
    227222/*}}}*/

  • issm/trunk-jpl/src/c/classes/Loads/Moulin.h

    r23020 r23588  
    7070                void  GetNodesLidList(int* lidlist);
    7171                int   GetNumberOfNodes(void);
    72                 bool  InAnalysis(int analysis_type);
    7372                bool  IsPenalty(void);
    7473                void  PenaltyCreateJacobianMatrix(Matrix<IssmDouble>* Jff,IssmDouble kmax){_error_("Not implemented yet");};

  • issm/trunk-jpl/src/c/classes/Loads/Neumannflux.cpp

    r23532 r23588  
    246246
    247247        return NUMNODES_BOUNDARY;
    248 }
    249 /*}}}*/
    250 bool  Neumannflux::InAnalysis(int in_analysis_type){/*{{{*/
    251         if (in_analysis_type==this->analysis_type) return true;
    252         else return false;
    253248}
    254249/*}}}*/

  • issm/trunk-jpl/src/c/classes/Loads/Neumannflux.h

    r23020 r23588  
    6363                void GetNodesSidList(int* sidlist);
    6464                int  GetNumberOfNodes(void);
    65                 bool InAnalysis(int analysis_type);
    6665                bool IsPenalty(void);
    6766                void PenaltyCreateJacobianMatrix(Matrix<IssmDouble>* Jff,IssmDouble kmax){_error_("Not implemented yet");};

  • issm/trunk-jpl/src/c/classes/Loads/Numericalflux.cpp

    r23532 r23588  
    345345        }
    346346
    347 }
    348 /*}}}*/
    349 bool  Numericalflux::InAnalysis(int in_analysis_type){/*{{{*/
    350         if (in_analysis_type==this->analysis_type) return true;
    351         else return false;
    352347}
    353348/*}}}*/

  • issm/trunk-jpl/src/c/classes/Loads/Numericalflux.h

    r20827 r23588  
    6464                void GetNodesSidList(int* sidlist);
    6565                int  GetNumberOfNodes(void);
    66                 bool InAnalysis(int analysis_type);
    6766                bool IsPenalty(void);
    6867                void PenaltyCreateJacobianMatrix(Matrix<IssmDouble>* Jff,IssmDouble kmax){_error_("Not implemented yet");};

  • issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp

    r23532 r23588  
    217217
    218218        return NUMVERTICES;
    219 }
    220 /*}}}*/
    221 bool  Pengrid::InAnalysis(int in_analysis_type){/*{{{*/
    222         if (in_analysis_type==this->analysis_type)return true;
    223         else return false;
    224219}
    225220/*}}}*/

  • issm/trunk-jpl/src/c/classes/Loads/Pengrid.h

    r23524 r23588  
    6565                void  GetNodesSidList(int* sidlist);
    6666                int   GetNumberOfNodes(void);
    67                 bool  InAnalysis(int analysis_type);
    6867                bool  IsPenalty(void);
    6968                void  PenaltyCreateJacobianMatrix(Matrix<IssmDouble>* Jff,IssmDouble kmax){_error_("Not implemented yet");};

  • issm/trunk-jpl/src/c/classes/Loads/Penpair.cpp

    r20827 r23588  
    165165
    166166        return NUMVERTICES;
    167 }
    168 /*}}}*/
    169 bool  Penpair::InAnalysis(int in_analysis_type){/*{{{*/
    170         if (in_analysis_type==this->analysis_type)return true;
    171         else return false;
    172167}
    173168/*}}}*/

  • issm/trunk-jpl/src/c/classes/Loads/Penpair.h

    r20827 r23588  
    5555                void  GetNodesSidList(int* sidlist);
    5656                int   GetNumberOfNodes(void);
    57                 bool  InAnalysis(int analysis_type);
    5857                bool  IsPenalty(void);
    5958                void  PenaltyCreateJacobianMatrix(Matrix<IssmDouble>* Jff,IssmDouble kmax);

  • issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp

    r23532 r23588  
    306306
    307307        return NUMVERTICES;
    308 }
    309 /*}}}*/
    310 bool  Riftfront::InAnalysis(int in_analysis_type){/*{{{*/
    311         if (in_analysis_type==this->analysis_type) return true;
    312         else return false;
    313308}
    314309/*}}}*/

  • issm/trunk-jpl/src/c/classes/Loads/Riftfront.h

    r20827 r23588  
    8181                void  GetNodesSidList(int* sidlist);
    8282                int   GetNumberOfNodes(void);
    83                 bool  InAnalysis(int analysis_type);
    8483                bool  IsPenalty(void);
    8584                void  PenaltyCreateJacobianMatrix(Matrix<IssmDouble>* Jff,IssmDouble kmax){_error_("Not implemented yet");};

  • issm/trunk-jpl/src/c/classes/Node.cpp

    r23585 r23588  
    395395}
    396396/*}}}*/
    397 bool Node::InAnalysis(int in_analysis_enum){/*{{{*/
    398         if (in_analysis_enum==this->analysis_enum) return true;
    399         else return false;
    400 }
    401 /*}}}*/
    402397int  Node::Lid(void){/*{{{*/
    403398        return lid;

  • issm/trunk-jpl/src/c/classes/Node.h

    r23585 r23588  
    6868                int   GetNumberOfDofs(int approximation_enum,int setenum);
    6969                void  HardDeactivate(void);
    70                 bool  InAnalysis(int analysis_enum);
    7170                bool  IsActive(void);
    7271                int   IsClone();

  • issm/trunk-jpl/src/c/modules/ConfigureObjectsx/ConfigureObjectsx.cpp

    r23490 r23588  
    3030        for(i=0;i<loads->Size();i++){
    3131                load=(Load*)loads->GetObjectByOffset(i);
    32                 _assert_(load->InAnalysis(configuration_type));
    3332                load->Configure(elements,loads,nodes,vertices,materials,parameters);
    3433        }

  • issm/trunk-jpl/src/c/modules/ConstraintsStatex/ConstraintsStatex.cpp

    r23490 r23588  
    3333        for(int i=0;i<femmodel->loads->Size();i++){
    3434                Load* load=(Load*)femmodel->loads->GetObjectByOffset(i);
    35                 _assert_(load->InAnalysis(analysis_type));
    3635                if(load->ObjectEnum()==PengridEnum){
    3736                        Pengrid* pengrid=(Pengrid*)load;

  • issm/trunk-jpl/src/c/modules/ConstraintsStatex/RiftConstraintsState.cpp

    r23490 r23588  
    2020        for (i=0;i<loads->Size();i++){
    2121                Load* load=(Load*)loads->GetObjectByOffset(i);
    22                 _assert_(load->InAnalysis(configuration_type));
    2322                if(RiftfrontEnum==loads->GetEnum(i)){
    2423                        found=1;
    … …  
    7675
    7776                        load=(Load*)loads->GetObjectByOffset(i);
    78                         _assert_(load->InAnalysis(configuration_type));
    7977                        riftfront=(Riftfront*)load;
    8078                        riftfront->Constrain(&unstable);
    … …  
    109107
    110108                        load=(Load*)loads->GetObjectByOffset(i);
    111                         _assert_(load->InAnalysis(configuration_type));
    112109                        riftfront=(Riftfront*)load;
    113110                        if (riftfront->IsFrozen()){
    … …  
    140137                if (RiftfrontEnum==loads->GetEnum(i)){
    141138                        load=(Load*)loads->GetObjectByOffset(i);
    142                         _assert_(load->InAnalysis(configuration_type));
    143139                        riftfront=(Riftfront*)load;
    144140                        riftfront->FreezeConstraints();

  • issm/trunk-jpl/src/c/modules/CreateJacobianMatrixx/CreateJacobianMatrixx.cpp

    r23490 r23588  
    3636        for (i=0;i<femmodel->loads->Size();i++){
    3737                load=(Load*)femmodel->loads->GetObjectByOffset(i);
    38                 _assert_(load->InAnalysis(configuration_type));
    3938                load->CreateJacobianMatrix(Jff);
    4039                load->PenaltyCreateJacobianMatrix(Jff,kmax);

  • issm/trunk-jpl/src/c/modules/ResetConstraintsx/ResetConstraintsx.cpp

    r23524 r23588  
    3131        for(int i=0;i<femmodel->loads->Size();i++){
    3232                Load* load=(Load*)femmodel->loads->GetObjectByOffset(i);
    33                 _assert_(load->InAnalysis(analysis_type));
    3433                if(load->ObjectEnum()==PengridEnum){
    3534                        Pengrid* pengrid=(Pengrid*)load;

  • issm/trunk-jpl/src/c/modules/UpdateDynamicConstraintsx/UpdateDynamicConstraintsx.cpp

    r14999 r23588  
    2020
    2121        for(int i=0;i<constraints->Size();i++){
     22                Constraint* constraint=(Constraint*)constraints->GetObjectByOffset(i);
     23                ((SpcDynamic*)constraint)->SetDynamicConstraint(nodes,yg_serial);
    2224
    23                 Constraint* constraint=(Constraint*)constraints->GetObjectByOffset(i);
    24 
    25                 /*Check this constraint belongs to this analysis: */
    26                 if(constraint->InAnalysis(configuration_type) && constraint->ObjectEnum()==SpcDynamicEnum){
    27 
    28                         ((SpcDynamic*)constraint)->SetDynamicConstraint(nodes,yg_serial);
    29 
    30                 }
    3125        }
    3226

  • issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp

    r23066 r23588  
    105105        for(int i=0;i<femmodel->constraints->Size();i++){
    106106                Constraint* constraint=(Constraint*)femmodel->constraints->GetObjectByOffset(i);
    107                 if(constraint->InAnalysis(configuration_type)){
    108                         constraint->PenaltyDofAndValue(&dof,&d,femmodel->nodes,femmodel->parameters);
    109                         if(dof!=-1){
    110                                 MatSetValue(LHS,dof,dof,dmax,INSERT_VALUES);
    111                         }
     107                constraint->PenaltyDofAndValue(&dof,&d,femmodel->nodes,femmodel->parameters);
     108                if(dof!=-1){
     109                        MatSetValue(LHS,dof,dof,dmax,INSERT_VALUES);
    112110                }
    113111        }
    … …  
    151149        for(int i=0;i<femmodel->constraints->Size();i++){
    152150                Constraint* constraint=(Constraint*)femmodel->constraints->GetObjectByOffset(i);
    153                 if(constraint->InAnalysis(configuration_type)){
    154                         constraint->PenaltyDofAndValue(&dof,&d,femmodel->nodes,femmodel->parameters);
    155                         d = d*dmax;
    156                         if(dof!=-1){
    157                                 VecSetValues(RHS,1,&dof,(const double*)&d,INSERT_VALUES);
    158                         }
     151                constraint->PenaltyDofAndValue(&dof,&d,femmodel->nodes,femmodel->parameters);
     152                d = d*dmax;
     153                if(dof!=-1){
     154                        VecSetValues(RHS,1,&dof,(const double*)&d,INSERT_VALUES);
    159155                }
    160156        }
Note: See TracChangeset for help on using the changeset viewer.