@@ -907,18 +907,18 @@ in its image or coimage.
907907 ` QRpos() ` , ` QL() ` , ` QLpos() ` , ` SVD() ` , ` SDD() ` , ` Polar() ` . ` QR() ` uses the underlying
908908 ` qr ` decomposition from ` LinearAlgebra ` , while ` QRpos() ` (the default) adds a correction
909909 to that to make sure that the diagonal elements of ` R ` are positive.
910- Both result in upper triangular ` R ` , which are square when ` t ` is full rank and
911- surjective otherwise. ` QL() ` and ` QLpos() ` similarly result
912- in a lower triangular matrices in ` R ` , but only work if all block matrices are
913- tall, i.e. ` blockdim(codomain(t), c) >= blockdim(domain(t), c) ` for all `c ∈
914- blocksectors(t)`.
915-
916- One can also use ` alg = SVD() ` or ` alg = SDD() ` , with extra
917- keywords to control the absolute (` atol ` ) or relative (` rtol ` ) tolerance. We then set
918- ` Q=U ` and ` R=Σ*Vʰ ` from the corresponding singular value decomposition, where only
919- these singular values ` σ >= max(atol, norm(t)*rtol) ` (and corresponding singular vectors
920- in ` U ` ) are kept. More finegrained control on the chosen singular values can be
921- obtained with ` tsvd ` and its ` trunc ` keyword.
910+ Both result in upper triangular ` R ` , which are square when ` codomain(t) ≾ domain(t) `
911+ and wide otherwise. ` QL() ` and ` QLpos() ` similarly result in a lower triangular
912+ matrices in ` R ` , but only work in the former case, i.e. ` codomain(t) ≾ domain(t) ` ,
913+ which amounts to ` blockdim(codomain(t), c) >= blockdim(domain(t), c) ` for all
914+ ` c ∈ blocksectors(t)` .
915+
916+ One can also use ` alg = SVD() ` or ` alg = SDD() ` , with extra keywords to control the
917+ absolute (` atol ` ) or relative (` rtol ` ) tolerance. We then set ` Q=U ` and ` R=Σ*Vʰ ` from
918+ the corresponding singular value decomposition, where only these singular values
919+ ` σ >= max(atol, norm(t)*rtol) ` (and corresponding singular vectors in ` U ` ) are kept.
920+ More finegrained control on the chosen singular values can be obtained with ` tsvd ` and
921+ its ` trunc ` keyword.
922922
923923 Finally, ` Polar() ` sets ` Q=U*Vʰ ` and ` R = (Vʰ)'*Σ*Vʰ ` , such that ` R ` is positive
924924 definite; in this case ` SDD() ` is used to actually compute the singular value
@@ -934,16 +934,17 @@ in its image or coimage.
934934 ` LQpos() ` , ` RQ() ` , ` RQpos() ` , ` SVD() ` , ` SDD() ` , ` Polar() ` . ` LQ() ` uses the underlying
935935 ` qr ` decomposition from ` LinearAlgebra ` on the transposed data, and leads to lower
936936 triangular matrices in ` L ` ; ` LQpos() ` makes sure the diagonal elements are
937- positive. ` RQ() ` and ` RQpos() ` similarly result in upper triangular matrices in
938- ` L ` , but only works for wide matrices, i.e. `blockdim(codomain(t), c) <=
939- blockdim(domain(t), c)` for all ` c ∈ blocksectors(t)`.
940-
941- One can also use ` alg = SVD() ` or ` alg = SDD() ` , with extra
942- keywords to control the absolute (` atol ` ) or relative (` rtol ` ) tolerance. We then set
943- ` L=U*Σ ` and ` Q=Vʰ ` from the corresponding singular value decomposition, where only these
944- singular values ` σ >= max(atol, norm(t)*rtol) ` (and corresponding singular vectors in
945- ` Vʰ ` ) are kept. More finegrained control on the chosen singular values can be obtained
946- with ` tsvd ` and its ` trunc ` keyword.
937+ positive. The matrices ` L ` are square when ` codomain(t) ≿ domain(t) ` and tall otherwise.
938+ Similarly, ` RQ() ` and ` RQpos() ` result in upper triangular matrices in ` L ` , but only
939+ works if ` codomain(t) ≿ domain(t) ` , i.e. when
940+ ` blockdim(codomain(t), c) <= blockdim(domain(t), c) ` for all ` c ∈ blocksectors(t) ` .
941+
942+ One can also use ` alg = SVD() ` or ` alg = SDD() ` , with extra keywords to control the
943+ absolute (` atol ` ) or relative (` rtol ` ) tolerance. We then set ` L=U*Σ ` and ` Q=Vʰ ` from
944+ the corresponding singular value decomposition, where only these singular values
945+ ` σ >= max(atol, norm(t)*rtol) ` (and corresponding singular vectors in ` Vʰ ` ) are kept.
946+ More finegrained control on the chosen singular values can be obtained with ` tsvd ` and
947+ its ` trunc ` keyword.
947948
948949 Finally, ` Polar() ` sets ` L = U*Σ*U' ` and ` Q=U*Vʰ ` , such that ` L ` is positive definite;
949950 in this case ` SDD() ` is used to actually compute the singular value decomposition and no
0 commit comments