mkl.fi interfaces are missing some dummy argument definitions
Some of the interfaces brought in by mkl.fi are missing some argument definitions. While not an problem, since the implicit typing takes over, it does generate compiler...
View ArticleConcurrency Issues with DSS
I have a very large program which does numerical modeling. Up until recently I've been debugging and everything has been fine. I'm now switching over to optimize areas that need improvements. In my...
View ArticleSIGSEGV with Poisson Solver
I have been running and debugging a code with the Intel visual fortran composer that uses the MKL poisson solver. It runs without issues in Windows. But when I run it in a UNIX system I get a SIGSEGV...
View ArticleRuntime error : cannot locate entry point on mkl_serv_print_verbose_info ?
Hello everyone, I used VS 2013 + IVF 2013 + MSMPI + PetSC + microsoft HPC those environment to build a fortran multithreaded console app.The OpenMP conditional compilation = yes, use intel math kernel...
View ArticleUnderflow exception when using vmsExp
The following code generates an exception, although I set the error mode to VML_ERRMODE_IGNORE. Is there a way to use high accuracy (VML_HA) computations without raising underflow exceptions when...
View ArticlePardiso always crashing on Linux
I consistently have Pardiso crash when running on Linux. I've trimmed it down to a relatively simple test case of building a Poisson matrix and solving it with Pardiso. The same program runs fine on...
View ArticleClik here to view.
Pardiso (direct/iterative)
Hi,I'm having a problem using pardiso solver in MKL library in direct or iterative mode. The manual says to choose between direct or iterative mode simply set the variable SOLVER (0=direct,...
View ArticleSOUP
Martin Gardner provided this rather interesting genetic program in his writing in Sci Amer in the mid 80's. I copied his program in basic and ran it on GWBASIC and it ran nicely without a problemAt...
View ArticleTrouble linking to BLAS library with fortran compiler
So I am trying to get accustomed to linking my f90 code to the MKL library and I decided to start with the $dot function. I have included my compile command and my code below; please note that I set...
View ArticleUnexpected EOF while reading module lapack95
Hello,I am trying to run examples in $(MKLROOT)/examples directory (specifically examples_f95.tgz). I have problem with lapack.As the user guide suggest if I want to use gfortran compiler I should...
View ArticleGfortran compiler linking help LAPACK95
Dear Mr. / Mrs. ,I'm having difficulty compiling code that uses the DGETRF and DGETRI Fortran functions. Below I placed a very small test example "test.f03", that illustrates my problem:program main...
View ArticleCustom Dynamic-link Library Builder and MKL_DIRECT_CALL
I'd like to set the MKL_DIRECT_CALL flag for the custom MKL dll that we are using from a C application and hopefully speed up the matrix operations. Will it be enough to set the flag for the C...
View ArticleClik here to view.
Pardiso problem
Hi,I'm having a problem using pardiso solver in MKL library in direct or iterative mode. The manual says to choose between direct or iterative mode simply set the variable SOLVER (0=direct,...
View ArticleBug in 64-bit mkl_dcsrmv
Hello,I need sparse matrix-vector multiplication in my code and want to use the MKL, as I want to also use the Xeon Phi cards in my system. However, the 64-bit version (the one activated by the...
View ArticleCluster pardiso gave wrong results when np>2
I'm trying to test the cluster pardiso to solver complex linear equations with Hermitian matrix. When I ran the parallel program using -np 2 option, everything looks fine. However, when I used more...
View ArticleHow Large is a Large Matrix ?
I was looking through the LAPACK documentation and came across various instances where there are different routines recommended for large and small matrices. For eg., SYEVD vs SYEVR (...
View ArticleErratic behavior of cluster_sparse_solver
The code attached below will either crash (access violation error) or produce the right result, or an approximate result depending on the number of MPI processes it runs on. This code is run on a Win 7...
View Articlelapack95 errors
I am trying to use an eigenvalue routine through the lapack95 interface of the mkl. But I end up with this:$ ifort Eigen.f90 -mkl...
View ArticlevdExp 32 bit vs 64 bit different results
Hi,I'm new to MKL so please bare with me. We're moving a project from 32 to 64 bit and we've encountered some inconsistencies between 32 & 64 bit vdExp.I have the following codedouble value =...
View ArticleLAPACKE_dgeev result differences on different CPU architectures
Hi, I am using MKL 11.1 to find eigen vectors with and am having issues with the results changing sign depending on which machine I run this. This gives completely different results in the end for an...
View Article