Crash in the Cluster Sparse Solver Windows-64
Hi, I am using the Cluster Sparse Solver on Windows 64 (Parallel Studio xe 2016.0.110) on a symmetric indefinite matrix and the numeric factorization crashes for 2 MPI process and 4 OpenMP threads....
View ArticleIs there a more efficient way for element-by-element multiplication for matrix?
Hi, everyone. I was writing some FORTRAN codes for computing the DOT of two matrix which idea is as same as dot of vectors. Firstly, it calculates the element-wise product of the two matrix. Secondly,...
View ArticlePardiso MKL_PROGRESS
I just upgraded MKL library from 10.1.1.022 to 11.3 (positive symmetric matrix)For 64-bit, I first tried to use static linking with sequential threading using the following lib...
View Articlemkl_cspblas_dcsrgemv gives wrong results
I was trying to compute the sparse matrix-vector multiplication using mkl_cspblas_dcsrgemv() in the MKL.The c compiler comes from composer_xe_2013.0.079 in the parallelstudio_2013 version with all the...
View Articlemkl_progress limited functionality
Hi, MKL team.We use MKL in multithreaded application. In each our thread we have heavy calculations: otherwise we wouldn't use such a giant like MKL. That is why we need to know where is MKL routine at...
View ArticleError with MKL and GAMESS
Hello guys!I've been trying to install GAMESS software and I'm mired on a particular step.I'm using the official guide provided with the software, and when I get to this step:./lked gamess 01 >&...
View Articlecopy for integer array
Hello,I have a naive question about that why there is only float array copy .Is the time cost of integer array copy much smaller than float number copy? Best Regards,Jiajun
View ArticleProposed GEMM Extension: Triangle Triangle Matrix Multiply
Hi,After speaking to one of the MKL developers recently I was wondering whether it might be beneficial to add some new functionality to MKL in order to compute triangular-triangular matrix products. As...
View ArticleSIGSEGV attempting to use mkl_dcoomm()
Hopefully this isn't a completely stupid question. I've been using BLAS dgemm() (with various different libraries) successfully, and now I want to optimize some of my sparse matrices and have been...
View ArticlePardiso Schur Complement Memory Allocation
Hello,I was just looking at the Schur Complement functionality of Pardiso and was studying the documentation for quite some time but there are still two things that I don't understand:- How much memory...
View ArticleUnexpected Results of 3D Not Inplace FFT
The following is my C code snipped illustrating a problem.I set dim1/dim2/dim3 to all be the next higher power of 2 based on the lengths of the actual input data. My test that works has all dims set...
View ArticleInsight as to why Valgrind shows memory leak for LAPACKE
This is the first time I've ever used Intel's MKL LAPACKE and Valgrind. Unfortunately, I get an error with something I have little to no experience with. I could use some advice on how to cure a...
View ArticleVS + MKL: dgemm_example.c - ERROR: mkl_intel_thread.dll is missing
Hi,I have a following problem. I have installed Intel® Parallel Studio XE Cluster Edition for Windows (with Intel MKL 11.3 Update 2), and during the installation I chose the option of integrating it...
View ArticleHow do I use mkl_sparse_d_spmmd()?
Continuing my journey into sparse storage, I need to do sparse * sparse. It seems like mkl_sparse_d_spmmd() fits the bill nicely, but the documentation at https://software.intel.com/en-us/node/590105...
View Articlezfeast_hrci initialization shows strange behaviour
Hi all,i am working my way through the example zexample_sparse_rci_c.c. I have it working now for my matrix, BUT found a strange behaviour for the input arrays Aq and Sq:when declaring them like...
View ArticleHow to do a 2d Fourier tranformation stored in column-major?
I have read the MKL manual but I can hardly understand the explanation below:MKL_LONG dims[] = { nd, …, n2, n1 }; DftiCreateDescriptor( &hand, precision, domain, d, dims ); // The above call...
View Article0x00007fffe006cbba in __vsldBRngWH () from...
STACK TRACE> x <- rxExec(rnorm, 100, 5, RNGkind="WH")Program received signal SIGILL, Illegal instruction. 0x00007fffe006cbba in __vsldBRngWH () from...
View ArticleWhat does mkl_sparse_optimize do in sparse BLAS?
I notice there is a routine called mkl_sparse_optimizeAlthough, the name suggest it is doing optimization. But What kind of optimization? The documentation didn't say about this.
View ArticleIntel® Parallel Studio XE 2017 Beta invitation – please register and provide...
Intel® Parallel Studio XE 2017 Beta is introducing a lot of new features and a couple of new products.Please register for the Beta program, checkout the new features introduced in Intel® Math Kernel...
View ArticleIntel Optimized Linpack Benchmark
Hello,I'm trying to use Intel's Linpack Benchmark to stress my cpu. What I need is being able to control the load it puts on the processor. Any of you have an idea how could I do that?Thanks in advance.
View Article