diff options
Diffstat (limited to 'lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c')
-rw-r--r-- | lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c | 22 |
1 files changed, 22 insertions, 0 deletions
diff --git a/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c b/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c index b9c2572..7a40991 100644 --- a/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c +++ b/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c | |||
@@ -1063,6 +1063,18 @@ int transR(KDMAT(x),DMAT(t)) { | |||
1063 | OK | 1063 | OK |
1064 | } | 1064 | } |
1065 | 1065 | ||
1066 | int transQ(KQMAT(x),QMAT(t)) { | ||
1067 | REQUIRES(xr==tc && xc==tr,BAD_SIZE); | ||
1068 | DEBUGMSG("transQ"); | ||
1069 | int i,j; | ||
1070 | for (i=0; i<tr; i++) { | ||
1071 | for (j=0; j<tc; j++) { | ||
1072 | ((complex*)tp)[i*tc+j] = ((complex*)xp)[j*xc+i]; | ||
1073 | } | ||
1074 | } | ||
1075 | OK | ||
1076 | } | ||
1077 | |||
1066 | int transC(KCMAT(x),CMAT(t)) { | 1078 | int transC(KCMAT(x),CMAT(t)) { |
1067 | REQUIRES(xr==tc && xc==tr,BAD_SIZE); | 1079 | REQUIRES(xr==tc && xc==tr,BAD_SIZE); |
1068 | DEBUGMSG("transC"); | 1080 | DEBUGMSG("transC"); |
@@ -1097,6 +1109,16 @@ int constantR(double * pval, DVEC(r)) { | |||
1097 | OK | 1109 | OK |
1098 | } | 1110 | } |
1099 | 1111 | ||
1112 | int constantQ(complex* pval, QVEC(r)) { | ||
1113 | DEBUGMSG("constantQ") | ||
1114 | int k; | ||
1115 | complex val = *pval; | ||
1116 | for(k=0;k<rn;k++) { | ||
1117 | ((complex*)rp)[k]=val; | ||
1118 | } | ||
1119 | OK | ||
1120 | } | ||
1121 | |||
1100 | int constantC(doublecomplex* pval, CVEC(r)) { | 1122 | int constantC(doublecomplex* pval, CVEC(r)) { |
1101 | DEBUGMSG("constantC") | 1123 | DEBUGMSG("constantC") |
1102 | int k; | 1124 | int k; |