C++ MKL BLAS wrappers vs expression templates
This is a conceptual question:Expression templates are a popular technique in C++ in order to implement Matrix and Array operations by avoiding unnecessary temporaries and loop unrolling. In other...
View ArticleHow can I reuse sparse factorizations in Pardiso
Hello, I have a serials of structurally identical matrixs such as {A1, A2, A3,....}and I need to solve A*X=Y, for A1,A2,A3......., note that rho vector Y changes as time goes while all matrixs are...
View ArticleUsing FEAST for large matrix
Hello, I am presently working with FEAST to find eigenvalues and eigenvectors for a symmetric matrix. I need to solve N X N matrix with N ~ 10^6- 10^8. Now I have few queries : 1. SInce the...
View ArticleIntel ODE Library on Mac OS
An application I would like to run on my macbook links to the Intel ODE library.Unfortunately, the library is available only for Windows and Linux: binaries but not source are available for...
View ArticleSVD produces wrong results in mkl=parallel (2013 sp1)
I have experienced a strange bug in MKL: zgesvd produces different results (some wrong) depending on the number of threads that MKL uses. Above 2, the singular values all become NaN, even though the...
View ArticleRunning multiple Pardiso solves concurrently
Hi,We're using MKL Pardiso inside an optimisation web service on Windows and Linux. Clients can spin up multiple optimisations in one call to the service, and so we have multiple runs occurring...
View ArticleClik here to view.
cluster_sparse_solver computes wrong solution
Hello,I'm trying to use cluster_sparse_solver and solve a system in-place (iparm(6) = 1), with a distributed format (iparm(40) = 1). I adapted the example cl_solver_unsym_distr_c.c as you can see...
View ArticleBuild Scipy With MKL failed
I use command to build numpy first:python %MYPWD%/%NUMPY_VER%/setup.py config --compiler=msvc build_clib --compiler=msvc build_extthe site.cfg content is:[mkl] library_dirs = C:\Program Files...
View ArticleClik here to view.
Help needed with bdsqr
Hello,I'm trying to compute a partial SVD of a rectangular matrix A. I tried to adapt an example of the MKL that uses gesvd. While I get the same singular values, I'm not able to compute the correct...
View ArticleParameters for ?stemr
I have a problem where I need to calculate a number of eigenvectors and a different number of eigenvalues. Instead of calling dsyevr twice I plan on calling dsytrd -> dstemr * 2 -> dormtr. (Or...
View Articleperformance of mkl in multithread application in numa memory architecture...
hii have a server with 80 logical core with NUMA memory architecture and windows server 2012. i want create one thread per each logical core and each thread execute fft, convolution,......
View Articleintelem not available as fcompiler
After install MKL 2015. I can only see "intelvem" is available. "intelem" is not available.What I need to do if I want to use intelem as fcompiler? IntelEM64VisualFCompiler instance properties:...
View Articledcsrmm is throwing 'integer division by zero' after upgrading to 11.0 update 5
After updating the installed MKL version when attempting to avoid another bug we discovered, previously working code has started throwing division by zero errors after sparse matrices exceed a certain...
View ArticlePardiso solver, error 157
Hi everyone,I'm using intel visual fortran composer xe 2011 on MVS 2012. I recently start having the error 157 when calling the pardiso solver. I keep receiving the following message : "forrtl: severe...
View ArticleCluster Pardiso error for iparm(24)=1
Hello all,Last week I have been trying out Cluster Pardiso, testing various options and so on, and I found a problem. Cluster Pardiso fails with various error messages when I enable two-level...
View ArticleHow to Control Pardiso Threading
Hi, dear engineers,Here is a little bit background: in our solver, we have a precondtioner made up of a block diagonal matrix and a sparse matrix. We use tbb+mkl to LU factorize the block diagonal...
View ArticleFortran FFT with MKL
I'm currently messing with some code that uses an old F77 looking FFT subroutine full of goto's and other nasty stuff. I have been looking around at using MKL's FFT routines, but I am having trouble...
View ArticleRank updates to pardiso factorization
Hi,Is it possible to do rank-1-updates of the factorization computed by pardiso? That is, given A = L D L^T, I need to be able to efficiently compute the decomposition of A + alpha w w^T (w is a...
View ArticleMultiple Xeon Phi's
Hello,Is it possible to use LEO directives to offload a ZGEMM call to multiple Xeon Phi coprocessors in a way such that they distribute the work among them automatically without me having to partition...
View ArticleProblem of multithreading mkl_?csrcsc
Hey GuysWhen I use mkl_?csrcsc to transpose my sparse matrix, I think the multithreading is supported internally: just like multithreading BLAS function cblas_dgemm(), call omp_set_num_threads() to set...
View Article