LU: workaround for panel threads in hetrf

icl-utk-edu · Oct 21, 2023 · 8354553 · 8354553
1 parent 2953ea6
commit 8354553
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/hetrf.cc b/src/hetrf.cc
@@ -50,7 +50,9 @@ int64_t hetrf(
         = get_option<double>( opts, Option::PivotThreshold, 1.0 );
     int64_t lookahead = get_option<int64_t>( opts, Option::Lookahead, 1 );
     int64_t ib = get_option<int64_t>( opts, Option::InnerBlocking, 16 );
-    int64_t max_panel_threads = std::max( omp_get_max_threads()/2, 1 );
+
+    // Using > 1 thread leads to hang, reason unclear.
+    int64_t max_panel_threads = 1;  //std::max( omp_get_max_threads()/2, 1 );
     max_panel_threads = get_option<int64_t>( opts, Option::MaxPanelThreads,
                                              max_panel_threads );