Nodalvalue.h File Reference

: header file for Nodalvalue object More...

#include "./Definition.h"
#include "./FemModel.h"

Go to the source code of this file.

Data Structures
class	Nodalvalue

Functions
void	NodalValuex (IssmDouble *pnodalvalue, int natureofdataenum, Elements *elements, Nodes *nodes, Vertices *vertices, Loads *loads, Materials *materials, Parameters *parameters)
IssmDouble	OutputDefinitionsResponsex (FemModel *femmodel, int output_enum)

Detailed Description

: header file for Nodalvalue object

Definition in file Nodalvalue.h.

Function Documentation

NodalValuex()

void NodalValuex	(	IssmDouble *	pnodalvalue,
		int	natureofdataenum,
		Elements *	elements,
		Nodes *	nodes,
		Vertices *	vertices,
		Loads *	loads,
		Materials *	materials,
		Parameters *	parameters
	)

Definition at line 10 of file NodalValuex.cpp.

                                                                                                                                                                                {
 
   IssmDouble value;
   int        index;
   int        found,sumfound,cpu_found,cpu;
 
   /*retrieve element we are interested in: */
   parameters->FindParam(&index,IndexEnum);
 
   /*This is the vertex id for which we want to collect the data. Go through elements, and for each 
    *element, figure out  if they hold the vertex, and the data. If so, return it: */
   cpu_found=-1;
   found=0;
 
   for(int i=0;i<elements->Size();i++){
      Element* element=xDynamicCast<Element*>(elements->GetObjectByOffset(i));
      found=element->NodalValue(&value,index,natureofdataenum);
      if(found){
         cpu_found=IssmComm::GetRank();
         break;
      }
   }
 
   /*Broadcast whether we found the element: */
   ISSM_MPI_Allreduce(&found,&sumfound,1,ISSM_MPI_INT,ISSM_MPI_SUM,IssmComm::GetComm());
   if(!sumfound)_error_("could not find element with vertex with id " << index << " to compute nodal value " << EnumToStringx(natureofdataenum));
 
   /*Broadcast and plug into response: */
   ISSM_MPI_Allreduce ( &cpu_found,&cpu,1,ISSM_MPI_INT,ISSM_MPI_MAX,IssmComm::GetComm());
   ISSM_MPI_Bcast(&value,1,ISSM_MPI_DOUBLE,cpu,IssmComm::GetComm()); 
 
   *pnodalvalue=value;
}

OutputDefinitionsResponsex()

IssmDouble OutputDefinitionsResponsex	(	FemModel *	femmodel,
		int	output_enum
	)

Definition at line 38 of file OutputDefinitionsResponsex.cpp.

                                                                         {
 
   /*Ok, go find the output definitions dataset in the parameters, where our responses are hiding: */
   DataSet* output_definitions=((DataSetParam*)femmodel->parameters->FindParamObject(OutputdefinitionEnum))->value;
 
   /*Now, go through the output definitions, and retrieve the object which corresponds to our requested response, output_enum: */
   for(int i=0;i<output_definitions->Size();i++){
 
      //Definition* definition=xDynamicCast<Definition*>(output_definitions->GetObjectByOffset(i));
      Definition* definition=dynamic_cast<Definition*>(output_definitions->GetObjectByOffset(i));
 
      int en = definition->DefinitionEnum();
      if(en==output_enum){
 
         /*This is the object that we have been chasing for. compute the response and return: */
         IssmDouble return_value=definition->Response(femmodel);
 
         /*return:*/
         return return_value;
      }
   }
 
   /*If we are here, did not find the definition for this response, not good!: */
   _error_("Could not find the response for output definition " << EnumToStringx(output_enum) 
            <<" ("<<output_enum<<")"
            << " because could not find the definition itself!");
}