VertexCoordinatesx.cpp File Reference

#include "./VertexCoordinatesx.h"
#include "../../shared/shared.h"
#include "../../toolkits/toolkits.h"

Go to the source code of this file.

Functions
void	VertexCoordinatesx (IssmDouble **px, IssmDouble **py, IssmDouble **pz, Vertices *vertices, bool spherical)

Function Documentation

VertexCoordinatesx()

void VertexCoordinatesx	(	IssmDouble **	px,
		IssmDouble **	py,
		IssmDouble **	pz,
		Vertices *	vertices,
		bool	spherical
	)

Definition at line 11 of file VertexCoordinatesx.cpp.

                                                                                                               {
 
   /*output: */
   IssmDouble* x=NULL;
   IssmDouble* y=NULL;
   IssmDouble* z=NULL;
 
   Vector<IssmDouble>* vx=NULL;
   Vector<IssmDouble>* vy=NULL;
   Vector<IssmDouble>* vz=NULL;
 
   /*intermediary: */
   int  numberofvertices;
   int  i;
 
   /*figure out how many vertices we have: */
   numberofvertices=vertices->NumberOfVertices();
 
   vx=new Vector<IssmDouble>(numberofvertices);
   vy=new Vector<IssmDouble>(numberofvertices);
   vz=new Vector<IssmDouble>(numberofvertices);
 
   /*march through our vertices: */
   for(i=0;i<vertices->Size();i++){
      Vertex* vertex=(Vertex*)vertices->GetObjectByOffset(i);
      vertex->VertexCoordinates(vx,vy,vz,spherical);
   }
 
   /*Assemble*/
   vx->Assemble();
   vy->Assemble();
   vz->Assemble();
 
   /*serialize: */
   x=vx->ToMPISerial();
   y=vy->ToMPISerial();
   z=vz->ToMPISerial();
 
   /*Free ressources: */
   delete vx;
   delete vy;
   delete vz;
 
   /*output: */
   if (px)*px=x;
   else xDelete<IssmDouble>(x);
   if (py)*py=y;
   else xDelete<IssmDouble>(y);
   if (pz)*pz=z;
   else xDelete<IssmDouble>(z);
}