VectorSlice.cpph

Go to the documentation of this file.

template<int SizeLhs, int SizeRhs, typename Scalar>
tarch::la::Vector<SizeLhs,Scalar> tarch::la::slice(const tarch::la::Vector<SizeRhs,Scalar>& vector, int fromIndex, int stride) {
  tarch::la::Vector<SizeLhs,Scalar> result;
 
  assertion5(fromIndex>=0,fromIndex, SizeLhs, SizeRhs, stride, ::toString(vector));
  assertion5(fromIndex+(SizeLhs-1)*stride<SizeRhs,fromIndex, SizeLhs, SizeRhs, stride, ::toString(vector));
 
  #ifdef CompilerICC
  #pragma ivdep
  #endif
  for (int i=0; i<SizeLhs; i++) {
    result(i) = vector(fromIndex+i*stride);
  }
 
  return result;
}
 
 
template<int SizeLhs, int SizeRhs, typename Scalar>
void tarch::la::slice(
  tarch::la::Vector<SizeLhs,Scalar>&        toVector,
  const tarch::la::Vector<SizeRhs,Scalar>&  fromVector,
  int                                       fromIndexInToVector,
  int                                       strideInToVector
) {
  assertion6(fromIndexInToVector>=0, SizeLhs, SizeRhs, ::toString(toVector), ::toString(fromVector), fromIndexInToVector, strideInToVector);
  assertion6(fromIndexInToVector+(SizeRhs-1)*strideInToVector<SizeLhs, SizeLhs, SizeRhs, ::toString(toVector), ::toString(fromVector), fromIndexInToVector, strideInToVector);
 
  #ifdef CompilerICC
  #pragma ivdep
  #endif
  for (int i=0; i<SizeRhs; i++) {
    toVector(fromIndexInToVector+i*strideInToVector) = fromVector(i);
  }
}
 
 
template<int SizeLhs, int SizeRhs, typename Scalar>
tarch::la::Vector<SizeLhs,Scalar> tarch::la::expandOrSlice(const Vector<SizeRhs,Scalar>& vector) {
  return slice<SizeLhs>(vector,0);
}
 
 
template<int SizeLhs, typename Scalar>
tarch::la::Vector<SizeLhs,Scalar> tarch::la::expandOrSlice(const Scalar&  scalar) {
  return Vector<SizeLhs,Scalar>(scalar);
}