Changeset 22739 for issm/trunk-jpl/src/c/classes/Inputs/ControlInput.cpp

Timestamp:

05/03/18 13:31:36 (7 years ago)

Author:

erobo

Message:

CHG: allow for transient ControlInputs

File:

• issm/trunk-jpl/src/c/classes/Inputs/ControlInput.cpp (modified) (11 diffs)

issm/trunk-jpl/src/c/classes/Inputs/ControlInput.cpp

@@ -26 +26 @@
         control_id=id;
         enum_type=in_enum_type;
+        layout_enum = input_layout_enum;
 
         _assert_(interp==P1Enum);
@@ -48 +49 @@
 }
 /*}}}*/
+ControlInput::ControlInput(int in_enum_type,int input_layout_enum,Input* input_pvalues,Input* input_pmin,Input* input_pmax,Input* input_pgrad,int interp,int id){/*{{{*/
+
+        this->control_id=id;
+        this->enum_type=in_enum_type;
+        this->layout_enum = input_layout_enum;
+
+        _assert_(interp==P1Enum);
+        if(input_layout_enum!=TransientInputEnum) _error_("Wrong type of layout_enum, needs to be a TransientInputEnum");
+
+                        this->values = input_pvalues;
+                        this->savedvalues= NULL;
+                        this->minvalues  = input_pmin;
+                        this->maxvalues  = input_pmax;
+                        this->gradient =  input_pgrad;
+        }
+/*}}}*/
 ControlInput::~ControlInput(){/*{{{*/
-        delete values;
-        delete savedvalues;
-        delete minvalues;
-        delete maxvalues;
-        delete gradient;
+        if(values)      delete values;
+        if(savedvalues) delete savedvalues;
+        if(minvalues)   delete minvalues;
+        if(maxvalues)   delete maxvalues;
+        if(gradient)    delete gradient;
 }
 /*}}}*/
@@ -65 +82 @@
         output->enum_type=this->enum_type;
         output->control_id=this->control_id;
+        output->layout_enum = this->control_id;
 
         if(values)      output->values      = xDynamicCast<Input*>(this->values->copy());
@@ -79 +97 @@
         _printf_("ControlInput:\n");
         _printf_(setw(15)<<"   ControlInput "<<setw(25)<<left<<EnumToStringx(this->enum_type)<<"\n");
+        _printf_(setw(15)<<"   ControlInput "<<setw(25)<<left<<EnumToStringx(this->layout_enum)<<"\n");
         _printf_("---values: \n");     if (values)      values->Echo();
         _printf_("---savedvalues: \n");if (savedvalues) savedvalues->Echo();
@@ -98 +117 @@
         MARSHALLING(enum_type);
         MARSHALLING(control_id);
+        MARSHALLING(layout_enum);
 
         if (marshall_direction == MARSHALLING_BACKWARD){
@@ -115 +135 @@
                                 gradient   =new PentaInput();
                                 break;
+                        case TransientInputEnum:
+                                values    =new TransientInput();
+                                savedvalues=new TransientInput();
+                                minvalues  =new TransientInput();
+                                maxvalues  =new TransientInput();
+                                gradient   =new TransientInput();
+                                break;
                         default:
                                 _error_("Input of Enum " << EnumToStringx(enum_type) << " not supported yet by ControlInput");
@@ -146 +173 @@
 }/*}}}*/
 void ControlInput::Configure(Parameters* parameters){/*{{{*/
+        if(this->values->ObjectEnum()==TransientInputEnum){
+                this->values->Configure(parameters);
+                this->minvalues->Configure(parameters);
+                this->maxvalues->Configure(parameters);
+                this->gradient->Configure(parameters);
+        }
         /*do nothing: */
 }
@@ -194 +227 @@
 }/*}}}*/
 void ControlInput::GetVectorFromInputs(Vector<IssmDouble>* vector,int* doflist,const char* data){/*{{{*/
-         if(strcmp(data,"value")==0){
+        if(strcmp(data,"value")==0){
                  _assert_(values);
                  values->GetVectorFromInputs(vector,doflist);
@@ -226 +259 @@
         this->savedvalues=xDynamicCast<Input*>(this->values->copy());
 }/*}}}*/
-void ControlInput::SetGradient(Input* gradient_in){/*{{{*/
+void ControlInput::SetGradient(Input* gradient_in,int timestep){/*{{{*/
+if(this->values->ObjectEnum()!=TransientInputEnum)_error_("you are in the wrong place, go home");
 
         /*Get enum for current gradient*/
@@ -243 +277 @@
         }
 
+        TransientInput* transient_input = xDynamicCast<TransientInput*>(this->gradient);
+        TransientInput* values_input = xDynamicCast<TransientInput*>(this->values);
+        if(values_input->numtimesteps==transient_input->numtimesteps){
+                TransientInput* new_trans_input = new TransientInput(ControlInputGradEnum);
+                IssmDouble time = transient_input->GetTimeByOffset(timestep);
+                for(int i=0;i<transient_input->numtimesteps;i++){
+                        if(transient_input->timesteps[i]==time) new_trans_input->AddTimeInput(xDynamicCast<TriaInput*>(gradient_in),time);
+                        else {
+                                Input* input = transient_input->GetTimeInput(transient_input->timesteps[i]);
+                                new_trans_input->AddTimeInput(xDynamicCast<TriaInput*>(input),transient_input->timesteps[i]);
+                        }
+                }
+                this->gradient=new_trans_input;
+        }
+        else{
+                IssmDouble time = values_input->GetTimeByOffset(timestep);
+                transient_input->AddTimeInput(gradient_in,time);
+        }
+}/*}}}*/
+void ControlInput::SetGradient(Input* gradient_in){/*{{{*/
+
+        /*Get enum for current gradient*/
+        switch(this->control_id){
+                case 1:
+                        gradient_in->ChangeEnum(Gradient1Enum);
+                        break;
+                case 2:
+                        gradient_in->ChangeEnum(Gradient2Enum);
+                        break;
+                case 3:
+                        gradient_in->ChangeEnum(Gradient3Enum);
+                        break;
+                default:
+                        _error_("more than 3 controls not implemented yet (Gradient " << this->control_id << " was requested). EnumDefinitions.h needs to be updated.");
+        }
+
         /*Delete old gradient and assign new gradient*/
         if(gradient) delete gradient;
@@ -250 +320 @@
 void ControlInput::SetInput(Input* in_input){/*{{{*/
 
+        if(layout_enum==TransientInputEnum)_error_("need two arguments in SetInput for TransientInput Controls");
         delete values; this->values=in_input;
         this->SaveValue(); //because this is what SpawnResult saves FIXME
+
+}/*}}}*/
+void ControlInput::SetInput(Input* in_input,int timeoffset){/*{{{*/
+        Input* input = this->values;
+        if(input->ObjectEnum()!=TransientInputEnum)_error_("cannot have timeoffset argument if not TransientInput Control");
+        TransientInput* transient_input = xDynamicCast<TransientInput*>(input);
+        IssmDouble time = transient_input->GetTimeByOffset(timeoffset);
+        TransientInput* new_trans_input = new TransientInput(ControlInputValuesEnum);
+        for(int i=0;i<transient_input->numtimesteps;i++){
+                if(transient_input->timesteps[i]==time) new_trans_input->AddTimeInput(xDynamicCast<TriaInput*>(in_input),time);
+                else {
+                        input = transient_input->GetTimeInput(transient_input->timesteps[i]);
+                        new_trans_input->AddTimeInput(xDynamicCast<TriaInput*>(input),transient_input->timesteps[i]);
+                }
+        }
+        this->values=new_trans_input;
+
+        //      this->values->Echo();
+        //this->values->Echo();
+        //new_trans_input->Echo();
 
 }/*}}}*/