1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This file is part of Eigen, a lightweight C++ template library 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for linear algebra. 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2006-2008 Benoit Jacob <jacob.benoit.1@gmail.com> 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2009 Ricard Marxer <email@ricardmarxer.com> 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This Source Code Form is subject to the terms of the Mozilla 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Public License v. 2.0. If a copy of the MPL was not distributed 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "main.h" 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <iostream> 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathusing namespace std; 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename MatrixType> void reverse(const MatrixType& m) 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef typename MatrixType::Index Index; 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef typename MatrixType::Scalar Scalar; 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, 1> VectorType; 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index rows = m.rows(); 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index cols = m.cols(); 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // this test relies a lot on Random.h, and there's not much more that we can do 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // to test it, hence I consider that we will have tested Random.h 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixType m1 = MatrixType::Random(rows, cols); 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VectorType v1 = VectorType::Random(rows); 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixType m1_r = m1.reverse(); 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Verify that MatrixBase::reverse() works 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int i = 0; i < rows; i++ ) { 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int j = 0; j < cols; j++ ) { 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(m1_r(i, j), m1(rows - 1 - i, cols - 1 - j)); 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Reverse<MatrixType> m1_rd(m1); 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Verify that a Reverse default (in both directions) of an expression works 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int i = 0; i < rows; i++ ) { 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int j = 0; j < cols; j++ ) { 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(m1_rd(i, j), m1(rows - 1 - i, cols - 1 - j)); 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Reverse<MatrixType, BothDirections> m1_rb(m1); 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Verify that a Reverse in both directions of an expression works 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int i = 0; i < rows; i++ ) { 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int j = 0; j < cols; j++ ) { 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(m1_rb(i, j), m1(rows - 1 - i, cols - 1 - j)); 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Reverse<MatrixType, Vertical> m1_rv(m1); 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Verify that a Reverse in the vertical directions of an expression works 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int i = 0; i < rows; i++ ) { 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int j = 0; j < cols; j++ ) { 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(m1_rv(i, j), m1(rows - 1 - i, j)); 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Reverse<MatrixType, Horizontal> m1_rh(m1); 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Verify that a Reverse in the horizontal directions of an expression works 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int i = 0; i < rows; i++ ) { 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int j = 0; j < cols; j++ ) { 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(m1_rh(i, j), m1(i, cols - 1 - j)); 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VectorType v1_r = v1.reverse(); 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Verify that a VectorType::reverse() of an expression works 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int i = 0; i < rows; i++ ) { 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v1_r(i), v1(rows - 1 - i)); 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixType m1_cr = m1.colwise().reverse(); 77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Verify that PartialRedux::reverse() works (for colwise()) 78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int i = 0; i < rows; i++ ) { 79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int j = 0; j < cols; j++ ) { 80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(m1_cr(i, j), m1(rows - 1 - i, j)); 81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixType m1_rr = m1.rowwise().reverse(); 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Verify that PartialRedux::reverse() works (for rowwise()) 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int i = 0; i < rows; i++ ) { 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for ( int j = 0; j < cols; j++ ) { 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(m1_rr(i, j), m1(i, cols - 1 - j)); 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Scalar x = internal::random<Scalar>(); 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index r = internal::random<Index>(0, rows-1), 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath c = internal::random<Index>(0, cols-1); 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m1.reverse()(r, c) = x; 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(x, m1(rows - 1 - r, cols - 1 - c)); 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /* 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m1.colwise().reverse()(r, c) = x; 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(x, m1(rows - 1 - r, c)); 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m1.rowwise().reverse()(r, c) = x; 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(x, m1(r, cols - 1 - c)); 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid test_array_reverse() 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < g_repeat; i++) { 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_1( reverse(Matrix<float, 1, 1>()) ); 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_2( reverse(Matrix2f()) ); 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_3( reverse(Matrix4f()) ); 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_4( reverse(Matrix4d()) ); 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_5( reverse(MatrixXcf(3, 3)) ); 117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_6( reverse(MatrixXi(6, 3)) ); 118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_7( reverse(MatrixXcd(20, 20)) ); 119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_8( reverse(Matrix<float, 100, 100>()) ); 120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_9( reverse(Matrix<float,Dynamic,Dynamic,RowMajor>(6,3)) ); 121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 122c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifdef EIGEN_TEST_PART_3 123c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Vector4f x; x << 1, 2, 3, 4; 124c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Vector4f y; y << 4, 3, 2, 1; 125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY(x.reverse()[1] == 3); 126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY(x.reverse() == y); 127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif 128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 129