How do I get the current number of MKL threads
Hi, I'm trying to make a section of code that runs with a certain number of MKL threads without affecting the number of MKL threads used in other sections of the code. To do this I want to store the...
View ArticleCan you use version 10.2 with an 11 license?
I'm using MKL 10.2 and and need a new license for another developer. I don't want to use version 11 yet for various reasons. If I buy a new license, can it be used for 10.2 while still allowing me to...
View ArticleMKL parallel execution confusion
I have some confusion regarding how MKL execute in parallel. The problem I have is that after making some changes to a program, calls to DGBTRS, DGETRS, DGETRF and DGBTRF are no longer executed in...
View ArticleCouple of Intel MKL automatic parallelisation/Phi coprocessor questions
Hi,I'm doing an evaluation of the MKL library with a Phi coprocessor for possible use in production at a bioinformatics institute. The production datasets are very large but I'm starting with some test...
View ArticleMKL and Visual Studio express 2013
Hello everyone,I'm am trying some of the given FFT examples from the MKL library in C++ and using VS Express 2013 on Windows 7 pro x64.I linked the my project to the correct libraries and added the...
View ArticleMKL+GCC+GOMP segfaults
Hi, I am experiencing the problem that I mention in the title. I have a program that we have built and linked in three different way: (i) ICC + MKL, (ii) GCC + MKL + GOMP and (iii) GCC + Atlas. Out of...
View ArticleTridigonal solver, gtsv or pardiso
Helle,I want to solve a tridigonal linear system. In mkl, there is gtsv for tridigonal solver and the direct solver Paridiso can also do it. My question is which one is better?Thanks.
View ArticleSummary statistics and padding
When using matrices in MKL, I'm careful to align my columns as recommended by Intel for maximum throughput (this seems to be particularly important on the Intel Xeon Phi). MKL's BLAS interface makes it...
View ArticleScaLapack / BLACS with MKL
I am investigating using ScaLapack/PBLAS for a molecular modelling code in C++ that I am working on. I have access to the MKL on the cluster computers that I work with, but I would also like to be...
View ArticleScalapack linear solver (memory problem)
Hi all,I am trying to use Scalapack in order to solve a distributed linear system. The C++ source code is reported in attachment. The code compiles and runs without problems and it gives the right...
View ArticleClik here to view.
Summary statistics correlation limited to 6 threads?
I've been experimenting with the summary statistics library for computing correlation matrices. I'm using a large sample matrix (10000x6000) of random double values. I've found that the...
View ArticleAre the trigonometric transforms thread-safe?
I'm using the trig transforms (MKL_STAGGERED_COSINE_TRANSFORM) within an OpenMP parallel region on Windows.Each OpenMP thread has it's own array of data to be transformed.When I use one OpenMP thread...
View ArticleDataFitting bug ?
Pls run the program bellow,and confirm the difference between NY==1 and NY==2 case. DFTaskPtr task; // int NY = 2; int NY = 1; double x[] = {1,2, 3, 4, 5}; double y[] =...
View ArticleProblem with Large Array Size in VML
Hi,This is probably a simple one, but I am new to MKL and can't figure it out. I have a large amount of time series data for which I am trying to construct a covariance matrix, and I can't seem to get...
View ArticleClik here to view.
Reduce overhead for repeated offload convolutions on Xeon Phi
I have a series of arrays that I need to convolve. They are not all the same size, but their size does have an upper bound which is ~50. The offload looks something like this (simplified): #pragma...
View ArticleWaht is wrong with cblas_dgemv
Hello, I am interested in the following operation:Y = A^T * Xusing the cblas_dgemv or cblas_dgemm routine in MKL. In the code I use small letters: y, a, and x.I have a test example as...
View ArticleClik here to view.
A'*B using mkl_dcscmm
I tried mkl_dcscmm to compute both A*B and A'*B using a Matlab mex file (64-bit Linux, Matlab 2013a and 2013b) similar to the code posted inhttp://software.intel.com/en-us/forums/topic/472320MKL is...
View ArticleCompiling with g++
How can I compile my program with mkl and g++?I already tried the mkl line advisor, but I got the "pthread_mutex_trylock" message.How`s it going ?Thanks.
View ArticleClik here to view.
FFT Failure Within SideFX Houdini
Hi guys,This is a long shot but I'm hoping you can help, as I'm totally out of ideas and don't know where to go with this. I'm trying to run a plugin inside SideFX's Houdini that uses MKL's FFT. I had...
View ArticleProblem linking MKL with pthread and gomp
A minimal test program:#include "mkl_lapacke.h"int main(int argc, char** argv) { float a = 0.5f; LAPACKE_slasrt('I', 1, &a); return 0; }I used the "Library Link Line Advisor" to get the...
View Article