void mitk::AnisotropicRegistrationCommon::PropagateMatrices( const MatrixList &src,
                                                             MatrixList &dst,
                                                             const Rotation &rotation )
{

 const vnl_matrix_fixed < double, 3, 3 > rotationT = rotation.GetTranspose();

#pragma omp parallel for
  for ( size_t i = 0; i < src.size(); ++i )
  {
    dst[i] = rotation * src[i] * rotationT;
    }
}