Fast

Bug in relational operations

9 år siden • opdateret af 9 år siden • 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

Stem

Under vurdering

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.

Fast

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.