Cholesky←{          ⍝ decomposition of a Hermitian positive-definite matrix.
    1≥n←≢⍵:⍵*0.5
    p←⌈n÷2
    q←⌊n÷2
    X(p,p)↑⍵⊣Y(p,-q)↑⍵⊣Z(-q,q)↑⍵
    L0←∇ X
    L1←∇ Z-(T(+⍉Y)+.×⌹X)+.×Y
    ((p,n)L0)(T+.×L0),L1
}
code_colours

test script

Back to: notes

Back to: Workspaces