topic dgetrf with mkl_progress in Intel® oneAPI Math Kernel Library

dgetrf with mkl_progress

BrianZHANG — Thu, 04 Jul 2024 05:00:10 GMT

I am currently calling the dgetrf function in Intel MKL in Fortran to solve a dense system of linear equations. The Developer Reference for Fortran metions that dgetrf supports the progress routine feature.

I customized the mkl_progress function and called mkl_set_progress to explicitly set it. However, dgetrf got the correct result but did not call mkl_progress to print any information. The pardiso function called elsewhere in my code prints progress information, so mkl_progress should be correct. But why is dgetrf not printing progress information?

My oneAPI component is the latest version 2024.2. Here is the code:

Integer Function mkl_progress( Thread, Step, Stage ) Integer*4 Thread, Step Character*(*) Stage Print*,'Thread:',Thread,',Stage:',Stage,',Step:',Step Mkl_Progress = 0 Return End Call mkl_set_progress(mkl_progress) Call dgetrf(......)

Re: dgetrf with mkl_progress

Gennady_F_Intel — Mon, 08 Jul 2024 06:52:37 GMT

The trivial example ( see _getrf.f90 is attached ) shows some progress. the test fails but the idea was to show how to see progress...

the part of the of progress log is attached also ( see _getrf.log ).

Thread: 0 ,Stage:DGETRF,Step: 8
Thread: 6 ,Stage:DGETRF,Step: 8
Thread: 29 ,Stage:DGETRF,Step: 8

.....................................

Thread: 9 ,Stage:DGETRF,Step: 32
Thread: 16 ,Stage:DGETRF,Step: 32
Thread: 13 ,Stage:DGETRF,Step: 32
Thread: 0 ,Stage:DGETRF,Step: 40
Thread: 0 ,Stage:DGETRF,Step: 48

--Gennady

Re:dgetrf with mkl_progress

Gennady_F_Intel — Tue, 09 Jul 2024 05:42:53 GMT

The topic is closing and will not longer be monitored.