1 ------------------------------------------------------------------------------
3 -- GNAT RUN-TIME COMPONENTS --
5 -- SYSTEM.GENERIC_COMPLEX_BLAS --
9 -- Copyright (C) 2006-2007, 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 2, 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. See the GNU General Public License --
17 -- for more details. You should have received a copy of the GNU General --
18 -- Public License distributed with GNAT; see file COPYING. If not, write --
19 -- to the Free Software Foundation, 51 Franklin Street, Fifth Floor, --
20 -- Boston, MA 02110-1301, USA. --
22 -- As a special exception, if other files instantiate generics from this --
23 -- unit, or you link this unit with other files to produce an executable, --
24 -- this unit does not by itself cause the resulting executable to be --
25 -- covered by the GNU General Public License. This exception does not --
26 -- however invalidate any other reasons why the executable file might be --
27 -- covered by the GNU Public License. --
29 -- GNAT was originally developed by the GNAT team at New York University. --
30 -- Extensive contributions were provided by Ada Core Technologies Inc. --
32 ------------------------------------------------------------------------------
34 -- Package comment required ???
36 with Ada
.Numerics
.Generic_Complex_Types
;
39 type Real
is digits <>;
40 with package Complex_Types
is new Ada
.Numerics
.Generic_Complex_Types
(Real
);
43 type Complex_Vector
is array (Integer range <>) of Complex
;
44 type Complex_Matrix
is array (Integer range <>, Integer range <>)
46 package System
.Generic_Complex_BLAS
is
49 -- Although BLAS support is only available for IEEE single and double
50 -- compatible floating-point types, this unit will accept any type
51 -- and apply conversions as necessary, with possible loss of
52 -- precision and range.
54 No_Trans
: aliased constant Character := 'N';
55 Trans
: aliased constant Character := 'T';
56 Conj_Trans
: aliased constant Character := 'C';
58 -- BLAS Level 1 Subprograms and Types
65 Inc_Y
: Integer := 1) return Complex
;
70 Inc_X
: Integer := 1) return Real
;
73 (Trans
: access constant Character;
76 Alpha
: Complex
:= (1.0, 0.0);
80 Inc_X
: Integer := 1; -- must be non-zero
81 Beta
: Complex
:= (0.0, 0.0);
82 Y
: in out Complex_Vector
;
83 Inc_Y
: Integer := 1); -- must be non-zero
87 -- gemm s, d, c, z Matrix-matrix product of general matrices
90 (Trans_A
: access constant Character;
91 Trans_B
: access constant Character;
95 Alpha
: Complex
:= (1.0, 0.0);
100 Beta
: Complex
:= (0.0, 0.0);
101 C
: in out Complex_Matrix
;
104 end System
.Generic_Complex_BLAS
;