Interfacing Ada to C for linear algebra

["vincent.diemunsch@gmail.com" <vincent.diemunsch@gmail.com>]

Hello everybody,

I would like to use in Ada some linear algebra libraries written in C,
especially the SuiteSparse :
 http://www.cise.ufl.edu/research/sparse/SuiteSparse/
These libraries handle Sparse Matrices that are in the format used
internally by MatLab : a column-
compressed version of the standard matrices.

Ada defines Generic_Real_Arrays to handle vectors and matrices. So I
would like to simply create a
Generic_Real_Sparse_Arrays package to define MatLab compatible
matrices. Since the internal
format of MatLab matrices is in C, in need to do a conversion between
my Sparse_Matrix and the C
arrays that represent the Sparse Matrix in MatLab.

Looking at the implementation of Generic_Real_Arrays in GNAT, I see
that they internally call the
BLAS and LAPACK packages, written in Fortran, through the use of two
packages :
System.Generic_Real_BLAS and System.Generic_Real_LAPACK. These package
do exactly the
same as what I want to do : they convert Ada arrays of Ada floating
point type in Fortran arrays
of Fortran Float type and back, without using copying, so that the
compiler produces no code
for this conversion (using pragma Inline_Always and so on). This is
important because if we
use optimized versions of linear algebra packages it is because we
manipulate huge arrays
and we don't want to loose time in copying values from on
representation to another !

I suppose that since C is very common, there must already exist a very
good implementation of
conversion of Ada arrays of floating point types in C arrays of
C.Float, idealy using the GNAT
implementation. This could save me a lot of time doing interfaces...

Thank you for your responses,

Vincent