determinant.cpp
Go to the documentation of this file.
1 // This file is part of Eigen, a lightweight C++ template library
2 // for linear algebra.
3 //
4 // Copyright (C) 2008 Benoit Jacob <jacob.benoit.1@gmail.com>
5 // Copyright (C) 2008 Gael Guennebaud <gael.guennebaud@inria.fr>
6 //
7 // This Source Code Form is subject to the terms of the Mozilla
8 // Public License v. 2.0. If a copy of the MPL was not distributed
9 // with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
10 
11 #include "main.h"
12 #include <Eigen/LU>
13 
14 template<typename MatrixType> void determinant(const MatrixType& m)
15 {
16  /* this test covers the following files:
17  Determinant.h
18  */
19  Index size = m.rows();
20 
21  MatrixType m1(size, size), m2(size, size);
22  m1.setRandom();
23  m2.setRandom();
24  typedef typename MatrixType::Scalar Scalar;
25  Scalar x = internal::random<Scalar>();
26  VERIFY_IS_APPROX(MatrixType::Identity(size, size).determinant(), Scalar(1));
27  VERIFY_IS_APPROX((m1*m2).eval().determinant(), m1.determinant() * m2.determinant());
28  if(size==1) return;
29  Index i = internal::random<Index>(0, size-1);
30  Index j;
31  do {
32  j = internal::random<Index>(0, size-1);
33  } while(j==i);
34  m2 = m1;
35  m2.row(i).swap(m2.row(j));
36  VERIFY_IS_APPROX(m2.determinant(), -m1.determinant());
37  m2 = m1;
38  m2.col(i).swap(m2.col(j));
39  VERIFY_IS_APPROX(m2.determinant(), -m1.determinant());
40  VERIFY_IS_APPROX(m2.determinant(), m2.transpose().determinant());
41  VERIFY_IS_APPROX(numext::conj(m2.determinant()), m2.adjoint().determinant());
42  m2 = m1;
43  m2.row(i) += x*m2.row(j);
44  VERIFY_IS_APPROX(m2.determinant(), m1.determinant());
45  m2 = m1;
46  m2.row(i) *= x;
47  VERIFY_IS_APPROX(m2.determinant(), m1.determinant() * x);
48 
49  // check empty matrix
50  VERIFY_IS_APPROX(m2.block(0,0,0,0).determinant(), Scalar(1));
51 }
52 
54 {
55  for(int i = 0; i < g_repeat; i++) {
56  int s = 0;
61  CALL_SUBTEST_5( determinant(Matrix<std::complex<double>, 10, 10>()) );
62  s = internal::random<int>(1,EIGEN_TEST_MAX_SIZE/4);
63  CALL_SUBTEST_6( determinant(MatrixXd(s, s)) );
65  }
66 }
Matrix3f m
SCALAR Scalar
Definition: bench_gemm.cpp:46
#define CALL_SUBTEST_6(FUNC)
#define CALL_SUBTEST_4(FUNC)
#define CALL_SUBTEST_3(FUNC)
void determinant(const MatrixType &m)
Definition: determinant.cpp:14
MatrixType m2(n_dims)
MatrixXf MatrixType
AnnoyingScalar conj(const AnnoyingScalar &x)
Scalar Scalar int size
Definition: benchVecAdd.cpp:17
#define VERIFY_IS_APPROX(a, b)
#define CALL_SUBTEST_1(FUNC)
Matrix3d m1
Definition: IOFormat.cpp:2
static int g_repeat
Definition: main.h:169
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
The Index type as used for the API.
Definition: Meta.h:74
RealScalar s
EIGEN_DECLARE_TEST(determinant)
Definition: determinant.cpp:53
#define TEST_SET_BUT_UNUSED_VARIABLE(X)
Definition: main.h:121
#define CALL_SUBTEST_5(FUNC)
#define EIGEN_TEST_MAX_SIZE
#define CALL_SUBTEST_2(FUNC)
internal::nested_eval< T, 1 >::type eval(const T &xpr)
The matrix class, also used for vectors and row-vectors.
set noclip points set clip one set noclip two set bar set border lt lw set xdata set ydata set zdata set x2data set y2data set boxwidth set dummy x
std::ptrdiff_t j


gtsam
Author(s):
autogenerated on Tue Jul 4 2023 02:34:10