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

Question

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

-Thanks

Answer 1

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

1. 删除C的缺失三角形
2. 执行pdgeadd（＆＃39; T＆＃39;，C，C）
3. 将C的给定三角形的所有元素除以

它有效，但我不知道是否有任何限制，因为它理论上禁止用A和C来调用pdgeadd来引用相同的内存......