Add shortcuts for (small) cases that do not need expensive buffer allocation

diff --git a/interface/trsv.c b/interface/trsv.c
index a054d8e..6a6e8f8 100644 (file)
--- a/interface/trsv.c
+++ b/interface/trsv.c
@@ -188,6 +188,12 @@ void CNAME(enum CBLAS_ORDER order, enum CBLAS_UPLO Uplo,
 
   if (n == 0) return;
 
+  if (incx == 1 && trans == 0 && n < 50) {
+    buffer = NULL;
+    (trsv[(trans<<2) | (uplo<<1) | unit])(n, a, lda, x, incx, buffer);
+    return;  
+  }
+
   IDEBUG_START;
 
   FUNCTION_PROFILE_START();

diff --git a/interface/ztrsv.c b/interface/ztrsv.c
index cbb7bba..cf750b0 100644 (file)
--- a/interface/ztrsv.c
+++ b/interface/ztrsv.c
@@ -199,6 +199,12 @@ void CNAME(enum CBLAS_ORDER order, enum CBLAS_UPLO Uplo,
 
   if (n == 0) return;
 
+  if (incx == 1 && trans == 0 && n < 50) {
+    buffer = NULL;
+    (trsv[(trans<<2) | (uplo<<1) | unit])(n, a, lda, x, incx, buffer);
+    return;
+  }
+
   IDEBUG_START;
 
   FUNCTION_PROFILE_START();