calls to pardiso from a DLL crash
I am a new user of MKL, and I'm using pardiso to factor a real nonsymmetric matrix. I've successfully run the sample program for this case, and I'm trying to do exactly the same thing, except my calls...
View Articlequestion about pardiso arguments at solving phase
I want to call pardiso with phase=33 to compute a solution to a given rhs. I'm sure I need to pass it the pt(:) array from the factoring phase, but do I also need to pass it any other arrays used in...
View Articlesystem crash when running dss on MKL_dss_indefinite
Hey there,my dss implementation freezes the whole system, probably due running out of memory, when switching from positive_definite to indefinite.The matrix provided to dss is 4x4, so not a big...
View ArticleAbout DFTI_BACKWARD_SCRAMBLED option usage
Hello Sir or Madam,Recently I am leaning and using MKL FFT in my work. I noticed that there is a option DFTI_BACKWARD_SCRAMBLED can applying an FFT to input data whose order is scrambled, or allowing...
View ArticleClik here to view.
May be a bug about Data Fitting module...
Hi all,I got a problem while using data fitting module in MKL (I am using version 11.3.0 Build 20150730), here is the sample code.#include <iostream> #ifdef __cplusplus extern "C" { #endif //...
View ArticleClik here to view.
Bug in xggsvd? "Parameter 5 was incorrect on entry to DORMQR"
When running the program in single threaded mode I get the error message "Parameter 5 was incorrect on entry to DORMQR". The message is not present if I either run the program with multiple threads or...
View ArticleSparse conversion routines; convert in place ?
I would like to convert from COO to CSR without making a new copy of my matrix. Is it OK to have actual arguments point to the same arrays, i.e. acsr and acoo are the same array and similar for the...
View ArticleProblem building numpy and scipy with MKL
After building numpy when I try to run the tests I get the error: undefined symbol: _intel_fast_memsetI'm not sure where that is coming from, could someone elaborate a bit? I followed the instructions...
View Articlewhat's the superb argument used for in the LAPACKE_dgesvd function?
I don't understand the documentation here (https://software.intel.com/en-us/node/521150); when I print it, it's all zeros. What's it useful for?
View Articlepardiso error in linear system solution compared to gaussian elimination
Dear IntelI have a mechanics problem which I solve with finite differences, and the matrix is sparse. I use pardiso to solve the linear system, and I regularly get solutions which do not match Gaussian...
View Articledtrnlsp_solve spinning/sleeping when called from multiple threads
I am using the trust region solver in MKL and having issues where dtrnlsp_solve takes significantly longer to complete. I have got many threads that all need to run an optimization using the trust...
View ArticleKronecker products
Hi all! I'm trying to find an efficient way to compute Kronecker products (of matrices) using Intel Fortran MKL. I was hoping to find a routine that directly does the Kronecker product, but couldn't...
View ArticleScipy 15.1 not linked properly when compiled against MKL
Hello,I have been trying to compile numpy and scipy against intel MKL. I have been following the following guide to the letter:https://software.intel.com/en-us/articles/numpyscipy-with-intel-mklI am...
View ArticleMKL and Inspector 2016 and uniniatialized reads in syrk
Using Inspector 2016 on Windows I can see there are uniniatialized reads dsyrk occasionally. When I replace the syrk call by the equivalent loops then it goes away.To be fair in both cases the results...
View ArticleCSR matrix handling and solution error
Hi, I built my own class for CSR3 matrices for Mkl Pardiso and I have to know how Mkl knows the length of arrays (values and columnIndex ): - Is it done by rowIndex[number_of_rows]? - Is it done by...
View Articlecan't figure out mkl_zcsrcoo
I have been able to use mkl_dscrcoo, but I can't figure out mkl_zcsrcoo. PROGRAM zcsrcoo_test IMPLICIT NONE include 'mkl_pardiso.f77' integer nx, ne parameter(nx=5,...
View ArticleWhat is the memory management strategy in sparse_matrix_t ?
Hi all,I have two questions about memory management of sparse_matrix_t :1. When I create a matrix handle using routine mkl_sparse_?_create_csr, the handle will copy all data to new memory blocks? Or it...
View ArticleDoes MKL support real*16 datatype?
I am trying to convert some old Cray routines with double precision to the Windows Visual Fortran platform, but somehow am having difficulties in making the routines work. Just wondering if MKL have...
View Articledifference result of mkl_sparse_s_export_csr in two subroutines
Hi all,I am using mkl_sparse_s_export_csr function for reading data in sparse_matrix_t object, but the return result is error while in other function. Here is the example:#include "mkl_spblas.h"...
View Articledifferent results for LAPACKE_dgels between MKL and LAPACKE
Hi guys,I wrote some code to do dgels originally using LAPACKE. I also tried to use MKL to do the dgels. From my understanding, with the exact same code base but just recompile/relink with MKL library...
View Article