core_sgetrip.c

Go to the documentation of this file.

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

#if defined(PLASMA_HAVE_WEAK)
#pragma weak CORE_sgetrip = PCORE_sgetrip
#define CORE_sgetrip PCORE_sgetrip
#endif
void CORE_sgetrip(int m, int n, float *A, float *W) {
    float t;
    int i, j;
    
    if( m != n ) {
        /* rectangular transposition (use workspace) */
        for (i=0; i<m; i++) {
            for (j=0; j<n; j++) {
                W[j+i*n] = A[i+j*m];
            }
        }
        memcpy(A, W, m*n*sizeof(float));
    }
    else {
        /* square transposition (swap pairwise) */
        for (i=0; i<m; i++) {
            for (j=i+1; j<n; j++) {
                t        = A[j+i*n];
                A[j+i*n] = A[i+j*m];
                A[i+j*m] = t;
            }
        }
    }
}

/***************************************************************************/
void QUARK_CORE_sgetrip(Quark *quark, Quark_Task_Flags *task_flags,
                        int m, int n, float *A, int szeA)
{
    DAG_CORE_GETRIP;
    QUARK_Insert_Task(quark, CORE_sgetrip_quark, task_flags,
        sizeof(int),                     &m,   VALUE,
        sizeof(int),                     &n,   VALUE,
        sizeof(float)*szeA, A,        INOUT,
        sizeof(float)*szeA, NULL,     SCRATCH,
        0);
}

/***************************************************************************/
#if defined(PLASMA_HAVE_WEAK)
#pragma weak CORE_sgetrip_quark = PCORE_sgetrip_quark
#define CORE_sgetrip_quark PCORE_sgetrip_quark
#endif
void CORE_sgetrip_quark(Quark *quark)
{
    int m;
    int n;
    float *A;
    float *W;

    quark_unpack_args_4(quark, m, n, A, W);
    CORE_sgetrip(m, n, A, W);
}

/***************************************************************************/
void QUARK_CORE_sgetrip_f1(Quark *quark, Quark_Task_Flags *task_flags,
                           int m, int n, 
                           float *A,    int szeA,
                           float *fake, int szeF, int paramF)
{
    DAG_CORE_GETRIP;
    QUARK_Insert_Task(
        quark, CORE_sgetrip_f1_quark, task_flags,
        sizeof(int),                     &m,   VALUE,
        sizeof(int),                     &n,   VALUE,
        sizeof(float)*szeA, A,        INOUT,
        sizeof(float)*szeA, NULL,     SCRATCH,
        sizeof(float)*szeF, fake,     paramF,
        0);
}

/***************************************************************************/
#if defined(PLASMA_HAVE_WEAK)
#pragma weak CORE_sgetrip_f1_quark = PCORE_sgetrip_f1_quark
#define CORE_sgetrip_f1_quark PCORE_sgetrip_f1_quark
#endif
void CORE_sgetrip_f1_quark(Quark *quark)
{
    int m;
    int n;
    float *A;
    float *W;
    float *fake;

    quark_unpack_args_5(quark, m, n, A, W, fake);
    CORE_sgetrip(m, n, A, W);
}

/***************************************************************************/
void QUARK_CORE_sgetrip_f2(Quark *quark, Quark_Task_Flags *task_flags,
                           int m, int n, 
                           float *A,    int szeA,
                           float *fake1, int szeF1, int paramF1,
                           float *fake2, int szeF2, int paramF2)
{
    DAG_CORE_GETRIP;
    QUARK_Insert_Task(
        quark, CORE_sgetrip_f2_quark, task_flags,
        sizeof(int),                     &m,   VALUE,
        sizeof(int),                     &n,   VALUE,
        sizeof(float)*szeA, A,        INOUT,
        sizeof(float)*szeA, NULL,     SCRATCH,
        sizeof(float)*szeF1, fake1,     paramF1,
        sizeof(float)*szeF2, fake2,     paramF2,
        0);
}

/***************************************************************************/
#if defined(PLASMA_HAVE_WEAK)
#pragma weak CORE_sgetrip_f2_quark = PCORE_sgetrip_f2_quark
#define CORE_sgetrip_f2_quark PCORE_sgetrip_f2_quark
#endif
void CORE_sgetrip_f2_quark(Quark *quark)
{
    int m;
    int n;
    float *A;
    float *W;
    float *fake1;
    float *fake2;

    quark_unpack_args_6(quark, m, n, A, W, fake1, fake2);
    CORE_sgetrip(m, n, A, W);
}