虽然这个问题是关于单个库的特定功能的，但我前段时间偶然发现了同样的问题。

据我所知，PETSc (3.8.3) 中没有现成的矩阵函数可以开箱即用。考虑到以下功能的存在，我很惊讶几年前和现在没有找到解决方案：

• VecAbs这将完全符合您的要求，但仅适用于向量
• MatImaginaryPart, MatRealPart, 和MatConjugate哪个会执行$A_{ij}=\Re(A_{ij})$,$A_{ij}=\Im(A_{ij})$， 和$A_{ij}=A^*_{ij}$， 分别。

考虑到支持矩阵上的元素 Re、Im 和共轭运算，我假设MatAb的所有必要构建块都在那里，但没有正式编码。

从MatRealPart子程序源代码中，以下摘录正在完成这项工作

(*mat->ops->realpart)(mat);

矩阵上支持和传播的操作列表存储在typedef的include/petscmat.h中。MatOperation从技术上讲，您可以编写自己的操作并将其传播到您需要的数据类型（也许，甚至通过使用自定义用户矩阵MatOperation和MatShellSetOperation）。然而，问题是这样的发展是否值得付出所需的时间和痛苦。