在当前时间点是否存在任何支持自动微分和对称矩阵的线性代数包？
在一系列计算过程中，某些矩阵保持对称通常很重要。A <- (A + transpose(A))/2一种可能的解决方案是在每次计算后将投影应用于对称部分。出于显而易见的原因，这充其量是次优的，在最坏的情况下速度非常慢。
我正在寻找的是一个线性代数包，它保证涉及对称矩阵的某些操作会产生对称矩阵。例如：
A + B, AB+BA, A^n, A', c*A, VAV', VV', V'V A^-1, exp(A), A∘B, A ⊗ B,
都应该再次对称，给定A，B对称，V任意，c标量。