summaryrefslogtreecommitdiff
path: root/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c
diff options
context:
space:
mode:
Diffstat (limited to 'lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c')
-rw-r--r--lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c22
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
1066int 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
1066int transC(KCMAT(x),CMAT(t)) { 1078int 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
1112int 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
1100int constantC(doublecomplex* pval, CVEC(r)) { 1122int constantC(doublecomplex* pval, CVEC(r)) {
1101 DEBUGMSG("constantC") 1123 DEBUGMSG("constantC")
1102 int k; 1124 int k;