From e2db4d1f08a8b9ac5bb500f0a248316784369b7d Mon Sep 17 00:00:00 2001
From: Daniel Lee <bearlee@alum.mit.edu>
Date: Mon, 5 Feb 2024 20:59:45 -0500
Subject: [PATCH] updating opencl multi_normal_cholesky_lpdf to match the other
 implementation

---
 stan/math/opencl/prim/multi_normal_cholesky_lpdf.hpp | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/stan/math/opencl/prim/multi_normal_cholesky_lpdf.hpp b/stan/math/opencl/prim/multi_normal_cholesky_lpdf.hpp
index 890b1e22958..728c8e65256 100644
--- a/stan/math/opencl/prim/multi_normal_cholesky_lpdf.hpp
+++ b/stan/math/opencl/prim/multi_normal_cholesky_lpdf.hpp
@@ -44,9 +44,11 @@ inline return_type_t<T_y_cl, T_loc_cl, T_covar_cl> multi_normal_cholesky_lpdf(
   static const char* function = "multi_normal_cholesky_lpdf(OpenCL)";
 
   check_consistent_sizes(function, "y", y, "mu", mu);
-  check_square(function, "covariance parameter", L);
   check_size_match(function, "Size of random variable", y.rows(),
                    "rows of covariance parameter", L.rows());
+  check_square(function, "Cholesky decomposition of a variance matrix", L);
+  check_cholesky_factor(function, "Cholesky decomposition of a variance matrix",
+                        L);
 
   if (max_size(y, mu, L) == 0) {
     return 0.0;