22#ifndef INCLUDE_SOLVERS_AND_PRECONDITIONERS_BLOCK_JACOBI_MATRICES_H_
23#define INCLUDE_SOLVERS_AND_PRECONDITIONERS_BLOCK_JACOBI_MATRICES_H_
30template<
typename Number>
32initialize_block_jacobi_matrices_with_zero(std::vector<dealii::LAPACKFullMatrix<Number>> & matrices)
35 for(
auto & m : matrices)
44template<
typename Number>
46calculate_lu_factorization_block_jacobi(std::vector<dealii::LAPACKFullMatrix<Number>> & matrices)
48 for(
auto & matrix : matrices)
52 matrix.compute_lu_factorization();
54 catch(std::exception & exc)
58 for(
unsigned int i = 0; i < matrix.m(); ++i)
60 for(
unsigned int j = 0; j < matrix.n(); ++j)
63 matrix(i, j) += 1.e-4;