From mboxrd@z Thu Jan  1 00:00:00 1970
X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on polar.synack.me
X-Spam-Level: 
X-Spam-Status: No, score=-1.3 required=5.0 tests=BAYES_00,INVALID_MSGID
	autolearn=no autolearn_force=no version=3.4.4
X-Google-Language: ENGLISH,ASCII-7-bit
X-Google-Thread: 103376,85917c9414146b3f
X-Google-Attributes: gid103376,public
From: Gautier <gautier.demontmollin@maths.unine.ch>
Subject: Re: Linear Algebra package for Ada ?
Date: 2000/04/25
Message-ID: <390559FA.B256DE5B@maths.unine.ch>#1/1
X-Deja-AN: 615386295
Content-Transfer-Encoding: 7bit
References: <8dv75n$qi1$1@news.uit.no>
To: Reinert.Korsnes@npolar.no
X-Accept-Language: en
Content-Type: text/plain; charset=us-ascii
X-Trace: 25 Apr 2000 10:40:26 +0100, mac13-32.unine.ch
Organization: Maths - Uni =?iso-8859-1?Q?Neuch=E2tel?=
MIME-Version: 1.0
Newsgroups: comp.lang.ada
Date: 2000-04-25T00:00:00+00:00
List-Id: <comp.lang.ada>

Reinert Korsnes wrote:

> Could anybody recommend me an Ada95 package for linear algebra
> (including manipulation of vectors/matrices ?)

There is a binding Lapack-Ada (Lapada ?); if you are interested
I have here

  - a generic package for band matrices (just extracted from a
    current program, see PS) and some bindings to LU, Cholesky, solvers
    (DGBTRF etc.)

  - a generic package for sparse matrices in CRS format with symmetric
    and asymmetric iterative solvers (stabilized conjuguate gradient) with
    an optional binding to BLAS.

BTW: a true, reference, Lapack in Ada (with generics for hitting
single/double/etc. in one move would be _very_ nice!
A thesis project maybe ?...)

_____________________________________________
Gautier  --  http://members.xoom.com/gdemont/


generic
  type real is digits <>;
  type index is range <>;
  type vector is array(index range <>) of real;
  type matrix is array(index range <>, index range <>) of real;

package BandMatr is

  --------------------------
  -- Define a band matrix --
  --------------------------

   type Band_matrix( max_size,    -- size of actual matrix
                     band,        -- bandwidth
                     multi_band   -- n-uple of bandwidth
                   : index ) is
   record
     val: matrix(1..max_size, 1..multi_band); -- transposed w/r Lapack due to Fortran 
     symmetric: boolean;
     size: index:= max_size;
   end record;
(...)