Using multiple DFTI DESCRIPTOR (FFT in MKL)
Is it possible to create and commit several different DFTI descriptor and re-use them later (the FFT of different sizes will be called many times, and creating the descriptor and free it for each call...
View ArticleScalapack raise error under certain circumstance
Dear All, I am using IntelMPI + ifort + MKL to compile Quantum-Espresso 6.1. Everthing works fine except invoking scalapack routines. Calls to PDPOTRF may exit with non-zero error code under...
View ArticleSVD hangs periodically
We are noticing that SVD, both dgesvd and dgesdd, will hang periodically. The call stack terminates with a call to either one of those functions. Killing the process and rerunning alleviates the...
View ArticleHow to add MKL with TBB threading to application that uses TBB under Visual...
I have installed IntelSWTools 2017.4: c:\Program Files (x86)\IntelSWTools\compilers_and_libraries\windowsMy application uses TBB: Release(/MD): tbb\lib\intel64\vc12\tbb.lib ->...
View ArticleAre DNN functions thread safe?
It seems that dnnExecute_F32 can not be called in multi-thread, does it?
View ArticleHow to store A to get fastest performance of AT*x using cblas_dgemv?
Hello,I am using cblas_dgemv to obtain AT*x. The size of the matrix A is about 10000 Rows x 20000 columns. I am storing A in row major format. Ai,j+1 is stored next to AijMy questions are as follows...
View ArticleWarnings from libiomp5.a when linking on Macintosh with Xcode 8
We are using MKL and linking with libiomp5.a. Since starting to use Xcode 8.3.3 to build, we have been getting a large number of warnings like this::-1: warning: pointer not aligned at address...
View ArticleRegarding cluster_sparse_solver
I am Mehdi and this is my first time using this forum.I need to used cluster_sparse_solver in my FORTRAN Finite Element program. Because the degree of freedom of my system is very high (1^6), the...
View ArticlePardiso and pardiso_64
Hello,I have a FORTRAN program for solving flow field (Stokes flow) problem using FEM. Firstly, I used PARDISO solver to solve the coupled problem. I was totally happy with PARDISO, because it is...
View Articlecluster_sparse_solve library and path setting in UBUNTU
p { margin-bottom: 0.1in; line-height: 120%; } Hello,I have a question regarding compiling a program containing cluster_sparse_solve.I have developed a program for 3D flow field calculation using...
View ArticlePardiso iparm(30) not returning equation number correctly
I am using Pardiso with the 2017 Update 2 Intel Fortran compiler in VS 2015 and I'm finding that when using mtype=2 (real and symmetric positive definite matrix) if my matrix has a singularity,...
View ArticleMKL Memory Allocator
Hello,I am looking for information regarding the Memory Allocator embedded in MKL.We are using intensively MKL_malloc/MKL_free in a project and are planning to add a memory manager on top of it. Our...
View Articlecluster_sparse_solver library and path setting in LINUX
p { margin-bottom: 0.1in; line-height: 120%; } Hello,I have a question regarding compiling a program containing cluster_sparse_solve.I have developed a program for 3D flow field calculation using...
View ArticleMKL_SINGLE_PATH_ENABLE
Hello,We have legacy code that call :mkl_enable_instructions(MKL_SINGLE_PATH_ENABLE);MKL_SINGLE_PATH_ENABLE is not in the documentation (anymore ?).In mkl_service.h, defines are :#define...
View ArticleCANNOT use cluster_sparse_solver in linux!
p { margin-bottom: 0.1in; line-height: 120%; }Hi everyone, I have developed a code which uses cluster_sparse_solver. I can compile it (static) in my PC, which has WINDOWS10 OS. But, I can not compile...
View ArticleAccess Violation Error while using dgesvd for c on Visual Studio
Hello, while running the following code on Visual Studio 2015: #include <thread> #include <mkl.h> #include <random> #include <ctime> const int MatrixLayout = LAPACK_COL_MAJOR;...
View ArticleClik here to view.
Pardiso generates wrong result for a small test matrix
Hi, I've written a test program to get familiar with the pardiso. However I find that comparing with the result from Eigen the pardiso spits out totally incorrect answer. I have totally no idea what...
View ArticleIs DFTI_NUMBER_OF_TRANSFORMS data-parallel?
If I set DFTI_NUMBER_OF_TRANSFORMS to 4 on a AVX computer, or 8 on a AVX-512 KNL, will MKL's DftiComputeForward/Backward compute the FFT's of similar but independant, non-overlapping arrays...
View Articlecluster_sparse_solver and cluster_sparse_solver_64
Hello all,I have developed a code for 3D fluid flow using FEM coupeled method. Once I have used cluster_sparse_solver. For compiling this file I have used the following terms, using intel link...
View Articlesoftmax
Does MKL support the popular DNN operation called softmax?I cannot find any suitable function.
View Article