Re: Interfacing Ada with C

[Simon Wright <simon@pushface.org>]

Simon Wright <simon@pushface.org> writes:

> Ada novice <posts@gmx.us> writes:
>> On Jul 26, 1:21 am, Simon Wright <si...@pushface.org> wrote:
>>
>>> I've encoded a general complex eigenvalues function, interfacing to the
>>> LAPACK procedure zgeev
>
>> Thank you very much for your commendable efforts. This is a very good
>> example to demonstrate how Ada can be binded with LAPACK.
>
> Well, it's a start!

I've taken the plunge and started a SourceForge project for this. It's
at http://sourceforge.net/projects/gnat-math-extn/ -- under "Ada 2005
Math Extensions", click on [Develop] then on [Code].

I've chosen to use Mercurial (Hg) as the VCS, mainly to get a real-world
feel for using a DVCS (Distributed Version Control System). To downoad
the code, you'll need to install Hg - http://mercurial.selenic.com/ -
because I haven't actually made a code release yet!

If anyone feels moved to join in, just say (of course you need a SF
account to update the SF repository, but with Hg it should be possible
to work via patchsets .. )

For interest, the test program now looks like

   with Ada.Text_IO; use Ada.Text_IO;
   with Ada.Numerics.Generic_Real_Arrays;
   with Ada.Numerics.Generic_Complex_Types;
   with Ada.Numerics.Generic_Complex_Arrays.Extensions;

   procedure Test_Extensions is

      package Real_Arrays
      is new Ada.Numerics.Generic_Real_Arrays (Long_Float);
      package Complex_Types
      is new Ada.Numerics.Generic_Complex_Types (Long_Float);
      package Complex_Arrays
      is new Ada.Numerics.Generic_Complex_Arrays (Real_Arrays,
                                                  Complex_Types);
      package Extensions
      is new Complex_Arrays.Extensions;

      use Complex_Arrays;

      Input : Complex_Matrix (1 .. 3, 1 .. 3);
      Result : Complex_Vector (1 .. Input'Length (1));

   begin

      --  Values in yc's example
      Input := (((8.0, 0.0), (-1.0, 0.0), (-5.0, 0.0)),
                ((-4.0, 0.0), (4.0, 0.0), (-2.0, 0.0)),
                ((18.0, 0.0), (-5.0, 0.0), (-7.0, 0.0)));

      Result := Extensions.Eigenvalues (Input);

      for J in Result'Range loop
         Put_Line (Long_Float'Image (Result (J).Re)
                     & " "
                     & Long_Float'Image (Result (J).Im));
      end loop;

   end Test_Extensions;

(results as before)