Wrong sgesvd answer
Hi,I want to compute the SVD of two fairly small square matrices (of size 32x32 or 16x16). I use the following code to compute SVD: int iRows = m1; //=32 int iCols = m2; //=16 int INFO; char JOBU =...
View ArticleVML vcAbs precision vs. AVX
Hi,I've been writing some vectorized C++ code to compute the magnitude of a vector of 8-byte complex samples using AVX instructions on an Intel Xeon SandyBridge CPU. I recently noticed some minor...
View ArticleSubstitution time consuming in Pardiso
Dear All,I have a reactive transport model that utilizing pardiso as the solver. The model is discretized by control volume method and solved by newton iteration method. I do sybmolic factorization at...
View ArticleSymbolic factorization in Pardiso
Hi All,I use Pardiso to solve reactive transport problem. I tried to do symbolic factorization at first and then in the following newton iterations, only do numerical factorization and substitution.For...
View Articledata fitting df?Interpolate1D
Hi,I am new in the data fitting land of mkl and have a couple of questions.a) the main interpolate function, is the cell input used if provided or is recalculated ( assuming that the pointer is not...
View ArticleMemory cost and efficiency of Pardiso
Hi All,I have implemented Pardiso to solve my problem and it works fine. But the memory cost and efficiency of Pardiso is not outstanding. We have another sparse solver (Solver A, not parallelized)...
View Articleusing all cores of a CPU with MKL
hii hope you are doing good. I am using intel MKL to develop a multi-CPU version of my linear system solver. The setup is as follows:I have say 8 nodes connected via infiniband. Each node is fitted...
View ArticleClik here to view.
How to call sygvd(a, b, w [,itype] [,jobz] [,uplo] [,info])
Dear all,I want to solve the Ax = vBx by calling subroutine of call sygvd(a,b,w[,itype][,jobz][,uplo][,info]) for g95.I refer to http://software.intel.com/en-us/articles/intel-mkl-link-line-advisorand...
View ArticleSVD multithreading bug in MKL
Hi,It seems like that mkl has a serious bug in ?gesvd.First, I tried to find singular values using the Armadillo library with mkl backends, but results are different according to MKL_NUM_THREADS.for...
View Articlehow to generate random numbers over [0.0, 1.0) with MKL VSL
Hi, everyone!The basic random number generators in MKL VSL can distribute the numbers within [0.0, 1.0], but how to generate random numbers over [0.0, 1.0)?Thanks in advance!
View ArticleClik here to view.
Pardiso sparse right hand side iparm(31) = 1 (or iparm[30]=1 in C++)
Dear AllI am testing an example with sparse right hand side and sparse solution iparm[30 ]=1 in my C++.What I have noticed thatthe solution of (perm[0] = 0, b[0] !=0) and (perm[0], b[0]=0) are...
View ArticleInconsistencies with declaration of 64-bit integers ( signed and unsigned )...
There are inconsistencies with declaration of 64-bit integers ( signed and unsigned ) in many MKL functions. In mkl_types.h there is a declaration:/* MKL integer types for LP64 and ILP64 */#if...
View ArticleClik here to view.
Sparse addition
Hello!My question is regarding the addition of a sparse matrix on to itself. For e.g in the dense sense, performing an operation like A = A + B, where A and B are dense matrices.Speaking in the sparse...
View ArticleIncorrect results or crashes found in MKL 11.0 Update 3 LAPACK; Patch...
Problems have been found in an LAPACK internal function in Intel MKL 11.0 update 3 (available in Intel® Composer XE 2013 update 3) which can cause many other LAPACK functions to crash or provide...
View ArticleParallel FGMRES
Hello!What paradigm of parallel processing is FGMRES designed for? SPMD or MPMD ?Is it possible to run FGMRES in parallel in the Single Program Multiple Data paradigm?I intend to partition my matrix...
View ArticleScaLAPACK pdpotrs
I want to use pdpotrs function to solves a system of linear equations with a Cholesky factored symmetric positivedefinite matrix. Can you explain the basic step of using the ScaLAPACK function for me,...
View ArticleMore problems with Forum software
At the top of the MKL forum there is a link with the text "Click here for more information about private thread." Clicking on the link leads to a "Page Not Found" message.Is is no longer possible to...
View ArticleA question on solving many similar matrices with PARDISO
Hello,I have a quick question about solving many structurely similar sparse matrices with PARDISO.I particular, if I have many linear systems A_i x_i = b_i , where i=1,...,N and A_i have the same...
View ArticleMKL_FREEBUFFERS error
I recently purchased Visual Fortran XE Composer 2013 and got MKL_FREEBUFFERS while compiling my code.The problem did not occur when I used the earlier Visual Fortran compiler.Any help would be greatly...
View ArticleFirst solve with dense right hand side and following solve with sparse rhs
Dear AllI am trying to using pardiso (sparse direct solver) as my subdomain (substructure) solver in my domain decomposition application.The first solve with dense right hand side (rhs) coming from...
View Article