Armadillo 12.8.3

Author: eddelbuettel

inst/include/armadillo_bits/arma_version.hpp

@@ -23,7 +23,7 @@
 
 #define ARMA_VERSION_MAJOR 12
 #define ARMA_VERSION_MINOR 8
-#define ARMA_VERSION_PATCH 2
+#define ARMA_VERSION_PATCH 3
 #define ARMA_VERSION_NAME  "Cortisol Injector"
 
 

inst/include/armadillo_bits/compiler_setup.hpp

@@ -377,14 +377,11 @@
 
 
 #if defined(__SUNPRO_CC)
-  
   // http://www.oracle.com/technetwork/server-storage/solarisstudio/training/index-jsp-141991.html
   // http://www.oracle.com/technetwork/server-storage/solarisstudio/documentation/cplusplus-faq-355066.html
-  
   #if (__SUNPRO_CC < 0x5140)
     #error "*** newer compiler required ***"
   #endif
-  
 #endif
 
 

inst/include/armadillo_bits/config.hpp

@@ -201,6 +201,11 @@
 
 #if defined(ARMA_EXTRA_DEBUG)
   #undef  ARMA_NO_DEBUG
+  #undef  ARMA_DONT_CHECK_NONFINITE
+  
+  #undef  ARMA_CHECK_NONFINITE
+  #define ARMA_CHECK_NONFINITE
+  
   #undef  ARMA_WARN_LEVEL
   #define ARMA_WARN_LEVEL 3
 #endif

inst/include/armadillo_bits/debug.hpp

@@ -579,7 +579,7 @@ arma_assert_same_size(const uword A_n_rows, const uword A_n_cols, const uword B_
 
 
 
-//! stop if given matrices have different sizes
+//! stop if given matrices do not have the same size
 template<typename eT1, typename eT2>
 arma_hot
 inline
@@ -600,7 +600,7 @@ arma_assert_same_size(const Mat<eT1>& A, const Mat<eT2>& B, const char* x)
 
 
 
-//! stop if given proxies have different sizes
+//! stop if given proxies do not have the same size
 template<typename eT1, typename eT2>
 arma_hot
 inline
@@ -804,7 +804,7 @@ arma_assert_same_size(const uword A_n_rows, const uword A_n_cols, const uword A_
 
 
 
-//! stop if given cubes have different sizes
+//! stop if given cubes do not have the same size
 template<typename eT1, typename eT2>
 arma_hot
 inline
@@ -883,7 +883,7 @@ arma_assert_same_size(const subview_cube<eT>& A, const ProxyCube<T1>& B, const c
 
 
 
-//! stop if given cube proxies have different sizes
+//! stop if given cube proxies do not have the same size
 template<typename eT1, typename eT2>
 arma_hot
 inline

inst/include/armadillo_bits/fft_engine_fftw3.hpp

@@ -23,6 +23,13 @@
 
 #if defined(ARMA_USE_FFTW3)
 
+struct fft_engine_fftw3_aux
+  {
+  #if (!defined(ARMA_DONT_USE_STD_MUTEX))
+  static inline std::mutex& get_plan_mutex() { static std::mutex plan_mutex; return plan_mutex; }
+  #endif
+  };
+
 template<typename cx_type, bool inverse>
 class fft_engine_fftw3
   {
@@ -74,7 +81,30 @@ class fft_engine_fftw3
     const int fftw3_flags = fftw3_flag_destroy | fftw3_flag_estimate;
 
     arma_extra_debug_print("fft_engine_fftw3::constructor: generating 1D plan");
-    fftw3_plan = fftw3::plan_dft_1d<cx_type>(N, X_work.memptr(), Y_work.memptr(), fftw3_sign, fftw3_flags);
+    
+    // only fftw3::execute() is thread safe, as per FFTW docs:
+    // https://www.fftw.org/fftw3_doc/Thread-safety.html
+    
+    #if defined(ARMA_USE_OPENMP)
+      {
+      #pragma omp critical (arma_fft_engine_fftw3)
+        {
+        fftw3_plan = fftw3::plan_dft_1d<cx_type>(N, X_work.memptr(), Y_work.memptr(), fftw3_sign, fftw3_flags);
+        }
+      }
+    #elif (!defined(ARMA_DONT_USE_STD_MUTEX))
+      {
+      std::mutex& plan_mutex = fft_engine_fftw3_aux::get_plan_mutex();
+      
+      const std::lock_guard<std::mutex> lock(plan_mutex);
+      
+      fftw3_plan = fftw3::plan_dft_1d<cx_type>(N, X_work.memptr(), Y_work.memptr(), fftw3_sign, fftw3_flags);
+      }
+    #else
+      {
+      fftw3_plan = fftw3::plan_dft_1d<cx_type>(N, X_work.memptr(), Y_work.memptr(), fftw3_sign, fftw3_flags);
+      }
+    #endif
 
     if(fftw3_plan == nullptr)  { arma_stop_runtime_error("fft_engine_fftw3::constructor: failed to create plan"); }
     }