1 ------------------------------------------------------------------------------
3 -- GNAT RUN-TIME COMPONENTS --
5 -- S Y S T E M . G E N E R I C _ C O M P L E X _ L A P A C K --
9 -- Copyright (C) 2006-2009, Free Software Foundation, Inc. --
11 -- GNAT is free software; you can redistribute it and/or modify it under --
12 -- terms of the GNU General Public License as published by the Free Soft- --
13 -- ware Foundation; either version 3, or (at your option) any later ver- --
14 -- sion. GNAT is distributed in the hope that it will be useful, but WITH- --
15 -- OUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY --
16 -- or FITNESS FOR A PARTICULAR PURPOSE. --
18 -- As a special exception under Section 7 of GPL version 3, you are granted --
19 -- additional permissions described in the GCC Runtime Library Exception, --
20 -- version 3.1, as published by the Free Software Foundation. --
22 -- You should have received a copy of the GNU General Public License and --
23 -- a copy of the GCC Runtime Library Exception along with this program; --
24 -- see the files COPYING3 and COPYING.RUNTIME respectively. If not, see --
25 -- <http://www.gnu.org/licenses/>. --
27 -- GNAT was originally developed by the GNAT team at New York University. --
28 -- Extensive contributions were provided by Ada Core Technologies Inc. --
30 ------------------------------------------------------------------------------
32 -- Package comment required ???
34 with Ada
.Numerics
.Generic_Complex_Types
;
36 type Real
is digits <>;
37 type Real_Vector
is array (Integer range <>) of Real
;
39 with package Complex_Types
is new Ada
.Numerics
.Generic_Complex_Types
(Real
);
42 type Complex_Vector
is array (Integer range <>) of Complex
;
43 type Complex_Matrix
is array (Integer range <>, Integer range <>)
45 package System
.Generic_Complex_LAPACK
is
48 type Integer_Vector
is array (Integer range <>) of Integer;
50 Upper
: aliased constant Character := 'U';
51 Lower
: aliased constant Character := 'L';
53 -- LAPACK Computational Routines
55 -- getrf computes LU factorization of a general m-by-n matrix
60 A
: in out Complex_Matrix
;
62 I_Piv
: out Integer_Vector
;
63 Info
: access Integer);
65 -- getri computes inverse of an LU-factored square matrix,
66 -- with multiple right-hand sides
70 A
: in out Complex_Matrix
;
72 I_Piv
: Integer_Vector
;
73 Work
: in out Complex_Vector
;
75 Info
: access Integer);
77 -- getrs solves a system of linear equations with an LU-factored
78 -- square matrix, with multiple right-hand sides
81 (Trans
: access constant Character;
86 I_Piv
: Integer_Vector
;
87 B
: in out Complex_Matrix
;
89 Info
: access Integer);
91 -- heevr computes selected eigenvalues and, optionally,
92 -- eigenvectors of a Hermitian matrix using the Relatively
93 -- Robust Representations
96 (Job_Z
: access constant Character;
97 Rng
: access constant Character;
98 Uplo
: access constant Character;
100 A
: in out Complex_Matrix
;
102 Vl
, Vu
: Real
:= 0.0;
103 Il
, Iu
: Integer := 1;
104 Abs_Tol
: Real
:= 0.0;
107 Z
: out Complex_Matrix
;
109 I_Supp_Z
: out Integer_Vector
;
110 Work
: out Complex_Vector
;
112 R_Work
: out Real_Vector
;
114 I_Work
: out Integer_Vector
;
116 Info
: access Integer);
118 -- steqr computes all eigenvalues and eigenvectors of a symmetric or
119 -- Hermitian matrix reduced to tridiagonal form (QR algorithm)
122 (Comp_Z
: access constant Character;
124 D
: in out Real_Vector
;
125 E
: in out Real_Vector
;
126 Z
: in out Complex_Matrix
;
128 Work
: out Real_Vector
;
129 Info
: access Integer);
131 end System
.Generic_Complex_LAPACK
;