diff options
Diffstat (limited to 'src/matrixOperations/powm/spowma.c')
-rw-r--r-- | src/matrixOperations/powm/spowma.c | 33 |
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 - + } |