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

["Nasser M. Abbasi" <nma@12000.org>]

On 7/24/2012 10:40 AM, Ada novice wrote:
> In http://12000.org/my_notes/ada/index.htm point 13 where a complete
>example is given to solve the system Ax=b, on running the code I get an error due to a line almost at the bottom of the code:
>
> for x of b loop
>
> The error message is "too few subscripts in array reference".
>
> The b vector is a 3 rows 1 column vector.
>
> YC
>

It works fine here. May be you copied the source from the page wrong.
(I'll make a text file of this later and add it later, for now
you can use the mouse to copy it, it is verbatim text)

--------------------------------
cd lapada/test  -- this is where you'll have the test program sit
                
>gnatmake -gnat2012 -I../ada  mysolve.adb
gcc -c -gnat2012 -I../ada mysolve.adb
gnatbind -I../ada -x mysolve.ali
gnatlink mysolve.ali

>ls -lrt
-rwxrwxrwx 1 me me 552678 Jul 24 11:35 mysolve

>./mysolve
-1.31250E+00
  3.50000E+00
  1.12500E+00

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

Make sure to copy the source correct. Here it is again

>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;
end mySolve;
--------------------------------------                     

--Nasser