stencil           325 src/secOrd_op_brusselator_impl.cpp   Dune::Fem::DiagonalStencil<FE_space, FE_space> stencil {df_space, df_space};
stencil           326 src/secOrd_op_brusselator_impl.cpp   matrix.reserve(stencil);