Subversion Repositories openrisc_me

------------------------------------------------------------------------------

--                                                                          --

--                         GNAT RUN-TIME COMPONENTS                         --

--                                                                          --

--               I N T E R F A C E S . F O R T R A N . B L A S              --

--                                                                          --

--                                 S p e c                                  --

--                                                                          --

--         Copyright (C) 2006-2009, Free Software Foundation, Inc.          --

--                                                                          --

-- GNAT is free software;  you can  redistribute it  and/or modify it under --

-- terms of the  GNU General Public License as published  by the Free Soft- --

-- ware  Foundation;  either version 3,  or (at your option) any later ver- --

-- sion.  GNAT is distributed in the hope that it will be useful, but WITH- --

-- OUT ANY WARRANTY;  without even the  implied warranty of MERCHANTABILITY --

-- or FITNESS FOR A PARTICULAR PURPOSE.                                     --

--                                                                          --

-- As a special exception under Section 7 of GPL version 3, you are granted --

-- additional permissions described in the GCC Runtime Library Exception,   --

-- version 3.1, as published by the Free Software Foundation.               --

--                                                                          --

-- You should have received a copy of the GNU General Public License and    --

-- a copy of the GCC Runtime Library Exception along with this program;     --

-- see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see    --

-- <http://www.gnu.org/licenses/>.                                          --

--                                                                          --

-- GNAT was originally developed  by the GNAT team at  New York University. --

-- Extensive contributions were provided by Ada Core Technologies Inc.      --

--                                                                          --

------------------------------------------------------------------------------

--  This package provides a thin binding to the standard Fortran BLAS library.

--  Documentation and a reference BLAS implementation is available from

--  ftp://ftp.netlib.org. The main purpose of this package is to facilitate

--  implementation of the Ada 2005 Ada.Numerics.Generic_Real_Arrays and

--  Ada.Numerics.Generic_Complex_Arrays packages. Bindings to other BLAS

--  routines may be added over time.

--  As actual linker arguments to link with the BLAS implementation differs

--  according to platform and chosen BLAS implementation, the linker arguments

--  are given in the body of this package. The body may need to be modified in

--  order to link with different BLAS implementations tuned to the specific

--  target.

package Interfaces.Fortran.BLAS is

   pragma Pure;

   pragma Elaborate_Body;

   No_Trans   : aliased constant Character := 'N';

   Trans      : aliased constant Character := 'T';

   Conj_Trans : aliased constant Character := 'C';

   --  Vector types

   type Real_Vector is array (Integer range <>) of Real;

   type Complex_Vector is array (Integer range <>) of Complex;

   type Double_Precision_Vector is array (Integer range <>)

     of Double_Precision;

   type Double_Complex_Vector is array (Integer range <>) of Double_Complex;

   --  Matrix types

   type Real_Matrix is array (Integer range <>, Integer range <>)

     of Real;

   type Double_Precision_Matrix is array (Integer range <>, Integer range <>)

     of Double_Precision;

   type Complex_Matrix is array (Integer range <>, Integer range <>)

     of Complex;

   type Double_Complex_Matrix is array (Integer range <>, Integer range <>)

     of Double_Complex;

   --  BLAS Level 1

   function sdot

     (N     : Positive;

      X     : Real_Vector;

      Inc_X : Integer := 1;

      Y     : Real_Vector;

      Inc_Y : Integer := 1) return Real;

   function ddot

     (N     : Positive;

      X     : Double_Precision_Vector;

      Inc_X : Integer := 1;

      Y     : Double_Precision_Vector;

      Inc_Y : Integer := 1) return Double_Precision;

   function cdotu

     (N     : Positive;

      X     : Complex_Vector;

      Inc_X : Integer := 1;

      Y     : Complex_Vector;

      Inc_Y : Integer := 1) return Complex;

   function zdotu

     (N     : Positive;

      X     : Double_Complex_Vector;

      Inc_X : Integer := 1;

      Y     : Double_Complex_Vector;

      Inc_Y : Integer := 1) return Double_Complex;

   function snrm2

     (N     : Natural;

      X     : Real_Vector;

      Inc_X : Integer := 1) return Real;

   function dnrm2

     (N     : Natural;

      X     : Double_Precision_Vector;

      Inc_X : Integer := 1) return Double_Precision;

   function scnrm2

     (N     : Natural;

      X     : Complex_Vector;

      Inc_X : Integer := 1) return Real;

   function dznrm2

     (N     : Natural;

      X     : Double_Complex_Vector;

      Inc_X : Integer := 1) return Double_Precision;

   --  BLAS Level 2

   procedure sgemv

     (Trans : access constant Character;

      M     : Natural := 0;

      N     : Natural := 0;

      Alpha : Real := 1.0;

      A     : Real_Matrix;

      Ld_A  : Positive;

      X     : Real_Vector;

      Inc_X : Integer := 1;  -- must be non-zero

      Beta  : Real := 0.0;

      Y     : in out Real_Vector;

      Inc_Y : Integer := 1); -- must be non-zero

   procedure dgemv

     (Trans : access constant Character;

      M     : Natural := 0;

      N     : Natural := 0;

      Alpha : Double_Precision := 1.0;

      A     : Double_Precision_Matrix;

      Ld_A  : Positive;

      X     : Double_Precision_Vector;

      Inc_X : Integer := 1;  -- must be non-zero

      Beta  : Double_Precision := 0.0;

      Y     : in out Double_Precision_Vector;

      Inc_Y : Integer := 1); -- must be non-zero

   procedure cgemv

     (Trans : access constant Character;

      M     : Natural := 0;

      N     : Natural := 0;

      Alpha : Complex := (1.0, 1.0);

      A     : Complex_Matrix;

      Ld_A  : Positive;

      X     : Complex_Vector;

      Inc_X : Integer := 1;  -- must be non-zero

      Beta  : Complex := (0.0, 0.0);

      Y     : in out Complex_Vector;

      Inc_Y : Integer := 1); -- must be non-zero

   procedure zgemv

     (Trans : access constant Character;

      M     : Natural := 0;

      N     : Natural := 0;

      Alpha : Double_Complex := (1.0, 1.0);

      A     : Double_Complex_Matrix;

      Ld_A  : Positive;

      X     : Double_Complex_Vector;

      Inc_X : Integer := 1;  -- must be non-zero

      Beta  : Double_Complex := (0.0, 0.0);

      Y     : in out Double_Complex_Vector;

      Inc_Y : Integer := 1); -- must be non-zero

   --  BLAS Level 3

   procedure sgemm

     (Trans_A : access constant Character;

      Trans_B : access constant Character;

      M       : Positive;

      N       : Positive;

      K       : Positive;

      Alpha   : Real := 1.0;

      A       : Real_Matrix;

      Ld_A    : Integer;

      B       : Real_Matrix;

      Ld_B    : Integer;

      Beta    : Real := 0.0;

      C       : in out Real_Matrix;

      Ld_C    : Integer);

   procedure dgemm

     (Trans_A : access constant Character;

      Trans_B : access constant Character;

      M       : Positive;

      N       : Positive;

      K       : Positive;

      Alpha   : Double_Precision := 1.0;

      A       : Double_Precision_Matrix;

      Ld_A    : Integer;

      B       : Double_Precision_Matrix;

      Ld_B    : Integer;

      Beta    : Double_Precision := 0.0;

      C       : in out Double_Precision_Matrix;

      Ld_C    : Integer);

   procedure cgemm

     (Trans_A : access constant Character;

      Trans_B : access constant Character;

      M       : Positive;

      N       : Positive;

      K       : Positive;

      Alpha   : Complex := (1.0, 1.0);

      A       : Complex_Matrix;

      Ld_A    : Integer;

      B       : Complex_Matrix;

      Ld_B    : Integer;

      Beta    : Complex := (0.0, 0.0);

      C       : in out Complex_Matrix;

      Ld_C    : Integer);

   procedure zgemm

     (Trans_A : access constant Character;

      Trans_B : access constant Character;

      M       : Positive;

      N       : Positive;

      K       : Positive;

      Alpha   : Double_Complex := (1.0, 1.0);

      A       : Double_Complex_Matrix;

      Ld_A    : Integer;

      B       : Double_Complex_Matrix;

      Ld_B    : Integer;

      Beta    : Double_Complex := (0.0, 0.0);

      C       : in out Double_Complex_Matrix;

      Ld_C    : Integer);

private

   pragma Import (Fortran, cdotu,  "cdotu_");

   pragma Import (Fortran, cgemm,  "cgemm_");

   pragma Import (Fortran, cgemv,  "cgemv_");

   pragma Import (Fortran, ddot,   "ddot_");

   pragma Import (Fortran, dgemm,  "dgemm_");

   pragma Import (Fortran, dgemv,  "dgemv_");

   pragma Import (Fortran, dnrm2,  "dnrm2_");

   pragma Import (Fortran, dznrm2, "dznrm2_");

   pragma Import (Fortran, scnrm2, "scnrm2_");

   pragma Import (Fortran, sdot,   "sdot_");

   pragma Import (Fortran, sgemm,  "sgemm_");

   pragma Import (Fortran, sgemv,  "sgemv_");

   pragma Import (Fortran, snrm2,  "snrm2_");

   pragma Import (Fortran, zdotu,  "zdotu_");

   pragma Import (Fortran, zgemm,  "zgemm_");

   pragma Import (Fortran, zgemv,  "zgemv_");

end Interfaces.Fortran.BLAS;