Random Number Generator
Hello, I am trying to generate a vector of random numbers using the following code. However, every time I execute the code, the output vector is the same as previous run! Can you please help? Thank you...
View Articlespffrt2 issue
I am trying to use mkl_dspffrt2, but it does not give a correct result. I used packed storage as recommended. The routine does not work for n=1. Let's say the matrix is A=[10]. Have anyone used this...
View ArticleHow to configure libs when trying to use Inspector-executor Sparse BLAS...
Here is my environment : Visual Studio Community 2017 + Parallel Studio XE 2019 update 1I`ve tried Link line advisor but it failed with 'No symbolic file loaded for mkl_avx2.dll'And the current libs...
View ArticleFinding inverse of a binary matrix by using LAPACKE_dgetrf and LAPACKE_dgetri
I'm trying to find out inverse of the following binary matrix by using following functions,lapack_int LAPACKE_dgetrf (int matrix_layout , lapack_int m , lapack_int n , double * a , lapack_int lda ,...
View ArticleClik here to view.
Different computation results on different processors using mkl_cbwr_set(int)
There are two desktop PC with the following processors: Intel Core i5 4570 (Haswell) и Intel Core i5 3330 (Ivy Bridge).Use MKL 2019.0.1, build 20180928.The app use PARDISO with the following settings:...
View ArticleClik here to view.
Pardiso Optimization
Hello,i am currently solving a relatively small system of equations (3000-10000 equations), but since i have to solve it millions of times, each time with slightly different system matrices and right...
View ArticleLAPACKE_dsyevr doesn't finish but LAPACKE_dsyev works well
I have a symmetrc matrix with some parameters to be changed. For some parameters, the LAPACKE_dsyevr does not return for a very long time, so I have to finish the program.If I change to LAPACKE_dsyev,...
View ArticleSignature of v?Fdim, v?Fmax, ...
Hello,why is the second array of the functions v?Fdim, v?Fmax, v?Fmin etc. not constant?The signature in the file mkl_vml_functions.h is:_Mkl_Api(void, vdFdim, (const MKL_INT n, const double a[],...
View ArticleQuestion about the allocate procedure
Hi!I`m wondering for what reason can it be that for an allocate statement likeif(.not. allocated(a)) thenallocate(a(size),stat=stat,errmsg=errmsg)else if (size(a) /= size)...
View ArticleClik here to view.
CPardiso phase 33 scaling
Hi,We want to use Cluster Pardiso for our finite element application. To get an estimate of performances, we used a simple code (attached file) to read a matrix from the Sparse Suite Collection (Matrix...
View ArticleBlock sparse matrix with sparse matrix
I have a sparse matrix M, and a sparse matrix K, both of the sparse matrices obtain the same non-zeros location, how to form a sparse matrix with following format:A= [ M, 0.2K; 0.1M, 0.1K];This was...
View ArticleThreading Problem in Pardiso
Hello,We meet a problem in using Pardiso with multiple threads. Our purpose is solving: AX=B, X is a non-symmetry sparse matrix, A & B are vector. What we do is initialize the Pardiso, LU...
View Articleundefined symbol: mkl_serv_check_ptr_and_warn when calling function from...
I'm calling a C function from within my Python code using the ctypes Python library. The C function calls a MKL function. This however gives me an error:symbol lookup error:...
View ArticleHelp on running GMRES with MKL in C++
Hi,I am new to MKL library. I tried to use to use it to solve a sparse system matrix using “GMRES” functions in MKL library.I tried the steps in the given in examples which come with the MKL library (...
View ArticleBatch normalization example
Hello,This is my first program using MKL library and i want to include a simple batch normalization call after convolution. This is my code for BN Does anyone have an idea what i did wrong and why the...
View ArticleProblem with random number generator ?
This isnt a BIG problem, but it might trip someone up if they are not aware of it.When I call RANDOM_NUMBER(XX) where xx is real(8),the first value is always a very small number, typically under...
View ArticleMKL and CMAKE
Is universal cmake script for finding and linking MKL exist?For example, I want to create solution which uses cmake and MKL:On the first machine I want to compile my application with Intel C++ Compiler...
View ArticleClik here to view.
Why VM library is so slow in the new mode of CPUs
Hi, I had reached a special problem and found that the calculation which using for ... loop is more efficiency than use the VM library of MKL. I test for several examples. It shows the same kind of...
View ArticleHow to implement numpy broadcast mechanism with mkl?
How to implement numpy broadcast mechanism with mkl? I have been confused, how to use mkl to efficiently implement the broadcast mechanism in numpy ((Element wise operator "+","-","*")?such as2-D...
View ArticleHow to implement numpy broadcast mechanism with mkl?
How to implement numpy broadcast mechanism with mkl? I have been confused, how to use mkl to efficiently implement the broadcast mechanism in numpy ((Element wise operator "+","-","*")?such as2-D...
View Article