分布式下/上三角矩阵的上/下部分的原位恢复

时间:2015-05-02 21:58:31

标签: matrix copy transpose triangular scalapack

我想知道是否存在任何Scalapack例程可以将三角矩阵转换回一般矩阵,即将下/上部分复制到上/下部分 - 我还没找到它。 我知道有一个简单的工作方式使用p?geadd操纵对角线,但是那个消耗的内存是理论上需要的两倍,我买不起...... 所以我觉得很奇怪,Scalapack / PBLAS没有提供解决这个相对简单问题的常规方法。如果有人比我更了解,请告诉我!

-Thanks

1 个答案:

答案 0 :(得分:0)

我找到了一个危险的解决方案:

  1. 删除C的缺失三角形
  2. 执行pdgeadd(' T',C,C)
  3. 将C的给定三角形的所有元素除以
  4. 它有效,但我不知道是否有任何限制,因为它理论上禁止用A和C来调用pdgeadd来引用相同的内存......