dss_solve_real takes more time to solve a linear system
I need to solve a system of linear equations with three righ-hand-side vectors. Initially, I was using the sequential version of MKL (compiling with "libmkl_sequential.a") and solving for each rhs...
View ArticleIntel® Math Kernel Library 11.1 is now available
Intel® Math Kernel Library (Intel® MKL) is a highly optimized, extensively threaded, and thread-safe library of mathematical functions for engineering, scientific, and financial applications that...
View Articlefrom CSR to CSC for non-square (m x n) matrix
Hello everyone. I need to convert a non-square sparse matrix (m x n) from CSR format to CSC. But MKL has routines only for square matrix (m x m): mkl_?csrcsc. What can I try to make to complete this my...
View ArticleDoes MKL and PARDISO work properly with more than 64 cores?
This is a continuation of my old question "Does PARDISO have core access limit" with new information. I'll put a quick summary and the new question:I run a code on Windows based HP DL980 server with 80...
View ArticleWhat performance I should expect from following code
Consider following two part of the codes:/* Perform LU factorization and store in DSS_handle */for(k = 0; k < N; k++){gettimeofday(&stTime, NULL);//DSS solver optionsMKL_INT solOpt =...
View ArticleClik here to view.
Segmentation fault trying to use mkl_dcscmm from a Matlab mex file
I'm trying to use mkl_dcscmm from inside a Matlab mex file (on 64-bit Linux, have tried Matlab versions 2009b and 2012a), but having trouble avoiding segfaults. I've tried static and dynamic linking to...
View ArticleMKL-DSS-DSS-Error, Out of memory
Hi all, I am using MKL DSS to solve a sequence of systems of equations in a finite element calculation. The structure of the systems is always the same, and so, DSS_CREATE, DSS_DEFINE_STRUCTURE...
View ArticleIndexing an array of size between 2^31 and 2^32-1 with LP64?
I am working with spblas. My matrix dimension is about 300k, nnz is between 2^31 and 2^32-1. To keep the memory consumption as small as possible. I would like to use 32bit unsigned integer to index my...
View ArticleHow to get LU sparse matrix decomposition?
Hello everyone,how could I get LU sparse matrix decomposition after PARDISO phase 22? May be, MKL has special routines to do it?Thanks.
View ArticleSparse FEAST
Hi,I ran into an issue with the sparse FEAST implementation of "dfeast_scsrev"& "dfeast_scsrgv".The actual matrix I need eigen-decomposed is of size 4096 x 4096. Here is a test case that presents...
View Article-mkl=parallel
I am using the Intel Fortran Composer with the MKL library. I am calling pardiso to solve a sparse system of equations. Everything was working well in OS X using the -mkl=parallel option with the...
View ArticleMKL source code
I'm evaluating MKL for internal use on both servers and potentially a Phi cluster. I've encountered a number of frustrating bugs, poor documentation, and performance problems. Does Intel ever make MKL...
View ArticleAvoiding unnecessary OpenMP synchronization
Say I wish to add a number of vectors:cblas_daxpy(n, 1.0, a, 1, b, 1);cblas_daxpy(n, 1.2, c, 1, d, 1);cblas_daxpy(n, 1.4, e, 1, f, 1);cblas_daxpy(n, 1.6, g, 1, h, 1);MKL will use OpenMP to parallelize...
View Articlemkl_dfti.h: No such file or directory
I installed ipp library successfully setting its environment paths as described in the manual and finally ran a program with:# gcc -o main main.c -I /opt/intel/ipp/include -L /opt/intel/ipp/lib/ia32...
View ArticleSolving non-linear Equations
Hi!I've recently discovered how great and optimal everything in the mkl library is but my knowledge of it is limited. I want to solve a set of 12 non-linear equations in c++ and was wondering how to do...
View ArticleWindows 7/8 64bits, intel mkl 32/64bits installed, 32bits win32 console...
Hello all,I have seen these two topicshttp://software.intel.com/en-us/forums/topic/283594http://software.intel.com/en-us/forums/topic/401425#comment-1764824strongly related to the problem I currently...
View ArticlePardiso - phase 11 errors
Thank you for taking the time to look into this problem with me. I am trying to implement the Pardiso solver for a nonlinear finite element code, particularly fluid flow, and I seem to be getting stuck...
View ArticleSherman Morrisson algorithm : almost tridiagonal matrices
Hello all,Let M be a real tridiagonal matrix of size n greater than or equal to 2. We perfectly know how to solve systems of the form MX = Y with dgttrf followed by dgttrs. Now imagine that I am not...
View ArticleDo I have Intel Compiler installed?
My company has a central license for this: w_mkl_11.0.0.089.exeWhen I installed, it installed Intel Composer XE 2013 (according to the screen and the installation folder, which is: C:\Program Files...
View ArticleHow to implement a DCT method in c++ using MKL?
I am trying to use MKL to implement a DCT method. I cannot find any methods or examples that handle the logic. I found some documentation on...
View Article