BACON outlier detection
I'm trying to run the following code:#include <iostream>#include "mkl.h"int main () { /* Define vector of obseravations (5 2D observation points) */ float pObservations [] = {1., 2., 3.,...
View Articlepardiso using a lot of memory
I'm using Pardiso to solve a large number of linear systems, and its using quite a lot of memory. My program takes the following stepsMake a large, sparse matrixCall Pardiso with phase=11 on that...
View Articlepardiso and pardiso_64 internal addressing
Hi, Using the LP64 interface, will the internal addressing of pardiso be 64-bit, or do I have to use the pardiso_64 version for that? I am using Intel Parallel studio 2016 with mkl 11.3 update 3.My...
View Articleextended eigensolver routines
Hello,I am trying to diagonalize a matrix using the feast algorithm routines. I have written this small program to testzfeast_heev. The following code tries to diagonalize a 4x4 complex Hermitian...
View ArticleRoutines mkl_?tppack and mkl_?tpunpack
The documentation reads the following names for these routines (C interface): LAPACKE_mkl_?tppack and LAPACKE_mkl_?tpunpack. However, these names are not defined in the mkl_lapacke.h header file....
View ArticleLAPACKE_dsyev on Xeon Phi - eigenvalues on xeon phi
Dear MKL Forum, I am testing a Xeon Phi x100 family with 5GB using automatic offload, MKL composer_xe_2013.1.117 and icc 13.0.1, a function as cblas_dgemm is working well with xeon phi. However...
View ArticleComplexity of functions ?potrs, ?potrf and cblas_dgemm
Dear MKL forum,I'm using the functions "?potrf" for Cholesky factorization of a matrix and "?potrs" for solving a linear equation system. Additionally I need the function "cblas_dgemm" (matrix...
View Articlecgelss: MKL 11.2.2 -> 11.3
Hi,Recently, the project I am working on switched from MKL version 11.2.2 to 11.3. However, after we switched to MKL 11.3, cgelss() started throwing an error: "Intel MKL ERROR: Parameter 4 was...
View ArticleCompiling R with MKL
Hi,I'm trying to compile R 3.3.1 on Ubuntu Server 16.04 with Intel MKL using the GNU Fortran and C compilers. Unfortunately, I'm unable to get the gcc to locate BLAS. I've tried several combinations...
View Articlevector array in mkl
Hi does anyone know what is the best mapping in MKL for ippmMul_vaca_64f.Thanks, Saar.Zone: Windows*Thread Topic: Question
View ArticleMulti-thread MKL cblas_sgemm with g++ problem
Here's an example of sgemm program.#include <mkl.h> #include <iostream> #include <cstdlib> #define ITERATION 1 int main() { int ra = 128; int lda = 75; int ldb = 55; float* left =...
View ArticleBest way to use MKL in Python
Hi,I want to install the Math Kernel Library for its use on Python, could you please advise me which one to install, Intel® System Studio 2016 or Intel® Parallel Studio XE?Thank you in advance and I...
View Articleask questions about mkl_csrmm function
Hi everyone,I used mkl_<>csrmm function in deep learning, but I met with a really strange problem. One parameter in mkl_<>csrmm is called pntrb (row pointer in compressed sparse row format)...
View ArticleCompiling Applications with MKL
Hi all,I have MKL installed on my Ubuntu Server 16.04. I've used the Intel-provided mklvars.sh script to provide the directories with the MKL libraries. Every time I try to use these directories with a...
View ArticleEfficient storage of 3D field data
Hi,in working on a numerical code, which solves a discretised equation on a three dimensional grid. I have multiple fields (around 80) I need to save on this grid and which are needed to compute my...
View ArticleBetter way to compute phi0 + sigma*vector?
Hello!I want to compute this quantity prob = phi0 + sigma*atilde, where phi0 and sigma are scalars and atilde a vector 1xind. I have computed it like this:for(i=0;i<ind;i++){ones[i] = 1.0;}...
View ArticleTroubles with undefined _MKLMPI_Get_wrappers
Hello all, recently, I successfully compiled and linked the MKL Pardiso solver within a Fortran 2008 code by including mkl_pardiso.f90 into the code and using the following compiler and linker flags:...
View ArticleClik here to view.
mkl_scscmm performance problem
Hello,I am building the attached program in OLCF's RHEA supercomputer (https://www.olcf.ornl.gov/computing-resources/rhea/) with Intel compiler icc (ICC) 14.0.4 20140805. Armadillo has a naive...
View ArticleInspector-Executor API for Sparse BLAS
Hi, I am trying to use the Inspector-Executor API for SpMV, but it seems like the mkl_sparse_optimize() routine is not performing any optimizations, as I am not seeing any performance difference for...
View ArticleShould VSL Leap Frog work with Wichmann Hill?
I tried vslLeapfrogStream and compared the numbers produced with a simple block of randoms.This worked for the congruential generators MCG31 and 59, but it did not put the numbers in the right places...
View Article