Skip to content

Commit 977f96e

Browse files
committed
Allow ad through cholesky_corr_free
1 parent df9d406 commit 977f96e

1 file changed

Lines changed: 2 additions & 1 deletion

File tree

stan/math/prim/constraint/cholesky_corr_free.hpp

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -14,6 +14,7 @@ template <typename T, require_eigen_t<T>* = nullptr>
1414
auto cholesky_corr_free(const T& x) {
1515
using Eigen::Dynamic;
1616
using Eigen::Matrix;
17+
using std::sqrt;
1718

1819
check_square("cholesky_corr_free", "x", x);
1920
// should validate lower-triangular, unit lengths
@@ -26,7 +27,7 @@ auto cholesky_corr_free(const T& x) {
2627
z.coeffRef(k++) = corr_free(x_ref.coeff(i, 0));
2728
auto sum_sqs = square(x_ref.coeff(i, 0));
2829
for (int j = 1; j < i; ++j) {
29-
z.coeffRef(k++) = corr_free(x_ref.coeff(i, j) / std::sqrt(1.0 - sum_sqs));
30+
z.coeffRef(k++) = corr_free(x_ref.coeff(i, j) / sqrt(1.0 - sum_sqs));
3031
sum_sqs += square(x_ref.coeff(i, j));
3132
}
3233
}

0 commit comments

Comments
 (0)