Experimenting with the symetric general eigen value solver dfeast_scsrgv.
Running it with a test matrix, I find that it overwrites the values for e and x on the stack?
dfeast_scsrgv (&uplo, &n, a, rows, columns,
b, rows, columns,
fpm, &epsout, &loop, &emin, &emax, &m0, eigs, eigv, &m, &res, &info);
exits with info=0, m=1, suggesting one eigenvalue was found but eigs and eigv have been changed and no longer point to the original arrays.
inputs are declared as
Hi Everyone,
I'm using the data fitting routines in MKL (w_mkl_11.0.3.171) for spline interpolation. The fortran routine where these are called is used in a C++ OpenMP multi-threaded application, inside a parallel for loop for completely independent datasets.
The first attempt produced random crashes and wrong results, which suggested a possible data race. When I added a OMP CRITICAL section around the task creation (dfdNewTask1D) and destruction (dfDeleteTask) the code started working perfectly, producing the same results as the single core invocation.
Hi,
we are using Intel MKL 11.02. since a few weeks and a first C++ example calling the includes works fine. But we have some trouble to import the mkl_rt.dll in C#. When calling the dll the program crashes without any notice. We have tried it for some functions (MKL_Set_Num_Threads, VDDIV, VDADD, DASUM, MKL_Get_Max_Threads) and always get the same result. Additionally we tested 32bit and 64bit versions of mkl_rt.dll and also the Intel C#-example "vddiv" - and result every time in a uncommented program stop.
namespace MKL_Test
public unsafe class MKL_Wrapper
Hello!
Any help will be highly appriciated)
Strange thing occures when i try to build a custom dll.
the make command i use looks like this: nmake ia32 buf_lib= export=func_list name=mkl_1
while func list contains only one function
DGEMM
And everything works perfect.
Then i add one more function to func_list so it looks like this:
DGEMM
DGEMV
, make another dll:
nmake ia32 buf_lib= export=func_list name=mkl_2
Matrix-Matrix product still works
but the dgemv function seems not working
Hi all,
MKL-pardiso gives back a totally wrong solution vector for one of my smaller test cases in C with a sparse unsymmetric 96x96 matrix. The solution clearly satisfies A*x != b.
Attached is the example file pardiso_unsym_c.c from the MKL example collection in which I replaced the data for ia[], ja[] and a[] with my matrix data. There are no other changes in the program, except that I print the solution vector at the end. For me this program gives a wrong result.
Hi,
I am attempting to add an OpenMP layer of parallelism into an MPI code writtten in Fortran. However, I am not doing this correctly as it does not appear to be thread safe. I can say this with some confidence as adding an !$OMP CRITICAL region around the library call results in a correct answer. I believe that this is a result of the FFTW wrappers.
Dear Intel support,
I am linking my program with the Single Dynamic Library (mkl_rt.dll). Could you please let me know exactly which dlls from the redistributable folder should we inlcude with our software? Currently it fails to run unless we include all of the dlls in the program folder.
Thanks,
Dear All
I am trying to using pardiso (sparse direct solver) as my subdomain (substructure) solver in my domain decomposition application.
The first solve with dense right hand side (rhs) coming from source terms and boundary conditions.
The following sloves using the factorized sparse matrices with sparse right hand side only assotiated with boundary conditions.
With pardiso, the first solve I set IPARM(31) = 0 and PERM(i)=0 with phase = 13.
The following solves, I set IPARM(31) = 1, and only some PERM(i) = 1.
