Reducevectorgtofx.h File Reference

reduce petsc vector from g set to f set (free dofs), using the nodeset partitioning vectors. More...

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

Go to the source code of this file.

Functions
void	Reducevectorgtofx (Vector< IssmDouble > **puf, Vector< IssmDouble > *ug, Nodes *nodes, Parameters *parameters)

Detailed Description

reduce petsc vector from g set to f set (free dofs), using the nodeset partitioning vectors.

Definition in file Reducevectorgtofx.h.

Function Documentation

Reducevectorgtofx()

void Reducevectorgtofx	(	Vector< IssmDouble > **	puf,
		Vector< IssmDouble > *	ug,
		Nodes *	nodes,
		Parameters *	parameters
	)

Definition at line 8 of file Reducevectorgtofx.cpp.

                                                                                                             {
 
   if(VerboseModule()) _printf0_("   Reduce vector from g to f set\n");
 
   /*first figure out fsize: */
   int fsize=nodes->NumberOfDofs(FsetEnum);
   int flocalsize = nodes->NumberOfDofsLocal(FsetEnum);
 
   /*If fsize is 0, return NULL vector*/
   if(fsize==0){
      *puf=NULL;
      return;
   }
 
   /*Get local vectors ug*/
   int        *indices_ug = NULL;
   IssmDouble *local_ug   = NULL;
   ug->GetLocalVector(&local_ug,&indices_ug);
 
   /*Allocate output*/
   Vector<IssmDouble>* uf=new Vector<IssmDouble>(flocalsize,fsize);
 
   /*Let nodes figure it out*/
   for(int i=0;i<nodes->Size();i++){
      Node* node=(Node*)nodes->GetObjectByOffset(i);
      node->VecReduce(uf,local_ug,indices_ug);
   }
 
   /*Assemble vector: */
   uf->Assemble();
 
   /*Cleanup and assing output pointer*/
   xDelete<int>(indices_ug);
   xDelete<IssmDouble>(local_ug);
   *puf=uf;
}