diff options
author | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2013-11-07 18:32:24 +0100 |
---|---|---|
committer | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2013-11-07 18:32:24 +0100 |
commit | ae83f5ede96cb1171fb5fd9303777a9035748ca9 (patch) | |
tree | 727b5ba43aee85da5f30728fefbc35642337be48 /unsupported/test/autodiff_scalar.cpp | |
parent | 76c230a84d4857722e4a72e9007302a504af0fb7 (diff) |
Fixed bug #702 and added unit test.
Thanks to Alexander Werner for the report.
Diffstat (limited to 'unsupported/test/autodiff_scalar.cpp')
-rw-r--r-- | unsupported/test/autodiff_scalar.cpp | 44 |
1 files changed, 44 insertions, 0 deletions
diff --git a/unsupported/test/autodiff_scalar.cpp b/unsupported/test/autodiff_scalar.cpp new file mode 100644 index 000000000..ba4b5aec4 --- /dev/null +++ b/unsupported/test/autodiff_scalar.cpp @@ -0,0 +1,44 @@ +// This file is part of Eigen, a lightweight C++ template library +// for linear algebra. +// +// Copyright (C) 2013 Christoph Hertzberg <chtz@informatik.uni-bremen.de> +// +// This Source Code Form is subject to the terms of the Mozilla +// Public License v. 2.0. If a copy of the MPL was not distributed +// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. + +#include "main.h" +#include <unsupported/Eigen/AutoDiff> + +/* + * In this file scalar derivations are tested for correctness. + * TODO add more tests! + */ + +template<typename Scalar> void check_atan2() +{ + typedef Matrix<Scalar, 1, 1> Deriv1; + typedef AutoDiffScalar<Deriv1> AD; + + AD x(internal::random<Scalar>(-3.0, 3.0), Deriv1::UnitX()); + + using std::exp; + Scalar r = exp(internal::random<Scalar>(-10, 10)); + + AD s = sin(x), c = cos(x); + AD res = atan2(r*s, r*c); + + VERIFY_IS_APPROX(res.value(), x.value()); + VERIFY_IS_APPROX(res.derivatives(), x.derivatives()); +} + + + + +void test_autodiff_scalar() +{ + for(int i = 0; i < g_repeat; i++) { + CALL_SUBTEST_1( check_atan2<float>() ); + CALL_SUBTEST_2( check_atan2<double>() ); + } +} |