0

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

-谢谢

4

1 回答 1

0

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

  1. 删除 C 的缺失三角形
  2. 执行 pdgeadd('T',C,C)
  3. 将 C 的给定三角形的所有元素除以二

它可以工作，但我不知道是否有任何限制，因为理论上禁止使用引用同一内存的 A 和 C 调用 pdgeadd ......

于 2015-05-17T20:45:00.277 回答