Question: A question about the memory requirement of generalized eigenvalue solver in Maple 10 and further versions

Dear everyone,

  I have a question about the memory requirement of generalized eigenvalue solver in Maple 10 and further versions.

Suppose I have a code like

***

with(LinearAlgebra):
Digits:=34;
n:=30;
A:=Matrix(n,n);
B:=Matrix(n,n);

for i from 1 to n do
  for j from 1 to n do
     if i=j then
       A[i,j]:=i*j;
       B[i,j]:=1;
     else
       A[i,j]:=0.01*i*j;
       B[i,j]:=0.001*i*j;
     end if;
  end do;
end do;

x,y:=Eigenvectors(A,B,output = [`vectors`, `values`]);
print(y[1]);
***

In Maple 14, the memory increases from 1M, 7M, 193M, 366M, to 681M for n=10,20,30,40, and 50. For larger n, it increases rapidly.

In Maple 10, the memory increases from 4M, 5M, 5M, 5M, to 6M for n=10,20,30,40,and 50.

I have some trial version of Maple 15, the memory cost goes like Maple 14.

How should I do to have similar memory requirement like Maple 10? Does this issue still similar in Maple 16?

Thank you very much in advance

with best regards