phynexis-deps/armadillo-code
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

clear output arrays before use

Browse Source
This commit is contained in:
conrad
2024-02-26 23:11:39 +10:00
parent 54c62da09b
commit ca8da8bdab
1 changed files with 13 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
include/armadillo_bits/translate_blas.hpp
+13 -3
View File
@@ -145,6 +145,9 @@ namespace blas
eT result[2]; // paranoia: using two elements instead of one
result[0] = eT(0);
result[1] = eT(0);
blas::gemv(&trans, &m, &n, &alpha, x, &m, y, &inc, &beta, &result[0], &inc);
return result[0];
@@ -186,11 +189,14 @@ namespace blas
eT result[2]; // paranoia: using two elements instead of one
result[0] = eT(0);
result[1] = eT(0);
blas::gemv(&trans, &m, &n, &alpha, x, &m, y, &inc, &beta, &result[0], &inc);
return result[0];
}
return eT(0);
}
@@ -205,11 +211,13 @@ namespace blas
if(is_float<eT>::value)
{
// WARNING: sasum() from Accelerate framework (macOS) may return 'double' instead of 'float'
blas_int n = blas_int(n_elem);
blas_int inc = 1;
typedef float T;
return arma_fortran(arma_sasum)(&n, (const T*)x, &inc);
return eT( arma_fortran(arma_sasum)(&n, (const T*)x, &inc) );
}
else
if(is_double<eT>::value)
@@ -218,7 +226,7 @@ namespace blas
blas_int inc = 1;
typedef double T;
return arma_fortran(arma_dasum)(&n, (const T*)x, &inc);
return eT( arma_fortran(arma_dasum)(&n, (const T*)x, &inc) );
}
return eT(0);
@@ -235,6 +243,8 @@ namespace blas
if(is_float<eT>::value)
{
// WARNING: snrm2() from Accelerate framework (macOS) may return 'double' instead of 'float'
blas_int n = blas_int(n_elem);
blas_int inc = 1;