mirror of
https://git.freebsd.org/ports.git
synced 2025-05-28 08:56:28 -04:00
789b5c5923
2. change e-mail address of maho. 3. honor CXX and CXXFLAGS. 4. modernize patch file names. 5. satisfy portlint. Reviewed by: knu (mentor) Reported by: bento
61 lines
1.5 KiB
C++
61 lines
1.5 KiB
C++
--- 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)
|