?trnlsq related routines for nonlinear fitting
I want to solve a nonlinear least square fitting problem to estimate a and b. g(x)=a(x^2)+bx. I have two vectors working as the training data. fx1={1.0, 2.0, 3.0, 4.0}; fx2={2.0, 4.0, 6.0, 8.0}. I...
View Articlestrnlsq related routines for nonlinear fitting
I want to solve a nonlinear least square fitting problem to estimate a and b. g(x)=a(x^2)+bx. I have two vectors working as the training data. fx1={1.0, 2.0, 3.0, 4.0}; fx2={2.0, 4.0, 6.0, 8.0}. I...
View ArticleForum poll: Intel MKL and threading
Intel MKL users,We would like to hear from you how you are using Intel MKL with threading. Do you use the parallel or sequential MKL? How do your multithreaded applications use MKL? We would appreciate...
View ArticleIs it possible to update Schur complement quickly with MKL/Pardiso?
I have a matrixA = [A11 A12;A21 A22];and I'd like to compute the Schur complement of A11 in A using pardiso: S = A22 - A21 * A11^-1 * A12. This seems well supported by the PARDISO-SCHUR routine...
View ArticleIntel MKL 11.1 java example build issue
Hello,I am trying to build Intel MKL 11.1 java example as discripted here (https://software.intel.com/en-us/articles/performance-tools-for-software...). I got the following error while running nmake....
View ArticleIs there any way to use VML on 2-dimensional array?
Suppose I have a 2-dimensional array a[m][n],and I want to get the result array of one row of a[m][n] adding another row of a[m][n].Is there any way to use VML in this case？I‘ve tried to directly use...
View Articleerror for R calling MKL with multithreading
Specifically, when using mkl functions (such as lapacke_dgesvd) in C++ with multi-threading, it works fine. But with the same function and R binding to call it under multi-threading, the result of SVD...
View ArticleMKL_malloc interface missing from Fortran module
Hi,It looks like the interface for MKL_malloc is missing from the mkl_service Fortran 90 module, at least in the versions I have. Is it supposed to be in there? The documentation seems to indicate so,...
View ArticleMKL ILU for BSR
I am the technical manager in a reservoir simulation project and working with various self made and existing math libraries for my speed comparisons. I use ILU (Incomplete LU) factorization on a...
View ArticlePardiso scaling inversely with number of threads for openmp
Hello. After a number of failed debugging attempts and tests, I'm hoping to get some input on using pardiso in parallel with openmp. The software in question uses the intel fortran compiler (ifort)...
View ArticlePardiso scaling inversely with number of threads for openmp
Hello. After a number of failed debugging attempts and tests, I'm hoping to get some input on using pardiso in parallel with openmp. The software in question uses the intel fortran compiler (ifort)...
View ArticlePardiso scaling inversely (openmp)
Hello. After a number of failed debugging attempts and tests, I'm hoping to get some input on using pardiso in parallel with openmp. The software in question uses the intel fortran compiler (ifort)...
View ArticleProblem using pdhseqr
Hi,I have some trouble using pdhseqr from ScaLAPACK (it is supposed to compute the eigen value and Schur factorization of a matrix already in Hessenberg form). It happens for both pshseqr and pdhseqr....
View ArticleProblems in running MKL fortran example
Dear all,After the example fortran code vdrnggaussian.f successfully compiled (in visual studio 2010+Intel Parallel Studio XE 2013), the resulting executable file returns the following error...
View Articlepardiso multiple right hand side with sparse matrix
Hi All,Pardiso can handle multiple right side, AX=B, meaning B is a matrix. Every column of B is a right hand side. If B is very sparse, is there any way to exploit this feature, instead of storing B...
View ArticleLinking BLA95 and LAPACK 95 to a shared library in Linux
I'm in the process of porting an application from Windows to Linux, and one component of my application is a shared library (a DLL on Windows) compiled with Intel Fortran and with heavy use of the BLAS...
View ArticleClik here to view.
Using Intel® C++ Composer XE for Multiple Simple Random Sampling without...
<!--break-->IntroductionRandom sampling is often used when pre- or post-processing of all records of the entire data set is expensive, as in the following examples. When the file of records or...
View ArticleLarge Scale Weighted Least Squares
Hello,I would like to use Intel MKL to solve a Large Scale Weighted Least Squares problem.The Matrices are sparse yet their elements number might be huge. This is a result of working on images and...
View ArticleThe PZGESV subroutine
Hello everybody, Please I need your help to understand the scalapck subroutine PZGESV. I am trying to use this subroutine to solve a linear equation system A*X = B where A is a (N by N) matrix...
View ArticleOctave 3.8.1 with MKL
I am trying to rebuild GNU Octave with the Intel MKL on Ubuntu 14.04uname -aLinux zoli-Precision-WorkStation-T3500 3.13.0-27-generic #50-Ubuntu SMP Thu May 15 18:06:16 UTC 2014 x86_64 x86_64 x86_64...
View Article