topic Re: Question about sparse matrix in Intel® oneAPI Math Kernel Library

Question about sparse matrix

kim__seongik — Fri, 10 Jul 2020 06:36:40 GMT

1. Is there an appropriate subroutine in mkl with respect to matrix assembly?
I made and used by myself, but it takes too long.
I want to something like below expression in fortran
(example-matlab case) global_matrix(dof,dof) = global_matrix(dof,dof) + element_matrix

2. Is there any subroutine in mkl to create a sparse matrix format?

3. Which subroutine should be used for sparse matrix multiplication?
I want to something like below expression in fortran

(example -matlab case) A*x=b (A : sparse matrix, x: vector, b:vector)

4. Is there a subroutine that can do in fortran mkl how to delete certain rows and columns?

I want to something like below expression in fortran

(example -matlab case) matrix(dof,dof) = [];

5. I want to use scalapack of intel mkl, but can't I use it except cluster??
When I used scalapack's dlasrt2 on the desktop, error lnk 2019 occurred.

Re: Question about sparse matrix

Gennady_F_Intel — Fri, 10 Jul 2020 09:40:57 GMT

1. there are no such routines

2. you can try to use conversion from/to dense - CSR, CSC, and others supported sparse formats. Please refer to use MKL Reference to see more details

3. please refer to the documentation follow with this link https://software.intel.com/content/www/us/en/develop/documentation/mkl-developer-reference-c/top/blas-and-sparse-blas-routines/sparse-blas-level-2-and-level-3-routines/sparse-blas-level-2-and-level-3-routines-1.html

You may see a lot of m*v routines with different sparse formal, e.x - tmkl_?csrgemv

4. There is no such routine.

5. you can build and run on the SMP system as well.

how did you link the case? how did you run it?

Re: Question about sparse matrix

Gennady_F_Intel — Fri, 10 Jul 2020 09:42:16 GMT

regarding the #5 - please check the mkl linker adviser how to properly link with the ScaLAPACK components

Re: Question about sparse matrix

kim__seongik — Mon, 13 Jul 2020 02:35:16 GMT

Thank you for your reply !!

In number 5 my setting is below figure::

Re: Question about sparse matrix

Gennady_F_Intel — Mon, 13 Jul 2020 09:12:46 GMT

Please check the user's guide on how to use the linker

https://software.intel.com/content/www/us/en/develop/documentation/mkl-linux-developer-guide/top/linking-your-application-with-the-intel-math-kernel-library/linking-quick-start/using-the-mkl-compiler-option.html

-mkl=cluster

to link with Intel® MKL cluster components (sequential) that use Intel MPI.

Re: Question about sparse matrix

Gennady_F_Intel — Mon, 13 Jul 2020 09:14:43 GMT

I missed that you use windows:

the linking would look like as follows: /Qmkl=cluster