3D convolution
I am using MKL's VSL multi-dimensional convolution routines to perform a 3D convolution. As a test, I am using the input function as a 3D gaussian as follows:f(x,y,z) = exp(- (x-2)**2) exp(-(y+2)**2)...
View ArticleSparse Matrix storage to compressed sparse row format
Hi :-)I am new here and I have a question about using "PARDISO - Parallel Direct Sparse Solver Interface". As you may know in order to use this subroutine for solving a set sparse linear equations as...
View Articleblacs for sgi mpt
I need to compile a software package which requires blacs, however since about composer_xe_2013.5.192 blacs seems to have disappeared from MKL (except blacs for intelmpi and MIC). Where do I get the...
View ArticleScalapack on Mac OS X
What are the chances of seeing scalapack in mkl for Mac OS X in near future?
View ArticleProblem with robust estimation of a covariance matrix.
Hi I have the matrix "x" and I want to compute the covariance matrix. The i column of the matrix stores the observations of the i variable.The matrix is 0.8147 0.9058 0.1270 0.9134...
View Article3D convolution
I am using MKL's VSL multi-dimensional convolution routines to perform a 3D convolution. As a test, I am using the input function as a 3D gaussian as follows:f(x,y,z) = exp(- (x-2)**2) exp(-(y+2)**2)...
View Articleproblem in computation of 2D FFT using two seperate real and imaginary arrays
Hi,i want to compute FFT of a complex 2D array. So first i have tried it in 2 ways.1) Directly using complex 2D array(real and imaginary interleaved) .2) 2 seperate arrays(where real and imaginary are...
View Articlevsrnguniform
Dear all,Can vsrnguniform generate random multidimensional arrays or it just works for 1D arrays? I guess that one solution could be to generate a 1D array and after reshape it. However, I would like...
View ArticleProblem computation of1D FFT
I want to compute a 1d fft with dfti.Here is my code #include "mkl_dfti.h" #include <stdio.h> int main() { float y[8]; for (int i=0; i<4; i++) { y[i] = i+1; } DFTI_DESCRIPTOR_HANDLE handle;...
View ArticleClik here to view.
Problem with dstebz and dormtr (dense eigenvalue solver)
Hello,I'm trying to solve a symmetric generalized eigenvalue problem with the BLAS Fortran bindings inside C and I'm facing two problem with the following piece of code:if I want all eigenvalues...
View ArticleProblem with scale factor in 2d fft
I am using mkl to do a 2d image fft. First, I used mkl to compute forward the 2d fft of the initial image and I computed backward. I think it should get a image where every pixel gray scale is scaled...
View ArticleClik here to view.
Different matrix multiplication results between C and Matlab.
Let covold be a matrix of size 72x72, theta be a vector of size 1x72 and meanold a vector of size 1x72. I want to compute [(w-1)/w ]* covold + [ 1/(1+w) ] * (theta - meanold)' * (theta - meanold)...
View ArticleMKL v11.1.2 link failure with scalapack
Hi,After updating from MKL 11.0.5 to 11.1.2, I keep getting link failures on mkl_scalapack_ilp64.lib. After some searching, I see that for 11.0.5, all mkl library files are in C:\Program Files...
View ArticleSpeed up the PARDISO solver
Hi,I have a problem while using the PARDISO in FORTRAN to solve a large sparse system. The processor of my computer is Intel i7-3770 CPU 3.40GHz. I tried to call the MKL_SET_NUM_THREADS() and the...
View ArticleBLACS context value and multiple MPI communicators
Hello everybody, I am trying to make some tests with BLACS/Scalapack (C interface, Intel MKL version) making use multiple MPI communicators: particularly, what I am trying to obtain is a set of BLACS...
View ArticleCopy a submatrix to another matrix.
Hi! I want to copy the first 2 columns of a matrix A to the last two columns of another matrix B. A and B is matrices of size 3x3. I found that this can be done with the function p?lacpy ( I don't...
View ArticleHow prevent termination due to Zero Pivot?
Hi,I use DSS interface to solve a sparse symmetric linear system. When a matrix is singular, the code is terminated at dss_factor_real() with a message "MKL-DSS-DSS-Error, Zero Pivot detected". It...
View ArticleParallel algorithm used by mkl_?csrmultcsr
I am wondering which algorithm and parallelization strategy is utilized in mkl_?csrmultcsr? Does this function scale well on Intel Xeon Phi architecture?
View ArticleSimple task from example failing to execute
Hi all,I' ve got a simple example from http://software.intel.com/sites/products/documentation/doclib/mkl_sa/11/mklman/index.htmwhich works fine.I wrote the pretty similar code #include <stdio.h>...
View ArticleClik here to view.
Excessive memory allocation in DftiCommitDescriptor()
Hi everyone!When designing your project, I was faced with the problem of excessive memory allocation. In the course of the analysis found that the function DftiCommitDescriptor in the implementation...
View Article