AverageOntoPartitionx.h File Reference

header file for averaging dakota responses onto a vertices partition More...

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

Go to the source code of this file.

Functions
void	AverageOntoPartitionx (double **paverage, Elements *elements, Nodes *nodes, Vertices *vertices, Loads *loads, Materials *materials, Parameters *parameters, double *vertex_response, double *qmu_part, int npart)

Detailed Description

header file for averaging dakota responses onto a vertices partition

Definition in file AverageOntoPartitionx.h.

Function Documentation

AverageOntoPartitionx()

void AverageOntoPartitionx	(	double **	paverage,
		Elements *	elements,
		Nodes *	nodes,
		Vertices *	vertices,
		Loads *	loads,
		Materials *	materials,
		Parameters *	parameters,
		double *	vertex_response,
		double *	qmu_part,
		int	npart
	)

Definition at line 16 of file AverageOntoPartitionx.cpp.

                                                                                                                                                                                                                  {
 
   /*output: */
   Vector<IssmDouble>* partition_contributions=NULL;
   Vector<IssmDouble>* partition_areas=NULL;
   Vector<IssmDouble>* vec_average=NULL;
   double* average=NULL;
 
   /*average onto the separate areas. The result will be a npart sized vector. */
 
   /*allocate: */
   partition_contributions=new Vector<IssmDouble>(npart);
   partition_areas=new Vector<IssmDouble>(npart);
   vec_average=new Vector<IssmDouble>(npart);
 
   /*loop on each element, and add contribution of the element to the partition (surface weighted average): */
   for(int i=0;i<elements->Size();i++){
      Element* element=xDynamicCast<Element*>(elements->GetObjectByOffset(i));
      element->AverageOntoPartition(partition_contributions,partition_areas,vertex_response,qmu_part);
   }
 
   /*Assemble: */
   partition_contributions->Assemble();
   partition_areas->Assemble();
 
   /*We have the partition_areas and the partition_contributions for each partition -> compute the surfae weighted average: */
   vec_average->PointwiseDivide(partition_contributions,partition_areas);
 
   /*serialize:*/
   average=vec_average->ToMPISerial();
 
   /*Free ressources:*/
   delete partition_contributions;
   delete partition_areas;
   delete vec_average;
 
   /*Assign output pointers:*/
   *paverage=average;
}