0
我有一個矩陣$ A $和矩陣$ B_i = L_i^\ mathsf {T} M_i L_i $。我想估計下面的量爲$ i的\ [1，n]的$，如何計算使用lapack/blas的線性系統，當矩陣是psd時？
\begin{eqnarray}
m_i^\mathsf{T}(A-B_i)^{-1}m_i
\end{eqnarray}
我使用LAPACK/BLAS獲得$ B_i $ S，但是當我想計算線性系統$（A- B_i）x = m_i $，使用lapack，例程dpotrf和dpotrs，他們要求矩陣處於下三角形式或上三角形式。如何使用lapack/blas來計算上述計算？ C或C++中是否有任何示例？從Full矩陣複製到更低三角形的代價是昂貴的，是否有任何方法可以直接在下三角形中執行$ B_i $的乘法，所以我們不需要複製它。