MKL Interpolation gives strange results (fortran)
With the help of the reference manual for MKL data fitting i have tried to put together a I have a piece of fortran code for interpolation using MKL. It seems to me to give somewhat strange results (I...
View ArticleUsing openmpi with the mkl's cluster library.
Hi all,I use the following tool chains to compile the vasp.5.3.5:1- ifort, icc, and mkl are the ones bundled in composer_xe_2011_sp1.11.339.2- openmpi version is openmpi-1.8.7.It's well known that...
View ArticleF90 SIMD/Parallel FFT
I found some related posts and links such as https://software.intel.com/en-us/node/507041#EEACAA31-C805-484F-B4AE-D8F...However nothing that jumped out.Most of what I do results in executables that run...
View ArticleGenerate Sobol Sequence
Hi,I would like to generate numbers from the Sobol sequence in dimension n (with n < 40). The MKL has a Sobol sequence generator, but reading the documentation only gives me headaches. I gave up...
View ArticleHPL benchmark compilation problem
Hello, I am trying to build from sources the Intel version of HPL benchmark which is included in MKL libraries. I am able to make the compilation using the 'non-offload' version. When I try to build...
View Articlemkl_?csrcoo for rectangular matrices?
Hello, I hope this is not too much of a beginner's question, I just started using sparse matrices with MKL a few days ago. This is work within a large codebase so working examples are very time...
View ArticleParameter double **mat to cblas_dgemm
Dear Intel Forum,I am developing a molecular dynamic system with a several MKL functions, the atoms position is a double **mat to the mkl function like cblas_dgemm. However, is necessary to convert...
View Articlesparse_matrix_checker error value constants
The MKL manual lists the following constants for the error value returned by the sparse_matrix_checker subroutine: MKL_SPARSE_CHECKER_SUCCESS, MKL_SPARSE_CHECKER_NON_MONOTONIC,...
View ArticlePARDISO Questions
Hi Pardiso experts,I have a couple of questions about PARDISO.First a simple one. I am using out of core PARDISO and specify PARDISO_OOC_KEEP_FILE=0 to keep the files.But then when I want to delete...
View ArticleFix arithmetic error
I am performing a distributed Cholesky factorization and inversion of a matrix. I assume that the serial versions of these routines give the correct result, thus I am trying to make my distributed...
View ArticleClik here to view.
[FORTRAN] MKL FFT issues
Hello,I am attempting to test my FFT routines with known solution transforms, and the results that I have been getting are not what is expected. In order to maintain generality for my final...
View ArticleIMSL and MKL
Dear All,Is there some documentation on how to live without ISML and make full use of the MKL library.I am particularly interested in CSVAL and CSINT.Many thanks for your help.
View ArticleHow to set the adaqute parameters for feast eigen solver
I'm using feast general Eigen solver to solve the Eigen value equations. However, I found it is difficult to follow a general rules to set the emin, emax and m0. It seems I need to know in advance...
View ArticleParallel cluster sparse direct solver on very large matrix
Hi there,We have some C code working for solving a big linear system that is quite sparse using the cluster sparse solver routines in MKL. We generally followed the `cl_solver_unsym_distr_c.c` example...
View ArticleFree technical webinar to learn new Intel MKL 11.3 and Intel IPP 9.0
Intel® MKL 11.3 introduces the new optimization for different functions, and provides MPI wrappers which support custom MPI* implementations such as Cray* MPI, and IBM* Platform MPI. It adds Sparse...
View ArticleUsing cluster_sparse_solver with 64 bit integers
I am trying to use `cluster_sparse_solver` with 64-bit integers (longs) for extended matrix sizes. The first thing that I have tried is compiling my code with `-DMKL_ILP64`, which changed the `MKL_INT`...
View ArticleMultithreading when called by C++, not when called by R
Hi everyone, I have been struggling with a problem for quite some time, and I would greatly appreciate your input.I have a program that uses Intel MKL's dgemm function many times. In fact, to...
View ArticleWhy no simple linear equation solver?
I looked in the INTEL MATH KERNEL library, and could not find a routine that just solvesA*X=Bwhere A is an N by N matrix, and X and B are vectors of length N.They throw a lot of confusing **** at you,...
View ArticleClik here to view.
floating point exception solving a non-symmetric matrix
The floating point exception occurs when solving the attached non-symmetric matrix.1 line: n2nd line: index base (i.e. 0 based)3rd line: nzfollowed by all nonzeros in coordinate (i.e. triplet)...
View ArticlePerforming DFTI along a single axis of 2D array
Hi there,I have a 2D array where the axes are positions (x, y). I wish to perform a Fourier transform along the x axis alone resulting in an array (Kx, y), where Kx is momentum, the Fourier transform...
View Article