interface/zgeadd.c

   1 /*********************************************************************/
   2 /* Copyright 2009, 2010 The University of Texas at Austin.           */
   3 /* All rights reserved.                                              */
   4 /*                                                                   */
   5 /* Redistribution and use in source and binary forms, with or        */
   6 /* without modification, are permitted provided that the following   */
   7 /* conditions are met:                                               */
   8 /*                                                                   */
   9 /*   1. Redistributions of source code must retain the above         */
  10 /*      copyright notice, this list of conditions and the following  */
  11 /*      disclaimer.                                                  */
  12 /*                                                                   */
  13 /*   2. Redistributions in binary form must reproduce the above      */
  14 /*      copyright notice, this list of conditions and the following  */
  15 /*      disclaimer in the documentation and/or other materials       */
  16 /*      provided with the distribution.                              */
  17 /*                                                                   */
  18 /*    THIS  SOFTWARE IS PROVIDED  BY THE  UNIVERSITY OF  TEXAS AT    */
  19 /*    AUSTIN  ``AS IS''  AND ANY  EXPRESS OR  IMPLIED WARRANTIES,    */
  20 /*    INCLUDING, BUT  NOT LIMITED  TO, THE IMPLIED  WARRANTIES OF    */
  21 /*    MERCHANTABILITY  AND FITNESS FOR  A PARTICULAR  PURPOSE ARE    */
  22 /*    DISCLAIMED.  IN  NO EVENT SHALL THE UNIVERSITY  OF TEXAS AT    */
  23 /*    AUSTIN OR CONTRIBUTORS BE  LIABLE FOR ANY DIRECT, INDIRECT,    */
  24 /*    INCIDENTAL,  SPECIAL, EXEMPLARY,  OR  CONSEQUENTIAL DAMAGES    */
  25 /*    (INCLUDING, BUT  NOT LIMITED TO,  PROCUREMENT OF SUBSTITUTE    */
  26 /*    GOODS  OR  SERVICES; LOSS  OF  USE,  DATA,  OR PROFITS;  OR    */
  27 /*    BUSINESS INTERRUPTION) HOWEVER CAUSED  AND ON ANY THEORY OF    */
  28 /*    LIABILITY, WHETHER  IN CONTRACT, STRICT  LIABILITY, OR TORT    */
  29 /*    (INCLUDING NEGLIGENCE OR OTHERWISE)  ARISING IN ANY WAY OUT    */
  30 /*    OF  THE  USE OF  THIS  SOFTWARE,  EVEN  IF ADVISED  OF  THE    */
  31 /*    POSSIBILITY OF SUCH DAMAGE.                                    */
  32 /*                                                                   */
  33 /* The views and conclusions contained in the software and           */
  34 /* documentation are those of the authors and should not be          */
  35 /* interpreted as representing official policies, either expressed   */
  36 /* or implied, of The University of Texas at Austin.                 */
  37 /*********************************************************************/
  38
  39 #include <stdio.h>
  40 #include "common.h"
  41 #ifdef FUNCTION_PROFILE
  42 #include "functable.h"
  43 #endif
  44
  45 #if defined(DOUBLE)
  46 #define ERROR_NAME "ZGEADD "
  47 #else
  48 #define ERROR_NAME "CGEADD "
  49 #endif
  50
  51 #ifndef CBLAS
  52
  53 void NAME(blasint *M, blasint *N, FLOAT *ALPHA, FLOAT *a, blasint *LDA,
  54                                   FLOAT *BETA,  FLOAT *c, blasint *LDC)
  55 {
  56
  57   blasint m = *M;
  58   blasint n = *N;
  59   blasint lda = *LDA;
  60   blasint ldc = *LDC;
  61
  62   blasint info;
  63
  64   PRINT_DEBUG_NAME;
  65
  66   info = 0;
  67
  68
  69   if (lda < MAX(1, m))  info = 6;
  70   if (ldc < MAX(1, m))  info = 8;
  71
  72   if (n < 0)            info = 2;
  73   if (m < 0)            info = 1;
  74
  75   if (info != 0){
  76     BLASFUNC(xerbla)(ERROR_NAME, &info, sizeof(ERROR_NAME));
  77     return;
  78   }
  79
  80 #else
  81 void CNAME( enum CBLAS_ORDER order, blasint m,  blasint n,  FLOAT *ALPHA, FLOAT *a,  blasint lda, FLOAT *BETA,
  82                   FLOAT *c,  blasint ldc)
  83 {
  84 /*
  85 void CNAME(enum CBLAS_ORDER order,
  86            blasint m, blasint n,
  87            FLOAT alpha,
  88            FLOAT  *a, blasint lda,
  89            FLOAT beta,
  90            FLOAT  *c, blasint ldc){ */
  91
  92   blasint info, t;
  93
  94   PRINT_DEBUG_CNAME;
  95
  96   info  =  0;
  97
  98   if (order == CblasColMajor) {
  99
 100     info = -1;
 101
 102     if (ldc < MAX(1, m))  info = 8;
 103     if (lda < MAX(1, m))  info = 5;
 104     if (n < 0)            info = 2;
 105     if (m < 0)            info = 1;
 106
 107   }
 108
 109   if (order == CblasRowMajor) {
 110     info = -1;
 111
 112     t = n;
 113     n = m;
 114     m = t;
 115
 116     if (ldc < MAX(1, m))  info = 8;
 117     if (lda < MAX(1, m))  info = 5;
 118     if (n < 0)            info = 2;
 119     if (m < 0)            info = 1;
 120   }
 121
 122   if (info >= 0) {
 123     BLASFUNC(xerbla)(ERROR_NAME, &info, sizeof(ERROR_NAME));
 124     return;
 125   }
 126
 127 #endif
 128
 129   if ((m==0) || (n==0)) return;
 130
 131
 132   IDEBUG_START;
 133
 134   FUNCTION_PROFILE_START();
 135
 136
 137   GEADD_K(m,n,ALPHA[0],ALPHA[1], a, lda, BETA[0], BETA[1], c, ldc);
 138
 139
 140   FUNCTION_PROFILE_END(1, 2* m * n ,  2 * m * n);
 141
 142   IDEBUG_END;
 143
 144   return;
 145
 146 }