SegInput2 Class Reference

#include <SegInput2.h>

Inheritance diagram for SegInput2:

Public Member Functions
	SegInput2 ()
	SegInput2 (int nbe_in, int nbv_in, int interp_in)
	~SegInput2 ()
Input2 *	copy ()
void	DeepEcho ()
void	Echo ()
int	Id ()
void	Marshall (char **pmarshalled_data, int *pmarshalled_data_size, int marshall_direction)
int	ObjectEnum ()
void	SetInput (int interp_in, int row, IssmDouble value_in)
void	SetInput (int interp_in, int numinds, int *rows, IssmDouble *values_in)
void	SetInput (int interp_in, int row, int numinds, IssmDouble *values_in)
int	GetInterpolation ()
void	GetInputDerivativeValue (IssmDouble *derivativevalues, IssmDouble *xyz_list, Gauss *gauss)
void	GetInputAverage (IssmDouble *pvalue)
IssmDouble	GetInputMin ()
IssmDouble	GetInputMax ()
IssmDouble	GetInputMaxAbs ()
SegInput2 *	GetSegInput ()
void	GetInputValue (IssmDouble *pvalue, Gauss *gauss)
void	Scale (IssmDouble scalar)
void	Pow (IssmDouble scalar)
void	AXPY (Input2 *xinput, IssmDouble scalar)
void	PointWiseMult (Input2 *xinput)
void	Serve (int numindices, int *indices)
void	Serve (int row, int numindices)
int	GetResultArraySize (void)
int	GetResultInterpolation (void)
int	GetResultNumberOfNodes (void)
void	Reset (int interp_in)
Public Member Functions inherited from ElementInput2
	ElementInput2 ()
	~ElementInput2 ()
int	GetInputInterpolationType ()
Public Member Functions inherited from Input2
int	InstanceEnum ()
void	ChangeEnum (int newenumtype)
virtual	~Input2 ()
virtual void	Configure (Parameters *parameters)
virtual TriaInput2 *	GetTriaInput ()
virtual PentaInput2 *	GetPentaInput ()
Public Member Functions inherited from Object
virtual	~Object ()
Public Member Functions inherited from SegRef
	SegRef ()
	~SegRef ()
void	GetInputDerivativeValue (IssmDouble *p, IssmDouble *plist, IssmDouble *xyz_list, GaussSeg *gauss, int finiteelement)
void	GetInputValue (IssmDouble *p, IssmDouble *plist, GaussSeg *gauss, int finiteelement)
void	GetJacobian (IssmDouble *J, IssmDouble *xyz_list, GaussSeg *gauss)
void	GetJacobianDeterminant (IssmDouble *Jdet, IssmDouble *xyz_list, GaussSeg *gauss)
void	GetJacobianInvert (IssmDouble *Jinv, IssmDouble *xyz_list, GaussSeg *gauss)
void	GetNodalFunctions (IssmDouble *basis, GaussSeg *gauss, int finiteelement)
void	GetNodalFunctionsDerivatives (IssmDouble *dbasis, IssmDouble *xyz_list, GaussSeg *gauss, int finiteelement)
void	GetNodalFunctionsDerivativesReference (IssmDouble *dbasis, GaussSeg *gauss, int finiteelement)
void	Marshall (char **pmarshalled_data, int *pmarshalled_data_size, int marshall_direction)
int	NumberofNodes (int finiteelement)

Additional Inherited Members
Data Fields inherited from ElementInput2
IssmDouble *	element_values
Protected Attributes inherited from ElementInput2
int	numberofelements_local
int	numberofvertices_local
int	interpolation
int	M
int	N
bool	isserved
IssmDouble *	values

Detailed Description

Definition at line 8 of file SegInput2.h.

Constructor & Destructor Documentation

SegInput2() [1/2]

SegInput2::SegInput2

(

void

)

Definition at line 16 of file SegInput2.cpp.

                        {/*{{{*/
 
   this->numberofelements_local = -1;
   this->numberofvertices_local = -1;
   this->isserved       = false;
   this->M = -1;
   this->N = -1;
   this->values         = NULL;
   this->element_values = NULL;
 
}/*}}}*/

SegInput2() [2/2]

SegInput2::SegInput2	(	int	nbe_in,
		int	nbv_in,
		int	interp_in
	)

Definition at line 27 of file SegInput2.cpp.

                                                       {/*{{{*/
 
   _assert_(nbe_in>0);
   _assert_(nbe_in<1e11);
   _assert_(nbv_in>0);
   _assert_(nbv_in<1e11);
   this->numberofelements_local = nbe_in;
   this->numberofvertices_local = nbv_in;
   this->isserved       = false;
 
   /*Reset takes care of the rest*/
   this->Reset(interp_in);
}/*}}}*/

~SegInput2()

SegInput2::~SegInput2

(

)

Definition at line 40 of file SegInput2.cpp.

                     {/*{{{*/
   if(this->element_values) xDelete<IssmDouble>(this->element_values);
   if(this->values)         xDelete<IssmDouble>(this->values);
}

Member Function Documentation

copy()

Input2 * SegInput2::copy ( void  )

virtual

Implements ElementInput2.

Definition at line 70 of file SegInput2.cpp.

                        {/*{{{*/
 
   SegInput2* output = new SegInput2(this->numberofelements_local,this->numberofvertices_local,this->interpolation);
 
   xMemCpy<IssmDouble>(output->values,this->values,this->M*this->N);
   xMemCpy<IssmDouble>(output->element_values,this->element_values,SegRef::NumberofNodes(this->interpolation));
 
   return output;
}

DeepEcho()

void SegInput2::DeepEcho ( void  )

virtual

Implements ElementInput2.

Definition at line 80 of file SegInput2.cpp.

                            {/*{{{*/
   _printf_("SegInput2 Echo:\n");
   _printf_("   interpolation: "<<EnumToStringx(this->interpolation)<<"\n");
   _printf_("   Size:          "<<M<<"x"<<N<<"\n");
   _printf_("   isserved:      "<<(isserved?"true":"false") << "\n");
   if(isserved){
      _printf_("   current values:      ");
      for(int i=0;i<3;i++) _printf_(" "<<this->element_values[i]);
      _printf_("] ("<<EnumToStringx(this->interpolation)<<")\n");
   }
   printarray(this->values,this->M,this->N);
   //_printf_(setw(15)<<"   SegInput2 "<<setw(25)<<left<<EnumToStringx(this->enum_type)<<" "<<(value?"true":"false") << "\n");
}

Echo()

void SegInput2::Echo ( void  )

virtual

Implements ElementInput2.

Definition at line 94 of file SegInput2.cpp.

                        {/*{{{*/
   _printf_("SegInput2 Echo:\n");
   _printf_("   interpolation: "<<EnumToStringx(this->interpolation)<<"\n");
   _printf_("   Size:          "<<M<<"x"<<N<<"\n");
   _printf_("   isserved:      "<<(isserved?"true":"false") << "\n");
   if(isserved){
      _printf_("   current values:      ");
      _printf_("[ ");
      for(int i=0;i<3;i++) _printf_(" "<<this->element_values[i]);
      _printf_("] ("<<EnumToStringx(this->interpolation)<<")\n");
   }
}

Id()

int SegInput2::Id ( void  )

virtual

Implements ElementInput2.

Definition at line 107 of file SegInput2.cpp.

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

Marshall()

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

virtual

Implements ElementInput2.

Definition at line 110 of file SegInput2.cpp.

                                                                                                  { /*{{{*/
 
   MARSHALLING_ENUM(SegInput2Enum);
   MARSHALLING(this->numberofelements_local);
   MARSHALLING(this->numberofvertices_local);
   MARSHALLING(this->interpolation);
   MARSHALLING(this->M);
   MARSHALLING(this->N);
   this->isserved = false;
   if(this->M*this->N){
      MARSHALLING_DYNAMIC(this->values,IssmDouble,this->M*this->N);
   }
   else this->values = NULL;
 
   if(marshall_direction == MARSHALLING_BACKWARD){
      this->element_values = xNewZeroInit<IssmDouble>(SegRef::NumberofNodes(this->interpolation));
   }
 
}

ObjectEnum()

int SegInput2::ObjectEnum ( void  )

virtual

Implements ElementInput2.

Definition at line 130 of file SegInput2.cpp.

                              {/*{{{*/
   return SegInput2Enum;
}

SetInput() [1/3]

void SegInput2::SetInput ( int  interp_in, int  row, IssmDouble  value_in  )

virtual

Implements ElementInput2.

Definition at line 136 of file SegInput2.cpp.

                                                                 {/*{{{*/
 
   _assert_(this);
   _assert_(row>=0);
   _assert_(row<this->M);
   _assert_(this->N==1);
 
   this->values[row] = value_in;
   this->isserved = false;
}

SetInput() [2/3]

void SegInput2::SetInput ( int  interp_in, int  numinds, int *  rows, IssmDouble *  values_in  )

virtual

Implements ElementInput2.

Definition at line 147 of file SegInput2.cpp.

                                                                                       {/*{{{*/
 
   _assert_(this);
   if(interp_in==P1Enum && this->interpolation==P1Enum){
      _assert_(this->N==1);
      for(int i=0;i<numindices;i++){
         int row = indices[i];
         _assert_(row>=0);
         _assert_(row<this->M);
         this->values[row] = values_in[i];
      }
   }
   else if(interp_in==P0Enum && this->interpolation==P0Enum){
      _assert_(this->N==1);
      for(int i=0;i<numindices;i++){
         int row = indices[i];
         _assert_(row>=0);
         _assert_(row<this->M);
         this->values[row] = values_in[i];
      }
   }
   else if(this->interpolation!=P1Enum && interp_in==P1Enum){
      this->Reset(interp_in);
      for(int i=0;i<numindices;i++){
         int row = indices[i];
         _assert_(row>=0);
         _assert_(row<this->M);
         this->values[row] = values_in[i];
      }
   }
   else{
      _error_("Cannot convert "<<EnumToStringx(this->interpolation)<<" to "<<EnumToStringx(interp_in));
   }
   this->isserved = false;
}

SetInput() [3/3]

void SegInput2::SetInput ( int  interp_in, int  row, int  numinds, IssmDouble *  values_in  )

virtual

Implements ElementInput2.

Definition at line 183 of file SegInput2.cpp.

                                                                                  {/*{{{*/
 
   _assert_(this);
   if(interp_in==this->interpolation){
      _assert_(this->N==numindices);
   }
   else{
      this->Reset(interp_in);
      _assert_(this->N==numindices);
   }
   for(int i=0;i<numindices;i++) this->values[row*this->N+i] = values_in[i];
   this->isserved = false;
}

GetInterpolation()

int SegInput2::GetInterpolation ( )

virtual

Implements ElementInput2.

Definition at line 227 of file SegInput2.cpp.

                                {/*{{{*/
   return this->interpolation;
}/*}}}*/

GetInputDerivativeValue()

void SegInput2::GetInputDerivativeValue ( IssmDouble *  derivativevalues, IssmDouble *  xyz_list, Gauss *  gauss  )

virtual

Implements ElementInput2.

Definition at line 279 of file SegInput2.cpp.

                                                                                                       {/*{{{*/
   _assert_(this);
   _assert_(this->isserved);
   _assert_(gauss->Enum()==GaussSegEnum);
   SegRef::GetInputDerivativeValue(derivativevalues,this->element_values,xyz_list,(GaussSeg*)gauss,this->interpolation);
}/*}}}*/

GetInputAverage()

void SegInput2::GetInputAverage ( IssmDouble *  pvalue )

virtual

Reimplemented from Input2.

Definition at line 230 of file SegInput2.cpp.

                                                 {/*{{{*/
   _assert_(this);
   _assert_(this->isserved);
 
   int        numnodes  = this->NumberofNodes(this->interpolation);
   IssmDouble numnodesd = reCast<int,IssmDouble>(numnodes);
   IssmDouble value     = 0.;
 
   for(int i=0;i<numnodes;i++) value+=this->element_values[i];
   value = value/numnodesd;
 
   *pvalue=value;
}/*}}}*/

GetInputMin()

IssmDouble SegInput2::GetInputMin ( void  )

virtual

Reimplemented from Input2.

Definition at line 243 of file SegInput2.cpp.

                                     {/*{{{*/
   _assert_(this);
   _assert_(this->isserved);
 
   int        numnodes  = this->NumberofNodes(this->interpolation);
   IssmDouble min=this->element_values[0];
 
   for(int i=1;i<numnodes;i++){
      if(this->element_values[i]<min) min=this->element_values[i];
   }
   return min;
}/*}}}*/

GetInputMax()

IssmDouble SegInput2::GetInputMax ( void  )

virtual

Reimplemented from Input2.

Definition at line 255 of file SegInput2.cpp.

                                     {/*{{{*/
   _assert_(this);
   _assert_(this->isserved);
 
   int        numnodes  = this->NumberofNodes(this->interpolation);
   IssmDouble max=this->element_values[0];
 
   for(int i=1;i<numnodes;i++){
      if(this->element_values[i]>max) max=this->element_values[i];
   }
   return max;
}/*}}}*/

GetInputMaxAbs()

IssmDouble SegInput2::GetInputMaxAbs ( void  )

virtual

Reimplemented from Input2.

Definition at line 267 of file SegInput2.cpp.

                                        {/*{{{*/
   _assert_(this);
   _assert_(this->isserved);
 
   int        numnodes  = this->NumberofNodes(this->interpolation);
   IssmDouble maxabs=fabs(this->element_values[0]);
 
   for(int i=1;i<numnodes;i++){
      if(fabs(this->element_values[i])>maxabs) maxabs=fabs(this->element_values[i]);
   }
   return maxabs;
}/*}}}*/

GetSegInput()

SegInput2* SegInput2::GetSegInput ( )

inlinevirtual

Reimplemented from Input2.

Definition at line 34 of file SegInput2.h.

{return this;};

GetInputValue()

void SegInput2::GetInputValue ( IssmDouble *  pvalue, Gauss *  gauss  )

virtual

Implements ElementInput2.

Definition at line 285 of file SegInput2.cpp.

                                                            {/*{{{*/
   _assert_(this);
   _assert_(this->isserved);
   _assert_(gauss->Enum()==GaussSegEnum);
   SegRef::GetInputValue(pvalue,this->element_values,(GaussSeg*)gauss,this->interpolation);
}/*}}}*/

Scale()

void SegInput2::Scale ( IssmDouble  scalar )

virtual

Reimplemented from Input2.

Definition at line 305 of file SegInput2.cpp.

                                     {/*{{{*/
 
   for(int i=0;i<this->M*this->N;i++) this->values[i] = alpha*this->values[i];
   for(int i=0;i<SegRef::NumberofNodes(this->interpolation);i++) this->element_values[i] = alpha*this->element_values[i];
}

Pow()

void SegInput2::Pow ( IssmDouble  scalar )

virtual

Reimplemented from Input2.

Definition at line 311 of file SegInput2.cpp.

                                   {/*{{{*/
 
   for(int i=0;i<this->M*this->N;i++) this->values[i] = pow(this->values[i],alpha);
   for(int i=0;i<SegRef::NumberofNodes(this->interpolation);i++) this->element_values[i] = pow(this->element_values[i],alpha);
}

AXPY()

void SegInput2::AXPY ( Input2 *  xinput, IssmDouble  scalar  )

virtual

Reimplemented from Input2.

Definition at line 317 of file SegInput2.cpp.

                                                   {/*{{{*/
 
   /*xinput is of the same type, so cast it: */
   if(xinput->ObjectEnum()!=SegInput2Enum) _error_("Operation not permitted because xinput is of type " << EnumToStringx(xinput->ObjectEnum()));
   SegInput2* xseginput=xDynamicCast<SegInput2*>(xinput);
   if(xseginput->GetInterpolation()!=this->interpolation) _error_("Operation not permitted because xinput is of type " << EnumToStringx(xinput->ObjectEnum()));
 
   /*Carry out the AXPY operation depending on type:*/
   for(int i=0;i<this->M*this->N;i++) this->values[i] = alpha*xseginput->values[i] + this->values[i];
   for(int i=0;i<SegRef::NumberofNodes(this->interpolation);i++) this->element_values[i] = alpha*xseginput->element_values[i] + this->element_values[i];
}

PointWiseMult()

void SegInput2::PointWiseMult ( Input2 *  xinput )

virtual

Reimplemented from Input2.

Definition at line 329 of file SegInput2.cpp.

                                           {/*{{{*/
 
   /*xinput is of the same type, so cast it: */
   if(xinput->ObjectEnum()!=SegInput2Enum) _error_("Operation not permitted because xinput is of type " << EnumToStringx(xinput->ObjectEnum()));
   SegInput2* xseginput=xDynamicCast<SegInput2*>(xinput);
   if(xseginput->GetInterpolation()!=this->interpolation) _error_("Operation not permitted because xinput is of type " << EnumToStringx(xinput->ObjectEnum()));
 
   /* we need to check that the vector sizes are identical*/
   if(xseginput->M!=this->M||xseginput->N!=this->N) _error_("Operation not permitted because the inputs have different sizes");
 
   /*Carry out the AXPY operation depending on type:*/
   for(int i=0;i<this->M*this->N;i++) this->values[i] = xseginput->values[i] * this->values[i];
   for(int i=0;i<SegRef::NumberofNodes(this->interpolation);i++) this->element_values[i] = xseginput->element_values[i] * this->element_values[i];
}

Serve() [1/2]

void SegInput2::Serve ( int  numindices, int *  indices  )

virtual

Implements ElementInput2.

Definition at line 197 of file SegInput2.cpp.

                                                {/*{{{*/
 
   _assert_(this);
   _assert_(this->N==1);
 
   for(int i=0;i<numindices;i++){
      int row = indices[i];
      _assert_(row>=0);
      _assert_(row<this->M);
      this->element_values[i] = this->values[row];
   }
 
   /*Set input as served*/
   this->isserved = true;
}

Serve() [2/2]

void SegInput2::Serve ( int  row, int  numindices  )

virtual

Implements ElementInput2.

Definition at line 213 of file SegInput2.cpp.

                                           {/*{{{*/
 
   _assert_(this);
   _assert_(this->N==numindices);
   _assert_(row<this->M);
   _assert_(row>=0);
 
   for(int i=0;i<numindices;i++){
      this->element_values[i] = this->values[row*this->N+i];
   }
 
   /*Set input as served*/
   this->isserved = true;
} /*}}}*/

GetResultArraySize()

int SegInput2::GetResultArraySize ( void  )

virtual

Implements ElementInput2.

Definition at line 291 of file SegInput2.cpp.

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

GetResultInterpolation()

int SegInput2::GetResultInterpolation ( void  )

virtual

Implements ElementInput2.

Definition at line 295 of file SegInput2.cpp.

                                          {/*{{{*/
   if(this->interpolation==P0Enum || this->interpolation==P0DGEnum){
      return P0Enum;
   }
   return P1Enum;
}/*}}}*/

GetResultNumberOfNodes()

int SegInput2::GetResultNumberOfNodes ( void  )

virtual

Implements ElementInput2.

Definition at line 301 of file SegInput2.cpp.

                                          {/*{{{*/
   return this->N;
}

Reset()

void SegInput2::Reset

(

int

interp_in

)

Definition at line 45 of file SegInput2.cpp.

                                  {/*{{{*/
 
   /*Clean up*/
   if(this->values)         xDelete<IssmDouble>(this->values);
   if(this->element_values) xDelete<IssmDouble>(this->element_values);
 
   /*Set interpolation*/
   this->interpolation  = interp_in;
 
   /*Create Sizes*/
   if(this->interpolation==P1Enum){
      this->M = this->numberofvertices_local;
      this->N = 1;
   }
   else{
      this->M = this->numberofelements_local;
      this->N = SegRef::NumberofNodes(interp_in);
   }
 
   /*Allocate Pointers*/
   this->values         = xNewZeroInit<IssmDouble>(this->M*this->N);
   this->element_values = xNewZeroInit<IssmDouble>(SegRef::NumberofNodes(interp_in));
}/*}}}*/

---

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

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