codelet_dzcpy.c

Go to the documentation of this file.

#include "../include/starpu_exageostat.h"

static void CORE_zcpy_starpu(void *buffers[], void *cl_arg){
        int m;
        double * A;
        int m0;
        double * r;

        A   = (double *)STARPU_MATRIX_GET_PTR(buffers[0]);
        starpu_codelet_unpack_args(cl_arg,  &m, &m0, &r);

        core_dzcpy(A,  m,  m0, r);
}

static struct starpu_codelet cl_zcpy =
{
                .where      = STARPU_CPU,
                .cpu_funcs  = {CORE_zcpy_starpu},
                .nbuffers   = 1,
                .modes      = {STARPU_W},
        .name       = "dzcpy"
};

/***************************************************************************/
int MORSE_MLE_zcpy_Tile_Async(MORSE_desc_t *descA, double * r,MORSE_sequence_t *sequence, MORSE_request_t  *request) {
        MORSE_context_t *morse;
        MORSE_option_t options;
        morse = morse_context_self();
        if (sequence->status != MORSE_SUCCESS)
                return -2;
        RUNTIME_options_init(&options, morse, sequence, request);

        int m, m0;
        int tempmm;
        MORSE_desc_t A = *descA;
        struct starpu_codelet *cl=&cl_zcpy;

        for (m = 0; m < A.mt; m++) {
                tempmm = m == A.mt-1 ? A.m - m * A.mb : A.mb;
                m0 = m * A.mb;

                starpu_insert_task(starpu_mpi_codelet(cl),
                                STARPU_VALUE, &tempmm, sizeof(int),
                                STARPU_VALUE, &m0, sizeof(int),
                                STARPU_VALUE, &r, sizeof(double),
                                STARPU_W, RTBLKADDR(descA, sizeof(double)*tempmm, m, 0),
                             #if defined(CHAMELEON_CODELETS_HAVE_NAME)
                                 STARPU_NAME, "dzcpy",
                                 #endif
                     0);

        }

    //MORSE_TASK_flush_desc( &options, MorseUpperLower, descA);
        RUNTIME_options_ws_free(&options);
        //MORSE_TASK_dataflush_all();
    return MORSE_SUCCESS;
}