Re: Real_Arrays on heap with overloaded operators and clean syntax

[Gautier write-only address <gautier_niouzes@hotmail.com>]

Note that Real_Arrays does not specify where things are allocated (heap or stack).
Only when you define "x : Real_Vector (1 .. n)", it is on stack. You can always write something like the snippet below.
Anyway, after a certain size, you may have to find compromises, like avoiding operators (they do too many allocations & deallocations in the background, even assuming elegant heap-allocated objects) and also give up plain matrices, against sparse matrices or band-stored matrices, typically for solving Partial Differential Equations.

with Ada.Numerics.Generic_Real_Arrays;

procedure Test_Large is

  type Float_15 is digits 15;

  package F15_R_A is new Ada.Numerics.Generic_Real_Arrays (Float_15);

  use F15_R_A;

  procedure Solve_it
    (x : in     Real_Vector;
     y :    out Real_Vector;
     A : in     Real_Matrix) is
  begin
    null;  --  Here, the big number-crunching
  end;
  
  n : constant := 10_000;

  type Vector_Access is access Real_Vector;
  type Matrix_Access is access Real_Matrix;

  x, y : Vector_Access := new Real_Vector (1 .. n);
  A    : Matrix_Access := new Real_Matrix (1 .. n, 1 .. n);
  
begin
  Solve_it (x.all, y.all, A.all);
  --  !! Deallocation here
end;