Re: Free AMD Core Math Library (BLAS/LAPACK) + Ada

[Ada novice <shai.lesh@gmx.com>]

First, I am using GNAT GPL 2011 and not the 2012 version. But I am adding the -gnat2012 switch as this is required or else I get:

$ gnatmake -I../ada mysolve.adbgcc -c -I../ada mysolve.adb
mysolve.adb:39:16: iterator is an Ada2012 feature
gnatmake: "mysolve.adb" compilation error

Anyway, back to using the -gnat2012 switch:

(my username here)~/work/ada/testlapack/lapada/test $ gnatmake -gnat2012 -I../ada mysolve.adb
gcc -c -gnat2012 -I../ada mysolve.adb
mysolve.adb:39:21: too few subscripts in array reference
gnatmake: "mysolve.adb" compilation error

and

(my username here)~/work/ada/testlapack/lapada/test $ cat mysolve.adb
with Ada.Text_IO;  use  Ada.Text_IO;

with Interfaces.Fortran; use Interfaces.Fortran;
with labase; -- from LAPACK binding
with ladrv;  -- from LAPACK binding

procedure mySolve is

   A:   labase.Fortran_Real_Matrix ( 1..3, 1..3 );
   b:   labase.Fortran_Real_Matrix ( 1..3, 1..1  );
   package Real_IO is new Ada.Text_IO.Float_IO( Real );
   INFO : Fortran_Integer;
   IPIV : labase.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 )
               );

         ladrv.SGESV (N    => A'Last(2),
                      NRHS => b'Last(2),
                      A    => A,
                      LDA  => A'Last(1),
                      IPIV => IPIV,
                      B    => b,
                      LDB  => b'Last(1),
                      INFO => INFO);

           if ( not(INFO = 0) ) then
           raise PROGRAM_ERROR;
        end if;

         for x of b loop  -- print solution
             real_io.PUT ( x ); new_line;
         end loop;


I am not convinced to think that using GNAT GPL 2011 instead of the 2012 version is the culprit here. If however I have to change to the 2012 version, I wonder if I can have both the 2011 and 202 on my machine. On windows it is simple as GNAT is stored in year-wise folders.

I am using the 2011 version since there is no in-built lapack/blas in GNAT 2012 and I want to be sure that I can still run the earlier math extension package from Simon Wright.

On Windows installing separate lapack/blas libraries is a nightmare. But I guess if I am using linux, then I can safely change to GPL 2012 and use the updated math extension package from Simon. Anyway, in this respect what is the way to uninstall GNAT GPL 2011 given that this cannot be done via the synaptic package manager since GNAT GPL 2011 does not come pre-packaged for Debian.

Thanks,
YC