summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat (limited to 'btl/accuracy/lapack/lapack_STEV.hh')
-rw-r--r--btl/accuracy/lapack/lapack_STEV.hh52
1 files changed, 0 insertions, 52 deletions
diff --git a/btl/accuracy/lapack/lapack_STEV.hh b/btl/accuracy/lapack/lapack_STEV.hh
deleted file mode 100644
index 3fc9147..0000000
--- a/btl/accuracy/lapack/lapack_STEV.hh
+++ /dev/null
@@ -1,52 +0,0 @@
-//=====================================================
-// Copyright (C) 2011 Andrea Arteaga <andyspiros@gmail.com>
-//=====================================================
-//
-// This program is free software; you can redistribute it and/or
-// modify it under the terms of the GNU General Public License
-// as published by the Free Software Foundation; either version 2
-// of the License, or (at your option) any later version.
-//
-// This program is distributed in the hope that it will be useful,
-// but WITHOUT ANY WARRANTY; without even the implied warranty of
-// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-// GNU General Public License for more details.
-// You should have received a copy of the GNU General Public License
-// along with this program; if not, write to the Free Software
-// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
-//
-#ifndef LAPACK_STEV_HH
-#define LAPACK_STEV_HH
-
-#include "LinearCongruential.hh"
-#include "diff.hh"
-
-double test_STEV(const int& N, const unsigned& seed = 0)
-{
- LinearCongruential lc(seed);
- vector<double> A(N*N), D(N), E(N-1), Z(N*N), DD(N*N), work(max(1, 2*N-2)), tmp(N*N);
-
- /* Fill D, E and A */
- for (int i = 0; i < N-1; ++i) {
- D[i] = lc.get_01(); A[i+i*N] = D[i];
- E[i] = lc.get_01(); A[(i+1)+i*N] = E[i]; A[i+(i+1)*N] = E[i];
- }
- D[N-1] = lc.get_01(); A[N*N-1] = D[N-1];
-
- /* Perform computation */
- int info;
- dstev_("V", &N, &D[0], &E[0], &Z[0], &N, &work[0], &info);
-
- /* Construct DD */
- for (int i = 0; i < N; ++i)
- DD[i+i*N] = D[i];
-
- /* A = V*D*V' */
- const double dONE = 1., dZERO = 0.;
- dgemm_("N", "N", &N, &N, &N, &dONE, &Z[0], &N, &DD[0], &N, &dZERO, &tmp[0], &N);
- dgemm_("N", "T", &N, &N, &N, &dONE, &tmp[0], &N, &Z[0], &N, &dZERO, &DD[0], &N);
-
- return diff(A, DD);
-}
-
-#endif /* LAPACK_STEV_HH */