Mergesolutionfromftogx.h File Reference

merge solution back from f set into g set More...

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

Go to the source code of this file.

Functions
void	Mergesolutionfromftogx (Vector< IssmDouble > **pug, Vector< IssmDouble > *uf, Vector< IssmDouble > *ys, Nodes *nodes, Parameters *parameters, bool flag_ys0=false)

Detailed Description

merge solution back from f set into g set

Definition in file Mergesolutionfromftogx.h.

Function Documentation

Mergesolutionfromftogx()

void Mergesolutionfromftogx	(	Vector< IssmDouble > **	pug,
		Vector< IssmDouble > *	uf,
		Vector< IssmDouble > *	ys,
		Nodes *	nodes,
		Parameters *	parameters,
		bool	flag_ys0 = false
	)

Definition at line 8 of file Mergesolutionfromftogx.cpp.

                                                                                                                                                            {
 
   /*Display message*/
   if(VerboseModule()) _printf0_("   Merging solution vector from fset to gset\n");
 
   /*first, get gsize, fsize and ssize: */
   int gsize=nodes->NumberOfDofs(GsetEnum);
   int gsize_local=nodes->NumberOfDofsLocal(GsetEnum);
   int fsize=nodes->NumberOfDofs(FsetEnum);
   int ssize=nodes->NumberOfDofs(SsetEnum);
 
   /*serialize uf and ys: those two vectors will be indexed by the nodes, who are the only ones 
    *that know which values should be plugged into ug and where: */
   if(ssize) if(flag_ys0) ys->Set(0.0);
 
   /*Get local vectors ys and uf*/
   int        *indices_ys = NULL;
   IssmDouble *local_ys   = NULL;
   ys->GetLocalVector(&local_ys,&indices_ys);
   int        *indices_uf = NULL;
   IssmDouble *local_uf   = NULL;
   uf->GetLocalVector(&local_uf,&indices_uf);
 
   /*initialize ug: */
   Vector<IssmDouble>* ug=new Vector<IssmDouble>(gsize_local,gsize);
 
   /*Let nodes figure it out*/
   for(int i=0;i<nodes->Size();i++){
      Node* node=(Node*)nodes->GetObjectByOffset(i);
      node->VecMerge(ug,local_uf,indices_uf,local_ys,indices_ys);
   }
 
   /*Assemble vector: */
   ug->Assemble();
 
   /*Cleanup and assign output pointer*/
   xDelete<int>(indices_uf);
   xDelete<int>(indices_ys);
   xDelete<IssmDouble>(local_uf);
   xDelete<IssmDouble>(local_ys);
   *pug=ug;
}