From: ESeNonFossiIo Date: Wed, 27 May 2015 09:58:52 +0000 (+0200) Subject: add template X-Git-Tag: v8.3.0-rc1~152^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=aff4d0ae9380d8f5d1d0abc33c2c1436b7a4428d;p=dealii.git add template --- diff --git a/include/deal.II/lac/linear_operator.h b/include/deal.II/lac/linear_operator.h index 73b0598c63..dc54741dea 100644 --- a/include/deal.II/lac/linear_operator.h +++ b/include/deal.II/lac/linear_operator.h @@ -2162,7 +2162,7 @@ lower_triangular_operator(const BlockMatrix &block_matrix) // And reinitialize every individual block with reinit_range_vectors: for (unsigned int i = 0; i < block_matrix.n_block_rows(); ++i) - linear_operator(block_matrix.block(i,0)).reinit_range_vector(v.block(i), fast); + linear_operator(block_matrix.block(i,0)).reinit_range_vector(v.block(i), fast); v.collect_sizes(); }; @@ -2175,7 +2175,7 @@ lower_triangular_operator(const BlockMatrix &block_matrix) // And reinitialize every individual block with reinit_domain_vectors: for (unsigned int i = 0; i < block_matrix.n_block_cols(); ++i) - linear_operator(block_matrix.block(0,i)).reinit_domain_vector(v.block(i), fast); + linear_operator(block_matrix.block(0,i)).reinit_domain_vector(v.block(i), fast); v.collect_sizes(); }; @@ -2277,7 +2277,7 @@ upper_triangular_operator(const BlockMatrix &block_matrix) // And reinitialize every individual block with reinit_range_vectors: for (unsigned int i = 0; i < block_matrix.n_block_rows(); ++i) - linear_operator(block_matrix.block(i,0)).reinit_range_vector(v.block(i), fast); + linear_operator(block_matrix.block(i,0)).reinit_range_vector(v.block(i), fast); v.collect_sizes(); }; @@ -2290,7 +2290,7 @@ upper_triangular_operator(const BlockMatrix &block_matrix) // And reinitialize every individual block with reinit_domain_vectors: for (unsigned int i = 0; i < block_matrix.n_block_cols(); ++i) - linear_operator(block_matrix.block(0,i)).reinit_domain_vector(v.block(i), fast); + linear_operator(block_matrix.block(0,i)).reinit_domain_vector(v.block(i), fast); v.collect_sizes(); }; @@ -2393,11 +2393,11 @@ block_triangular_inverse( const BlockMatrix &block_matrix, if (lower) { - op_a = lower_triangular_operator(block_matrix); + op_a = lower_triangular_operator(block_matrix); } else { - op_a = upper_triangular_operator(block_matrix); + op_a = upper_triangular_operator(block_matrix); } auto id = identity_operator(op_a.reinit_range_vector);