core_dgeadd.c

Go to the documentation of this file.

#include "common.h"

/***************************************************************************/
#if defined(PLASMA_HAVE_WEAK)
#pragma weak CORE_dgeadd = PCORE_dgeadd
#define CORE_dgeadd PCORE_dgeadd
#endif
void CORE_dgeadd(int M, int N,  double alpha,
                double *A, int LDA,
                double *B, int LDB)
{
    int j;

    if (M == LDA && M == LDB)
        cblas_daxpy(M*N, (alpha), A, 1, B, 1);
    else {
        for (j = 0; j < N; j++)
            cblas_daxpy(M, (alpha), &A[j*LDA], 1, &B[j*LDB], 1);
    }
}

/***************************************************************************/
void QUARK_CORE_dgeadd(Quark *quark, Quark_Task_Flags *task_flags,
                      int m, int n, int nb, double alpha,
                      double *A, int lda,
                      double *B, int ldb)
{
    DAG_CORE_GEADD;
    QUARK_Insert_Task(quark, CORE_dgeadd_quark, task_flags,
        sizeof(int),                        &m,     VALUE,
        sizeof(int),                        &n,     VALUE,
        sizeof(double),         &alpha, VALUE,
        sizeof(double)*nb*nb,    A,             INPUT,
        sizeof(int),                        &lda,   VALUE,
        sizeof(double)*nb*nb,    B,             INOUT,
        sizeof(int),                        &ldb,   VALUE,
        0);
}

/***************************************************************************/
#if defined(PLASMA_HAVE_WEAK)
#pragma weak CORE_dgeadd_quark = PCORE_dgeadd_quark
#define CORE_dgeadd_quark PCORE_dgeadd_quark
#endif
void CORE_dgeadd_quark(Quark *quark)
{
    int M;
    int N;
    double alpha;
    double *A;
    int LDA;
    double *B;
    int LDB;

    int j;

    quark_unpack_args_7(quark, M, N, alpha, A, LDA, B, LDB);
    if (M == LDA && M == LDB)
        cblas_daxpy(M*N, (alpha), A, 1, B, 1);
    else {
        for (j = 0; j < N; j++)
            cblas_daxpy(M, (alpha), &A[j*LDA], 1, &B[j*LDB], 1);
    }
}