Clik here to view.
Rewrite subroutine to use MKL - matrix multiplication
p { margin-bottom: 0.1in; line-height: 120%; }Hi,My first post on the forum!I'm working with an older code written i Fortran at uni. The program consists of several subroutines and I would like to...
View ArticlePARDISO b array
In most problems involving a matrix inversion, one only looks at a single b array in the Ax=b problem.PARDISO has however an arrangement, where b is actually b[*,nhrs] in both the Geneva and the MKL...
View ArticleCall objective function in Parallel in jacobi matrix calculation function...
I'm learning to use the Non-linear optimization in the Intel Math Kernel Library. My starting point is the the example code "ex_nlsqp_bc_c_x.c" in the directory "examples_core\solverc\source". In its...
View ArticleInteresting Observations
I have had a lot of fun and some heart ache playing with the MKL Library and the PARDISO solver.The PARDISO solver is very fast. The problem I am solving is the hydraulic network problem, in this case...
View ArticleClik here to view.
CSCMM for Armadillo Sparse dense multiplications
Environment : armadillo 4.320.0 and 4.400 Compiler : Intel CPP compiler OS : Ubuntu 12.04I am trying to replace the Armadillo's native sparse dense multiplication with Intel MKL's CSCMM call. I wrote...
View Articlexerbla overrides VML callback when mixing with BLAS, LAPACK, VSL
I have encountered an issue with the MKL error handling functions when mixing calls to BLAS, LAPACK, VML, and VSL in C. Essentially, the VML callback, which is being set using vmlSetErrorCallBack, is...
View ArticleDifferent Eigenvalues on Different Machines
Hello, I am computing eigenvalues on different machines. With the exact same binary and the same shared libraries, I obtain different eigenvalues for near zero entries when I run on different machines....
View ArticleMKL random number stream equivalent to Matlab default RandStream
Is it possible to generate random numbers with MKL that are equivalent to Matlab random numbers? I use the following Matlab/ Fortran codes but the results are differentMatlab...
View ArticleWant to link MKL to existing code
Hi, I am replacing certain operations specifically vector dot product,matrix vector multiplication and solve AX=b in existing code by calls to lapack subroutines. Do I need to link the MKL in the...
View ArticleGetting very low efficiency for mkl_dcsrmv function
Hi,I am trying to benchmark the intel mkl library function mkl_dcsrmv. I tried running the function for single, three,five and seven diagonals nnz entries and i repeated it multiple times with...
View ArticleWarming up strategy for MIC dgemm call
In my computation, I manually offload some computation to MIC using offload pragmas. Offloaded computation also involves a call to MKL's Double precision general matrix-matrix multiplication (dgemm)....
View ArticlePerformance of offloaded MKL FFTs on the MIC, anyone?
My initial experiments offloading MKL FFTs into the MIC (using C language in Linux) give me approximately 9.3 GFLOPS of performance, judging by the reported [MIC Time] numbers when I set the...
View ArticleLinking against both the sequential and threaded mkl
I have two dlls that link against the static mkl libraries. One of the dlls links against the sequential version and the other against the multi-threaded version. Those two dlls are then loaded in to...
View Articledjacobix only uses 4 threads on a 16 CPUs virtual machine
Hello,I'm using djacobix in Intel MKL. My testing machine is a virtual Windows Server 2012 with 16 CPUs, . I'm use the following statements in my code: mkl_set_dynamic(0);...
View Articlephase 13 does not work in cluster_sparse_solver
Hi: I take mkl/examplescluster_sparse_solver/cl_solver_unsym_distr_c.c and compile it with mpic++ -cxx=icpc -mkl -xHost and it works well. From the mkl document, I know that I can combine phase 11,...
View Articlehow to split matrix in cluster_sparse_solver, similar number of rows or...
Hi: I need to use cluster_sparse_solver to solve a larger complex symmetric matrix. So I need to split matrix into two parts and put them into two different compute nodes respectively. The first part...
View Articledgemm - Intel inspector: "uninitialized memory acces"
Hi,I launched an Memory Error Analysis with Intel Inspector on simple test of the dgemm function and i got this warning:"uninitialized memory access"located on the dgemm call.Here is my test (really...
View ArticleProblem building 64 bit numpy using MKL and vc11 (Windows)
Hi, I’ve built numpy (1.9) 64 bit using vc11, the Intel Fortran compiler and the MKL ‘mkl_rt’ library.*why? (see end of message for the reason, if interested) Any advice or assistance would be greatly...
View Articlepardiso_getdiag
HelloI am calling "pardiso_getdiag( pt, df, da, mnum, error)" after PARDISO call (phase = 22) and it returns diagonal terms of initial and final (factorized) matrix.As noticed, these results returned...
View ArticleXcode config scheme for composer 2015
I'm successfully using composer 2015 Fortran via command line but I'm attempting to integrate into Xcode 5 (Xcode 6.x is not working yet)I'm using one of the demo programs to test compilation and...
View Article