summaryrefslogtreecommitdiff
path: root/src/matrixOperations/powm/spowma.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/matrixOperations/powm/spowma.c')
-rw-r--r--src/matrixOperations/powm/spowma.c33
1 files changed, 2 insertions, 31 deletions
diff --git a/src/matrixOperations/powm/spowma.c b/src/matrixOperations/powm/spowma.c
index 65bfb74f..aa0f7043 100644
--- a/src/matrixOperations/powm/spowma.c
+++ b/src/matrixOperations/powm/spowma.c
@@ -11,39 +11,13 @@
*/
#include "matrixPow.h"
-#include "lapack.h"
#include "eye.h"
#include "matrixMultiplication.h"
void spowma(float* in, int size, float expand, float* out){
-#ifndef WITHOUT_BLAS
- switch ((int)expand){
- case 0 :
- seyea(out,size,size);
- break;
- case 1 :
- {
- int i;
- for (i=0;i<size*size;i++) out[i]=in[i];
- }
- break;
- default :
- {
- int i=0,j=0,k=0;
- int One=1;
- for (i=1; i<expand; i++){
- for (j=0;j<size;j++) /*column*/ {
- for (k=0;k<size;k++)/*row*/{
- out[k+j*size]=ddot_(&size,in+k,&size,in+j*size,&One);
- }
- }
- }
- }
- break;
- }
-#else
+
switch ((int)expand){
case 0 :
seyea(out,size,size);
@@ -67,8 +41,5 @@ void spowma(float* in, int size, float expand, float* out){
}
break;
}
-
-
-#endif
-
+
}