Re: heap size exceeded for large matrices

[Brian Drummond <brian_drummond@btconnect.com>]

On Sun, 29 Aug 2010 12:27:37 +0200, Pascal Obry <pascal@obry.net> wrote:

>John,
>
>> I wish to process both real and complex large matrices eg (1..10000,
>> 1..10000)
>> I have no difficulty in increasing the storage size of a task.
>> I cannot do the same for matrices within a procedure which reports
>> that the heap is exceeded.
>
>You probably mean stack size and not heap size. In that case you need to
>add the proper linker option to increase the stack size for the application.
>
>Pascal.

I have had trouble with this in the past - the documented linker options
apparently didn't work (around the 2008 timeframe - may not still be true) and
Pragma Storage_Size only worked for new tasks, not the main program..

I resorted to dynamic allocation on the heap, and renaming to hide the pointer
dereference. So, given
my_array_type is array (positive range <>, positive range <>) of float;

I replaced

my_array: my_array_type (1..10000,  1..10000);

with

my_array_ptr : access my_array_type := new my_array_type (1..10000,  1..10000);
my_array : my_array_type renames my_array_ptr.all;

The rename meant I didn't have to rewrite anything else in the application.

- Brian