From: kanschat Date: Sun, 1 Oct 2006 23:02:08 +0000 (+0000) Subject: test for LU decomposition X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=8b5968befe370fb1f20c4527d5cfad38f711e4fc;p=dealii-svn.git test for LU decomposition git-svn-id: https://svn.dealii.org/trunk@13974 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/lapack/full_matrix_02.cc b/tests/lapack/full_matrix_02.cc new file mode 100644 index 0000000000..b4e478a039 --- /dev/null +++ b/tests/lapack/full_matrix_02.cc @@ -0,0 +1,88 @@ +//-------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2005, 2006 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//-------------------------------------------------------------------- + +// LAPACKFullMatrix::compute_lu_factorization +// LAPACKFullMatrix::apply_lu_factorization + +#include "../tests.h" +#include +#include +#include + +#include +#include + + +// Fill a matrix with the values of the Hilbert matrix +template +void +hilbert (LAPACKFullMatrix& M, + const bool nonsymmetric) +{ + const unsigned int n = M.n_rows(); + for (unsigned int i=0;i M(size, size); + hilbert(M, nonsymmetric); + + Vector u(size); + Vector v(size); + Vector x(size); + Vector y(size); + + for (unsigned int i=0;i