codelet_zplgsy.c

Go to the documentation of this file.

#include "morse_starpu.h"

/*
 * Codelet CPU
 */
static void cl_zplgsy_cpu_func(void *descr[], void *cl_arg)
{
    PLASMA_Complex64_t bump;
    int m;
    int n;
    PLASMA_Complex64_t *A;
    int lda;
    int bigM;
    int m0;
    int n0;
    unsigned long long int seed;

    starpu_unpack_cl_args( cl_arg, &bump, &m, &n, &lda, &bigM, &m0, &n0, &seed );

    A = (PLASMA_Complex64_t *)STARPU_MATRIX_GET_PTR(descr[0]);

    CORE_zplgsy( bump, m, n, A, lda, bigM, m0, n0, seed );
}

/*
 * Codelet Multi-cores
 */
#if defined(MORSE_USE_MULTICORE) && 0
static void cl_zplgsy_mc_func(void *descr[], void *cl_arg)
{
}
#endif

/*
 * Codelet GPU
 */
#if (defined MORSE_USE_CUDA) && 0
static void cl_zplgsy_cuda_func(void *descr[], void *cl_arg)
{
    cudaThreadSynchronize();
}
#endif

/*
 * Codelet definition
 */
CODELETS_CPU(zplgsy, 1, cl_zplgsy_cpu_func)

/*
 * Wrapper
 */
void MORSE_zplgsy( MorseOption_t *options, 
                   PLASMA_Complex64_t bump, int m, int n, magma_desc_t *A, int Am, int An,
                   int bigM, int m0, int n0, unsigned long long int seed )
{
    starpu_codelet *zplgsy_codelet;
    void (*callback)(void*) = options->profiling ? cl_zplgsy_callback : NULL;
    int lda = BLKLDD( A, Am );

#ifdef MORSE_USE_MULTICORE
    zplgsy_codelet = options->parallel ? &cl_zplgsy_mc : &cl_zplgsy;
#else
    zplgsy_codelet = &cl_zplgsy;
#endif
    
    starpu_Insert_Task(
            &cl_zplgsy,
            VALUE,  &bump,  sizeof(PLASMA_Complex64_t),
            VALUE,  &m,     sizeof(int),
            VALUE,  &n,     sizeof(int), 
            OUTPUT,  BLKADDR( A, PLASMA_Complex64_t, Am, An ),
            VALUE,  &lda,   sizeof(int),
            VALUE,  &bigM,  sizeof(int),
            VALUE,  &m0,    sizeof(int), 
            VALUE,  &n0,    sizeof(int),
            VALUE,  &seed,  sizeof(unsigned long long int), 
            PRIORITY, options->priority,
            CALLBACK, callback, NULL,
            0);
}