From 7cf0686e9c6f81ddf6b9b6c6d01ddd2612b56668 Mon Sep 17 00:00:00 2001 From: fred Date: Sun, 24 Mar 2002 19:50:04 +0000 Subject: [PATCH] lilypond-0.1.1 --- flower/test/mat-test.cc | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/flower/test/mat-test.cc b/flower/test/mat-test.cc index f7ce01bfc5..8fdc7d794e 100644 --- a/flower/test/mat-test.cc +++ b/flower/test/mat-test.cc @@ -47,7 +47,7 @@ matrix() Matrix hilbert(N,N), h2(hilbert); for (int i=0; i < N; i++) { for (int j=0; j < N; j++) { - hilbert(i,j) = 1/(i+j+1); + hilbert(i,j) = 1/Real(i+j+1); h2 (i,j) = (abs(i-j) > 3) ?0 : hilbert(i,j); } } @@ -55,6 +55,10 @@ matrix() Choleski_decomposition ch(h2); cout << "red Hilbert " << h2; cout << "choleski " << ch.L; + Matrix T =ch.L.transposed(); + cout << "L^T " << T; + cout << "L * L^T" << ch.L * T; + cout << "H2^{-1} * H2" << h2 * ch.inverse(); } ADD_TEST(matrix); -- 2.39.5