summaryrefslogtreecommitdiffstats
path: root/math/lapack++/files/patch-src:eigslv.cc
diff options
context:
space:
mode:
Diffstat (limited to 'math/lapack++/files/patch-src:eigslv.cc')
-rw-r--r--math/lapack++/files/patch-src:eigslv.cc61
1 files changed, 0 insertions, 61 deletions
diff --git a/math/lapack++/files/patch-src:eigslv.cc b/math/lapack++/files/patch-src:eigslv.cc
deleted file mode 100644
index c087bce..0000000
--- a/math/lapack++/files/patch-src:eigslv.cc
+++ /dev/null
@@ -1,61 +0,0 @@
---- src/eigslv.cc~ Wed Mar 1 02:52:47 2000
-+++ src/eigslv.cc Wed May 22 21:49:45 2002
-@@ -58,14 +58,14 @@
- // char *fname = NULL; // HP C++ does not support string initalization!
- //#endif
-
-- // view symmetric matrix S as a LaGenMatDouble
-- // assumes S is contiguous...
-- LaSymmMatDouble tmp(S);
-- LaEigSolveVecIP(tmp, eigvals);
--
-- LaGenMatDouble tmp2(&tmp(0,0), S.size(0), S.size(1));
-- eigvec.ref(tmp2);
--
-+ long int N = S.size(0);
-+ long int i,j;
-+
-+ for(j=0;j<N;j++){
-+ for(i=j;i<N;i++){
-+ eigvec(i,j)=S(i,j);
-+ }}
-+ LaEigSolveVecIP(eigvec, eigvals);
-
- }
-
-@@ -100,18 +100,18 @@
-
- }
-
--void LaEigSolveVecIP(LaSymmMatDouble &S, LaVectorDouble &eigvals)
-+void LaEigSolveVecIP(LaGenMatDouble &A, LaVectorDouble &eigvals)
- {
- #ifndef HPPA
- const char fname[] = "LaEigSolveVecIP(LaGenMatDouble &A, &eigvals)";
- #else
- char *fname = NULL; // HP C++ does not support string initalization!
- #endif
-- long int N = S.size(0);
-+ long int N = A.size(0);
- char jobz = 'V';
- char uplo = 'L';
- long int info;
-- long int lda = S.gdim(0);
-+ long int lda = A.gdim(0);
-
- if (eigvals.size() < N)
- {
-@@ -119,11 +119,11 @@
- }
-
-
-- long int w = (LaEnvBlockSize("SSYTRD", S) +2) * N;
-+ long int w = (LaEnvBlockSize("SSYTRD", A) +2) * N;
- LaVectorDouble Work(w);
-
-
-- F77NAME(dsyev)(&jobz, &uplo, &N, S.addr(), &lda, &eigvals(0), &Work(0),
-+ F77NAME(dsyev)(&jobz, &uplo, &N, A.addr(), &lda, &eigvals(0), &Work(0),
- &w, &info);
-
- if (info != 0)
OpenPOWER on IntegriCloud