0
Fixed

error LU function (the devil is in the details)

john doe 5 years ago updated by Pavel Holoborodko 5 years ago 7
%% test_lu
clear all %#ok<*CLSCR>
clc
mem = memory %#ok<*NASGU,*NOPTS>

tic, [l,u,p] = lu(9), toc %#ok<*ASGLU>
% now do this...
tic, [l,u,p] = lu(mp(9)), toc
... and look at differences ...
mem = memory
% ... wait for while ( for me about 10..40 sec, randomly !)

% ...
% and then I get out of memory error window

... here is my output ...

mem =

MaxPossibleArrayBytes: 3.3001e+09
MemAvailableAllArrays: 3.3001e+09
MemUsedMATLAB: 735469568


l =

1


u =

9


p =

1

Elapsed time is 0.024654 seconds.

l =

{3x1 cell}
{3x1 cell}
[ 0]


u =

{3x1 cell}
{3x1 cell}
[ 0]


p =

0

Elapsed time is 2.856406 seconds.

mem =

MaxPossibleArrayBytes: 3.2757e+09
MemAvailableAllArrays: 3.2757e+09
MemUsedMATLAB: 760942592

>> mp.Info
-------------------------------------------------------------------------------------------------------------
Multiprecision Computing Toolbox, (c) 2008-2015 Advanpix LLC.
Version : 3.8.4 Build 8915
Platform: 64-bit (Win64)
Release : 2015-07-22

P.S.
I reduced the virtual machine's memory to 1G and get much faster now out of memory

P.S. #2
L = lu(mp(9)) % Matlab crash
Under review
Yes, this is hilarious bug! Fixing it now. Thank you as always.

We can offer you nice discount for such great help.
Contact me by e-mail (pavel@advanpix.com) if you are interested.
;)

I used mail address :)

Thanks !
Discount is sent to your mail. Thank you, John Doe :)!
Fixed version has been released. Please update your installation.
[...] = lu(A,'vector') still not work; please take a look at my workaround (with 'vector' the third output MUST be VECTOR)
Certainly, my code is valid only for full matrices.
+1
I have just added that option to LU function (it wasn't implemented before). Please use newest toolbox.