如何查看Intel MKL中dsyev（或通用LAPACK）函数的运行进度？

针对20000×20000实对称矩阵的dsyev对角化任务，以下几种方法可以实现进度或剩余时间的监控：

1. 利用MKL_VERBOSE环境变量获取阶段耗时 #

MKL提供了verbose模式，会输出每个底层计算步骤的耗时信息，能帮你大致判断整体进度。

• 使用方式：在运行程序前设置环境变量：

  export MKL_VERBOSE=1
  

  程序运行时，终端会打印MKL各函数的调用细节，包括每个子步骤的执行时间。你可以通过这些阶段耗时，结合已完成步骤的总时间，对比同规模矩阵的历史运行数据，估算剩余时间。

2. 自定义计时+多线程进度监控 #

由于dsyev是阻塞式调用，单线程下无法在函数执行中实时获取进度，可通过多线程实现：

• 步骤：
  1. 先使用小规模矩阵（比如2000×2000）测试dsyev的性能，计算出单位规模的平均耗时，预估20000×20000矩阵的总运行时间。
  2. 在主线程调用dsyev前，启动一个独立的监控线程：
     • 监控线程每隔固定时间（比如10秒）读取主线程的已运行时间，用已用时间/预估总时间计算当前进度，同时输出剩余时间（预估总时间 - 已用时间）。
  3. 注意：监控线程需要避免干扰主线程的计算资源，可设置较低的优先级。

3. MKL进度回调函数（需确认函数支持性） #

MKL提供了mkl_set_progress_function接口，可注册回调函数来获取部分计算任务的进度更新。

• 实现步骤：
  1. 定义进度回调函数，示例：

     int progress_callback(double progress, void *user_data) {
         printf("当前进度：%.2f%%\n", progress * 100);
         // 返回0继续执行，返回1终止计算
         return 0;
     }
     

  2. 在调用dsyev前注册回调：

     mkl_set_progress_function(progress_callback, NULL);
     

  注意：并非所有MKL LAPACK函数都支持进度回调，需确认你的MKL版本中dsyev是否触发该回调。如果dsyev不支持，可考虑改用MKL的分块对角化函数（如dsyevr），这类函数通常更易实现进度监控。

内容的提问来源于stack exchange，提问作者velenos14