Re: Ann: Mathpaqs, release Feb. 2011

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

On 2/21/2011 10:43 AM, Gautier write-only wrote:
> Hello,
>
> There is a new release of Mathpaqs @ http://sf.net/projects/mathpaqs/
> .
>
> What's new:
>
> - *new* Discrete_random_simulation package (this is for simulating any
> discrete random distribution)
> - updated the Finite_distributed_random function (this is for
> simulating a random distribution for an enumerated type)
> - cleanup of ConjGrad (Conjugate gradient iterative methods for
> solving the matrix equation Ax=b)
>

Thanks Gautier,

I do not use Ada now, but I was looking at your ConjGrad solver,
and was just wondering, why new Matrix and Vectors are
defined in that package, since with Ada now, there exist
such types now?

For example, in  ConjGrad.ads it says

type Vector is array(Index range <>) of Real

type Any_matrix (<>) is private;
   -- NB: 2 syntaxes for instanciating that as unconstrained type :
   -- [Ada 95+] type Any_matrix (<>) is private;
   -- [Ada 83]  type Any_matrix is private;

But when I look at Ada 2005,

http://www.adaic.org/resources/add_content/standards/05rat/html/Rat-7-6.html

"Ada 2005 includes two new packages which are Ada.Numerics.Generic_Real_Arrays and Ada.Numerics.Generic_Complex_Arrays.

type Real_Vector is array (Integer range <>) of Real'Base;
type Real_Matrix is array (Integer range <>, Integer range <>) of Real'Base;
"

So, was just wondering the reasoning, to learn something, that's all.

I wrote a conjugate gradient solver myself for a HW, with preconditioning
in Matlab, and I saw your post, so wanted to see how it would 'look' in Ada.

thanks
--Nasser