Your comments
Dear Ginette,
First of all, I am very happy to here from you!
Hope you are doing fine and everything is good.
As for octuple precision, if we setup mp.Digits(71), toolbox will provide computations with ~99% compatibility with the standard. The only major difference would be exponent size - IEEE requires 19 bits, whereas toolbox uses 32 bits (for better compatibility with CPU registers).
Basically, if numbers stay inside IEEE exponent range (=most common case) - than computations with toolbox will be fully compatible with the standard. And also reproducible on different platforms.
In general, toolbox follows IEEE guidelines and requirements in implementation of arbitrary precision computations.
***
Maybe I am missing something here, do you have some particular compatibility test in mind?
Well, mp.Test() includes only the very basic tests - just to make sure toolbox is installed properly.
We run full set of quality tests (with corner cases, etc.) in-house before every release. Usually it takes 3-6 hours.....
I've made the question public, as it might be helpful for the others.
Hello!
Associated Legendre functions? Yes, of course we have them, but without the 'mp' prefix:
>> X = rand(2,4,5,'mp');
>> n = 2;
>> P = legendre(n,X) >> P(:,1,2,3)
ans =
0.441388679511820695786134758975951
-1.450338406517212223598006351398264
1.117222640976358608427730482048098
Also toolbox provides 'legendreP' - function to compute Legendre polynomials. See Function Reference for the complete list of the functions implemented in toolbox.
It is not included in toolbox yet, but it is in our ToDo list with medium priority.
You might speed things up if you share some details on the problem you solve and why it is important to have 'mp' version of it.
Thank you,
Hi Stefan,
Thank you for using the toolbox!
This bug was fixed long time ago - please use newer version of toolbox.
One year has passed - 3.9.4 is an ancient fossil by now!!!
There is another reason to switch to newest version of toolbox - now it has full-featured EIGS with the support of all modes ('SM', etc.) and variants (including when matrix is supplied as operator). Generalized problems are supported as well.
The only (temporary) downside - now EIGS is written in M-language and quite slow. We are working on this and it will be much faster soon.
I highly recommend you to test new version.
Hi Denis,
To get machine epsilon for particular number 'x' please use:
before the 4.3.3.12177:
>> mp.eps(x)
starting from 4.3.3.12177, the prefix 'mp.' is no longer needed:
>> eps(x)
The 'rank(A,tol)' has been included in recent versions of toolbox as well.
I highly suggest you to use the newest version of toolbox - 4.3.3.12177.
The 'accumarray' has been added to new version of toolbox - 4.3.2.12172.
Please download and re-install toolbox on your computer.
Currently 'accumarray' is implemented in MATLAB language and thus can be slow for large arrays.
Other than that - it supports all the features, including sparse matrices, fill value, custom functions, etc.
Let me know if you encounter any issues.
Customer support service by UserEcho
By the way, please update toolbox to the latest version (if you use older one).
Now it has EIGS and many other improvements....