C examples in MKL's ScaLAPACK
I am aware of these examples, but I do not know any Fortran, thus I can not understand much. So finding at least an example with MKL and ScaLAPACK in C would be critical for me. I know there is a C...
View ArticleSparse Matrix mkl_?csrmultd problem
Hello,I want to use the mkl_?csrmultd to do 2 matrix product and the output is a dense matrix.But i am confused when i read the manual, The ldc (leading dimension of dense matrix C) is a output...
View ArticlePARDISO consistent crash for INCORE RUN
Hi, I am trying to solve several big 3d solid FE models with PARDISO 11.2Although the out-of-core run is successful I am consistently getting segmentation fault errors for the in core runs.This also...
View ArticlePARDISO for complex symmetric matrix
Dear all,I used PARDISO for complex symmetric matrix. The code is as follows:void PardisoSolver(int n, int nz, int *rowind, int *colind, MKL_Complex16 *a, int nrhs, MKL_Complex16 *x, MKL_Complex16 *b)...
View ArticlePARDISO for complex symmetric matrix
Dear all,I used PARDISO for complex symmetric matrix. The code is as follows:void PardisoSolver(int n, int nz, int *rowind, int *colind, MKL_Complex16 *a, int nrhs, MKL_Complex16 *x, MKL_Complex16 *b)...
View ArticleBug report: Pardiso get wrong result when solving complex symmetric matrix
Dear all,I use the following code to solve complex symmetric matrix by Pardiso and get a wrong result, could anyone help me to take a look at it? Thanks #include <stdio.h> #include...
View ArticlePARDISO segfault in 64 bit mode on Mac OS X
Hi, I am using Intel MKL version 11.2 update 3 (latest I believe) on Mac OS X Yosemite (10.10.3). Specifically I am noticing that the PARDISO direct sparse solver is crashing in the 64 bit mode...
View ArticleClik here to view.
How to use both: FFTW3 and MKL
Hey!I've got the the following problem:I want to keep the flexibility within my code to switch between FFTW3 and MKL-FFT (with the purpose to test performance issues of both FFTs and thus, decide which...
View ArticleIs iterative solver in Parallel with OPENMP?
Hi, all, Right now I have to solve a linear equation using gmres with ILU precontioner on SMP system. I believe MKL is able to do so. I am just wondering which parts of the iterative solver are...
View Articlesee lots of warnings when compile Lib with MKL used
Hi, I recently create and build a lib on window host with mkl lib function used, however, I see lots of weird warnings , can you tell me how to get rid of them, thanks a lot. I included below...
View ArticleCompling error: undefined reference to `__isoc99_sscanf'
Hello, everyone, I tried to apply MKL to compile several fortran code. However, I got some errors below. Also, I attached the makefile. It seemed something was not linked. Is there anyone who can give...
View ArticleFEAST Eigensolver not returning all eigenvalues in specified range
Hi,I have came across a problem with FEAST in MKL 11.2, whereby it is not returning all the eigenvalues in the specified search range. The problem type is Generalized Sparse (feast_scsrgv).Default...
View ArticleFailure during collective
Hi , I have compiled espresso with intel mpi and MKL library but getting error Failure during collective error when ever it is working fine with openmpi. is there problem with intel mpi Fatal error in...
View ArticleSame numbers differente signal
Dear MKL forum, I am rewriting some code using the MKL library from a fortran + open source math library, in the validation process i realize the function LAPACKE_dgesv giving the exact solution...
View Articleabout in-place VML operations
Hi all, I am writing to check if MKL VML library supports in-place operations like vzmul(100, x, y, x) or vzmul(100, x, y, y); I saw a post in this forum and there it said VML support in-place...
View ArticleWhy mkl 11.2.3 can not use all the cpu in Intel(R) Core(TM) i7-4700HQ CPU @...
Hi: I use mkl_dss to solve a problem. I already use mkl_set_num_threads(8) to set the maximum threads of my computer. But when I run the program, I use top command and just can see only 4 cpu are...
View ArticleMKL 3D FFT of 1D function (fortran)
Dear all,I want to perform a 3D Fourier transformation (using fortran) of a complex function which is in a 10^9 array, where each element corresponds to a 3D vector k. So each element of the array is...
View ArticleProblem with ?gemm and AMD processor
My program is compiled with the latest revision of IVF 2015 Composer edition. It makes use of a number of mkl functions including CGEMM and ZGEMM. The code involving the mkl functions is unchanged and...
View ArticleMKL sparse BLAS
Hi all, Please, are threaded the following MKL sparse BLAS functions? 'mkl_cspblas_?coosymv', 'mkl_?coosymv' and 'mkl_sparse_?_mv1 (11.3 beta) I'm not able to run with more than one thread these...
View Article3D FFT in MKL with data larger than cache
Hi, I am working on a 3D numerical integrator for a non-linear PDE using the parallel FFT library included in MKL. My arrays consist of 2^30 data points which is much much larger than the cache. This...
View Article