Eigan Value Solution (LAPACKE_sggev)
LAPACKE_sggev gives the Eigan Value in my codeThe function is not giving correct answer when used on multiple nodes. When nodes = 1 it gives a correct answerAlso, if LAPACKE_sggev is not compatible...
View ArticleVery large sparse matrix solving problem
Hi dear Intel Now I'm using 'MKL2017 update 1' and 'MPICH3.1.4'. And I have 2 machines with 512GB memory in each.When I tired to solve the SPD sparse matrix having 488 and 1500 million elements with...
View ArticleTwo questions about 'iparm[1]=10'&'cluster_sparse_solver' speed
Hi dear Intel I'm the user of 'MKL2017 update 1' and 'MPICH3.1.4'. Now a days, I tried to solve the large SPD sparse matrix containing on the order 10^8 rows. Therefore, I'm troubled with reducing the...
View ArticleProblem creating sparse matrix handle coordinate format
I have a problem creating a sparse matrix handle in coordinate format. I have a problem to make this function work in my code. I made a very simple example. I define a 10x10 identity matrix. My rows...
View Articleruntime dlls required
Hi,I would like to know which runtime dlls are required to use the vslNewStream, vdRngUniform and vslDeleteStream functions of the MKL. My application (dll) loads some of them manually (using the...
View ArticleLocating Pardiso Example files
As mentioned in MKL Pardiso manual there are example files located in "examples/solverf/source".But for me I cannot locate those files. I would appreciate if you could tell me why they are missing?I...
View ArticleIntel MKL Extended Eigensolvers: SIGSEGV, segmentation fault occurred
Hi All,I am trying to use Intel MKL Extended Eigensolver with intel/17.0.1.132 compiler on nersc/cori supercomputer. I am compiling with: ftn -qopenmp -mkl diagonalise_sparse.f90 read_write_mtx.f90 -o...
View ArticleProblem with nested parallelism using OpenMP and MKL fft
Hi folks.I had a problem with nested parallelism using OpenMP and MKl fft.The part of my code that uses MKL calls within omp parallel region looks like this:...
View ArticleSparse index base for export csr
I do not understand how is working the sparse_index_base in mkl_sparse_d_export_csr. It looks to me that the exported arrays are independant of what is chosen for sparse_index_base parameter. To...
View ArticleAll eigevalues and eigenvectors of large sparse symmetric matrix
Hi Everyone,In advance I am sorry for my noobish question. I am a physics PHD student and basically I use python for my math/physics problems. But now I have a problem which requires more computing...
View ArticleCan I free original matrix after Pardiso factorization
Hi, I am wondering when can I release the original matrix memory (*a, *ia, *ja)? Since pardiso will keep the factor matrix internally, I think the original matrix memory can be release after...
View ArticleNewbie questions on Pardiso : allocate and initialize right-hand-side B and...
Hello,Newbie to MKL/Pardiso. I was running my first program yesterday. Some quick questions on the right hand side B and solution X: - is it correct that both B and X have to be allocated in every MPI...
View Articlepardiso performance vs Mumps
Hello,Running my first Pardiso (cluster) program, and benchmark Mumps. I was expecting that pardiso would be faster or at least close enough but the result is not very encouraging.My environment-...
View ArticleInterpolate1D "in place" allowed ?
Hi everyone,I notice that df?Interpolate1D and df?InterpolateEx1D have different behaviour when working InPlace (&site == &r) : InterpolateEx1D works, Interpolate1D does not.Here is my sample...
View ArticleNested parallelisation problem OMP + MKL
I am attempting to parallelise calls to mkl within a parallel omp region to test whether or not the code executes faster. Simply parallelising part of the code does not yield linear increase in...
View ArticleSORMRQ fails erroneously when LDA < M
I've encountered a bug after updating from a fairly old (~2013) version of MKL to the latest version that SORMRQ now fails if LDA < M, even though the only stated constraint is LDA >= K.It seems...
View ArticleHow improve efficiency using Pardiso in an iterative problem?
I am trying to improve the efficiency of my code but I am new to Pardiso and there are many things left as a mystery to me.I have to solve many simultaneous equations that change a little bit in every...
View ArticlePardiso only uses 2 cpu cores out of 4
I am running a Pardiso solver but it only uses two cores out of 4 cores.I read the documentations but I cannot find any solution.I set windows environmental variables MKL_NUM_THREADS = 4 and...
View ArticleProblem getting mkl_scsrsv function to work.
Hi, I am just trying to solve a simple linear equation with this test program. I am using MKL 2017 with Microsoft Visual Studio and have tested other functions with the array matdescra with them...
View ArticleIs there a parallel implementation of mkl_sparse_d_trsv()?
Hello, I'm using function mkl_sparse_d_trsv(). No matter how many threads I use (environment variable MKL_NUM_THREADS), the performance is always the same. So it seems to me that this function has no...
View Article