How to link Calculix with Pardiso using Intel MKL
Hi,I was wondering if it is possible to link Calculix v2.7 (an open source application) with Pardiso using Intel MKL somehow. I am a complete beginner in this field so I was wondering if it is possible...
View Article1D convolution of a 3D array using Intel MKL
I have a 3D array which is stored in a columnwise fashion.for( int k = 0; k < nTop; k++ ) // Loop through the tops. for( int j = 0; j < nCol; j++ ) // Loop through the columns. for( int i = 0; i...
View Articlematrix multiplication speedup
Hi, I'm using cblas_dgemm to calculate matrix multiplication. For random generated matrix X of size N * N (N could be 100), I calculate Y = X^T * X. (X^T is the tranpose of X). I can do it in two...
View Articleexample code of Direct Sparse Solver (DSS) Interface gives wrong result
I am trying to use DSS routine to solver linear equation with a sparse matrix. I found the example code under the intel compiler directory nameddss_sym_f90.f90 and compiled it as ifort...
View ArticleLinking error
Greetings,I'm trying to compile and link a use-subroutine for a commercial FE code. In my code, I have following lapack calls: use lapack95 call potrf(cmMinusCfInv) call potri(cmMinusCfInv)To compile,...
View ArticleUnitialized issue in cblas_dsyrk
I have created a very small test program that illustrates an an issue in syrk. In fact when I run valgrind on my program I get.==44033====44033== Conditional jump or move depends on uninitialised...
View ArticleAfter making static linking using \MT my Project still have dependency on...
I am using IntelMKL in my project. I use 'PARDISO' API from IntelMKL.For parallel processing, I have made changes in the project settings as follows, "Configuration Properties => Intel...
View ArticleDGESVD/DGESDD computation complexity
I hope to know that the real computation complexity of DGESVD and DGESDD functions of mkl for an N by N matrix. If the complexity can be written as T(N) = C1*N^3 + C2*N^2, I hope to know the values of...
View ArticleMKL PARDISO pivot function fault
HI all,I tried to implement mkl_pardiso_pivot() in my program.Sometime, it will cause a segmentation fault in a function in mkl_intel_thread.so.Does anyone experience the same situation?
View ArticleClik here to view.
Issue during replacing ipp DCT function with MKL DCT function
Hi,I want to replace my IPP based DCT function with mkl based DCT function .I am getting different output data when I will cross check with the ipp DCT vs mkl DCT function output.I used below functions...
View ArticleC++ 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 Article