diff --git a/MathLib/LinAlg/Sparse/CRSMatrixOpenMP.h b/MathLib/LinAlg/Sparse/CRSMatrixOpenMP.h index 6ef94b58b5883bf2e998b1f157fa404232660f9f..70fa7e399d69e27051794dc2d87d57dd421980d8 100644 --- a/MathLib/LinAlg/Sparse/CRSMatrixOpenMP.h +++ b/MathLib/LinAlg/Sparse/CRSMatrixOpenMP.h @@ -39,8 +39,8 @@ public: virtual ~CRSMatrixOpenMP() {} - - virtual void amux(FP_TYPE const d, FP_TYPE const * const x, FP_TYPE *y) const + virtual void amux(FP_TYPE const d, FP_TYPE const* const __restrict__ x, + FP_TYPE* __restrict__ y) const { amuxCRSParallelOpenMP(d, this->_n_rows, this->_row_ptr, this->_col_idx, this->_data, x, y); } diff --git a/MathLib/LinAlg/Sparse/SparseMatrixBase.h b/MathLib/LinAlg/Sparse/SparseMatrixBase.h index 5ac0ded8974099f1fcbb25ddc3f17f963d2d92a7..579d1ad208248d00601e52e79e86fce925085c59 100644 --- a/MathLib/LinAlg/Sparse/SparseMatrixBase.h +++ b/MathLib/LinAlg/Sparse/SparseMatrixBase.h @@ -18,7 +18,8 @@ public: * @param x vector to multiply with * @param y result vector */ - virtual void amux(FP_TYPE d, FP_TYPE const * const __restrict__ x, FP_TYPE * __restrict__ y) const = 0; + virtual void amux(FP_TYPE const d, FP_TYPE const* const __restrict__ x, + FP_TYPE* __restrict__ y) const = 0; virtual ~SparseMatrixBase() {} /** * get the number of rows