/usr/include/atlas/atlas_ssyr2.h is in libatlas-dev 3.10.2-7.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 | /*
* This file generated on line 445 of /tmp/buildd/atlas-3.10.2/build/atlas-base/../..//tune/blas/ger/r2hgen.c
*/
#ifndef ATLAS_SSYR2_H
#define ATLAS_SSYR2_H
#include "atlas_sr2.h"
#define ATL_s2U_NU 1
#define ATL_s2L_NU 1
#define ATL_GetPartS2 ATL_GetPartR2
#define ATL_MIN_RESTRICTED_M 2
#define ATL_URGERK ATL_sger2k__900001
static void ATL_GENGERK(ATL_CINT M, ATL_CINT N, const TYPE *X,
const TYPE *Y, const TYPE *W, const TYPE *Z,
TYPE *A, ATL_CINT lda)
{
int nu, minM, minN, i, FNU, aX, aX2A, aY, aW, aZ;
ATL_INT CEL;
ATL_r2kern_t gerk;
gerk = ATL_GetR2Kern(M, N, A, lda, &i, &nu, &minM, &minN,
&aX, &aX2A, &aY, &FNU, &CEL);
if (aX2A)
{
aX = ((size_t)A) % ATL_Cachelen == ((size_t)X) % ATL_Cachelen;
aW = ((size_t)A) % ATL_Cachelen == ((size_t)W) % ATL_Cachelen;
} else
{
aW = (aX) ? (((size_t)W)/aX)*aX == (size_t)W : 1;
aX = (aX) ? (((size_t)X)/aX)*aX == (size_t)X : 1;
}
aZ = (aY) ? (((size_t)Z)/aY)*aY == (size_t)Z : 1;
aY = (aY) ? (((size_t)Y)/aY)*aY == (size_t)Y : 1;
if (M >= minM && N >= minN && aX && aY && aW && aZ)
{
if (FNU)
{
ATL_CINT n = (N/nu)*nu, nr=N-n;
gerk(M, n, X, Y, W, Z, A, lda);
if (nr)
ATL_sger2k_Nlt8(M, nr, ATL_rone, X, 1, Y+n, 1, ATL_rone, W, 1, Z+n, 1, A+n*lda, lda);
} /* end if (FNU) */
else
gerk(M, N, X, Y, W, Z, A, lda);
} /* end if can call optimized kernel */
else
ATL_sger2k_Mlt16(M, N, ATL_rone, X, 1, Y, 1, ATL_rone, W, 1, Z, 1, A, lda);
}
#define ATL_SYR2U_nu(A_, lda_, x_, y_) \
{ \
TYPE *aa=(A_); \
ATL_CINT lda0_ = 0; \
const TYPE x0_=*(x_); \
const TYPE y0_=*(y_); \
aa[lda0_+0] += x0_*y0_ + y0_*x0_; \
}
#define ATL_SYR2L_nu(A_, lda_, x_, y_) \
{ \
TYPE *aa=(A_); \
ATL_CINT lda0_ = 0; \
const TYPE x0_=*(x_); \
const TYPE y0_=*(y_); \
aa[lda0_+0] += x0_*y0_ + y0_*x0_; \
}
#endif
|