- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi
I am having very serious precision issues by using intel mkl-lapack for matrix inversion:
Steps:
(1) I inverted a matrix using Matlab/Octave
(2) I use dgetrf and dgetri to invert the same matrix on two processors (a) Intel(R) Core(TM) i7-2600 CPU 3.40GHz for the test code/16GB of RAM on a windows machine using Intel Parallel Composer XE 2013, and (b) Intel(R) Xeon(R) CPU X5660 2.80GHz on a linux machine by using Composer xe 2011
(3) The problem is that the difference between the inverse obtained using Matlab/Octave and by using dgetrf and dgetri is different. There are differences is not an issue but the differences are based on processors is creating problems in large simulations. The answers received by using Intel Xeon Processors and Intel Composer XE 2011 Machine are more accurate than what is obtained by using the same code on windows machine
At this moment I think I am overlooking something i.e. creating a big mistake. An advice on solving this issue would be greatly appreciated. I have attached a sample code to highlight this issue. I have included the sample code but I was not able to upload input binary files on the forum (It was taking long long time)
Many thanks
Vineet
Link Copied
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Here are the command lines you requested. Send me the C/C++ files and I will execute them to verify precision control
For linux (Intel Xeon processor)
ifort source1.f90 –heap-arrays -openmp -L /share/apps/intel/composer_xe_2011_sp1.7.256/mkl/lib/intel64/ -I /share/apps/intel/composer_xe_2011_sp1.7.256/mkl/include/ -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -liomp5 –lpthread –o source1.exe
For Windows
Compiling with Intel(R) Visual Fortran Compiler XE 13.1.0.149 [Intel(R) 64]...
ifort /nologo /debug:full /O2 /I"C:\Program Files (x86)\Intel\Composer XE 2013\mkl\include" /warn:interfaces /module:"x64\Debug\\" /object:"x64\Debug\\" /Fd"x64\Debug\vc100.pdb" /traceback /check:none /libs:static /threads /dbglibs /Qmkl:parallel /c -heap-arrays /Qvc10 /Qlocation,link,"C:\Program Files (x86)\Microsoft Visual Studio 10.0\VC\\bin\amd64" "C:\Users\Vineet Work\Documents\Visual Studio 2010\Projects\source1.f90\source1.f90\Source1.f90"
Linking...
Link /OUT:"x64\Debug\source1.f90.exe" /INCREMENTAL:NO /NOLOGO /LIBPATH:"C:\Program Files (x86)\Intel\Composer XE 2013\mkl\lib\intel64" /MANIFEST /MANIFESTFILE:"C:\Users\Vineet Work\Documents\Visual Studio 2010\Projects\source1.f90\source1.f90\x64\Debug\source1.f90.exe.intermediate.manifest" /MANIFESTUAC:"level='asInvoker' uiAccess='false'" /DEBUG /PDB:"C:\Users\Vineet Work\Documents\Visual Studio 2010\Projects\source1.f90\source1.f90\x64\Debug\source1.f90.pdb" /SUBSYSTEM:CONSOLE /IMPLIB:"C:\Users\Vineet Work\Documents\Visual Studio 2010\Projects\source1.f90\source1.f90\x64\Debug\source1.f90.lib" mkl_intel_lp64.lib mkl_intel_thread.lib mkl_core.lib libiomp5md.lib mkl_lapack95_lp64.lib "x64\Debug\Source1.obj" "x64\Debug\Source2.obj"
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page