PARDISO with METIS re-ordering
Hello:I am using (INTEL_MKL) PARDISO solver in a visco-elasto-plastic Stokes flow code. It is a finite difference code which uses marker-in-cell technique and has a mixed formulation with linear...
View ArticleJacobian matrix - incorrect result calculation
Hi,The 'sjacobi_rci_c.c' contained in the examples MKL library returns the error results: x f1 f20.0000 0.0000 0.00000.0000 1.#INF 1.#INF1.#INF 1.#INF 1.#INF650751168283084220000000000000.0000 1.#INF...
View ArticleMKL FFT - Real input data to comples result
Hello everyone,I am currently working on doing fft on my wav audio signal. My goal is to recreate Matlab's function called "spectrogram", which does a STFT with moving Hammington window.My first...
View ArticleRAM Performance Question
Hello AllI'm building a new workstation for solving a lot of banded matrices. I only have to solve for a single solution, and I use LU decomposition of the matrix using the functino 'cgbtrf'. The...
View ArticleClik here to view.
How to setup environment to let PARDISO use all CPU resource?
Dear all,I am trying to use PARDISO in machines with 2 X5675 CPUs. I have enabled hyperthread in CMOS so I can see 24 cores in windows task manager. I run PARDISO with parallel paramerter settings:...
View Article“NUMROC” and “DESCINIT”
HelloI am new in this Community and I have same problems with the Math Kernel Library.I hat loaded all libs and header Files from the Composer XE from Intel. The most Functions work but the Problem is...
View ArticleLearn about Intel MKL PARDISO solutions for Finite Element Analysis
Intel MKL is in the feature article of Parallel Universe Magazine this quarter!Find out the Intel MKL Pardiso solution for NEi's NASTRAN software here...
View ArticleClik here to view.
SVD changes makes least squares implementation much more fragile
I am using the linear least squares in MKL through Numpy and found that a program I have been using for quite some time is now misbehaving. I first reported the bug to Numpy (...
View ArticleUsing mkl_set_xerbla
I'm trying to set a custom xerbla function in MKL 11 update 2. I'm calling the library from c# using the code below. However, it does not seem to replace the error function and still calls the inbuilt...
View Articlecsrcoo garbles unsorted complex16 matrices
While using csrcoo to convert an unsorted (rows and colums) complex16 (double precision complex) matrix from COO to CSR using the option to sort the columns in the output, I discovered the array of...
View ArticleIntel MKL PARDISO seg. faults within MPI codes
Dear forum users,my MPI codes are linked with Intel MKL, and intensively use PARDISO sparse direct solver. These MPI codes are actually MPMD in the sense that the last MPI task in the communicator does...
View Articlecsrcoo garbles unsorted complex16 matrices
While using csrcoo to convert an unsorted (rows and colums) complex16 (double precision complex) matrix from COO to CSR using the option to sort the columns in the output, I discovered the array of...
View ArticleFast Small Dense Matrix Solver
I have a general square dense matrix A (not symmetric) which is formed by A=PTBP where B was in a compressed storage scheme and P is a rectangular matrix. The size of A ranges from 10x10 to 500x500,...
View ArticleSubroutine Call and LAPACK Utility Function
Hi,I'm trying to run a very simple code that finds machine epsilon with DLAMCH from LAPACK utility function and calls some random subroutine afterward. When I try to run this code, the error I get...
View ArticleVML Rounding Functions v?Modf v?Round questions
MKL reference Manual says 'Halfway values, such as 0.5, -1.5, are rounded off towards even values. An inexact result exception is raised for each changed value.' for v?Modf (VML) function.I tried it...
View Articlerandom number generators
Hi,Why the first element of the vector is generated from the distribution with the mean value equal to 0.0. How should I correct this code:void RNorm(int n ,float *x ,float mu=100.0f, float...
View Articleneed the source code or exe binary.
I want to get those source code ro binary on cpu and mic. show as in http://software.intel.com/en-us/intel-mklBenchmarks--> intel xeon phi Cooprocessor.I want to benchmark them on mic.Who can tell...
View ArticlePardiso performance 10.1.1 vs 11.0.2
We have recently upgraded to version 11.0.2 from version 10.1.1. We are finding that pardiso is typically around 10% slower in 11.0.2 than it was previously. The matrix is not too large for this...
View ArticleSparse matrix - vector product
I use the MKL_ZCSRSYMV routine in a FORTRAN code. The sparse matrix (a) is very large and then I need that indexing vectors (ia and ja) where INTEGER*8 instead of INTEGER.How can I do?
View ArticleMultiple Core Usage Issues
I recompiled my installation of numpy and scipy agaisnt Intel MKL. I am trying to speed up a script that fits tensor for DT MRI. The script bottlenecks during the svd operation, in particular the call...
View Article