GetVectorFromControlInputsx.h File Reference

#include "../../classes/classes.h"

Go to the source code of this file.

Functions
void	GetVectorFromControlInputsx (Vector< IssmDouble > **pvector, Elements *elements, Nodes *nodes, Vertices *vertices, Loads *loads, Materials *materials, Parameters *parameters, const char *data="value")
void	GetVectorFromControlInputsx (IssmDouble **pvector, int *pN, Elements *elements, Nodes *nodes, Vertices *vertices, Loads *loads, Materials *materials, Parameters *parameters, const char *data="value")
void	GetPassiveVectorFromControlInputsx (double **pvector, int *pN, Elements *elements, Nodes *nodes, Vertices *vertices, Loads *loads, Materials *materials, Parameters *parameters, const char *data="value")

Function Documentation

GetVectorFromControlInputsx() [1/2]

void GetVectorFromControlInputsx	(	Vector< IssmDouble > **	pvector,
		Elements *	elements,
		Nodes *	nodes,
		Vertices *	vertices,
		Loads *	loads,
		Materials *	materials,
		Parameters *	parameters,
		const char *	data = "value"
	)

Definition at line 9 of file GetVectorFromControlInputsx.cpp.

                                                                                                                                                                                                {/*{{{*/
 
   bool isautodiff;
   parameters->FindParam(&isautodiff,AutodiffIsautodiffEnum);
   if(isautodiff){
      int*  N = NULL;
      int*  M = NULL;
      int  num_controls;
      int* control_type = NULL;
      Vector<IssmDouble>*  vector=NULL;
 
      /*Retrieve some parameters*/
      parameters->FindParam(&num_controls,InversionNumControlParametersEnum);
      parameters->FindParam(&control_type,NULL,InversionControlParametersEnum);
      parameters->FindParam(&N,NULL,ControlInputSizeNEnum);
      parameters->FindParam(&M,NULL,ControlInputSizeMEnum);
 
      /*1. Get vector size*/
      int size = 0;
      for(int i=0;i<num_controls;i++) size+=M[i]*N[i];
 
      /*2. Allocate vector*/
      vector=new Vector<IssmDouble>(size);
 
      /*3. Populate vector*/
      int offset = 0;
      for(int i=0;i<num_controls;i++){
         for(int j=0;j<elements->Size();j++){
            Element* element=(Element*)elements->GetObjectByOffset(j);
            element->GetVectorFromControlInputs(vector,control_type[i],i,data,offset);
         }
         offset += M[i]*N[i];
      }
 
      vector->Assemble();
 
      /*Assign output pointers:*/
      xDelete<int>(control_type);
      xDelete<int>(M);
      xDelete<int>(N);
 
      *pvector=vector;
   }
   else{
      int  num_controls;
      int* control_type = NULL;
      Vector<IssmDouble>*  vector=NULL;
 
      /*Retrieve some parameters*/
      parameters->FindParam(&num_controls,InversionNumControlParametersEnum);
      parameters->FindParam(&control_type,NULL,InversionControlParametersEnum);
 
      /*2. Allocate vector*/
      vector=new Vector<IssmDouble>(num_controls*vertices->NumberOfVertices());
 
      /*3. Populate vector*/
      int offset = 0;
      for(int i=0;i<num_controls;i++){
         for(int j=0;j<elements->Size();j++){
            Element* element=(Element*)elements->GetObjectByOffset(j);
            element->GetVectorFromControlInputs(vector,control_type[i],i,data);
         }
      }
      vector->Assemble();
 
      /*Assign output pointers:*/
      xDelete<int>(control_type);
      *pvector=vector;
   }
 
}/*}}}*/

GetVectorFromControlInputsx() [2/2]

void GetVectorFromControlInputsx	(	IssmDouble **	pvector,
		int *	pN,
		Elements *	elements,
		Nodes *	nodes,
		Vertices *	vertices,
		Loads *	loads,
		Materials *	materials,
		Parameters *	parameters,
		const char *	data = "value"
	)

Definition at line 80 of file GetVectorFromControlInputsx.cpp.

                                                                                                                                                                                                  {/*{{{*/
 
   /*intermediary: */
   int                 N;
   Vector<IssmDouble>* vec_vector=NULL;
 
   /*Get PETSc vector*/
   GetVectorFromControlInputsx( &vec_vector, elements,nodes, vertices, loads, materials, parameters,data);
 
   /*Serialize*/
   vec_vector->GetSize(&N);
   IssmDouble* vector=vec_vector->ToMPISerial();
   delete vec_vector;
 
   /*Assign output pointers:*/
   *pvector=vector;
   if(pN) *pN=N;
}/*}}}*/

GetPassiveVectorFromControlInputsx()

void GetPassiveVectorFromControlInputsx	(	double **	pvector,
		int *	pN,
		Elements *	elements,
		Nodes *	nodes,
		Vertices *	vertices,
		Loads *	loads,
		Materials *	materials,
		Parameters *	parameters,
		const char *	data = "value"
	)

Definition at line 127 of file GetVectorFromControlInputsx.cpp.

                                                                                                                                                                                                         {/*{{{*/
 
   GetVectorFromControlInputsx(pvector,pN,elements,nodes, vertices, loads, materials, parameters,data);
}/*}}}*/