workaround missing blasfeo_dtrcpsc_l symbol on skylake

PREDICT-EPFL · Nov 20, 2024 · 3225493 · 3225493
1 parent 1a4b6ac
commit 3225493
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/include/piqp/sparse/blocksparse_stage_kkt.hpp b/include/piqp/sparse/blocksparse_stage_kkt.hpp
@@ -972,7 +972,14 @@ class BlocksparseStageKKT : public KKTSystem<T>
                         blasfeo_dgead(m, m, delta_inv, AtA.D[i]->ref(), 0, 0, kkt_mat.D[i]->ref(), 0, 0);
                     } else {
                         // D_i = delta^{-1} * AtA.D_i, lower triangular
+#ifdef TARGET_X64_INTEL_SKYLAKE_X
+                        // blasfeo_dtrcpsc_l not implemented on Skylake yet
+                        // and reference implementation not exported ...
+                        kkt_mat.D[i]->setZero();
+                        blasfeo_dgead(m, m, delta_inv, AtA.D[i]->ref(), 0, 0, kkt_mat.D[i]->ref(), 0, 0);
+#else
                         blasfeo_dtrcpsc_l(m, delta_inv, AtA.D[i]->ref(), 0, 0, kkt_mat.D[i]->ref(), 0, 0);
+#endif
                         mat_set = true;
                     }
                 }