numpi / hm-toolbox

Toolbox for HODLR and HSS matrices in MATLAB
GNU General Public License v3.0
42 stars 8 forks source link

Issue with mldivide and 0-rank blocks #4

Closed bonevbs closed 4 years ago

bonevbs commented 4 years ago

Hi again,

I suspect there is an issue with 0-rank blocks and hss_ulv_solve. I get

Index in position 1 is invalid. Array indices must be positive integers or logical values.

Error in hss_ulv_solve>hss_ulv_solve_rec (line 108)
    U = U(end-k+1:end,:);

when I try to use mldivide with the matrix attached in the file. It would be great if you guys could have a look!

Thanks a lot, Boris inversion_problem.mat.zip