MatrixScalarOperations.cpph

Go to the documentation of this file.

template<int Rows, int Cols, typename Scalar>
tarch::la::Matrix<Rows,Cols,Scalar>& tarch::la::operator*=(
  Matrix<Rows,Cols,Scalar>&  matrix,
  const Scalar&              scalar
) {
  #ifdef CompilerICC
  #pragma ivdep
  #endif
  for (int i=0; i < Rows; i++) {
    for (int j=0; j < Cols; j++) {
      matrix(i,j) *= scalar;
    }
  }
  return matrix;
}
 
 
 
template<int Rows, int Cols, typename Scalar>
tarch::la::Matrix<Rows,Cols,Scalar> tarch::la::operator* (
  const tarch::la::Matrix<Rows,Cols,Scalar>&  matrix,
  const Scalar&                               scalar
) {
  Matrix<Rows,Cols,Scalar> result(matrix);
  #ifdef CompilerICC
  #pragma ivdep
  #endif
  for (int i=0; i < Rows; i++) {
    for (int j=0; j < Cols; j++) {
      result(i,j) = matrix(i,j)*scalar;
    }
  }
  return result;
}
 
 
template<int Rows, int Cols, typename Scalar>
tarch::la::Matrix<Rows,Cols,Scalar> tarch::la::operator* (
  const Scalar&                               scalar,
  const tarch::la::Matrix<Rows,Cols,Scalar>&  matrix
) {
  return matrix * scalar;
}
 
 
template<int Rows, int Cols>
tarch::la::Matrix<Rows,Cols,std::complex<double> > tarch::la::operator* (
  const double&                                              scalar,
  const tarch::la::Matrix<Rows,Cols,std::complex<double> >&  matrix
) {
  return std::complex<double>( scalar ) * matrix;
}