36 #ifndef NEKTAR_LIB_UTILITIES_COMMUNICATION_XXT_HPP
37 #define NEKTAR_LIB_UTILITIES_COMMUNICATION_XXT_HPP
70 unsigned int n, *Lrp, *Lj;
76 unsigned int n, *Arp, *Aj;
131 double *vl, *vc, *
vx, *combuf;
136 unsigned int nz,
const unsigned int *Ai,
137 const unsigned int *Aj,
const double *
A,
138 unsigned int null_space,
169 #ifdef NEKTAR_USE_MPI
170 unsigned int nz = pAr.num_elements();
173 ASSERTL1(vCommMpi,
"Failed to cast MPI Comm object.");
175 MPI_Comm_dup(vCommMpi->GetComm(), &vComm.
c);
176 vComm.
id = vCommMpi->GetRank();
177 vComm.
np = vCommMpi->GetSize();
180 MPI_Comm_free(&vComm.
c);
194 #ifdef NEKTAR_USE_MPI
208 #ifdef NEKTAR_USE_MPI
210 MPI_Finalized(&finalized);
211 if (pCrs && !finalized)
static void Solve(Nektar::Array< OneD, NekDouble > pX, struct crs_data *pCrs, Nektar::Array< OneD, NekDouble > pB)
Solve the matrix system for a given input vector b.
void nektar_crs_free(struct crs_data *data)
boost::shared_ptr< Comm > CommSharedPtr
Pointer to a Communicator object.
boost::shared_ptr< CommMpi > CommMpiSharedPtr
Pointer to a Communicator object.
static struct crs_data * Init(unsigned int pRank, const Nektar::Array< OneD, unsigned long > pId, const Nektar::Array< OneD, unsigned int > pAi, const Nektar::Array< OneD, unsigned int > pAj, const Nektar::Array< OneD, NekDouble > pAr, const LibUtilities::CommSharedPtr &pComm)
Initialise the matrix-solve.
struct crs_data * nektar_crs_setup(unsigned int n, const unsigned long *id, unsigned int nz, const unsigned int *Ai, const unsigned int *Aj, const double *A, unsigned int null_space, const struct comm *comm)
static void Finalise(crs_data *pCrs)
Deallocates the crs mapping data.
void nektar_crs_stats(struct crs_data *data)
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
void nektar_crs_solve(double *x, struct crs_data *data, double *b)