Changeset 198db5a in tol

Timestamp:

Dec 29, 2008, 4:44:56 PM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

1c5053a

Parents:

f743e2a

Message:

Avoiding internal cblas exit calling when some dimension is zero

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@354 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol/bmath/blinalg/tol_blas.cpp (modified) (4 diffs)

tol/bmath/blinalg/tol_blas.cpp

@@ -88 +88 @@
     C.SetAllValuesTo(0);
   }
-  const double* A_ = A.Data().Buffer();
-  const double* B_ = B.Data().Buffer();
-  double* C_ = C.GetData().GetBuffer();
-  cblas_dgemm
-  (
-    CblasRowMajor, TransA, TransB,
-    M, N, K,
-    alpha,
-    A_, Ac,
-    B_, Bc,
-    beta,
-    C_, N
-  );
+  if(N&&K&&M)
+  {
+    const double* A_ = A.Data().Buffer();
+    const double* B_ = B.Data().Buffer();
+    double* C_ = C.GetData().GetBuffer();
+    cblas_dgemm
+    (
+      CblasRowMajor, TransA, TransB,
+      M, N, K,
+      alpha,
+      A_, Ac,
+      B_, Bc,
+      beta,
+      C_, N
+    );
+  }
 //VBR: Comprobar si hay error
   return(0);
@@ -154 +157 @@
     C.SetAllValuesTo(0);
   }
-  cblas_dsyrk
-  (
-    CblasRowMajor, Uplo, Trans,
-    N, K,
-    alpha,
-    A.Data().Buffer(), Ac,
-    beta,
-    C.GetData().GetBuffer(), Cc
-  );
+  if(N&&K)
+  {
+    cblas_dsyrk
+    (
+      CblasRowMajor, Uplo, Trans,
+      N, K,
+      alpha,
+      A.Data().Buffer(), Ac,
+      beta,
+      C.GetData().GetBuffer(), Cc
+    );
+  }
 //VBR: Comprobar si hay error
   return(0);
@@ -205 +211 @@
   int lda = (Side==CblasLeft)?M:N;
   int ldb = N;
-  cblas_dtrmm(CblasRowMajor,Side,Uplo,TransA,Diag,M,N,alpha, 
-              A.Data   ().Buffer   (),lda,
-              X.GetData().GetBuffer(),ldb);
+  if(N&&M)
+  {
+    cblas_dtrmm(CblasRowMajor,Side,Uplo,TransA,Diag,M,N,alpha, 
+                A.Data   ().Buffer   (),lda,
+                X.GetData().GetBuffer(),ldb);
+  }
 //VBR: Comprobar si hay error
   return(0);
@@ -255 +264 @@
   int lda = (Side==CblasLeft)?M:N;
   int ldb = N;
-  cblas_dtrsm(CblasRowMajor,Side,Uplo,TransA,Diag,M,N,alpha, 
-              A.Data   ().Buffer   (),lda,
-              X.GetData().GetBuffer(),ldb);
+  if(N&&M)
+  {
+    cblas_dtrsm(CblasRowMajor,Side,Uplo,TransA,Diag,M,N,alpha, 
+                A.Data   ().Buffer   (),lda,
+                X.GetData().GetBuffer(),ldb);
+  }
 //VBR: Comprobar si hay error
   return(0);