Fast Poisson solver with inner boundary conditions
Dear all,I would like to calculate the potential due to a point charge in the proximities of a conducting cylinder. For this I started by calculating the potential due to the point charge alone using...
View ArticlePARDISO Complex
HiI am currently working with PARDISO to solve an iterative process that uses complex numbers; the iterative process requires that the sparse matrix be factorized only once (i.e. phase 12). My concern...
View ArticleCompilation Errors with MKL_DIRECT_CALL_SEQ
Hi,I'm compiling a cpp file that includes mkl.h and contains a single cblas_dgemm call with the following:icpc -DMKL_DIRECT_CALL_SEQ -DMKL_ILP64 -I/opt/intel/mkl/include -c -o main.o...
View ArticleMatrix product with MKL
Hi,I want to do a complicated matrix product which is defined asC[i,j] = \sum_{n,m} A[i,n]*A[j,n]* B[n,m] *A[i,m]*A[j,m] Where A and B are two matrix with size of NxN, (N~1000).To accomplish this...
View ArticleLooking for MKL FFT complementary functions
Hey I migrated my DFT code from IPP to FFT code using MKL FFT functions, I'm looking forward for performance improvement.1st, I want to do phase correlation in spectrum domain. In my algorithm I have 2...
View ArticleIntel® MKL 11.3.3 patch
There are two listed below limitations with Intel® Math Kernel Library (Intel® MKL) 11.3 Update 3 which were discovered recently. The official fix of these issues will be available the nearest update...
View ArticleCluster Sparse Solver -- Matrix input format & MPI problem
Hello,I am implementing MKL Cluster Sparse Solver (11.3.3.210 on Linux). I've been able to succesfully run a routine that works with centralized input format (nonsymetric real matrix). I am trying to...
View Articlesparse matrix multiplication mkl_?csrmultcsr example and limits
Are there any examples of mkl_?csrmultcsr functions? I have MKL 11.3 but have not found examples for these functions. I am also interested in knowing the maximum dimension sizes for rows/cols of the...
View ArticleFeast can't find eigenvalues in the specified interval
Hi,I am trying to solve eigenvalue problem using zfeast_hcsrev routine, here is the excerpt from my codeMKL_INT fpm[128]; feastinit(fpm); fpm[0] = 1; fpm[1] = 20; char uplo = 'U'; MKL_INT n = 85;...
View ArticleSolving an ill-conditioned, sparce unsymmetric linear system
Hi,I am trying to solve a set of linear equations for a material model through c++. The system was already solved in Mathematica using Krylov method (BiCGStab).Through Mathematica, I found that the...
View ArticleHow to use MKL Routines for Matrix Inversion
Hi everyone!I never used libraries in my programs, usually I do all my routines on my own, but now I'm having efficiency problems and I decided to use libraries to improve the efficiency of my...
View ArticlePardiso 11.3 does not work for iparm[30]=2 for partial solve
I recently tested the partial solve option with iparm[30]=1. I found setting only one element of perm to be 1 works, but when I set more than one element of perm of be 1 does not work. I also tried to...
View ArticleMKL Inspector-executor
Hi,Is MKL Inspector-executor API written specifically for use on multi core processors, Phi processors, or both?Moreover, if it is for use on Phi, and at this point in time, is it for both KNC and...
View ArticleEntry point mkl_serv_check_ptr_and_warn not located in mkl_core.dll
We have a program that uses 32 bit version of MKL v10.2.6.037. The program had been tested and released for more than one year. Recently one user had this error dialog poping up saying "The procedure...
View ArticleFortran PDGETRF + PDGETRI for Matrix Inversion - parameter with an illegal value
Hi everyone,This is the first time that I am using libraries. I need to calculate the inverse matrix of a square matrix. I found the routine PDGETRF + PDGETRI from ScaLAPACK in the MKL documentation...
View ArticleEigen+MKL vs. MKL-Only
Currently, I have a code that uses Eigen (a C++ template library for linear algebra) to save a square general dense matrix in the following wayZMatrix = new...
View ArticleDeal with huge matrix
I have written some codes with MKL library, but recently these codes do not work. They can be complied well and do not report any error when I run it. However, they did not work for some huge matrix...
View ArticleAre VSL Routines ThreadSafe?
Hi,I have written a simple code to generate Random Numbers using VSL Routines. It is an OpenMP code written in Fortran:INCLUDE "mkl_vsl.f90" program vsl_ex use omp_lib use mkl_vsl_type use mkl_vsl...
View ArticleSetting the MKL dist dlls path, for the single dynamic library
Hi,In Windows, the mkl_rt.dll searches for the MKL dlls it needs in the PATH env variable. Is there other way besides setting the PATH to tell the mkl_rt.dll where the MKL dlls resides.The motivation...
View Article