how to set the number of right hand side in phase 12 and -1 when I use...
Hi:I need to use pardiso to solve a big matrix with big right hand side and the physical memory is not big enough to store the matrix and right hand side at the same time. So I split the right hand...
View ArticleIntel Math Kernel Library for Free – requires registration, includes forum...
Intel MKL is a popular math library used by many to create fast and reliable applications in science, engineering, and finance. Do you know it is now available for free (at no cost)? The community...
View ArticlePardiso error
Hi,I am using Pardiso, from the MKL that ships with icpc 16.0.0 on Mac OSX.I have already computed a LU factorization, and I want to use Pardiso as an iterative solver using the previous LU...
View ArticleBuild R from source with Inetl MKL and ICC
OS: Fedora 22 parallel_studio_xe_2016 Hardware : 8 Thread(s) per core: 2 Vendor ID: GenuineIntel Model name: Intel(R) Core(TM) i7-2600K CPU @ 3.40GHz - Sandybridge R-3.2.2Here is my build...
View ArticleQR factorisation from LAPACK
Hi,I would like to do a QR factorization using LAPACK. From the documentation available here https://software.intel.com/en-us/node/521003#E832D468-0891-40EC-9468-686... , I've decided to use geqrf for...
View Articlecalls 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 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 Article