G = size 3 by 3 Matrix of double:
[ 1 1 1 ]
[ 1 2 3 ]
[ 1 3 6 ]

A = syrk(G) = size 3 by 3 LowerMatrix of double:
[ 3 * * ]
[ 6 14 * ]
[ 10 25 46 ]

X = size 3 by 2 Matrix of double:
[ 4 6 ]
[ 3 2 ]
[ 5 1 ]

B = symm(A,X) = size 3 by 2 Matrix of double:
[ 80 40 ]
[ 191 89 ]
[ 345 156 ]

L = potrf(A) = size 3 by 3 LowerMatrix of double:
[ 1.73205 * * ]
[ 3.4641 1.41421 * ]
[ 5.7735 3.53553 0.408248 ]

A\B = L\L'\B = size 3 by 2 Matrix of double:
[ 4 6 ]
[ 3 2 ]
[ 5 1 ]

frobenius(A\B-X) = 1.42081e-12

