pardiso Memory leak?
I feel very puzzled with the Win32 Pardiso version.My integrated enviroment is Intel(R) Visual Fortran Compiler XE 12.1.7.371.My matrix has 5 million non zeros ,and n =200,000 , real and symmetric...
View Articlepardiso intermediate result output
Hi,can pardiso save intermediate results in hard disk?For example ,after phase=22,I save the intermediate results in the hard disk.Few days later ,I call pardiso to finish the calculation of phase=33....
View ArticleMemory Leak
Hi,I have an issue with dgemm() and other mkl functions not returning the memory that is used internal to the function.How can I release the memory used up. If this example is ran long enough the...
View ArticleClik here to view.
non-linear optimization: jacobi_solve question
I have tried to use the same _JACOBIMATRIX_HANDLE_t to make evaluations of the jacobian at successive iteration points of the optimizer,but the output matrix (fjac) is not getting updated.I created an...
View ArticleMissing FFTW libraries
I'm trying to compile a code written in fortran90/95 and I'm getting this error:user:~> makeifort -O3 -lmpi -L/opt/local/intel/fftw/lib -I/opt/local/intel/fftw/include -...
View ArticleClik here to view.
Data race in Pardiso Solver?
Hi All,I ran into a tricky problem when test data race for my codes. The Intel Inspector XE 2013 shows that there is data race in the following calling:282 !C.. Factorization.283 phase = 22...
View ArticleError in feast eigenvalue solver for sparse matrices
Hi,We are trying to use the MKL FEAST 11.0.3.1 solver for symmetric sparse CSR matrices of doubles. The call we are using is something like:--------------------------------------------#include...
View ArticleMinimum working example for mkl_ddnscsr
Hi. Can anyone provide me with a minimum working example for mkl_ddnscsr? I have tried this so far#include <stdio.h> #include <stdlib.h> #include <mkl.h> int main(int argc, char...
View ArticleIssues with declarations of MKL functions: 'remark #424: extra ";" ignored'...
When Intel C++ compiler option /W5 is turned on the compiler shows lots of 'remark #424: extra ";" ignored' messages related declaration of some functions in MKL headers. Here is a small...
View ArticleClik here to view.
How to solve this symmetric indefinite matrix with pardiso?
Hello, Greetings! I would like to know how to effectively solve this(attached file) symmetric indefinite system using PARDISO. This particular matrix has super diagonal dominancy and almost zero...
View ArticleMKL vs Microsoft exp() function
I have a client that is migrating a large C++ software base from 32- to 64-bit code in MS Visual Studio. One of the problems they are having is that the 32- and 64-bit versions of the C library exp()...
View Articleefficiently solving least squares problems iteratively
I'm performing an iterative routine, where at each iteration I am solving a least squares problem using ?gels. At each iteration step I update a column of the matrix A, which will ulitmately converge,...
View ArticleClik here to view.
Program crash in mkl_avx.dll on Windows with Intel MKL 11.0
One of our customers has encountered a problem with our library (NumPy) when linked against MKL version 11.0.3 on the Windows platform. The program dies during an eigenvalue decomposition. The...
View ArticleClik here to view.
Pardiso result problem
Hi All,I have implemented pardiso solver in our model for flow simulation. Pardiso can solve most of our problem well and produce correct results. But recently I have a case that pardiso generates...
View ArticleClik here to view.
Program crash in mkl_avx.dll on Windows with Intel MKL 11.0
One of our customers has encountered a problem with our library (NumPy) when linked against MKL version 11.0.3 on the Windows platform. The program dies during an eigenvalue decomposition. The...
View ArticleClik here to view.
Pardiso memory leak
Hi,I have the following subroutine works perfectly. However, I am using it in a nonlinear solution. I call it many times in a program. After each call, the memory usage increases, though I am using...
View Articlepardiso iparm(2) parameter
Dear AllFrom the talbe http://software.intel.com/en-us/articles/pardiso-parameter-table#table2It states thatiparm(2)=0 is the MD algorithm and iparm(2)=2 is metis package.But from the former thread...
View ArticlePerformance gets worse over time for the same instructions
First, I'm not if this is the right forum for this question. I don't know what is as the reason could be due to hardware or MKL or .NET or some other hidden factors.I have a neural network code in C#...
View ArticleClik here to view.
sparse right hand side reordering problem
Dear AllI am trying to use the sparse right hand feature of pardiso ( iparm(31)=1 ).I tested with an identity matrix A.When I set perm[i] = 1 for last several entries except the last one, I got the...
View Articleifort: error #10104: unable to open '--start-group' (MPI with FGMRES)
Hello!Can someone help me get past this error please...I get the following error while generating my executable:ifort: error #10104: unable to open '--start-group'I'm trying to run FGMRES sequentially...
View Article