1 files changed, 0 insertions, 78 deletions

diff --git a/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c b/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c
index 61cb002..310f6ee 100644
--- a/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c
+++ b/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c
@@ -814,81 +814,3 @@ int lu_l_C(KCMAT(a), DVEC(ipiv), CMAT(r)) {
     free(auxipiv);
     OK
 }
-
-////////////////////////////////////////////////////////////
-
-// int multiplyR(KDMAT(a),KDMAT(b),DMAT(r)) {
-//     REQUIRES(ac==br && ar==rr && bc==rc,BAD_SIZE);
-//     int i,j,k;
-//     for (i=0;i<ar;i++) {
-//         for(j=0;j<bc;j++) {
-//             double temp = 0;
-//             for(k=0;k<ac;k++) {
-//                 temp += ap[i*ac+k]*bp[k*bc+j];
-//             }
-//             rp[i*rc+j] = temp;
-//         }
-//     }
-//     OK
-// }
-
-
-int multiplyC(KCMAT(a),KCMAT(b),CMAT(r)) {
-    REQUIRES(ac==br && ar==rr && bc==rc,BAD_SIZE);
-    int i,j,k;
-    for (i=0;i<ar;i++) {
-        for(j=0;j<bc;j++) {
-            doublecomplex temp = {0,0};
-            for(k=0;k<ac;k++) {
-                doublecomplex A = ((doublecomplex*)ap)[i*ac+k];
-                doublecomplex B = ((doublecomplex*)bp)[k*bc+j];
-                double w = A.r * B.r - A.i * B.i;
-                double w1 = A.r * B.r;
-                double w2 = A.i * B.i;
-                if(w != w) {
-                    printf("at : %d %d %d\n",i,j,k);
-                    printf("%f %f %f\n",A.i,B.i, A.i * B.i);
-                    printf("%f %f %f\n",A.i,B.i, w2);
-                    exit(1);
-                }
-                temp.r += (w + w1-w2)/2;
-                //temp.r += w;
-                temp.i += A.r * B.i + A.i * B.r;
-            }
-            ((doublecomplex*)rp)[i*rc+j] = temp;
-        }
-    }
-    OK
-}
-
-// void dgemm_(char *, char *, integer *, integer *, integer *,
-//            double *, const double *, integer *, const double *,
-//            integer *, double *, double *, integer *);
-// 
-// void zgemm_(char *, char *, integer *, integer *, integer *,
-//            doublecomplex *, const doublecomplex *, integer *, const doublecomplex *,
-//            integer *, doublecomplex *, doublecomplex *, integer *);
-// 
-// 
-// int multiplyR2(KDMAT(a),KDMAT(b),DMAT(r)) {
-//     REQUIRES(ac==br && ar==rr && bc==rc,BAD_SIZE);
-//     double alpha = 1;
-//     double beta = 0;
-//     integer m = ar;
-//     integer n = bc;
-//     integer k = ac;
-//     int i,j;
-//     dgemm_("N","N",&m,&n,&k,&alpha,ap,&m,bp,&k,&beta,rp,&m);
-//     OK
-// }
-// 
-// int multiplyC2(KCMAT(a),KCMAT(b),CMAT(r)) {
-//     REQUIRES(ac==br && ar==rr && bc==rc,BAD_SIZE);
-//     integer m = ar;
-//     integer n = bc;
-//     integer k = ac;
-//     doublecomplex alpha = {1,0};
-//     doublecomplex beta = {0,0};
-//     zgemm_("N","N",&m,&n,&k,&alpha,(doublecomplex*)ap,&m,(doublecomplex*)bp,&k,&beta,(doublecomplex*)rp,&m);
-//     OK
-// }