From 72b8ce30e268f7eefc86a5fea6a731b6eef4f707 Mon Sep 17 00:00:00 2001
From: Jiri Pittner <jiri@pittnerovi.com>
Date: Mon, 20 Oct 2025 17:27:43 +0200
Subject: [PATCH] bugfix in cublas diagmultr

---
 mat.cc | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/mat.cc b/mat.cc
index d6fe756..df652fa 100644
--- a/mat.cc
+++ b/mat.cc
@@ -2031,7 +2031,7 @@ void NRMat<double>::diagmultr(const NRVec<double> &rhs) {
 		for(register int i=0; i<mm; i++){ cblas_dscal(nn, rhs[i], &(*this)(0, i), mm); }
 #ifdef CUDALA
 	}else{
-		for(register int i=0; i<mm; i++){ cublasDscal(mm, rhs[i], v + i, mm); }
+		for(register int i=0; i<mm; i++){ cublasDscal(nn, rhs[i], v + i, mm); }
 	}
 #endif
 }