(Arne Babenhauserheide)
2016-12-19: note the reason for using the mostly-exact-sqrt note the reason for using the mostly-exact-sqrt
diff --git a/examples/cholesky.w b/examples/cholesky.w --- a/examples/cholesky.w +++ b/examples/cholesky.w @@ -72,7 +72,7 @@ a is represented as list of lists." matrix-ref a j j : > sum 0 ; diagonal element matrix-set! a i i - mostly-exact-sqrt sum + mostly-exact-sqrt sum ; preserves the exactness, since the result is an exact number, though not an exact result else throw 'matrix-numerically-not-symmetric-positive-definite a do-ec (: i n)