TransientInput2 Class Reference

#include <TransientInput2.h>

Inheritance diagram for TransientInput2:

Public Member Functions
	TransientInput2 ()
	TransientInput2 (int in_enum_type, int nbe, int nbv, IssmDouble *times, int N)
	~TransientInput2 ()
void	AddTimeInput (Input2 *input, IssmDouble time)
void	AddTriaTimeInput (IssmDouble time, int numindices, int *indices, IssmDouble *values_in, int interp_in)
void	AddPentaTimeInput (IssmDouble time, int numindices, int *indices, IssmDouble *values_in, int interp_in)
void	AddTriaTimeInput (int step, int numindices, int *indices, IssmDouble *values_in, int interp_in)
void	AddPentaTimeInput (int step, int numindices, int *indices, IssmDouble *values_in, int interp_in)
Input2 *	copy ()
void	Configure (Parameters *params)
void	DeepEcho ()
void	Echo ()
int	Id ()
void	Marshall (char **pmarshalled_data, int *pmarshalled_data_size, int marshall_direction)
int	ObjectEnum ()
void	GetAllTimes (IssmDouble **ptimesteps, int *pnumtimesteps)
TriaInput2 *	GetTriaInput ()
TriaInput2 *	GetTriaInput (IssmDouble time)
TriaInput2 *	GetTriaInput (IssmDouble start_time, IssmDouble end_time, int averaging_method)
TriaInput2 *	GetTriaInput (int offset)
PentaInput2 *	GetPentaInput ()
PentaInput2 *	GetPentaInput (IssmDouble time)
PentaInput2 *	GetPentaInput (int offset)
PentaInput2 *	GetPentaInput (IssmDouble start_time, IssmDouble end_time, int averaging_method)
Input2 *	GetTimeInput (IssmDouble time)
IssmDouble	GetTimeByOffset (int offset)
int	GetTimeInputOffset (IssmDouble time)
void	SetCurrentTimeInput (IssmDouble time)
void	SetAverageAsCurrentTimeInput (IssmDouble start_time, IssmDouble end_time, int averaging_method)
Public Member Functions inherited from Input2
int	InstanceEnum ()
void	ChangeEnum (int newenumtype)
virtual	~Input2 ()
virtual void	GetInputAverage (IssmDouble *pvalue)
virtual IssmDouble	GetInputMax (void)
virtual IssmDouble	GetInputMaxAbs (void)
virtual IssmDouble	GetInputMin (void)
virtual void	GetInputDerivativeValue (IssmDouble *derivativevalues, IssmDouble *xyz_list, Gauss *gauss)
virtual void	GetInputValue (IssmDouble *pvalue, Gauss *gauss)
virtual int	GetInputInterpolationType ()
virtual SegInput2 *	GetSegInput ()
virtual void	AXPY (Input2 *xinput, IssmDouble scalar)
virtual void	PointWiseMult (Input2 *xinput)
virtual void	Pow (IssmDouble scale_factor)
virtual void	Scale (IssmDouble scale_factor)
virtual int	GetResultArraySize (void)
virtual int	GetResultInterpolation (void)
virtual int	GetResultNumberOfNodes (void)
Public Member Functions inherited from Object
virtual	~Object ()

Data Fields
int	enum_type
int	numtimesteps
Input2 **	inputs
IssmDouble *	timesteps
Parameters *	parameters
IssmDouble	current_step
Input2 *	current_input

Private Attributes
int	numberofelements_local
int	numberofvertices_local

Detailed Description

Definition at line 13 of file TransientInput2.h.

Constructor & Destructor Documentation

TransientInput2() [1/2]

TransientInput2::TransientInput2

(

)

Definition at line 18 of file TransientInput2.cpp.

                                {/*{{{*/
 
   enum_type=UNDEF;
   inputs=NULL;
   this->numtimesteps=0;
   this->parameters=NULL;
   this->timesteps=NULL;
 
   this->current_input=NULL;
   this->current_step=-1;
 
}

TransientInput2() [2/2]

TransientInput2::TransientInput2	(	int	in_enum_type,
		int	nbe,
		int	nbv,
		IssmDouble *	times,
		int	N
	)

Definition at line 31 of file TransientInput2.cpp.

                                                                                          {/*{{{*/
 
   /*Set Enum*/
   this->enum_type=in_enum_type;
   this->numberofelements_local = nbe;
   this->numberofvertices_local = nbv;
 
   /*Allocate values and timesteps, and copy: */
   _assert_(N>=0 && N<1e6);
   this->numtimesteps=N;
   if(N>0){
      this->timesteps=xNew<IssmDouble>(N);
      xMemCpy(this->timesteps,timesin,N);
 
      this->inputs     = xNew<Input2*>(N);
      for(int i=0;i<N;i++) this->inputs[i] = NULL;
   }
   else{
      this->timesteps=0;
      this->inputs   =0;
   }
   this->parameters = NULL;
   this->current_input=NULL;
   this->current_step=-1;
}

~TransientInput2()

TransientInput2::~TransientInput2

(

)

Definition at line 57 of file TransientInput2.cpp.

                                 {/*{{{*/
 
   for(int i=0;i<this->numtimesteps;i++){
      delete this->inputs[i];
   }
   xDelete<Input2*>(this->inputs);
   xDelete<IssmDouble>(this->timesteps);
 
   if(this->current_input) delete this->current_input;
}

Member Function Documentation

AddTimeInput()

void TransientInput2::AddTimeInput	(	Input2 *	input,
		IssmDouble	time
	)

AddTriaTimeInput() [1/2]

void TransientInput2::AddTriaTimeInput	(	IssmDouble	time,
		int	numindices,
		int *	indices,
		IssmDouble *	values_in,
		int	interp_in
	)

Definition at line 138 of file TransientInput2.cpp.

                                                                                                                     {/*{{{*/
 
   /*Check whether this is the last time step that we have*/
   if(this->numtimesteps){
      if(fabs(this->timesteps[this->numtimesteps-1]-time)<1.0e-5){
         this->AddTriaTimeInput(this->numtimesteps-1,numindices,indices,values_in,interp_in);
         return;
      }
   }
 
   /*This is a new time step! we need to add it to the list*/
   if(this->numtimesteps>0 && time<this->timesteps[this->numtimesteps-1]) _error_("timestep values must increase sequentially");
 
   IssmDouble *old_timesteps = NULL;
   Input2    **old_inputs    = NULL;
   if (this->numtimesteps > 0){
      old_timesteps=xNew<IssmDouble>(this->numtimesteps);
      xMemCpy(old_timesteps,this->timesteps,this->numtimesteps);
      xDelete<IssmDouble>(this->timesteps);
      old_inputs=xNew<Input2*>(this->numtimesteps);
      xMemCpy(old_inputs,this->inputs,this->numtimesteps);
      xDelete<Input2*>(this->inputs);
   }
 
   this->numtimesteps=this->numtimesteps+1;
   this->timesteps=xNew<IssmDouble>(this->numtimesteps);
   this->inputs   = xNew<Input2*>(this->numtimesteps);
 
   if (this->numtimesteps > 1){
      xMemCpy(this->inputs,old_inputs,this->numtimesteps-1);
      xMemCpy(this->timesteps,old_timesteps,this->numtimesteps-1);
      xDelete(old_timesteps);
      xDelete<Input2*>(old_inputs);
   }
 
   /*go ahead and plug: */
   this->timesteps[this->numtimesteps-1] = time;
   this->inputs[this->numtimesteps-1]    = NULL;
   this->AddTriaTimeInput(this->numtimesteps-1,numindices,indices,values_in,interp_in);
 
}

AddPentaTimeInput() [1/2]

void TransientInput2::AddPentaTimeInput	(	IssmDouble	time,
		int	numindices,
		int *	indices,
		IssmDouble *	values_in,
		int	interp_in
	)

Definition at line 180 of file TransientInput2.cpp.

                                                                                                                      {/*{{{*/
 
   /*Check whether this is the last time step that we have*/
   if(this->numtimesteps){
      if(fabs(this->timesteps[this->numtimesteps-1]-time)<1.0e-5){
         this->AddPentaTimeInput(this->numtimesteps-1,numindices,indices,values_in,interp_in);
         return;
      }
   }
 
   /*This is a new time step! we need to add it to the list*/
   if(this->numtimesteps>0 && time<this->timesteps[this->numtimesteps-1]) _error_("timestep values must increase sequentially");
 
   IssmDouble *old_timesteps = NULL;
   Input2    **old_inputs    = NULL;
   if (this->numtimesteps > 0){
      old_timesteps=xNew<IssmDouble>(this->numtimesteps);
      xMemCpy(old_timesteps,this->timesteps,this->numtimesteps);
      xDelete<IssmDouble>(this->timesteps);
      old_inputs=xNew<Input2*>(this->numtimesteps);
      xMemCpy(old_inputs,this->inputs,this->numtimesteps);
      xDelete<Input2*>(this->inputs);
   }
 
   this->numtimesteps=this->numtimesteps+1;
   this->timesteps=xNew<IssmDouble>(this->numtimesteps);
   this->inputs   = xNew<Input2*>(this->numtimesteps);
 
   if (this->numtimesteps > 1){
      xMemCpy(this->inputs,old_inputs,this->numtimesteps-1);
      xMemCpy(this->timesteps,old_timesteps,this->numtimesteps-1);
      xDelete(old_timesteps);
      xDelete<Input2*>(old_inputs);
   }
 
   /*go ahead and plug: */
   this->timesteps[this->numtimesteps-1] = time;
   this->inputs[this->numtimesteps-1]    = NULL;
   this->AddPentaTimeInput(this->numtimesteps-1,numindices,indices,values_in,interp_in);
 
}

AddTriaTimeInput() [2/2]

void TransientInput2::AddTriaTimeInput	(	int	step,
		int	numindices,
		int *	indices,
		IssmDouble *	values_in,
		int	interp_in
	)

Definition at line 222 of file TransientInput2.cpp.

                                                                                                              {/*{{{*/
 
   _assert_(step>=0 && step<this->numtimesteps);
 
   /*Create it if necessary*/
   if(this->inputs[step]){
      if(this->inputs[step]->ObjectEnum()!=TriaInput2Enum) _error_("cannot add Element values to a "<<EnumToStringx(this->inputs[step]->ObjectEnum()));
   }
   else{
      this->inputs[step] = new TriaInput2(this->numberofelements_local,this->numberofvertices_local,interp_in);
   }
 
   /*Set input*/
   TriaInput2* input = xDynamicCast<TriaInput2*>(this->inputs[step]);
   input->SetInput(interp_in,numindices,indices,values_in);
 
}

AddPentaTimeInput() [2/2]

void TransientInput2::AddPentaTimeInput	(	int	step,
		int	numindices,
		int *	indices,
		IssmDouble *	values_in,
		int	interp_in
	)

Definition at line 240 of file TransientInput2.cpp.

                                                                                                               {/*{{{*/
 
   _assert_(step>=0 && step<this->numtimesteps);
 
   /*Create it if necessary*/
   if(this->inputs[step]){
      if(this->inputs[step]->ObjectEnum()!=PentaInput2Enum) _error_("cannot add Element values to a "<<EnumToStringx(this->inputs[step]->ObjectEnum()));
   }
   else{
      this->inputs[step] = new PentaInput2(this->numberofelements_local,this->numberofvertices_local,interp_in);
   }
 
   /*Set input*/
   PentaInput2* input = xDynamicCast<PentaInput2*>(this->inputs[step]);
   input->SetInput(interp_in,numindices,indices,values_in);
 
}

copy()

Input2 * TransientInput2::copy ( void  )

virtual

Implements Input2.

Definition at line 70 of file TransientInput2.cpp.

                              {/*{{{*/
 
   TransientInput2* output=NULL;
 
   output = new TransientInput2();
   output->enum_type=this->enum_type;
   output->numtimesteps=this->numtimesteps;
   if(this->numtimesteps>0){
      output->timesteps=xNew<IssmDouble>(this->numtimesteps);
      xMemCpy(output->timesteps,this->timesteps,this->numtimesteps);
      output->inputs = xNew<Input2*>(this->numtimesteps);
      for(int i=0;i<this->numtimesteps;i++){
         output->inputs[i] = this->inputs[i]->copy();
      }
   }
   output->parameters=this->parameters;
 
   return output;
}/*}}}*/

Configure()

void TransientInput2::Configure ( Parameters *  params )

virtual

Reimplemented from Input2.

Definition at line 104 of file TransientInput2.cpp.

                                                 {/*{{{*/
   this->parameters=params;
}

DeepEcho()

void TransientInput2::DeepEcho ( void  )

virtual

Implements Object.

Definition at line 89 of file TransientInput2.cpp.

                                  {/*{{{*/
 
   int i;
 
   _printf_("TransientInput2:\n");
   _printf_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")\n");
   _printf_("   numtimesteps: " << this->numtimesteps << "\n");
   _printf_("---inputs: \n");
   for(i=0;i<this->numtimesteps;i++){
      _printf_("   time: " << this->timesteps[i]<<"  ");
      if(this->inputs[i]) this->inputs[i]->Echo();
      else                _printf_(" NOT SET! \n");
   }
}

Echo()

void TransientInput2::Echo ( void  )

virtual

Implements Object.

Definition at line 108 of file TransientInput2.cpp.

                              {/*{{{*/
   this->DeepEcho();
}

Id()

int TransientInput2::Id ( void  )

virtual

Implements Object.

Definition at line 112 of file TransientInput2.cpp.

{ return -1; }/*{{{*/

Marshall()

void TransientInput2::Marshall ( char **  pmarshalled_data, int *  pmarshalled_data_size, int  marshall_direction  )

virtual

Implements Object.

Definition at line 114 of file TransientInput2.cpp.

                                                                                                        { /*{{{*/
 
   if (marshall_direction == MARSHALLING_BACKWARD){
      _error_("not implmented");
      //inputs = new Inputs();
   }
 
   MARSHALLING_ENUM(TransientInput2Enum);
 
   MARSHALLING(enum_type);
   MARSHALLING(numtimesteps);
   MARSHALLING_DYNAMIC(this->timesteps,IssmDouble,numtimesteps);
   //inputs->Marshall(pmarshalled_data,pmarshalled_data_size,marshall_direction);
   _error_("not implemented");
}

ObjectEnum()

int TransientInput2::ObjectEnum ( void  )

virtual

Implements Object.

Definition at line 130 of file TransientInput2.cpp.

                                    {/*{{{*/
 
   return TransientInput2Enum;
 
}

GetAllTimes()

void TransientInput2::GetAllTimes	(	IssmDouble **	ptimesteps,
		int *	pnumtimesteps
	)

Definition at line 258 of file TransientInput2.cpp.

                                                                           {/*{{{*/
 
   if(ptimesteps){
      *ptimesteps=xNew<IssmDouble>(this->numtimesteps);
      xMemCpy(*ptimesteps,this->timesteps,this->numtimesteps);
   }
   if(pnumtimesteps){
      *pnumtimesteps = this->numtimesteps;
   }
 
}

GetTriaInput() [1/4]

TriaInput2 * TransientInput2::GetTriaInput ( )

virtual

Reimplemented from Input2.

Definition at line 270 of file TransientInput2.cpp.

                                         {/*{{{*/
 
   IssmDouble time;
   this->parameters->FindParam(&time,TimeEnum);
   return this->GetTriaInput(time);
 
}

GetTriaInput() [2/4]

TriaInput2 * TransientInput2::GetTriaInput

(

IssmDouble

time

)

Definition at line 278 of file TransientInput2.cpp.

                                                        {/*{{{*/
 
   /*Set current time input*/
   this->SetCurrentTimeInput(time);
   _assert_(this->current_input);
 
   /*Cast and return*/
   if(this->current_input->ObjectEnum()!=TriaInput2Enum){
      _error_("Cannot return a TriaInput2");
   }
   return xDynamicCast<TriaInput2*>(this->current_input);
 
}

GetTriaInput() [3/4]

TriaInput2 * TransientInput2::GetTriaInput	(	IssmDouble	start_time,
		IssmDouble	end_time,
		int	averaging_method
	)

Definition at line 292 of file TransientInput2.cpp.

                                                                                                         {/*{{{*/
 
   /*Set current time input*/
   this->SetAverageAsCurrentTimeInput(start_time,end_time,averaging_method);
   _assert_(this->current_input);
 
   /*Cast and return*/
   if(this->current_input->ObjectEnum()!=TriaInput2Enum){
      _error_("Cannot return a TriaInput2");
   }
   return xDynamicCast<TriaInput2*>(this->current_input);
 
}

GetTriaInput() [4/4]

TriaInput2 * TransientInput2::GetTriaInput

(

int

offset

)

Definition at line 306 of file TransientInput2.cpp.

                                                   {/*{{{*/
 
   /*Check offset*/
   if(offset<0 || offset>this->numtimesteps-1){
      _error_("Cannot return input for offset "<<offset);
   }
   Input2* input = this->inputs[offset];
 
   /*Cast and return*/
   _assert_(input);
   if(input->ObjectEnum()!=TriaInput2Enum) _error_("Cannot return a TriaInput2");
   return xDynamicCast<TriaInput2*>(input);
 
}

GetPentaInput() [1/4]

PentaInput2 * TransientInput2::GetPentaInput ( )

virtual

Reimplemented from Input2.

Definition at line 321 of file TransientInput2.cpp.

                                           {/*{{{*/
 
   IssmDouble time;
   this->parameters->FindParam(&time,TimeEnum);
   return this->GetPentaInput(time);
}

GetPentaInput() [2/4]

PentaInput2 * TransientInput2::GetPentaInput

(

IssmDouble

time

)

Definition at line 328 of file TransientInput2.cpp.

                                                          {/*{{{*/
 
   /*Set current time input*/
   this->SetCurrentTimeInput(time);
   _assert_(this->current_input);
 
   /*Cast and return*/
   if(this->current_input->ObjectEnum()!=PentaInput2Enum){
      _error_("Cannot return a PentaInput2");
   }
   return xDynamicCast<PentaInput2*>(this->current_input);
 
}

GetPentaInput() [3/4]

PentaInput2 * TransientInput2::GetPentaInput

(

int

offset

)

Definition at line 342 of file TransientInput2.cpp.

                                                     {/*{{{*/
 
 
   /*Check offset*/
   if(offset<0 || offset>this->numtimesteps-1){
      _error_("Cannot return input for offset "<<offset);
   }
   Input2* input = this->inputs[offset];
 
   /*Cast and return*/
   if(input->ObjectEnum()!=PentaInput2Enum) _error_("Cannot return a PentaInput2");
   return xDynamicCast<PentaInput2*>(input);
 
}

GetPentaInput() [4/4]

PentaInput2 * TransientInput2::GetPentaInput	(	IssmDouble	start_time,
		IssmDouble	end_time,
		int	averaging_method
	)

Definition at line 357 of file TransientInput2.cpp.

                                                                                                           {/*{{{*/
 
   /*Set current time input*/
   this->SetAverageAsCurrentTimeInput(start_time,end_time,averaging_method);
   _assert_(this->current_input);
 
   /*Cast and return*/
   if(this->current_input->ObjectEnum()!=PentaInput2Enum){
      _error_("Cannot return a PentaInput2");
   }
   return xDynamicCast<PentaInput2*>(this->current_input);
 
}

GetTimeInput()

Input2* TransientInput2::GetTimeInput ( IssmDouble  time )

inline

Definition at line 58 of file TransientInput2.h.

{_error_("This should not happen!");};

GetTimeByOffset()

IssmDouble TransientInput2::GetTimeByOffset

(

int

offset

)

Definition at line 539 of file TransientInput2.cpp.

                                                      {/*{{{*/
   if(offset<0) offset=0;
   _assert_(offset<this->numtimesteps);
   return this->timesteps[offset];
}

GetTimeInputOffset()

int TransientInput2::GetTimeInputOffset

(

IssmDouble

time

)

Definition at line 545 of file TransientInput2.cpp.

                                                       {/*{{{*/
 
   int offset;
 
   /*go through the timesteps, and figure out which interval we
    *     *fall within. Then interpolate the values on this interval: */
   int found=binary_search(&offset,time,this->timesteps,this->numtimesteps);
   if(!found) _error_("Input not found (is TransientInput sorted ?)");
 
   return offset;
}

SetCurrentTimeInput()

void TransientInput2::SetCurrentTimeInput

(

IssmDouble

time

)

Definition at line 372 of file TransientInput2.cpp.

                                                        {/*{{{*/
 
   /*First, recover current time from parameters: */
   bool linear_interp;
   this->parameters->FindParam(&linear_interp,TimesteppingInterpForcingsEnum);
 
   /*Figure step out*/
   int offset;
   if(!binary_search(&offset,time,this->timesteps,this->numtimesteps)){
      _error_("Input not found (is TransientInput sorted ?)");
   }
 
   if (offset==-1){
 
      /*get values for the first time: */
      _assert_(time<this->timesteps[0]);
 
      /*If already processed return*/
      if(this->current_step==0.) return;
 
      /*Prepare input*/
      if(this->current_input) delete this->current_input;
      this->current_step = 0.;
      this->current_input = this->inputs[0]->copy();
 
   }
   else if(offset==(this->numtimesteps-1) || !linear_interp){
 
      /*get values for the last time: */
      _assert_(time>=this->timesteps[offset]);
 
      /*If already processed return*/
      if(this->current_step==reCast<IssmDouble>(offset)) return;
 
      /*Prepare input*/
      if(this->current_input) delete this->current_input;
      this->current_step  = reCast<IssmDouble>(offset);
      this->current_input = this->inputs[offset]->copy();
   }
   else {
 
      /*Interpolate */
      _assert_(time>=this->timesteps[offset] && time<this->timesteps[offset+1]);
 
      /*get values between two times [offset:offset+1[, Interpolate linearly*/
      IssmDouble deltat=this->timesteps[offset+1]-this->timesteps[offset];
      IssmDouble this_step = reCast<IssmDouble>(offset) + (time - this->timesteps[offset])/deltat;
 
      /*If already processed return*/
      if(fabs(this->current_step-this_step)<1.e-5) return;
      //    if(this->current_step>this_step-1.e-5 && this->current_step<this_step+1.e-5) return;
 
      /*Prepare input*/
      if(this->current_input) delete this->current_input;
      this->current_step = this_step;
      IssmDouble alpha2=(time-this->timesteps[offset])/deltat;
      IssmDouble alpha1=(1.0-alpha2);
 
      Input2* input1=this->inputs[offset];
      Input2* input2=this->inputs[offset+1];
 
      this->current_input = input1->copy();
      this->current_input->Scale(alpha1);
      this->current_input->AXPY(input2,alpha2);
   }
 
}/*}}}*/

SetAverageAsCurrentTimeInput()

void TransientInput2::SetAverageAsCurrentTimeInput	(	IssmDouble	start_time,
		IssmDouble	end_time,
		int	averaging_method
	)

Definition at line 439 of file TransientInput2.cpp.

                                                                                                                 {/*{{{*/
 
   IssmDouble  dt,durinv;
   IssmDouble  dtsum=0;
   IssmDouble  timespan,mid_step;
   int         found,start_offset,end_offset,input_offset;
 
   /*go through the timesteps, and grab offset for start and end*/
   found=binary_search(&start_offset,start_time,this->timesteps,this->numtimesteps);
   if(!found) _error_("Input not found (is TransientInput sorted ?)");
   found=binary_search(&end_offset,end_time,this->timesteps,this->numtimesteps);
   if(!found) _error_("Input not found (is TransientInput sorted ?)");
 
   if(start_offset==-1){
      timespan=this->timesteps[end_offset]-start_time;
   }
   else{
      timespan=this->timesteps[end_offset]-this->timesteps[start_offset];
   }
   mid_step=reCast<IssmDouble>(start_offset)+0.5*timespan;
 
   /*If already processed return, we set step in the middle of the interval*/
   if(fabs(this->current_step-mid_step)<1.e-5) return;
   /*If not processed set current_step*/
   if(this->current_input) delete this->current_input;
   this->current_step = mid_step;
 
   int offset=start_offset;
   while(offset < end_offset){
      if(offset==start_offset){
         dt=this->timesteps[offset+1]-start_time;
         _assert_(dt>0.);
         if(offset==end_offset-1){
            dt=end_time-start_time;
            _assert_(dt>0.);
         }
      }
      else if(offset==end_offset-1){
         dt=end_time-this->timesteps[offset];
         _assert_(dt>0.);
      }
      else{
         dt=this->timesteps[offset+1]-this->timesteps[offset];
         _assert_(dt>0.);
      }
      Input2* stepinput=this->inputs[offset+1];
 
      switch(averaging_method){
         case 0: /*Arithmetic mean*/
            if(offset==start_offset){
               this->current_input=stepinput->copy();
               this->current_input->Scale(dt);
            }
            else{
               this->current_input->AXPY(stepinput,dt);
            }
            break;
         case 1: /*Geometric mean*/
            if(offset==start_offset){
               this->current_input = stepinput->copy();
               this->current_input->Scale(dt);
            }
            else{
               stepinput->Scale(dt);
               this->current_input->PointWiseMult(stepinput);
            }
            break;
         case 2: /*Harmonic mean*/
            if(offset==start_offset){
               this->current_input = stepinput->copy();
               this->current_input->Pow(-1);
               this->current_input->Scale(dt);
            }
            else{
               stepinput->Pow(-1);
               this->current_input->AXPY(stepinput,dt);
            }
         default:
            _error_("averaging method is not recognised");
      }
      dtsum+=dt;
      offset+=1;
   }
   _assert_(dtsum>0);
   durinv=1./dtsum;
   /*Integration done, now normalize*/
   switch(averaging_method){
      case 0: //Arithmetic mean
         this->current_input->Scale(durinv);
         break;
      case 1: /*Geometric mean*/
         this->current_input->Pow(durinv);
         break;
      case 2: /*Harmonic mean*/
         this->current_input->Scale(durinv);
         this->current_input->Pow(-1);
      default:
         _error_("averaging method is not recognised");
   }
}/*}}}*/

Field Documentation

numberofelements_local

int TransientInput2::numberofelements_local

private

Definition at line 16 of file TransientInput2.h.

numberofvertices_local

int TransientInput2::numberofvertices_local

private

Definition at line 17 of file TransientInput2.h.

enum_type

int TransientInput2::enum_type

Definition at line 20 of file TransientInput2.h.

numtimesteps

int TransientInput2::numtimesteps

Definition at line 21 of file TransientInput2.h.

inputs

Input2** TransientInput2::inputs

Definition at line 22 of file TransientInput2.h.

timesteps

IssmDouble* TransientInput2::timesteps

Definition at line 23 of file TransientInput2.h.

parameters

Parameters* TransientInput2::parameters

Definition at line 24 of file TransientInput2.h.

current_step

IssmDouble TransientInput2::current_step

Definition at line 26 of file TransientInput2.h.

current_input

Input2* TransientInput2::current_input

Definition at line 27 of file TransientInput2.h.

---

The documentation for this class was generated from the following files:

• src/c/classes/Inputs2/TransientInput2.h
• src/c/classes/Inputs2/TransientInput2.cpp