ext/nmatrix/math/rotg.h in nmatrix-0.1.0 vs ext/nmatrix/math/rotg.h in nmatrix-0.2.0
- old
+ new
@@ -83,33 +83,24 @@
*b = z;
}
}
template <>
-inline void rotg(float* a, float* b, float* c, float* s) {
- cblas_srotg(a, b, c, s);
-}
-
-template <>
-inline void rotg(double* a, double* b, double* c, double* s) {
- cblas_drotg(a, b, c, s);
-}
-
-template <>
inline void rotg(Complex64* a, Complex64* b, Complex64* c, Complex64* s) {
- cblas_crotg(reinterpret_cast<void*>(a), reinterpret_cast<void*>(b), reinterpret_cast<void*>(c), reinterpret_cast<void*>(s));
+ rb_raise(rb_eNotImpError, "BLAS not available, and existing template requires modification for complex");
}
template <>
inline void rotg(Complex128* a, Complex128* b, Complex128* c, Complex128* s) {
- cblas_zrotg(reinterpret_cast<void*>(a), reinterpret_cast<void*>(b), reinterpret_cast<void*>(c), reinterpret_cast<void*>(s));
+ rb_raise(rb_eNotImpError, "BLAS not available, and existing template requires modification for complex");
}
+
template <typename DType>
inline void cblas_rotg(void* a, void* b, void* c, void* s) {
rotg<DType>(reinterpret_cast<DType*>(a), reinterpret_cast<DType*>(b), reinterpret_cast<DType*>(c), reinterpret_cast<DType*>(s));
}
} } //nm::math
-#endif // ROTG_H
\ No newline at end of file
+#endif // ROTG_H