Fix

Author: michel2323

deps/src/onemkl_dft.cpp

@@ -186,44 +186,111 @@ int onemklDftGetValueConfigValue(onemklDftDescriptor_t desc, onemklDftConfigPara
         } \
     } while (0)
 
+// Pointer (USM) dispatch with proper element typing rather than using void* directly.
+// Using void* caused instantiation of compute_forward/backward with <void> template
+// parameters on some oneMKL versions, leading to unresolved symbols at runtime.
 int onemklDftComputeForward(onemklDftDescriptor_t desc, void *inout) {
     if (!desc || !inout) return -2;
     try {
-        ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, inout));
+        if (desc->dom == domain::REAL) {
+            if (desc->prec == precision::SINGLE) {
+                auto *p = static_cast<float*>(inout);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, p));
+            } else {
+                auto *p = static_cast<double*>(inout);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, p));
+            }
+        } else { // COMPLEX
+            if (desc->prec == precision::SINGLE) {
+                auto *p = static_cast<std::complex<float>*>(inout);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, p));
+            } else {
+                auto *p = static_cast<std::complex<double>*>(inout);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, p));
+            }
+        }
         return 0;
-    } catch (...) {
-        return -1;
-    }
+    } catch (...) { return -1; }
 }
 
 int onemklDftComputeForwardOutOfPlace(onemklDftDescriptor_t desc, void *in, void *out) {
     if (!desc || !in || !out) return -2;
     try {
-        ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, in, out));
+        if (desc->dom == domain::REAL) {
+            if (desc->prec == precision::SINGLE) {
+                auto *pi = static_cast<float*>(in);
+                auto *po = static_cast<float*>(out);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, pi, po));
+            } else {
+                auto *pi = static_cast<double*>(in);
+                auto *po = static_cast<double*>(out);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, pi, po));
+            }
+        } else { // COMPLEX
+            if (desc->prec == precision::SINGLE) {
+                auto *pi = static_cast<std::complex<float>*>(in);
+                auto *po = static_cast<std::complex<float>*>(out);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, pi, po));
+            } else {
+                auto *pi = static_cast<std::complex<double>*>(in);
+                auto *po = static_cast<std::complex<double>*>(out);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_forward(*d, pi, po));
+            }
+        }
         return 0;
-    } catch (...) {
-        return -1;
-    }
+    } catch (...) { return -1; }
 }
 
 int onemklDftComputeBackward(onemklDftDescriptor_t desc, void *inout) {
     if (!desc || !inout) return -2;
     try {
-        ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, inout));
+        if (desc->dom == domain::REAL) {
+            if (desc->prec == precision::SINGLE) {
+                auto *p = static_cast<float*>(inout);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, p));
+            } else {
+                auto *p = static_cast<double*>(inout);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, p));
+            }
+        } else { // COMPLEX
+            if (desc->prec == precision::SINGLE) {
+                auto *p = static_cast<std::complex<float>*>(inout);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, p));
+            } else {
+                auto *p = static_cast<std::complex<double>*>(inout);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, p));
+            }
+        }
         return 0;
-    } catch (...) {
-        return -1;
-    }
+    } catch (...) { return -1; }
 }
 
 int onemklDftComputeBackwardOutOfPlace(onemklDftDescriptor_t desc, void *in, void *out) {
     if (!desc || !in || !out) return -2;
     try {
-        ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, in, out));
+        if (desc->dom == domain::REAL) {
+            if (desc->prec == precision::SINGLE) {
+                auto *pi = static_cast<float*>(in);
+                auto *po = static_cast<float*>(out);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, pi, po));
+            } else {
+                auto *pi = static_cast<double*>(in);
+                auto *po = static_cast<double*>(out);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, pi, po));
+            }
+        } else { // COMPLEX
+            if (desc->prec == precision::SINGLE) {
+                auto *pi = static_cast<std::complex<float>*>(in);
+                auto *po = static_cast<std::complex<float>*>(out);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, pi, po));
+            } else {
+                auto *pi = static_cast<std::complex<double>*>(in);
+                auto *po = static_cast<std::complex<double>*>(out);
+                ONEMKL_DFT_DISPATCH(desc->ptr, compute_backward(*d, pi, po));
+            }
+        }
         return 0;
-    } catch (...) {
-        return -1;
-    }
+    } catch (...) { return -1; }
 }
 
 // Keep dispatch macros defined for buffer variants below; undef at end of file.