File tree Expand file tree Collapse file tree 1 file changed +10
-7
lines changed Expand file tree Collapse file tree 1 file changed +10
-7
lines changed Original file line number Diff line number Diff line change @@ -53,6 +53,7 @@ impl_lu!(c32, lapack_sys::cgetrf_);
5353impl_lu ! ( f64 , lapack_sys:: dgetrf_) ;
5454impl_lu ! ( f32 , lapack_sys:: sgetrf_) ;
5555
56+ #[ cfg_attr( doc, katexit:: katexit) ]
5657/// Helper trait to abstract `*getrs` LAPACK routines for implementing [Lapack::solve]
5758///
5859/// If the array has C layout, then it needs to be handled
@@ -63,13 +64,15 @@ impl_lu!(f32, lapack_sys::sgetrf_);
6364/// or "no transpose", respectively. For the "Hermite" case, we
6465/// can take advantage of the following:
6566///
66- /// ```text
67- /// A^H x = b
68- /// ⟺ conj(A^T) x = b
69- /// ⟺ conj(conj(A^T) x) = conj(b)
70- /// ⟺ conj(conj(A^T)) conj(x) = conj(b)
71- /// ⟺ A^T conj(x) = conj(b)
72- /// ```
67+ /// $$
68+ /// \begin{align*}
69+ /// A^H x &= b \\\\
70+ /// \Leftrightarrow \overline{A^T} x &= b \\\\
71+ /// \Leftrightarrow \overline{\overline{A^T} x} &= \overline{b} \\\\
72+ /// \Leftrightarrow \overline{\overline{A^T}} \overline{x} &= \overline{b} \\\\
73+ /// \Leftrightarrow A^T \overline{x} &= \overline{b}
74+ /// \end{align*}
75+ /// $$
7376///
7477/// So, we can handle this case by switching to "no transpose"
7578/// (which is equivalent to transposing the array since it will
You can’t perform that action at this time.
0 commit comments