Gallery

Collection of example matrices. More...

Functions
template<typename Method , typename MatrixType >
void	cusp::gallery::diffusion (MatrixType &matrix, size_t m, size_t n, const double eps=1e-5, const double theta=M_PI/4.0)
template<typename MatrixType >
void	cusp::gallery::grid2d (MatrixType &matrix, size_t m, size_t n)
template<typename MatrixType >
void	cusp::gallery::grid3d (MatrixType &matrix, size_t m, size_t n, size_t l)
template<typename MatrixType >
void	cusp::gallery::poisson5pt (MatrixType &matrix, size_t m, size_t n)
template<typename MatrixType >
void	cusp::gallery::poisson9pt (MatrixType &matrix, size_t m, size_t n)
template<typename MatrixType >
void	cusp::gallery::poisson7pt (MatrixType &matrix, size_t m, size_t n, size_t k)
template<typename MatrixType >
void	cusp::gallery::poisson27pt (MatrixType &matrix, size_t m, size_t n, size_t l)
template<typename MatrixType >
void	cusp::gallery::random (MatrixType &matrix, size_t m, size_t n, size_t num_samples)

Function Documentation

template<typename Method , typename MatrixType >

void cusp::gallery::diffusion	(	MatrixType &	matrix,
		size_t	m,
		size_t	n,
		const double	eps = 1e-5,
		const double	theta = M_PI/4.0
	)

diffusion: Create a matrix representing an anisotropic Poisson problem discretized on an m by n grid with the a given direction.

Parameters

matrix	output
m	number of grid rows
n	number of grid columns
eps	magnitude of anisotropy
theta	angle of anisotropy in radians

Template Parameters

MatrixType

matrix container

template<typename MatrixType >

void cusp::gallery::grid2d	(	MatrixType &	matrix,
		size_t	m,
		size_t	n
	)

grid2d: Create a matrix representing a 2d m by n grid.

Parameters

matrix	output
m	number of grid rows
n	number of grid columns

Template Parameters

MatrixType

matrix container

#include <cusp/gallery/grid.h>
#include <cusp/coo_matrix.h>
#include <cusp/print.h>

int main(void)
{
    cusp::coo_matrix<int, float, cusp::device_memory> A;

    // create a matrix for a 4x4 grid
    cusp::gallery::grid2d(A, 4, 4);

    // print matrix
    cusp::print(A);

    return 0;
}

template<typename MatrixType >

void cusp::gallery::grid3d	(	MatrixType &	matrix,
		size_t	m,
		size_t	n,
		size_t	l
	)

grid3d: Create a matrix representing a 3d m by n by l grid.

Parameters

matrix	output
m	number of grid rows
n	number of grid columns
l	number of grid layers

Template Parameters

MatrixType

matrix container

#include <cusp/gallery/grid.h>
#include <cusp/coo_matrix.h>
#include <cusp/print.h>

int main(void)
{
    cusp::coo_matrix<int, float, cusp::device_memory> A;

    // create a matrix for a 4x4x4 grid
    cusp::gallery::grid3d(A, 4, 4, 4);

    // print matrix
    cusp::print(A);

    return 0;
}

template<typename MatrixType >

void cusp::gallery::poisson27pt	(	MatrixType &	matrix,
		size_t	m,
		size_t	n,
		size_t	l
	)

poisson27pt: Create a matrix representing a 27pt Poisson problem discretized on an m by n by l grid with the standard 3D 27-point finite-difference stencil.

Parameters

matrix	output
m	number of grid rows
n	number of grid columns
l	number of grid layers

Template Parameters

MatrixType

matrix container

#include <cusp/gallery/poisson.h>
#include <cusp/coo_matrix.h>
#include <cusp/print.h>

int main(void)
{
    cusp::coo_matrix<int, float, cusp::device_memory> A;

    // create a matrix for a Poisson problem on a 4x4x4 grid
    cusp::gallery::poisson27pt(A, 4, 4, 4);

    // print matrix
    cusp::print(A);

    return 0;
}

template<typename MatrixType >

void cusp::gallery::poisson5pt	(	MatrixType &	matrix,
		size_t	m,
		size_t	n
	)

poisson5pt: Create a matrix representing a 5pt Poisson problem discretized on an m by n grid with the standard 2D 5-point finite-difference stencil.

Parameters

matrix	output
m	number of grid rows
n	number of grid columns

Template Parameters

MatrixType

matrix container

#include <cusp/gallery/poisson.h>
#include <cusp/coo_matrix.h>
#include <cusp/print.h>

int main(void)
{
    cusp::coo_matrix<int, float, cusp::device_memory> A;

    // create a matrix for a Poisson problem on a 4x4 grid
    cusp::gallery::poisson5pt(A, 4, 4);

    // print matrix
    cusp::print(A);

    return 0;
}

Examples:

ainv.cu, bicgstab.cu, cg.cu, cg_m.cu, cr.cu, gmres.cu, maximal_independent_set.cu, monitor.cu, poisson.cu, smoothed_aggregation.cu, and verbose_monitor.cu.

template<typename MatrixType >

void cusp::gallery::poisson7pt	(	MatrixType &	matrix,
		size_t	m,
		size_t	n,
		size_t	k
	)

poisson7pt: Create a matrix representing a 7pt Poisson problem discretized on an m by n by k grid with the standard 3D 7-point finite-difference stencil.

Parameters

matrix	output
m	number of grid rows
n	number of grid columns
k	number of grid layers

Template Parameters

MatrixType

matrix container

#include <cusp/gallery/poisson.h>
#include <cusp/coo_matrix.h>
#include <cusp/print.h>

int main(void)
{
    cusp::coo_matrix<int, float, cusp::device_memory> A;

    // create a matrix for a Poisson problem on a 4x4x4 grid
    cusp::gallery::poisson7pt(A, 4, 4, 4);

    // print matrix
    cusp::print(A);

    return 0;
}

template<typename MatrixType >

void cusp::gallery::poisson9pt	(	MatrixType &	matrix,
		size_t	m,
		size_t	n
	)

poisson9pt: Create a matrix representing a 9pt Poisson problem discretized on an m by n grid with the standard 2D 9-point finite-difference stencil.

Parameters

matrix	output
m	number of grid rows
n	number of grid columns

Template Parameters

MatrixType

matrix container

#include <cusp/gallery/poisson.h>
#include <cusp/coo_matrix.h>
#include <cusp/print.h>

int main(void)
{
    cusp::coo_matrix<int, float, cusp::device_memory> A;

    // create a matrix for a Poisson problem on a 4x4 grid
    cusp::gallery::poisson9pt(A, 4, 4);

    // print matrix
    cusp::print(A);

    return 0;
}

template<typename MatrixType >

void cusp::gallery::random	(	MatrixType &	matrix,
		size_t	m,
		size_t	n,
		size_t	num_samples
	)

random: Create a matrix with random connections

Parameters

matrix	output
m	number of grid rows
n	number of grid columns
num_samples	number of random edges

Template Parameters

MatrixType

matrix container

#include <cusp/gallery/random.h>
#include <cusp/coo_matrix.h>
#include <cusp/print.h>

int main(void)
{
    cusp::coo_matrix<int, float, cusp::device_memory> A;

    // create a matrix for a Poisson problem on a 4x4 grid
    cusp::gallery::random(A, 4, 4, 12);

    // print matrix
    cusp::print(A);

    return 0;
}