Merge commit '4cf25de56d29ef39c93eb2595d13dcfd65f97818' into HEAD

Author: yashikno

.github/ISSUE_TEMPLATE.md

@@ -25,4 +25,4 @@ If this is a **feature request**, show what you expect to happen if the feature
 
 
 #### Current Version:
-v4.6.2
+v4.7.0

.github/workflows/header_checks.yml

@@ -26,7 +26,7 @@ jobs:
     runs-on: windows-latest
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
     - uses: actions/setup-python@v4
       with:
         python-version: '3.x'
@@ -52,7 +52,7 @@ jobs:
     runs-on: ubuntu-latest
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
 
     - name: Run header tests
       run: |
@@ -63,7 +63,7 @@ jobs:
     runs-on: ubuntu-latest
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
 
     - name: Run header tests
       run: |

.github/workflows/main.yml

@@ -24,7 +24,7 @@ jobs:
     runs-on: windows-latest
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
     - uses: actions/setup-python@v4
       with:
         python-version: '3.x'
@@ -74,7 +74,7 @@ jobs:
     runs-on: windows-latest
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
     - uses: actions/setup-python@v4
       with:
         python-version: '3.x'
@@ -128,7 +128,7 @@ jobs:
     runs-on: windows-latest
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
     - uses: actions/setup-python@v4
       with:
         python-version: '3.x'
@@ -177,7 +177,7 @@ jobs:
     runs-on: windows-latest
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
     - uses: actions/setup-python@v4
       with:
         python-version: '3.x'

RELEASE-NOTES.txt

@@ -1,5 +1,24 @@
 Stan Math Library Release Notes
 
+======================================================================
+v4.7.0 (5 September 2023)
+======================================================================
+
+- Improved accuracy of `matrix_exp_multiply`. (#2619)
+- Added framework for vectorising ternary functions (#2642)
+- Added `select()` function for vectorised ternary operations, as well as the `any()` and `all()` boolean reduction functions (#2853)
+- Added a minimum number of iterations (5) to the `grad_2f1` function to avoid early convergence (#2858)
+- Fixed some spelling errors in documentation (#2908)
+- Added new function in `stan::math::internal`, `finite_diff_hessian_vector_product_auto` (#2914)
+- Fixed an issue with the return type of `eigenvalues` being incorrect (#2915)
+- Delete jquery files from vendored Boost documentation (#2916)
+- Fixed fvar/higher-order autodiff compatibility with GLM distributions (#2917)
+- Fixed a type issue in the tests for inv_wishart_cholesky. (#2920)
+- The constraint tolerance can now be adjusted by defining the `STAN_MATH_CONSTRAINT_TOLERANCE` macro before including Stan headers. The default value remains `1E-8`. (#2921)
+- Use -dumpfullversion if available in makefiles (#2922)
+- Vectorised binary signature for `log_sum_exp` (#2930)
+- Added new functions `qr_thin`, `eigendecompose_sym`, `eigendecompose`, `complex_schur_decompose`, `svd`, and `csr_extract`. Each of these is equivalent to calling several existing functions, but should be more efficient due to sharing work. For example, `svd(m)` is equivalent to `(svd_U(m), singular_values(m), svd_V(m))`. (#2931)
+
 ======================================================================
 v4.6.2 (15 May 2023)
 ======================================================================

doxygen/doxygen.cfg

@@ -48,7 +48,7 @@ PROJECT_NAME           = "Stan Math Library"
 # could be handy for archiving the generated documentation or if some version
 # control system is used.
 
-PROJECT_NUMBER         = 4.6.2
+PROJECT_NUMBER         = 4.7.0
 
 # Using the PROJECT_BRIEF tag one can provide an optional one line description
 # for a project that appears at the top of each page and should give viewer a

stan/math/fwd/fun/mdivide_left.hpp

@@ -60,9 +60,6 @@ template <typename T1, typename T2,
 inline Eigen::Matrix<value_type_t<T2>, T1::RowsAtCompileTime,
                      T2::ColsAtCompileTime>
 mdivide_left(const T1& A, const T2& b) {
-  constexpr int S1 = T1::RowsAtCompileTime;
-  constexpr int C2 = T2::ColsAtCompileTime;
-
   check_square("mdivide_left", "A", A);
   check_multiplicable("mdivide_left", "A", A, "b", b);
   if (A.size() == 0) {

stan/math/fwd/fun/mdivide_left_tri_low.hpp

@@ -56,7 +56,6 @@ inline Eigen::Matrix<value_type_t<T2>, T1::RowsAtCompileTime,
                      T2::ColsAtCompileTime>
 mdivide_left_tri_low(const T1& A, const T2& b) {
   constexpr int S1 = T1::RowsAtCompileTime;
-  constexpr int C2 = T2::ColsAtCompileTime;
 
   check_square("mdivide_left_tri_low", "A", A);
   check_multiplicable("mdivide_left_tri_low", "A", A, "b", b);

stan/math/fwd/fun/mdivide_right.hpp

@@ -97,7 +97,6 @@ inline Eigen::Matrix<value_type_t<EigMat2>, EigMat1::RowsAtCompileTime,
 mdivide_right(const EigMat1& A, const EigMat2& b) {
   using T = typename value_type_t<EigMat2>::Scalar;
   constexpr int R1 = EigMat1::RowsAtCompileTime;
-  constexpr int C1 = EigMat1::ColsAtCompileTime;
   constexpr int R2 = EigMat2::RowsAtCompileTime;
   constexpr int C2 = EigMat2::ColsAtCompileTime;
 

stan/math/fwd/fun/mdivide_right_tri_low.hpp

@@ -98,7 +98,6 @@ inline Eigen::Matrix<value_type_t<EigMat2>, EigMat1::RowsAtCompileTime,
 mdivide_right_tri_low(const EigMat1& A, const EigMat2& b) {
   using T = typename value_type_t<EigMat2>::Scalar;
   constexpr int R1 = EigMat1::RowsAtCompileTime;
-  constexpr int C1 = EigMat1::ColsAtCompileTime;
   constexpr int R2 = EigMat2::RowsAtCompileTime;
   constexpr int C2 = EigMat2::ColsAtCompileTime;
   check_square("mdivide_right_tri_low", "b", b);

stan/math/opencl/kernel_generator/broadcast.hpp

@@ -98,7 +98,6 @@ class broadcast_
    * @return pair of indices - bottom and top diagonal
    */
   inline std::pair<int, int> extreme_diagonals() const {
-    int bottom, top;
     std::pair<int, int> arg_diags
         = this->template get_arg<0>().extreme_diagonals();
     return {Colwise ? std::numeric_limits<int>::min() : arg_diags.first,