Re: GNAT vs Matlab - operation on multidimensional complex matrices

[Shark8 <onewingedshark@gmail.com>]

IIRC, one method to ensure proper cleanup w/o using Unchecked Deallocation is to use Generics:

    Generic
        Type T is private;
        Type IndexA is (<>);
        Type IndexB is (<>);
        Type IndexC is (<>);
        Type Cubic is Array(IndexA, IndexB, indexC) of T;
        Default : in T;
    Package Test_Data is
        Type Access_Cubic is not null access all Cubic;
        
        Access_Data : Constant Access_Cubic := new Cubic'(others => (others => (others => Default)));
        Data        : Cubic renames Access_Data.all;
    End Test_Data;

Now you can do something like:

    COMPLEX_CUBE_TEST:
    Declare
        Package Test_Set is new Test_Data(
           T       => Complex,
           IndexA  => Channel,
           IndexB  => Angle,
           IndexC  => Range_T,
           Cubic   => tCubeComplex,
           Default => (Re => 1.0, Im => 1.0)
          );
        Complex_Cube : tCubeComplex renames Test_Set.Data;
        SComplex     : Complex;
    Begin
        TIO.Put_Line("Complex cube");
        TIO.Put_Line("Complex type size is:" & Integer'Image(Complex'Size));
        
        SpeedSumComplexCube(
           Mtx          => Complex_Cube,
           S            => SComplex,
           Iterations   => NumIteration
          );
        
        TIO.Put("Sum is:"); CTIO.Put(SComplex); TIO.New_Line;
    End COMPLEX_CUBE_TEST;

and upon exit of the COMPLEX_CUBE_TEST block your test-data is cleaned up when Test_Set is destroyed.