Skip to content

Commit

Permalink
[Jenkins] auto-formatting by clang-format version 10.0.0-4ubuntu1
Browse files Browse the repository at this point in the history
  • Loading branch information
stan-buildbot committed Aug 16, 2023
1 parent 24c8f0f commit 4a71d4f
Show file tree
Hide file tree
Showing 7 changed files with 12 additions and 15 deletions.
2 changes: 1 addition & 1 deletion stan/math/prim/fun/complex_schur_decompose.hpp
Original file line number Diff line number Diff line change
Expand Up @@ -82,7 +82,7 @@ inline std::tuple<Eigen::Matrix<complex_return_t<scalar_type_t<M>>, -1, -1>,
Eigen::Matrix<complex_return_t<scalar_type_t<M>>, -1, -1>>
complex_schur_decompose(const M& m) {
if (m.size() == 0)
return std::make_tuple(m,m);
return std::make_tuple(m, m);
check_square("complex_schur_decompose", "m", m);
using MatType = Eigen::Matrix<scalar_type_t<M>, -1, -1>;
// copy because ComplexSchur requires Eigen::Matrix type
Expand Down
13 changes: 7 additions & 6 deletions stan/math/rev/fun/eigendecompose_sym.hpp
Original file line number Diff line number Diff line change
Expand Up @@ -39,14 +39,15 @@ inline auto eigendecompose_sym(const T& m) {

reverse_pass_callback([eigenvals, arena_m, eigenvecs]() mutable {
// eigenvalue reverse calculation
arena_m.adj()
+= eigenvecs.val_op() * eigenvals.adj().asDiagonal() * eigenvecs.val_op().transpose();
arena_m.adj() += eigenvecs.val_op() * eigenvals.adj().asDiagonal()
* eigenvecs.val_op().transpose();
// eigenvector reverse calculation
const int p = arena_m.val().cols();
Eigen::MatrixXd f = (1
/ (eigenvals.val_op().rowwise().replicate(p).transpose()
- eigenvals.val_op().rowwise().replicate(p))
.array());
Eigen::MatrixXd f
= (1
/ (eigenvals.val_op().rowwise().replicate(p).transpose()
- eigenvals.val_op().rowwise().replicate(p))
.array());
f.diagonal().setZero();
arena_m.adj()
+= eigenvecs.val_op()
Expand Down
5 changes: 2 additions & 3 deletions test/unit/math/mix/fun/complex_schur_decompose_test.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -51,9 +51,9 @@ TEST(mathMixFun, complexSchurDecompose) {

template <typename V>
void test_complex_schur_decompose(const Eigen::MatrixXd& x) {
using stan::math::complex_schur_decompose;
using stan::math::complex_schur_decompose_t;
using stan::math::complex_schur_decompose_u;
using stan::math::complex_schur_decompose;
using stan::math::get_real;
using stan::math::value_of_rec;
Eigen::Matrix<V, -1, -1> X = x;
Expand All @@ -67,14 +67,13 @@ void test_complex_schur_decompose(const Eigen::MatrixXd& x) {
std::tie(U, T) = complex_schur_decompose(X);
auto X3 = U * T * U.adjoint();
EXPECT_MATRIX_NEAR(x, value_of_rec(get_real(X3)), 1e-8);

}

template <typename V>
void test_complex_schur_decompose_complex(const Eigen::MatrixXd& x) {
using stan::math::complex_schur_decompose;
using stan::math::complex_schur_decompose_t;
using stan::math::complex_schur_decompose_u;
using stan::math::complex_schur_decompose;
using stan::math::value_of_rec;
Eigen::Matrix<std::complex<V>, -1, -1> X(x.rows(), x.cols());
for (int i = 0; i < x.size(); ++i)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -27,7 +27,7 @@ void expectComplexEigenvectorsId() {

std::tie(eigenvectors, eigenvalues) = stan::math::eigendecompose(c22);
auto I2 = (eigenvectors.inverse() * c22 * eigenvectors
* eigenvalues.asDiagonal().inverse());
* eigenvalues.asDiagonal().inverse());

expect_identity_matrix_complex(I2);
}
Expand Down
2 changes: 1 addition & 1 deletion test/unit/math/mix/fun/eigendecompose_identity_test.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -19,7 +19,7 @@ void expectEigenvectorsId() {
auto I = (vecs.inverse() * m * vecs * vals.asDiagonal().inverse()).real();
expect_identity_matrix(I);

std::tie(vecs,vals) = eigendecompose(m);
std::tie(vecs, vals) = eigendecompose(m);
auto I2 = (vecs.inverse() * m * vecs * vals.asDiagonal().inverse()).real();
expect_identity_matrix(I2);
}
Expand Down
2 changes: 0 additions & 2 deletions test/unit/math/prim/fun/complex_schur_decompose_test.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -26,8 +26,6 @@ TEST(primFun, complex_schur_decompose_ut) {
EXPECT_MATRIX_COMPLEX_NEAR(B, B_recovered, 1e-8);
}



TEST(primFun, complex_schur_decompose) {
using c_t = std::complex<double>;
using stan::math::complex_schur_decompose;
Expand Down
1 change: 0 additions & 1 deletion test/unit/math/prim/fun/qr_thin_test.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -34,5 +34,4 @@ TEST(MathMatrixPrimMat, qr_thin) {
EXPECT_NEAR(m1(i, j), m3(j, i), 1e-12);
}
}

}

0 comments on commit 4a71d4f

Please sign in to comment.