From c171b8ad13054518869cdc54db5af5cf6b886089 Mon Sep 17 00:00:00 2001
From: Martin Kroeker <martin@ruby.chemie.uni-freiburg.de>
Date: Tue, 13 Nov 2018 13:57:18 +0100
Subject: [PATCH] Handle special case INCX=0,INCY=0 in the axpy interface

---
 interface/axpy.c  | 5 +++++
 interface/zaxpy.c | 6 ++++++
 2 files changed, 11 insertions(+)

diff --git a/interface/axpy.c b/interface/axpy.c
index 39edea6af..9032946d2 100644
--- a/interface/axpy.c
+++ b/interface/axpy.c
@@ -75,6 +75,11 @@ void CNAME(blasint n, FLOAT alpha, FLOAT *x, blasint incx, FLOAT *y, blasint inc
 
   if (alpha == ZERO) return;
 
+  if (incx == 0 && incy == 0) {
+    *y += n * alpha *(*x);
+    return;
+  }
+    
   IDEBUG_START;
 
   FUNCTION_PROFILE_START();
diff --git a/interface/zaxpy.c b/interface/zaxpy.c
index 1a0259c96..dbd559628 100644
--- a/interface/zaxpy.c
+++ b/interface/zaxpy.c
@@ -82,6 +82,12 @@ void CNAME(blasint n, FLOAT *ALPHA, FLOAT *x, blasint incx, FLOAT *y, blasint in
 
   if ((alpha_r == ZERO) && (alpha_i == ZERO)) return;
 
+  if (incx == 0 && incy == 0) {
+  *y += n * (alpha_r * (*x) - alpha_i* (*(x+1)) );
+  *(y+1) += n * (alpha_i * (*x) + alpha_r * (*(x +1)) );
+  return;
+  }
+  
   IDEBUG_START;
 
   FUNCTION_PROFILE_START();