aboutsummaryrefslogtreecommitdiff
path: root/redist/dclhelpers.c
diff options
context:
space:
mode:
authorcnlohr <lohr85@gmail.com>2018-03-17 01:00:07 -0400
committercnlohr <lohr85@gmail.com>2018-03-17 01:00:07 -0400
commit783a88f8904b5758195f96c65c7e0b6e87f8a51e (patch)
tree8c78896e796ad33a2d2b8cb9e7766a164413a351 /redist/dclhelpers.c
parent1fe0debd6474c0368f81159cd6a0451edbc6d381 (diff)
downloadlibsurvive-783a88f8904b5758195f96c65c7e0b6e87f8a51e.tar.gz
libsurvive-783a88f8904b5758195f96c65c7e0b6e87f8a51e.tar.bz2
Update DCLapack with the C helpers to cize stuff.
Diffstat (limited to 'redist/dclhelpers.c')
-rw-r--r--redist/dclhelpers.c63
1 files changed, 63 insertions, 0 deletions
diff --git a/redist/dclhelpers.c b/redist/dclhelpers.c
new file mode 100644
index 0000000..9d3d0eb
--- /dev/null
+++ b/redist/dclhelpers.c
@@ -0,0 +1,63 @@
+#include "dclhelpers.h"
+#define FLOAT DCL_FLOAT
+#define DYNAMIC_INDEX
+#include <stdio.h>
+#include "dclapack.h"
+
+
+void dclPrint( const DCL_FLOAT * A, int n, int m )
+{
+ PRINT( A, n, m );
+}
+
+void dclIdentity( DCL_FLOAT * A, int n )
+{
+ IDENTITY( A, n );
+}
+
+void dclTransp( const DCL_FLOAT * A, DCL_FLOAT * B, int n, int m)
+{
+ TRANSP(A,B,n,m);
+}
+
+void dclLU( const DCL_FLOAT * A, DCL_FLOAT * L, DCL_FLOAT * U, int * Piv, int n )
+{
+ LU(A,L,U,Piv,n);
+}
+
+void dclPivot( const DCL_FLOAT * A, DCL_FLOAT * B, int * Piv, int n, int m )
+{
+ PIVOT(A,B,Piv,n,m);
+}
+
+void dclLSub( const DCL_FLOAT * L, DCL_FLOAT * X, const DCL_FLOAT * B, int n, int m )
+{
+ L_SUB(L,X,B,n,m);
+}
+
+void dclUSub( const DCL_FLOAT * U, DCL_FLOAT * X, const DCL_FLOAT * B, int n, int m )
+{
+ U_SUB(U,X,B,n,m);
+}
+
+void dclInv( const DCL_FLOAT * A, DCL_FLOAT * Ainv, int n )
+{
+ INV(A,Ainv,n,n);
+}
+
+void dclMul( const DCL_FLOAT * A, const DCL_FLOAT * B, DCL_FLOAT * C, int n, int m, int p )
+{
+ MUL(A,B,C,n,m,p);
+}
+
+void dclMulAdd( const DCL_FLOAT * A, const DCL_FLOAT * B, const DCL_FLOAT * C, DCL_FLOAT * D, int n, int m, int p )
+{
+ MULADD(A,B,C,D,n,m,p);
+}
+
+void dclGMulAdd( const DCL_FLOAT * A, const DCL_FLOAT * B, const DCL_FLOAT * C, DCL_FLOAT * D, DCL_FLOAT alpha, DCL_FLOAT beta, int n, int m, int p )
+{
+ GMULADD(A,B,C,D,alpha,beta,n,m,p);
+}
+
+