Equivalent function NCONF/NLQPL
hello folks!I am looking for an equivalent function to NCONF/DNCONF which has been defined in the IMSL library. Here is the description of this function:Solve a general nonlinear programming problem...
View ArticleMKL 11.2 + PARDISO
HelloJust noticed that the latest version of MKL (11.2 Update 2) has these comments in the release notes:"Added ability to free up memory used by the input matrix after the factorization step. This...
View Articlemulti-threaded matrix multiplication
Dear Intel MKL developers,I am integrating the MKL subroutine mkl_zcsrmultcsr in my MPI code. I tested a case with 16 processors, and mkl_zcrsmultcsr is called in every processor in parallel. Once it...
View ArticlePARDISO + OOC Memory Consumption
HelloI recently noticed this and wanted to post it here: when I using PARDISO solver with IC or OOC, I noticed that the solver consumes the same amount of memory in both cases (except I noticed that...
View ArticleA bug in zgelsd in MKL 15.0
Hi there, Thank you for reading this post. I got these error messages when calling zgelsd in MKL 15.0 to solve a fairly large matrix,Intel MKL INTERNAL ERROR: Condition 1 detected in function...
View Articlecan I use routine dfeast_scsrgv when the matrix B is positive semidefinite?
Hi,I want to use MKL to solve the vibration problem:Kx=λMxwhere the matrix M is positive semidefinite.can I use routine dfeast_scsrgv to solve the problem? think you
View ArticleApplication Shutdown on call to DftiCreateDescriptor for Intel MKL c++
Hello all, I have built a C++plugin to that employs MKL's fft libraries for filtering purposes. This plugin is used by a C++ compute system that generates data on demand by any other application. This...
View ArticleDo I need Pardiso_64? iparm(18) reports negative value
I am attempting to solve a large complex, symmetric matrix using PARDISO. The matrix has approximately 1,000,000 equations, and 91,000,000 non-zeros. I suspect that during Pardiso's factorization...
View ArticleClik here to view.
Inplace permutation
Hello,I need to permute a vector according to an index array jpvt returned by ?geqp3. In the example lapack/source/dgeqpfx.f it is done using an auxiliary array. I would like to do the transformation...
View ArticleBack-substitution of Pardiso not executing in parallel
All,In an iterative algorithm I'm reusing the factorization of a symmetric indefinite matrix to solve for different right-hand sides. The back-substitution is performed in a loop with a single...
View ArticleCan I use mkl lp64 vs ilp64 together in a singe program?
In my program, it has used my other libraries.I want to use some blas rountines with integer in 64bit. But I don't want to change the other parts of the program.So, I have to use blas in 64version and...
View Articledggev deterministic behaviour
Hello, I am using dggev (mkl versions 11.2.1 & 11.2.2) and i notice some non deterministic behaviour concerning eigenvalues near zero and infinite eigenvalues and eigenvectors. I use the parallel...
View ArticleDoes Intel MKL Pardiso come with selected Matrix Inversion
Hi, I am trying to compute the inverse of a sparse matrix. I saw in the pardiso website that manipulating certain parameters may result in computing the inverse of selected rows and columns of the...
View ArticleClik here to view.
ZGETRS memory corruption(?) with denormal numbers
We recently encountered a rather strange issue when passing denormal numbers in the right-hand-side matrix to ZGETRS. Attached is a small C++ file that reproduces this error.In this file, we set up a...
View Articlebug of sgemm, 2015 mkl gives wrong result, but 2013 gives correct
The compile option of 2015:/work1/soft/intel2015/composer_xe_2015.0.090/bin/intel64/ifort -i8 -openmp main.f90 -o i64_2015 -L /work1/soft/intel2015/composer_xe_2015.0.090/mkl/lib/intel64...
View Articleconversion from COO format to CSR
Dear all,I am using Intel MKL subroutine mkl_zcoocsr to convert a sparse matrix in COO form to CSR, but the results do not seem right. Here is my...
View ArticleClik here to view.
A bug in vsldConvExec
I use 64-bit AMD Linux machine. When I run the test file (attached) I get wrong results with the current MKL 11.2 update 2:EXAMPLE executing a convolution task result: z[0]=1343342.8 result:...
View ArticleClik here to view.
QR factorization without pivoting
Hello,I'm trying to compute the solution of a singular square system using a QR factorization. When I'm using a QR factorization with pivoting (i.e. ?geqp3), everything works OK. However, I'm concerned...
View ArticleClik here to view.
Strange (and seemingly wrong) result from Pardiso
I obtained strange (and seemingly wrong) result from Pardiso. I used Pardiso to solve a very ill-conditioned complex linear system arising from a circuital analysis of a three-dimensional conducting...
View Articleproblems setting MKL environment variables
Hello,I have searched the forums for a solution to my problem but haven't seen a similar topic.I have been unable to get the scripts that set the mkl environment varaibles (described here:...
View Article