Xxt Namespace Reference

Classes
struct	comm
struct	sparse_cholesky
struct	csr_mat
struct	crs_data

Typedefs
typedef int	comm_ext
typedef int	comm_req

Functions
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)
void	nektar_crs_solve (double *x, struct crs_data *data, double *b)
void	nektar_crs_stats (struct crs_data *data)
void	nektar_crs_free (struct crs_data *data)
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.
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.
static void	Finalise (crs_data *pCrs)
	Deallocates the crs mapping data.

Typedef Documentation

typedef int Xxt::comm_ext

Definition at line 57 of file Xxt.hpp.

typedef int Xxt::comm_req

Definition at line 58 of file Xxt.hpp.

Function Documentation

static void Xxt::Finalise ( crs_data *  pCrs )

inlinestatic

Deallocates the crs mapping data.

Definition at line 198 of file Xxt.hpp.

    {
#ifdef NEKTAR_USE_MPI
        if (pCrs)
        {
            //nektar_crs_free(pCrs);
        }
#endif
    }

static struct crs_data* Xxt::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  )

staticread

Initialise the matrix-solve.

On each process an array of IDs for each global degree of freedom is supplied which corresponds to a unique numbering of universal degrees of freedom. Three vectors describing the matrix are also provided. The parallel matrix solve is then set up.

Parameters

pId	Array of integers providing universal IDs for each global DOF on the process.
pAi	Row indices of matrix entries
pAj	Column indices of matrix entries
pAr	Values of matrix entries
pComm	Communication object used for inter-process communication.

Returns

crs_data structure

Definition at line 157 of file Xxt.hpp.

References ASSERTL1, Xxt::comm::c, Xxt::comm::id, nektar_crs_setup(), and Xxt::comm::np.

    {
#ifdef NEKTAR_USE_MPI
        unsigned int nz = pAr.num_elements();
        LibUtilities::CommMpiSharedPtr vCommMpi = boost::dynamic_pointer_cast<LibUtilities::CommMpi> (pComm);
        ASSERTL1(vCommMpi, "Failed to cast MPI Comm object.");
        comm vComm;
        MPI_Comm_dup(vCommMpi->GetComm(), &vComm.c);
        vComm.id = vCommMpi->GetRank();
        vComm.np = vCommMpi->GetSize();
        return nektar_crs_setup(pRank, &pId[0], nz, &pAi[0], &pAj[0], &pAr[0], 0, &vComm);
#else
        return 0;
#endif
    }

void Xxt::nektar_crs_free

(

struct crs_data *

data

)

struct crs_data* Xxt::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  )

read

Referenced by Init().

void Xxt::nektar_crs_solve	(	double *	x,
		struct crs_data *	data,
		double *	b
	)

Referenced by Solve().

void Xxt::nektar_crs_stats

(

struct crs_data *

data

)

Referenced by Nektar::MultiRegions::GlobalLinSysXxtFull::AssembleMatrixArrays(), and Nektar::MultiRegions::GlobalLinSysXxtStaticCond::v_AssembleSchurComplement().

static void Xxt::Solve ( Nektar::Array< OneD, NekDouble >  pX, struct crs_data *  pCrs, Nektar::Array< OneD, NekDouble >  pB  )

inlinestatic

Solve the matrix system for a given input vector b.

Definition at line 182 of file Xxt.hpp.

References nektar_crs_solve().

    {
#ifdef NEKTAR_USE_MPI
        if (!pCrs) {
            return;
        }
        nektar_crs_solve(&pX[0], pCrs, &pB[0]);
#endif
    }