Problem in subsequent calls to PARDISO MKL
I am using INTEL MKL PARDISO to solve Laplace equation. However, the Matrix becomes more sparse every time I call PARDISO due to my physical problem. PARDISO works perfectly fine until it fails at the...
View ArticleSstrange behavior MKL-FEAST
At the moment I am trying to implement a minimal working example of the <<dfeast_scsrgv>> routine in c++ of the FEAST package to solve the generalized eigenvalue problem Ax=lBx. The...
View ArticleProblem with FGMRES method
Dear All,I have been programming in fortran using RCI FGMRES to solve large sparse linear system. However, when the system size is very large, say, 10^5, the program hold on the dfgmres_init and never...
View ArticleHow to simulate several logical nodes on one phyisical node (MPI FFT)?
I'm trying to run my fftw3 mpi program compiled with icc compiler (wrapped within mpicc wrapper) on my single computer. When I was running regular (openmpi) mpirun I was able to run it with arbitrary...
View ArticleProblem in calling MKL Trust Region routine from Abaqus user subrotuine
Hello,I want to use the Trust Region solver of MKL in Abaqus user subroutine. I have a problem in using dtrnlsp_solve. As a test, I took an example from the manual (Example 14-1. dtrnlsp Usage in...
View ArticleClik here to view.
daxpby doesn't scale
I've discovered a strange performance problem: replacing a daxpy() call with daxpby() significantly degrades program scalability. I've written a simple test case (hopefully attached) which shows that,...
View ArticleHow to shift zero frequency with 2d fftw data
Hi,gays! I do SFFTW_PLAN_DFT_R2C_2D(FFTW_PLAN,nx,ny,density,fft,FFTW_ESTIMATE) and then get a complex array fft(ny/2+1,nx)! I try to shift the zero frequency to the center but failed! Who can tell me...
View Articledgemm: slow performance
I am using dgemm function from Math Kernel Library in my FORTRAN source code to do the following matrix calculation: X = A - transpose(B) *CA is 400*400 dense matrix, B is 10000*400 sparse matrix and...
View ArticleTrust Region solver memory corruption
Hi there,I'm trying to use MKL to replace existing Levenberg-Marquardt fitting using MINPACK or MPFit. As a test case, I'm fitting a Gaussian function. Things work pretty well, except that my program...
View ArticleClik here to view.
Scalapack BUG in PZUNGQR (Again)
Hi,Previously, we reported a possible Scalapack bug in the PZUNGQR function http://software.intel.com/en-us/forums/topic/473803That issue has still not been resolved, but it was stated that...
View ArticlePardiso slow performance with multiple RHS
I am solving K*D=F equation with pardiso sparse solver. K is a 11000*11000 sparse symmetric matrix and F is also sparse matrix with 11000*363 dimensions. On a dell computer with Intel Core 2 Duo CPU...
View ArticleLower triangular matrix of a symmetric sparse matrix
Hi,Is there any way to get lower (or upper) triangular matrix of a symmetric sparse matrix from pardiso solver? If not, what is the best way to do that using MKL intel in FORTRAN? I know I can get...
View ArticleGlibc error for dgesvd
Hi, I am using MKL to compute SVD of a real matrix. I use C++ code with std::vectors. Here is a portion of my code related to SVDvoid Decompose_InterpolMatrix_To_SVD(int rows, int columns,...
View ArticleProblem with MKL dpotrf function
Can anyone see why I am getting all zeroes back in [A] when calling the MKL dpotrf function below? Info comes back as 0 which indicates success, but A is just full of zeroes. What am I missing? I have...
View ArticleProblem with compiling 64 bit Cluster MKL FFT program
I've been trying to compile cluster version of MKL FFT (MPI version). I've followed link adviser for 64-bit architecture and used this line:mpicc -o fftmkl_mpi.out fftmkl_mpi.c -Wl,--start-group...
View Articlenlsq solver with linear (bound) constraints generates out of boundary guesses
Hello everybody,This is my second post with a set of questions about using the nonlinear least square solver from mkl (the first one is about the OMP parallelization is here...
View ArticleInterrupt a thread with a MKL function
Hi,How can I interrupt a thread that is performing a MKL subroutine? My threads are created with boost::thread, and normally I use interruptions to stop my threads...
View ArticleProblem with Intel Fortran error #6633
I am trying to use the Pardiso matrix routines but I keep getting an error #6633 “The type of the actual argument differs from the type of the dummy” with the PT variable when I compile. The relevant...
View ArticleMKL: Partial factorization and Schur Complement?
Hi,I have a large symmetric sparse matrix (10000*10000) and I want to reduce the matrix to condense form of last N elements. In order to do that I need to calculate the schur complement form of the...
View ArticleBug in *GESVD when using MPI
Hi,Attached is a test case and screenshots where ZGESVD is sometimes computing the wrong lwork value when using MPI. This is the regular Lapack ZGESVD and not the Scalapack PZGESVD. I have not...
View Article