TransferFromPETSc.cpp File Reference

Copy from PETSc vector to HyPar array. More...

#include <stdlib.h>
#include <vector>
#include <basic.h>
#include <arrayfunctions.h>
#include <simulation_object.h>
#include <petscinterface_struct.h>

Go to the source code of this file.

Macros
#define	__FUNCT__   "TransferVecFromPETSc"

Functions
int	TransferVecFromPETSc (double *const u, const Vec Y, void *ctxt, const int sim_idx, const int offset)

Detailed Description

Copy from PETSc vector to HyPar array.

Author

Debojyoti Ghosh

Definition in file TransferFromPETSc.cpp.

Macro Definition Documentation

#define __FUNCT__   "TransferVecFromPETSc"

Definition at line 16 of file TransferFromPETSc.cpp.

Function Documentation

int TransferVecFromPETSc	(	double *const	u,
		const Vec	Y,
		void *	ctxt,
		const int	sim_idx,
		const int	offset
	)

Copy data from a PETSc vector (used by PETSc time integrators, and with no ghost points) to a HyPar::u array (with ghost points).

See Also

TransferVecToPETSc()

Parameters

u	HyPar::u type array (with ghost points)
Y	PETSc vector
ctxt	Object of type PETScContext
sim_idx	Simulation object index
offset	Offset

Definition at line 23 of file TransferFromPETSc.cpp.

{
  PETScContext* context = (PETScContext*) ctxt;
  SimulationObject* sim = (SimulationObject*) context->simobj;
  const double* Yarr;

  PetscFunctionBegin;
  VecGetArrayRead(Y,&Yarr);
  std::vector<int> index(sim[sim_idx].solver.ndims,0);
  ArrayCopynD(  sim[sim_idx].solver.ndims,
                (Yarr+offset),
                u,
                sim[sim_idx].solver.dim_local,
                0,
                sim[sim_idx].solver.ghosts,
                index.data(),
                sim[sim_idx].solver.nvars );
  VecRestoreArrayRead(Y,&Yarr);

  PetscFunctionReturn(0);
}