0
Fixed

Bug in relational operations

Michael_ 8 years ago updated by Pavel Holoborodko 8 years ago 3

Hi.

I think there is a bug in the relational operations "equal" and "unequal"


test(:,1)=(1:4)';

testmp=mp(test);


This works:

test==1
ans =
1
0
0
0
testmp==1
ans =
1
0
0
0

This doesn't

1==test
ans =
1
0
0
0
1==testmp
ans =
1


This has worked in previous versions though.

Best regards,

Michael

Under review

Dear Michael,


Thank you for your report, will fix it asap.

(Recently I have added support for sparse matrices in relational operations - probably something went wrong).

Bug has been fixed, updated version will be released in a few hours (after compilation finishes).

All relational operations were affected, and also element-wise power x.^y.

Fixed

Updated version has been released, download link: http://goo.gl/pMXV3


Now I am working on speed optimization of Bessel functions for special cases.

For example new version already has much faster bessely(n,x) when n - integer, x - non-negative real in quadruple precision.


Please try it (if you use it in your work). I am finishing besselj - will be available very soon as well.