Question about parameter setting of PARDISO
Dear all,I use the MKL under WIndows 7 64bit + VS 2012 + MKL 11.0 or 11.2.I need to solve very large FEM matrix by PARDISO. Currently I set parameters as follows: iparm(1) = 1 ! no solver default...
View ArticleIntel MKL's call from mex-files are crashing MATLAB
Hello,this is an issue already discussed but I have not been able to find a solution. I'm using Intel MKL 10.2 and MATLAB R2011b. My C function, which has many MKL BLAS and MKL Lapack calls, runs well...
View ArticleCompiling R 3.1.1 with MKL - undefined reference to symbol '__kmpc_end@@VERSION'
Hello,I need to compile R with shared library (--enable-R-shlib) in order to use a specific IDE (rstudio) and I am running into trouble. Hopefully I'll be able to get some help from you.I previously...
View ArticleSparse Matrix
Steve:I have stumbled into this interesting problem - I need to invert a large (5000 by 5000) sparse matrix.Clearly I am going to use MKL but I have not played previously with sparse matrices. The...
View ArticleNumpy install with intel fortran on windows
Hi ! I'm trying to achieve in using f2py on windows 8, x64, with Intel Fortran Compiler.I've tried two things :1) install via windows 64x installer. http://www.lfd.uci.edu/~gohlke/pythonlibs/ NumPy is...
View ArticleLAPACKE_zgesvd has a bug
I think LAPACKE_zgesvd has a bug.When I use this function to calculate the SVD of a 2*1 matrix [1;3], the result seems to be wrong.Thanks.
View ArticlePotential error in examples/cblas/common_func.c
When compiling CBLAS examples on Intel64 I get a number of related warnings, which indicate a benign error present in the example code:common_func.c(863): warning #810: conversion from...
View ArticleClik here to view.
Calling MKL Library from C# - Intel broken links
Dear Intel:All of the links on the Intel site that provide sample of how to call MKL from C# are broken - I cannot download the files. Yes I have tried all of them.I am trying to call the PARDISO...
View Articlemkl_cspblas_dcsrgemv: heap corruption when transposing
When I'm using the mkl_cspblas_dcsrgemv function, if I request the input matrix to be transposed, the program heap gets corrupted.This happens in threaded and sequential MKL, this example program uses...
View ArticleMKL Linking
Dear Fellow Programmers:I am lost - the MKL library includes a set of routines held in the REDIST folders. Intel supplied a nice make file to compile PARDISO as called from C#. The make file I can run...
View ArticleClik 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 Article