3D Convolution by Using Intel MKL
There is an example for using Intel MKL to perform FFT in 3D, as posted here. It would be absolutely helpful to know how MKL FFT libraries can be used to compute 3D convolution (in a serial manner) for...
View Article1D FFT of a 3D array
My goal is to compute 1D FFT of a 3D array along all its dimensions. This 3D array is stored as a 1D array in a columnwise fashion. For example,for( int k = 0; k < nk; k++ ) // Loop through the...
View ArticleExtended Eigensolver Routines: Strange eigenvectors
Hello,I am currently trying to solve the eigenproblem Av = λv where A is a complex hermitian matrix using the MKL feast implementation. For testing purposes I constructed the following example (using...
View Articleerror #6284: There is no matching specific function for this generic...
Hi! I have the questions:I want to use the FFT, but now the errors appeared, I don't know how to deal whit them .Can you tell me the reasons of the errors,and tell me how to deal with them.Thank you !
View ArticleFast Discrete Fourier Transform with MKL
Hi all,In R programming, there is the "fft" function: https://stat.ethz.ch/R-manual/R-devel/library/stats/html/fft.html> x <- c(102, 55, 89, 12, 3, 45, 9)> fft(x) [1] 315.00000+ 0.00000i...
View Articlemkldcsrbsr gives core dump when trying to fill all arrays
Hi I am using the mkldcsrbsr routine. I ma using MKL version 10.3 (i have tried with version 11.0 as well). It seems i can call the routine perfectly for the job type = -1. I get the number of blocks...
View ArticleIncosistent in-place and out-of-place results in mkl_dfti
Hi,I am trying to transfer from out-of-place calculation to in-place using MKL 10.3.12. I don't see any problem when doing forward FFT. However in backward FFT, for dimensions larger than 4, I get...
View ArticleHow to configure properly for cluster PARDISO in WinServer 2012?
We are running an application with (currently) the non-cluster PARDISO under Windows Server 2012 SP1, the machine is a HP DL980 with 160 logical cores (80 cores with HT turned on) and 2TB RAMs. We are...
View Article-mkl gives different results from -lfftw3
I am developing a time-stepping code that calls fft routines in every step. While writing and testing the code, I used the -lfftw3 flag to link to the fftw3 library. Now that the code is functional, I...
View Articleabout parallelism on BLAS level-1 routines and VML
Hi all,I am running BLAS routines in MKL with intel compiler (icpc). Following the example given in the compiler, I try to set the numbers of threads from 1 to 10 while running dgemm routine for...
View Articleoccured linking errors while using CLUSTER_SPARSE_SOLVER
Hello Everyone, While using 'CLUSTER_SPARSE_SOLVER' for solving sparse matrix, I got linking errors. I have included "mkl_cluster_sparse_solver.h" as well as "mpi.h" files also. What should I...
View Article[REQUEST] Looking for old MKL Versions
Hi,I'm searching for some older MKL versions, but am having trouble locating them. Does anyone have (or know where I can find) some/any of these? Thanks in advancel_mkl_p_8.0.019.tgz...
View ArticleFFT-Based 3D Convolution With Zero Padding
I have been trying to figure out how I can use Intel MKL to perform a FFT-based 3D convolution with zero-padding. I have been searching and posting in online forums (including Intel MKL forum),...
View Articletrsm memory leak
Hello,When runing valgrind on "source/cblas_dtrsmx.out" compiled from within "mkl/examples/cblas", I end up with a log similar to:==9740== 3,131,264 bytes in 1 blocks are still reachable in loss record...
View Article[Scalapack] Please Help with using pdgesv
Hello all: I'm trying to solve a linear system (9 by 9 full matrix) by pdgesv in c. I use the example code (http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=2&t=1683&sid=26b4f253...) and...
View Articlecomplex system GETRF+GETRS
Dear all,I would like to solve a linear complex system with MKL libraries. As I have done with real system I use GETRF with GETRS. The MKL reference says that I can also use getrs also for comple...
View ArticleSparse-Sparse Matrix Multiplication
HiI have used mkl_dcsrmultcsr in my research. However, it is performing double pass to compute sparse*sparse matrix product. For small size problems, this is not a problem, however for large size...
View ArticleDifferent results using 11.2.1 vs 11.1.0
I have the identical SPD matrix and rhs but got the different solution using mkl 11.2.1 on win vs 11.1.0 on linux.on winMajor version: 11 Minor version: 2 Update version: 1 Product status: Product...
View ArticleDiverge in Newton Method
Hi,I encounter a diverge problem when using Intel MKL PARDISO to solve a transient simulation by Newton method.This case is a unsymmetrical matrix with size is 64,000 and 1,721,082 non-zeros (L+U).I...
View ArticleBasic Code of Using MKL FFT on MIC
Hello,In our local cluster we have a bunch of MIC's, which are (almost) never esed. I would like to give it the try, but I have no experience iwth the MKL library or MIC. My programs are very simple...
View Article