How to compile the sample code in minGW complier?
Hi all,I am using minGW as the compiler and our code has a lot of FFT. In the past we use fftw and we want to try Intel's fft to improve it. I download the trial version of the MKL and install it on my...
View ArticleData Fitting Functions - dfdInterpolateEx1d troubles
Hi,I am trying to use dfdInterpolateEx1d to approximate a function by a linear interpolation and i get some troubles.The function dfdInterpolateEx1d seems to be unable to do the interpolation if the...
View ArticlePARDISO gives wrong answers when given too many rhs's
When I call PARDISO with too many rhs's I get wrong answers.I first call PARDISO for phase 11 and phase 22. I keep the factored matrix in memory, and, subsequently leave the subroutine that calls...
View ArticleCompilation Interface issue with dfeast_scsrev
Hi,I recently upgrades to XE2013 to get access to the extended eigensolvers in the MKL11 libraries.I am trying to use dfeast_scsrev but am getting some compilation errors.Here is my code: SUBROUTINE...
View Articleadvice for getting rid of temporary creation and fortran95 interface for dsecnd
I am using intel fortran compiler and intel mkl for a performance check. I am passing some array sections to Fortran 77 interface with calls like<code>call dgemm( transa,transb,sz_s,P,P,& a,...
View Articlemixing fortran77 and fortran95 interfaces
Dear all,I programmed my code with the fortran 77 interfaces. However I was wondering if that is possible to mix fortran77 and fortran95 routines. For instance, looking at the sparse blas level 2...
View Articlefortran95 interface syevr return -1001 as info
Dear all,As the subject line suggests, fortran95 interface of dsyevr, namely, syevr returns -1001 as info parameter. I could not understand what this means, could you please help me with this...
View Article2D FFT
Hello,I have been assigned the task of converting a matlab script to C++, and am currently working on the FFT part.Given this input:octave:67> r57 r57 = 0.00000 0.20000 0.30000 0.30000...
View Articleinconsistent results from mkl_dcscmv and mkl_dcsrmv
Hello,I tried to calculate the (sparse) matrix-vector product using mkl_dcscmv and mkl_dcsrmv. However, sometimes they gave me different result. Here is the example I was using:Let M : = [ 1 1 0 0...
View Articlehow to tell mkl the underlying microarchitecture
Hi: I read the gcc document as well as intel compiler document, they both said that the default behavior will not detect the underlying microarchitecture in linux x86-64 (default is -xsse2). As a...
View ArticleRun time problems with eigenvalue decomposition
Dear all, I have the following problem with the LAPACK driver routines used for solving symmetric eigenvalue problems in C++. As the input, I have complex Hermitian matrices with the dimension ~30X30....
View ArticleWhere to find the mkl_lapack lib
Hi, I'm trying to compile numpy for python 3.4 with Intel MKL on Mac OS X. In this process I have solved many problems myself, but I'm stuck on this one:icc -m64 -O3 -g -fPIC -fp-model strict...
View ArticleHow to write ILP64-version vdRngUniform? (Namely, very long random number array)
...... MKL_LONG status; VSLStreamStatePtr stream; ...... int main(){ vslNewStream(&stream, VSL_BRNG_SFMT19937, 777); rnd_data = (double*)malloc(total * sizeof(double)); status =...
View Articleusing gcc compile this sample program
I am testing a program which diagnalize a hermitian matrix, using LAPACKE_zheev function describe in the MKL manual.When I try to compile this program, using gcc myprogram.c command. I got an error...
View ArticleErrors with using FGMRES solving large scale problems
I'm using FGMRES to solve some large scale problems.And now I have encountered with some errors when I allocate the memory for the tmp array.The MKL references say that,the FGMRES subroutine must have...
View ArticleHow to use complex data type in MKL C program?
Dear all:I am facing a problem on how to use lapack_complex_double data type, Including multiply, divide,conjugate(I can do this using z_imag=-z_imag, is this the best way?), and print the...
View ArticleODE Solver Support for C#
Does ODE Solver support calling solver API from .Net/C# managed environment?
View ArticleDistributing MKL libraries
I know this topic has come up before, but I can't see that anyone has found a satisfactory resolution to it.Is there any progress on finding out exactly which MKL libraries I need to distribute with my...
View ArticleClik here to view.
Pardiso now fails after upgrading from 11.0.5.1 to 11.1.3.1
Hi,I just upgraded my Mkl library from 11.0.5.1 to 11.1.3.1 and now Pardiso crashes. All I'm doing is running a solve on a 4x4 matrix for 3 vectors, it used to work fine but now gets a weird crash....
View ArticleQuestions with pardiso
Hello everyone,I met a tricky problem when solving a sparse unsymmetric complex matrics with fotran pardiso MKL within Windows environment:I installed Mircosoft Visual Studio 2012 and Intel Parallel...
View Article