From: fred Date: Sun, 24 Mar 2002 19:50:04 +0000 (+0000) Subject: lilypond-0.1.1 X-Git-Tag: release/1.5.59~4336 X-Git-Url: https://git.donarmstrong.com/?a=commitdiff_plain;h=7cf0686e9c6f81ddf6b9b6c6d01ddd2612b56668;p=lilypond.git lilypond-0.1.1 --- 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);