OpenMP* Threading for FFT Computation
Dear Sir/Madam,As I understood from an example on https://software.intel.com/en-us/mkl-developer-reference-c-examples-of-u..., you calculate with 4 threads 4 different 2D-FFT (50x100 points). My...
View Articlecorrect paths for mkl pardiso dll to install ipopt on windows using bash
Hi, can someone tell me what path should exist to the dlls at configure time to use the ia32 version of pardiso on windows.I am trying to build ipopt using a bash shell but using the Microsoft C++...
View ArticleProblem with Fortran version Cluster FFT Compiling
Hi all,I get serve DftiCreateDescriptorDM initialization error in fortran version Cluster FFT functions (while I have no problem with c version compiling with similar code), could someone help me find...
View ArticleProblem compiling Fortran code with MKL
Hi all,I'm trying to replace my own LU decomposition code with the fortran subroutine gbtrf from MKL. Please bear with me as I am somewhat new to Fortran.I successfully compiled, run and validated the...
View ArticleClik here to view.
Inaccurate results while using PARDISO
Hello, I am trying to utilize the PARDISO sparse matrix solver to test out a solution with my sparse matrix. The goal of this experiment was to see if I can replicate the same solution that I had while...
View ArticleIMSL Library not included in INTEL PARALLEL STUDIO XE COMPOSER for WINDOWS 2018
Hi , IMSL library is not included in Intel studio XE 2018 for windows . I have used IMSL functions in my earlier code. How Can I build that program in the parallel studio 2018 without IMSL library....
View Articlepardiso with parallel MKL and TBB on large systems
Hi,I'm solving a large sparse system with pardiso. I'm using the Parallel MKL library with TBB enabled on Windows 64bit system. My A matrix is 57k x 57k and my right hand side B matrix is 57k x...
View ArticleMKL LLC usage policy
Hello, experts,Assume that I have 4 cores machine, each core has 2MB of LLC slice and LLC includes L2. 1) If I use single-threaded MKL, the MKL instance will use 2MB of LLC or use 8MB LLC?2) If I use...
View ArticleProblem of cblas_zgemmt
I would like to use cblas_zgemmt() to calculate C=C-A*B. The final C is symmetrical.The input matrices is verified by cblas_zgemm() and the result is correct. However, the output C is unchanged when...
View ArticleUsing Armadillo with MKL - Beginner
Hi all. I had a few questions regarding the usage of the Armadillo library with MKL. Apologies in advance if this is not the right place for them. I'm currently on Ubuntu 16.04, using gcc 5.4.0. When...
View ArticleIn my program, a thread is stop at the intel mkl API DftiComputeForward...
In my program, a thread is stop at the intel mkl API DftiComputeForward randomly, Please give me a mind to solve the problem
View Article[Help] Example works with extention .c but fails with .cpp
Hi,I have a big (n~100,000) complex matrix which I would like to diagonlize.After googling extensively, I found the following link provides excellent...
View Article[help]Data type in zheev.c example.
Hi,In the example ofhttps://software.intel.com/sites/products/documentation/doclib/mkl_sa/11...The to-be-diagonalized matrix is declared to be a dcomplex type(struct _dcomplex { double re, im; };...
View ArticleCalling Fortran Blas routine from C
Hello,I am trying to call Fortran subroutine that calls zdotu from C.I am using linking flags from Advisor https://software.intel.com/en-us/articles/intel-mkl-link-line-advisor/These linking flags for...
View ArticleFailure of DGESVXX
Hi all, I want to report a problem with DGESVXX. When I compile my program that uses this routine with the option "initialize stack variables to an unusual value" it fails when run with this message:...
View ArticleQuestion about mkl_?csrgemv
Hi, when I read reference manual of mkl, the mkl_?csrgemv routine calculates y=Ax (or y=A'x) where A is m-by-m matrix in 3-array CSR format.However, when I tested, I figured out that this routine also...
View Articlegetrf hangs when using intel threads in calling application
I am calling getrf from one thread in my application, the other threads are waiting at a taskwait. Task manager looks like all threads are maxed out and getrf never returns.The matrix size is quite...
View ArticlePossible Heap Corruption after calling DftiCreateDescriptor
Hello,Our application crashes some time after calling DftiCreateDescriptor; we believe because the heap has become corrupt. If we delete the call, the application no longer crashes, or if we reduce the...
View ArticleMemory usage during solve in function of number of cores
During the solve step of Pardiso of the system A*X=B where B is multiple right-hand-side we see that the memory consumed by Pardiso is highly dependent on the number of cores. So we observe mem(total)...
View ArticleSparse Matrix Multiplication Failed
Hello,I am trying to develop a function that multiplies matrix A and B which are in general format but sparse in nature. These matrices contain complex numbers. My issue is, when I do not use the...
View Article