Epetra Package Browser (Single Doxygen Collection)  Development
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Epetra_LAPACK.h
Go to the documentation of this file.
1 /*
2 //@HEADER
3 // ************************************************************************
4 //
5 // Epetra: Linear Algebra Services Package
6 // Copyright 2011 Sandia Corporation
7 //
8 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9 // the U.S. Government retains certain rights in this software.
10 //
11 // Redistribution and use in source and binary forms, with or without
12 // modification, are permitted provided that the following conditions are
13 // met:
14 //
15 // 1. Redistributions of source code must retain the above copyright
16 // notice, this list of conditions and the following disclaimer.
17 //
18 // 2. Redistributions in binary form must reproduce the above copyright
19 // notice, this list of conditions and the following disclaimer in the
20 // documentation and/or other materials provided with the distribution.
21 //
22 // 3. Neither the name of the Corporation nor the names of the
23 // contributors may be used to endorse or promote products derived from
24 // this software without specific prior written permission.
25 //
26 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37 //
38 // Questions? Contact Michael A. Heroux (maherou@sandia.gov)
39 //
40 // ************************************************************************
41 //@HEADER
42 */
43 
44 #ifndef EPETRA_LAPACK_H
45 #define EPETRA_LAPACK_H
46 
47 #if defined(Epetra_SHOW_DEPRECATED_WARNINGS)
48 #ifdef __GNUC__
49 #warning "The Epetra package is deprecated"
50 #endif
51 #endif
52 
53 
54 
56 
72 #include "Epetra_ConfigDefs.h"
73 #include "Epetra_Object.h"
74 
75 class EPETRA_LIB_DLL_EXPORT Epetra_LAPACK {
76 
77  public:
79 
80 
83  Epetra_LAPACK(void);
84 
85 
87 
89  Epetra_LAPACK(const Epetra_LAPACK& LAPACK);
90 
92  virtual ~Epetra_LAPACK(void);
94 
95 
97 
98 
100  void POTRF( const char UPLO, const int N, float * A, const int LDA, int * INFO) const;
102  void POTRF( const char UPLO, const int N, double * A, const int LDA, int * INFO) const;
103 
105  void POTRS( const char UPLO, const int N, const int NRHS, const float * A, const int LDA, float * X, const int LDX, int * INFO) const;
107  void POTRS( const char UPLO, const int N, const int NRHS, const double * A, const int LDA, double * X, const int LDX, int * INFO) const;
108 
110  void POTRI( const char UPLO, const int N, float * A, const int LDA, int * INFO) const;
112  void POTRI( const char UPLO, const int N, double * A, const int LDA, int * INFO) const;
113 
115  void POCON( const char UPLO, const int N, const float * A, const int LDA, const float ANORM,
116  float * RCOND, float * WORK, int * IWORK, int * INFO) const;
118  void POCON( const char UPLO, const int N, const double * A, const int LDA, const double ANORM,
119  double * RCOND, double * WORK, int * IWORK, int * INFO) const;
120 
122  void POSV( const char UPLO, const int N, const int NRHS, float * A, const int LDA, float * X, const int LDX, int * INFO) const;
124  void POSV( const char UPLO, const int N, const int NRHS, double * A, const int LDA, double * X, const int LDX, int * INFO) const;
125 
127  void POEQU(const int N, const float * A, const int LDA, float * S, float * SCOND, float * AMAX, int * INFO) const;
129  void POEQU(const int N, const double * A, const int LDA, double * S, double * SCOND, double * AMAX, int * INFO) const;
130 
132  void PORFS(const char UPLO, const int N, const int NRHS, const float * A, const int LDA, const float * AF, const int LDAF,
133  const float * B, const int LDB, float * X, const int LDX,
134  float * FERR, float * BERR, float * WORK, int * IWORK, int * INFO) const;
136  void PORFS(const char UPLO, const int N, const int NRHS, const double * A, const int LDA, const double * AF, const int LDAF,
137  const double * B, const int LDB, double * X, const int LDX,
138  double * FERR, double * BERR, double * WORK, int * IWORK, int * INFO) const;
139 
141  void POSVX(const char FACT, const char UPLO, const int N, const int NRHS, float * A, const int LDA, float * AF, const int LDAF,
142  const char EQUED, float * S, float * B, const int LDB, float * X, const int LDX, float * RCOND,
143  float * FERR, float * BERR, float * WORK, int * IWORK, int * INFO) const;
145  void POSVX(const char FACT, const char UPLO, const int N, const int NRHS, double * A, const int LDA, double * AF, const int LDAF,
146  const char EQUED, double * S, double * B, const int LDB, double * X, const int LDX, double * RCOND,
147  double * FERR, double * BERR, double * WORK, int * IWORK, int * INFO) const;
149 
151 
152 
154  void GELS( const char TRANS, const int M, const int N, const int NRHS, double* A, const int LDA,
155  double* B, const int LDB, double* WORK, const int LWORK, int * INFO) const;
157  void GETRF( const int M, const int N, float * A, const int LDA, int * IPIV, int * INFO) const;
159  void GETRF( const int M, const int N, double * A, const int LDA, int * IPIV, int * INFO) const;
160 
162  void GEQRF( const int M, const int N, float * A, const int LDA, float * TAU, float * WORK, const int lwork, int * INFO) const;
164  void GEQRF( const int M, const int N, double * A, const int LDA, double * TAU, double * WORK, const int lwork, int * INFO) const;
165 
167  void GETRS( const char TRANS, const int N, const int NRHS, const float * A, const int LDA, const int * IPIV, float * X, const int LDX, int * INFO) const;
169  void GETRS( const char TRANS, const int N, const int NRHS, const double * A, const int LDA, const int * IPIV, double * X, const int LDX, int * INFO) const;
170 
172  void GETRI( const int N, float * A, const int LDA, int * IPIV, float * WORK, const int * LWORK, int * INFO) const;
174  void GETRI( const int N, double * A, const int LDA, int * IPIV, double * WORK, const int * LWORK, int * INFO) const;
175 
177  void GECON( const char NORM, const int N, const float * A, const int LDA, const float ANORM,
178  float * RCOND, float * WORK, int * IWORK, int * INFO) const;
180  void GECON( const char NORM, const int N, const double * A, const int LDA, const double ANORM,
181  double * RCOND, double * WORK, int * IWORK, int * INFO) const;
182 
184  void GESV( const int N, const int NRHS, float * A, const int LDA, int * IPIV, float * X, const int LDX, int * INFO) const;
186  void GESV( const int N, const int NRHS, double * A, const int LDA, int * IPIV, double * X, const int LDX, int * INFO) const;
187 
189  void GEEQU(const int M, const int N, const float * A, const int LDA, float * R, float * C, float * ROWCND, float * COLCND, float * AMAX, int * INFO) const;
191  void GEEQU(const int M, const int N, const double * A, const int LDA, double * R, double * C, double * ROWCND, double * COLCND, double * AMAX, int * INFO) const;
192 
194  void GERFS(const char TRANS, const int N, const int NRHS, const float * A, const int LDA, const float * AF, const int LDAF,
195  const int * IPIV, const float * B, const int LDB, float * X, const int LDX,
196  float * FERR, float * BERR, float * WORK, int * IWORK, int * INFO) const;
198  void GERFS(const char TRANS, const int N, const int NRHS, const double * A, const int LDA, const double * AF, const int LDAF,
199  const int * IPIV, const double * B, const int LDB, double * X, const int LDX,
200  double * FERR, double * BERR, double * WORK, int * IWORK, int * INFO) const;
201 
203  void GESVX(const char FACT, const char TRANS, const int N, const int NRHS, float * A, const int LDA, float * AF, const int LDAF, int * IPIV,
204  const char EQUED, float * R, float * C, float * B, const int LDB, float * X, const int LDX, float * RCOND,
205  float * FERR, float * BERR, float * WORK, int * IWORK, int * INFO) const;
207  void GESVX(const char FACT, const char TRANS, const int N, const int NRHS, double * A, const int LDA, double * AF, const int LDAF, int * IPIV,
208  const char EQUED, double * R, double * C, double * B, const int LDB, double * X, const int LDX, double * RCOND,
209  double * FERR, double * BERR, double * WORK, int * IWORK, int * INFO) const;
210 
211 
213  void GEHRD(const int N, const int ILO, const int IHI, float * A, const int LDA, float * TAU, float * WORK, const int LWORK, int * INFO) const;
215  void GEHRD(const int N, const int ILO, const int IHI, double * A, const int LDA, double * TAU, double * WORK, const int LWORK, int * INFO) const;
217 
219 
220  void HSEQR( const char JOB, const char COMPZ, const int N, const int ILO, const int IHI, float * H, const int LDH, float * WR, float * WI,
222  float * Z, const int LDZ, float * WORK, const int LWORK, int * INFO) const;
224  void HSEQR( const char JOB, const char COMPZ, const int N, const int ILO, const int IHI, double * H, const int LDH, double * WR, double * WI,
225  double * Z, const int LDZ, double * WORK, const int LWORK, int * INFO) const;
227 
229 
230  void ORGQR( const int M, const int N, const int K, float * A, const int LDA, float * TAU, float * WORK, const int LWORK, int * INFO) const;
233  void ORGQR( const int M, const int N, const int K, double * A, const int LDA, double * TAU, double * WORK, const int LWORK, int * INFO) const;
234 
236  void ORGHR( const int N, const int ILO, const int IHI, float * A, const int LDA, float * TAU, float * WORK, const int LWORK, int * INFO) const;
238  void ORGHR( const int N, const int ILO, const int IHI, double * A, const int LDA, double * TAU, double * WORK, const int LWORK, int * INFO) const;
239 
241  void ORMHR( const char SIDE, const char TRANS, const int M, const int N, const int ILO, const int IHI, const float * A, const int LDA,
242  const float * TAU, float * C,
243  const int LDC, float * WORK, const int LWORK, int * INFO) const;
245  void ORMHR( const char SIDE, const char TRANS, const int M, const int N, const int ILO, const int IHI, const double * A, const int LDA,
246  const double * TAU, double * C,
247  const int LDC, double * WORK, const int LWORK, int * INFO) const;
249  void LARFT( const char DIRECT, const char STOREV, const int N, const int K, double * V, const int LDV, double * TAU, double * T, const int LDT) const;
251  void LARFT( const char DIRECT, const char STOREV, const int N, const int K, float * V, const int LDV, float * TAU, float * T, const int LDT) const;
253 
255 
256 
258 
260  void TREVC( const char SIDE, const char HOWMNY, int * SELECT, const int N, const float * T, const int LDT, float *VL, const int LDVL,
261  float * VR, const int LDVR, const int MM, int * M, float * WORK, int * INFO) const;
263 
265  void TREVC( const char SIDE, const char HOWMNY, int * SELECT, const int N, const double * T, const int LDT, double *VL, const int LDVL,
266  double * VR, const int LDVR, const int MM, int *M, double * WORK, int * INFO) const;
267 
269  void TREXC( const char COMPQ, const int N, float * T, const int LDT, float * Q, const int LDQ, int IFST, int ILST,
270  float * WORK, int * INFO) const;
272  void TREXC( const char COMPQ, const int N, double * T, const int LDT, double * Q, const int LDQ, int IFST, int ILST,
273  double * WORK, int * INFO) const;
275 
277 
278 
280  void GESVD( const char JOBU, const char JOBVT, const int M, const int N, float * A, const int LDA, float * S, float * U,
281  const int LDU, float * VT, const int LDVT, float * WORK, const int * LWORK, int * INFO) const;
283  void GESVD( const char JOBU, const char JOBVT, const int M, const int N, double * A, const int LDA, double * S, double * U,
284  const int LDU, double * VT, const int LDVT, double * WORK, const int * LWORK, int * INFO) const;
285 
287  void GGSVD(const char JOBU, const char JOBV, const char JOBQ, const int M, const int N, const int P, int * K, int * L, double* A, const int LDA, double* B, const int LDB,
288  double* ALPHA, double* BETA, double* U, const int LDU, double* V, const int LDV, double* Q, const int LDQ, double* WORK,
289  #ifdef HAVE_EPETRA_LAPACK_GSSVD3
290  const int LWORK,
291  #endif
292  int* IWORK, int* INFO) const;
294  void GGSVD(const char JOBU, const char JOBV, const char JOBQ, const int M, const int N, const int P, int * K, int * L, float* A, const int LDA, float* B, const int LDB,
295  float* ALPHA, float* BETA, float* U, const int LDU, float* V, const int LDV, float* Q, const int LDQ, float* WORK,
296  #ifdef HAVE_EPETRA_LAPACK_GSSVD3
297  const int LWORK,
298  #endif
299  int* IWORK, int* INFO) const;
301 
303 
304  void GEEV(const char JOBVL, const char JOBVR, const int N, double* A, const int LDA, double* WR, double* WI,
306  double* VL, const int LDVL, double* VR, const int LDVR, double* WORK, const int LWORK, int* INFO) const;
308  void GEEV(const char JOBVL, const char JOBVR, const int N, float* A, const int LDA, float* WR, float* WI,
309  float* VL, const int LDVL, float* VR, const int LDVR, float* WORK, const int LWORK, int* INFO) const;
310 
312  void SPEV(const char JOBZ, const char UPLO, const int N, double* AP, double* W, double* Z, int LDZ, double* WORK, int* INFO) const;
314  void SPEV(const char JOBZ, const char UPLO, const int N, float* AP, float* W, float* Z, int LDZ, float* WORK, int* INFO) const;
315 
317  void SPGV(const int ITYPE, const char JOBZ, const char UPLO, const int N, double* AP, double* BP, double* W, double* Z, const int LDZ, double* WORK, int* INFO) const;
319  void SPGV(const int ITYPE, const char JOBZ, const char UPLO, const int N, float* AP, float* BP, float* W, float* Z, const int LDZ, float* WORK, int* INFO) const;
320 
322  void SYEV(const char JOBZ, const char UPLO, const int N, double* A, const int LDA, double* W, double* WORK, const int LWORK, int* INFO) const;
324  void SYEV(const char JOBZ, const char UPLO, const int N, float* A, const int LDA, float* W, float* WORK, const int LWORK, int* INFO) const;
325 
327  void SYEVD(const char JOBZ, const char UPLO, const int N, double* A, const int LDA, double* W,
328  double* WORK, const int LWORK, int* IWORK, const int LIWORK, int* INFO) const;
330  void SYEVD(const char JOBZ, const char UPLO, const int N, float* A, const int LDA, float* W,
331  float* WORK, const int LWORK, int* IWORK, const int LIWORK, int* INFO) const;
332 
334  void SYEVX(const char JOBZ, const char RANGE, const char UPLO, const int N, double* A, const int LDA,
335  const double* VL, const double* VU, const int* IL, const int* IU,
336  const double ABSTOL, int * M, double* W, double* Z, const int LDZ, double* WORK,
337  const int LWORK, int* IWORK, int* IFAIL,
338  int* INFO) const;
340  void SYEVX(const char JOBZ, const char RANGE, const char UPLO, const int N, float* A, const int LDA,
341  const float* VL, const float* VU, const int* IL, const int* IU,
342  const float ABSTOL, int * M, float* W, float* Z, const int LDZ, float* WORK,
343  const int LWORK, int* IWORK, int* IFAIL,
344  int* INFO) const;
345 
347  void SYGV(const int ITYPE, const char JOBZ, const char UPLO, const int N, double* A, const int LDA, double* B,
348  const int LDB, double* W, double* WORK, const int LWORK, int* INFO) const;
350  void SYGV(const int ITYPE, const char JOBZ, const char UPLO, const int N, float* A, const int LDA, float* B,
351  const int LDB, float* W, float* WORK, const int LWORK, int* INFO) const;
352 
354  void SYGVX(const int ITYPE, const char JOBZ, const char RANGE, const char UPLO, const int N,
355  double* A, const int LDA, double* B, const int LDB, const double* VL, const double* VU,
356  const int* IL, const int* IU, const double ABSTOL, int* M, double* W, double* Z,
357  const int LDZ, double* WORK, const int LWORK, int* IWORK,
358  int* IFAIL, int* INFO) const;
360  void SYGVX(const int ITYPE, const char JOBZ, const char RANGE, const char UPLO, const int N,
361  float* A, const int LDA, float* B, const int LDB, const float* VL, const float* VU,
362  const int* IL, const int* IU, const float ABSTOL, int* M, float* W, float* Z,
363  const int LDZ, float* WORK, const int LWORK, int* IWORK,
364  int* IFAIL, int* INFO) const;
365 
367  void SYEVR(const char JOBZ, const char RANGE, const char UPLO, const int N, double* A, const int LDA, const double* VL, const double* VU, const int *IL, const int *IU,
368  const double ABSTOL, int* M, double* W, double* Z, const int LDZ, int* ISUPPZ, double* WORK, const int LWORK, int* IWORK,
369  const int LIWORK, int* INFO) const;
371  void SYEVR(const char JOBZ, const char RANGE, const char UPLO, const int N, float* A, const int LDA,
372  const float* VL, const float* VU, const int *IL, const int *IU,
373  const float ABSTOL, int* M, float* W, float* Z, const int LDZ, int* ISUPPZ,
374  float* WORK, const int LWORK, int* IWORK,
375  const int LIWORK, int* INFO) const;
376 
378  void GEEVX(const char BALANC, const char JOBVL, const char JOBVR, const char SENSE, const int N, double* A, const int LDA, double* WR, double* WI, double* VL,
379  const int LDVL, double* VR, const int LDVR, int* ILO, int* IHI, double* SCALE, double* ABNRM, double* RCONDE,
380  double* RCONDV, double* WORK, const int LWORK, int* IWORK, int* INFO) const;
382  void GEEVX(const char BALANC, const char JOBVL, const char JOBVR, const char SENSE, const int N, float* A, const int LDA, float* WR, float* WI, float* VL,
383  const int LDVL, float* VR, const int LDVR, int* ILO, int* IHI, float* SCALE, float* ABNRM, float* RCONDE,
384  float* RCONDV, float* WORK, const int LWORK, int* IWORK, int* INFO) const;
385 
387  void GESDD(const char JOBZ, const int M, const int N, double* A, const int LDA, double* S, double* U, const int LDU, double* VT, const int LDVT, double* WORK,
388  const int LWORK, int* IWORK, int* INFO) const;
390  void GESDD(const char JOBZ, const int M, const int N, float* A, const int LDA, float* S, float* U, const int LDU, float* VT, const int LDVT, float* WORK,
391  const int LWORK, int* IWORK, int* INFO) const;
393 
394  void GGEV(const char JOBVL, const char JOBVR, const int N, double* A, const int LDA, double* B, const int LDB, double* ALPHAR, double* ALPHAI,
395  double* BETA, double* VL, const int LDVL, double* VR, const int LDVR, double* WORK, const int LWORK, int* INFO) const;
397  void GGEV(const char JOBVL, const char JOBVR, const int N, float* A, const int LDA, float* B, const int LDB, float* ALPHAR, float* ALPHAI,
398  float* BETA, float* VL, const int LDVL, float* VR, const int LDVR, float* WORK, const int LWORK, int* INFO) const;
399 
401 
403 
404  void GGLSE(const int M, const int N, const int P, double* A, const int LDA, double* B, const int LDB,
406  double* C, double* D, double* X, double* WORK, const int LWORK, int* INFO) const;
408  void GGLSE(const int M, const int N, const int P, float* A, const int LDA, float* B, const int LDB,
409  float* C, float* D, float* X, float* WORK, const int LWORK, int* INFO) const;
411 
413 
414  void LAMCH ( const char CMACH, float & T) const;
417  void LAMCH ( const char CMACH, double & T) const;
419 
421 
423 
424  void TRTRS(const char UPLO, const char TRANS, const char DIAG, const int N, const int NRHS, const float *A,
426  const int LDA, float *B, const int LDB, int *INFO) const;
428  void TRTRS(const char UPLO, const char TRANS, const char DIAG, const int N, const int NRHS, const double *A,
429  const int LDA, double *B, const int LDB, int *INFO) const;
430 };
431 
432 // Epetra_LAPACK constructor
434 // Epetra_LAPACK constructor
435 inline Epetra_LAPACK::Epetra_LAPACK(const Epetra_LAPACK& LAPACK){(void)LAPACK;}
436 // Epetra_LAPACK destructor
438 
439 #endif /* EPETRA_LAPACK_H */
Epetra_LAPACK: The Epetra LAPACK Wrapper Class.
Definition: Epetra_LAPACK.h:75
virtual ~Epetra_LAPACK(void)
Epetra_LAPACK Destructor.
Epetra_LAPACK(void)
Epetra_LAPACK Constructor.