#include <stdlib.h>
#include "common.h"

Include dependency graph for descriptor.c:

Functions
magma_desc_t	magma_desc_init (MAGMA_enum dtyp, int mb, int nb, int bsiz, int lm, int ln, int i, int j, int m, int n)
magma_desc_t	magma_desc_submatrix (magma_desc_t descA, int i, int j, int m, int n)
int	magma_desc_check (magma_desc_t *desc)
int	magma_desc_mat_alloc (magma_desc_t *desc)
int	magma_desc_mat_free (magma_desc_t *desc)
int	MAGMA_Desc_Create (magma_desc_t *desc, void mat, MAGMA_enum dtyp, int mb, int nb, int bsiz, int lm, int ln, int i, int j, int m, int n)
int	MAGMA_Desc_Destroy (magma_desc_t **desc)

Function Documentation

int magma_desc_check ( magma_desc_t * desc )

Check for descriptor correctness

Definition at line 88 of file descriptor.c.

References magma_desc_s::desc, MAGMA_ERR_ILLEGAL_VALUE, MAGMA_ERR_UNALLOCATED, magma_error, and MAGMA_SUCCESS.

{
    if (desc->desc.mat == NULL) {
        magma_error("magma_desc_check", "NULL matrix pointer");
        return MAGMA_ERR_UNALLOCATED;
    }
    if (desc->desc.dtyp != PlasmaRealFloat &&
        desc->desc.dtyp != PlasmaRealDouble &&
        desc->desc.dtyp != PlasmaComplexFloat &&
        desc->desc.dtyp != PlasmaComplexDouble  ) {
        magma_error("magma_desc_check", "invalid matrix type");
        return MAGMA_ERR_ILLEGAL_VALUE;
    }
    if (desc->desc.mb <= 0 || desc->desc.nb <= 0) {
        magma_error("magma_desc_check", "negative tile dimension");
        return MAGMA_ERR_ILLEGAL_VALUE;
    }
    if (desc->desc.bsiz < desc->desc.mb*desc->desc.nb) {
        magma_error("magma_desc_check", "tile memory size smaller than the product of dimensions");
        return MAGMA_ERR_ILLEGAL_VALUE;
    }
    if (desc->desc.lm <= 0 || desc->desc.ln <= 0) {
        magma_error("magma_desc_check", "negative matrix dimension");
        return MAGMA_ERR_ILLEGAL_VALUE;
    }
    if (desc->desc.i >= desc->desc.lm || desc->desc.j >= desc->desc.ln) {
        magma_error("magma_desc_check", "beginning of the matrix out of scope");
        return MAGMA_ERR_ILLEGAL_VALUE;
    }
    if (desc->desc.i+desc->desc.m > desc->desc.lm || desc->desc.j+desc->desc.n > desc->desc.ln) {
        magma_error("magma_desc_check", "submatrix out of scope");
        return MAGMA_ERR_ILLEGAL_VALUE;
    }
    return MAGMA_SUCCESS;
}

Here is the caller graph for this function:

int MAGMA_Desc_Create	(	magma_desc_t **	desc,
		void *	mat,
		MAGMA_enum	dtyp,
		int	mb,
		int	nb,
		int	bsiz,
		int	lm,
		int	ln,
		int	i,
		int	j,
		int	m,
		int	n
	)

MAGMA_Desc_Create - Create matrix descriptor.

Parameters:

[out]	desc	On exit, descriptor of the matrix.
[in]	mat	Memory location of the matrix.
[in]	dtyp	Data type of the matrix: MagmaRealFloat: single precision real (S), MagmaRealDouble: double precision real (D), MagmaComplexFloat: single precision complex (C), MagmaComplexDouble: double precision complex (Z).
[in]	mb	Number of rows in a tile.
[in]	nb	Number of columns in a tile.
[in]	bsiz	Size in bytes including padding.
[in]	lm	Number of rows of the entire matrix.
[in]	ln	Number of columns of the entire matrix.
[in]	i	Row index to the beginning of the submatrix.
[in]	j	Column indes to the beginning of the submatrix.
[in]	m	Number of rows of the submatrix.
[in]	n	Number of columns of the submatrix.

Returns:

Return values:

MAGMA_SUCCESS successful exit

Definition at line 211 of file descriptor.c.

References magma_context_self(), magma_desc_check(), magma_desc_init(), MAGMA_ERR_NOT_INITIALIZED, MAGMA_ERR_OUT_OF_RESOURCES, magma_error, MAGMA_SUCCESS, and morse_desc_create().

{
    magma_context_t *magma;
    int status;
    magma = magma_context_self();
    if (magma == NULL) {
        magma_error("MAGMA_Desc_Create", "MAGMA not initialized");
        return MAGMA_ERR_NOT_INITIALIZED;
    }
    /* Allocate memory and initialize the descriptor */
    *desc = (magma_desc_t*)malloc(sizeof(magma_desc_t));
    if (*desc == NULL) {
        magma_error("MAGMA_Desc_Create", "malloc() failed");
        return MAGMA_ERR_OUT_OF_RESOURCES;
    }
    **desc = magma_desc_init(dtyp, mb, nb, bsiz, lm, ln, i, j, m, n);
    (**desc).desc.mat = mat;
    /* Create scheduler structure like registering data */
    morse_desc_create( *desc );
    status = magma_desc_check(*desc);
    if (status != MAGMA_SUCCESS) {
        magma_error("MAGMA_Desc_Create", "invalid descriptor");
        return status;
    }
   
    return MAGMA_SUCCESS;
}

Here is the call graph for this function:

Here is the caller graph for this function:

int MAGMA_Desc_Destroy ( magma_desc_t ** desc )

MAGMA_Desc_Destroy - Destroys matrix descriptor.

Parameters:

[in] desc Matrix descriptor.

Returns:

Return values:

MAGMA_SUCCESS successful exit

Definition at line 260 of file descriptor.c.

References magma_context_self(), MAGMA_ERR_NOT_INITIALIZED, MAGMA_ERR_UNALLOCATED, magma_error, MAGMA_SUCCESS, and morse_desc_destroy().

{
    magma_context_t *magma;
    magma = magma_context_self();
    if (magma == NULL) {
        magma_error("MAGMA_Desc_Destroy", "MAGMA not initialized");
        return MAGMA_ERR_NOT_INITIALIZED;
    }
    if (*desc == NULL) {
        magma_error("MAGMA_Desc_Destroy", "attempting to destroy a NULL descriptor");
        return MAGMA_ERR_UNALLOCATED;
    }
    /* Clean at scheduler level like unregistering data */
    morse_desc_destroy( *desc );
    free(*desc);
    *desc = NULL;
    return MAGMA_SUCCESS;
}

Here is the call graph for this function:

Here is the caller graph for this function:

magma_desc_t magma_desc_init	(	MAGMA_enum	dtyp,
		int	mb,
		int	nb,
		int	bsiz,
		int	lm,
		int	ln,
		int	i,
		int	j,
		int	m,
		int	n
	)

Internal static descriptor initializer

Definition at line 20 of file descriptor.c.

References magma_desc_s::desc, morse_desc_init(), and magma_desc_s::occurences.

{
    magma_desc_t desc;
    PLASMA_desc *pdesc = &(desc.desc);
    // Matrix address
    pdesc->mat = NULL;
    pdesc->A21 = (lm - lm%mb)*(ln - ln%nb);
    pdesc->A12 = (     lm%mb)*(ln - ln%nb) + pdesc->A21;
    pdesc->A22 = (lm - lm%mb)*(     ln%nb) + pdesc->A12;
    // Matrix properties
    pdesc->dtyp = dtyp;
    pdesc->mb = mb;
    pdesc->nb = nb;
    pdesc->bsiz = bsiz;
    // Large matrix parameters
    pdesc->lm = lm;
    pdesc->ln = ln;
    // Large matrix derived parameters
    pdesc->lm1 = (lm/mb);
    pdesc->ln1 = (ln/nb);
    pdesc->lmt = (lm%mb==0) ? (lm/mb) : (lm/mb+1);
    pdesc->lnt = (ln%nb==0) ? (ln/nb) : (ln/nb+1);
    // Submatrix parameters
    pdesc->i = i;
    pdesc->j = j;
    pdesc->m = m;
    pdesc->n = n;
    // Submatrix derived parameters
    pdesc->mt = (i+m-1)/mb - i/mb + 1;
    pdesc->nt = (j+n-1)/nb - j/nb + 1;
    desc.occurences = 0;
    morse_desc_init( &desc );
    return desc;
}

Here is the call graph for this function:

Here is the caller graph for this function:

int magma_desc_mat_alloc ( magma_desc_t * desc )

Definition at line 127 of file descriptor.c.

References magma_desc_s::desc, MAGMA_ERR_OUT_OF_RESOURCES, magma_error, MAGMA_SUCCESS, morse_desc_create(), and plasma_element_size().

{
    size_t size;
    size = (size_t)desc->desc.lm * (size_t)desc->desc.ln * (size_t)plasma_element_size(desc->desc.dtyp);
    if ((desc->desc.mat = malloc(size)) == NULL) {
        magma_error("magma_desc_mat_alloc", "malloc() failed");
        return MAGMA_ERR_OUT_OF_RESOURCES;
    }
    morse_desc_create( desc );
    return MAGMA_SUCCESS;
}

Here is the call graph for this function:

int magma_desc_mat_free ( magma_desc_t * desc )

Definition at line 145 of file descriptor.c.

References magma_desc_s::desc, MAGMA_SUCCESS, and morse_desc_destroy().

{
    morse_desc_destroy( desc );
    if (desc->desc.mat != NULL) {
        free(desc->desc.mat);
        desc->desc.mat = NULL;
    }
    return MAGMA_SUCCESS;
}

Here is the call graph for this function:

Here is the caller graph for this function:

magma_desc_t magma_desc_submatrix	(	magma_desc_t	descA,
		int	i,
		int	j,
		int	m,
		int	n
	)

Internal static descriptor initializer for submatrices

Definition at line 62 of file descriptor.c.

References magma_desc_s::desc, and morse_desc_submatrix().

{
    magma_desc_t descB;
    PLASMA_desc *pdescB = &(descB.desc);
    int mb, nb;
    descB = descA;
    mb = descA.desc.mb;
    nb = descA.desc.nb;
    // Submatrix parameters
    pdescB->i = i;
    pdescB->j = j;
    pdescB->m = m;
    pdescB->n = n;
    // Submatrix derived parameters
    pdescB->mt = (i+m-1)/mb - i/mb + 1;
    pdescB->nt = (j+n-1)/nb - j/nb + 1;
    morse_desc_submatrix( &descB );
    return descB;
}

Here is the call graph for this function:

Functions

Function Documentation