Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- Quad precision for MKL functions

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

S__MPay

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-01-2017
12:10 AM

147 Views

Quad precision for MKL functions

Is it possible to use quad precision for MKL functions and maintain the accuracy?

Link Copied

5 Replies

TimP

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-01-2017
01:07 AM

147 Views

Is your question answered by previous threads about MKL and quad precision or real*16 ?

S__MPay

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-01-2017
03:12 AM

147 Views

Tim P. wrote:

Is your question answered by previous threads about MKL and quad precision or real*16 ?

According to a similar threat:

Tim P. wrote:

You would have to compile the public source code (netlib.org) to make the quad precision versions of the functions you want. You would be giving up on the "performance" aspect of MKL, as there is no support for vectorization in quad precision.

But my confusion is for some inaccuracies that are causing trouble are gone when I use quad precision even though all I did was to simply call double precision version of MKL.

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-01-2017
06:25 AM

147 Views

MKL doesn't support quad precision. MKL's API supports single & double precision only.

TimP

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-01-2017
08:01 AM

147 Views

S__MPay

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-01-2017
09:27 AM

147 Views

When I try to calculate the determinant of a matrix there is error in double precision while the error of quad precision is exactly zero:

print *, 'zero determinant double precision' , dprMklDet(reshape((/ 6.d0, 3.d0, 5.d0, 4.d0, 4.d0, 7.d0, 12.d0, 6.d0, 10.d0 /),(/3, 3/))) print *, 'zero determinant quad precision ' , qprMklDet(reshape((/ 6.q0, 3.q0, 5.q0, 4.q0, 4.q0, 7.q0, 12.q0, 6.q0, 10.q0 /),(/3, 3/)))

Output is:

zero determinant double precision 1.065814103640150E-014

zero determinant quad precision 0.000000000000000000000000000000000E+0000

I will post more detail in this topic https://software.intel.com/en-us/forums/intel-math-kernel-library/topic/733238

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

For more complete information about compiler optimizations, see our Optimization Notice.