Reducevectorgtofx.cpp

Go to the documentation of this file.

 
#include "./Reducevectorgtofx.h"
 
void Reducevectorgtofx(Vector<IssmDouble>** puf, Vector<IssmDouble>* ug, Nodes* nodes,Parameters* parameters){
 
   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;
}