diff options
Diffstat (limited to 'math/lapack++/files/patch-src:eigslv.cc')
-rw-r--r-- | math/lapack++/files/patch-src:eigslv.cc | 61 |
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) |