dgetrf with mkl_progress

[BrianZHANG]

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(......)

[Gennady_F_Intel]

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

[Gennady_F_Intel]

The topic is closing and will not longer be monitored.