aboutsummaryrefslogtreecommitdiff
path: root/redist/dclhelpers.c
blob: 9d3d0eb227b269535260a860f1590c512fa695a8 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
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);
}