topic LAPACKE_zgesv Array Limitations in Intel® oneAPI Math Kernel Library

LAPACKE_zgesv Array Limitations

Gianluca_G_ — Wed, 21 Oct 2015 10:26:55 GMT

I'm using LAPACKE_zgesv (MKL 11.3.01) in c# with arrays of Complex.
I have problems when the arrays exceeded about 12000x12000 elements (2.145GB)

.NET Have limitations but I overcomed it with <gcAllowVeryLargeObjects enabled="true" /> directive.

By the way it seems that the problem is with this function LAPACKE_zgesv().
this is the message: "The array size exceed the adresses limits"

Is there any array size limitation?

Thank you very much

Gianluca

There is not such limitation.

Gennady_F_Intel — Wed, 21 Oct 2015 10:45:42 GMT

There is not such limitation. That's the problem with .NET. You may check how it works with the C example.

ok, thankyou I will check

Gianluca_G_ — Wed, 21 Oct 2015 10:48:45 GMT

ok, thankyou I will check further

Gianluca_G_ — Mon, 26 Oct 2015 09:40:00 GMT

Dear Gennady,

the problem seems to be with LAPACKE_zgesv() or at least something linked with it.

I have a Matrix Complex[12145,12145] and a vector of Complex[12145 ] if I apply the LU decomposition (implemented by us) to resolve that system it works. It takes a lot of time but it works.

If we pass these parameters to LAPACKE_zgesv function it doesn't work. Maybe it depends of C# declaration, but with smaller System it works.

I'm using a notebook with windows 10 and 10GB of RAM memory.

I attach the simple example I used to test it.

If you want I have also the files with data but are very large (http://we.tl/ypYNEAk2ZU).

Thank you for any help you can give

Gianluca