NodalValuex.h File Reference

header file for NodalValuex More...

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

Go to the source code of this file.

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

Detailed Description

header file for NodalValuex

Definition in file NodalValuex.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;
}