Re: fyi, small update to Ada LAPACK and BLAS binding

[Niklas Holsti <niklas.holsti@tidorum.invalid>]

On 12-07-31 09:37 , Nasser M. Abbasi wrote:

> Here is a complete example using LAPACK from Ada, with the gnatmake
> command to build it:
> 
> gnatmake -Iada_lapack/binding mysolve.adb -largs -L/usr/lib -lblas -llapack
> 
> ------ mysolve.adb -------------------------------------
> with Ada.Text_IO; use Ada.Text_IO;
> 
> with Interfaces.Fortran; use Interfaces.Fortran;
> with lapack;
> 
> procedure mysolve is
> 
>    A    : lapack.Fortran_Real_Matrix (1 .. 3, 1 .. 3);
>    b    : lapack.Fortran_Real_Matrix (1 .. 3, 1 .. 1);
>    info : Fortran_Integer;
>    ipiv : lapack.Fortran_Integer_Vector (1 .. A'Last (2));
> 
> begin -- solve A x=b
> 
>    A := ((2.0, 3.0, 1.0), (2.0, 1.0, 1.0), (4.0, -1.0, 6.0));
>    b := ((1 => 9.0), (1 => 2.0), (1 => -2.0));
> 
>    lapack.SGESV
>      (N    => A'Last (2),

Wouldn't A'Length(2) be more to the point? And the same for the other
uses of 'Last below. Of course, when 'First = 1 then 'Length = 'Last,
but just as a general point of principle...

>       NRHS => b'Last (2),
>       A    => A,
>       LDA  => A'Last (1),
>       IPIV => ipiv,
>       B    => b,
>       LDB  => b'Last (1),
>       INFO => info);


-- 
Niklas Holsti
Tidorum Ltd
niklas holsti tidorum fi
      .      @       .