From 725629ca73eae311c62aa076c71efb6cdb117e1e Mon Sep 17 00:00:00 2001
From: conrad <conrad@localhost.localdomain>
Date: Fri, 25 Mar 2022 13:23:53 +1000
Subject: [PATCH] detect NaN

---
 include/armadillo_bits/glue_quantile_meat.hpp | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/include/armadillo_bits/glue_quantile_meat.hpp b/include/armadillo_bits/glue_quantile_meat.hpp
index 6ad73d64..6f456c62 100644
--- a/include/armadillo_bits/glue_quantile_meat.hpp
+++ b/include/armadillo_bits/glue_quantile_meat.hpp
@@ -178,6 +178,8 @@ glue_quantile::apply(Mat<typename T2::elem_type>& out, const mtGlue<typename T2:
   const quasi_unwrap<T1> UA(expr.A);
   const quasi_unwrap<T2> UB(expr.B);
   
+  arma_debug_check((UA.M.has_nan() || UB.M.has_nan()), "quantile(): detected NaN");
+  
   if(UA.is_alias(out) || UB.is_alias(out))
     {
     Mat<eTb> tmp;
@@ -208,6 +210,8 @@ glue_quantile_default::apply(Mat<typename T2::elem_type>& out, const mtGlue<type
   
   const uword dim = (T1::is_xvec) ? uword(UA.M.is_rowvec() ? 1 : 0) : uword((T1::is_row) ? 1 : 0);
   
+  arma_debug_check((UA.M.has_nan() || UB.M.has_nan()), "quantile(): detected NaN");
+  
   if(UA.is_alias(out) || UB.is_alias(out))
     {
     Mat<eTb> tmp;