From 765a32e41cc7985ec8ecbec2698abb4262dbc8df Mon Sep 17 00:00:00 2001 From: Matthias Maier Date: Thu, 27 Aug 2015 23:29:01 -0500 Subject: [PATCH] Testsuite: another iteration... --- .../memory_consumption_01.cc | 0 ...nsumption_01.with_64bit_indices=off.output | 0 ...onsumption_01.with_64bit_indices=on.output | 0 tests/deal.II/helmholtz.h | 26 ------- tests/deal.II/helmholtz1.th | 69 ------------------- tests/deal.II/helmholtz1mg.th | 65 ----------------- tests/{deal.II => dofs}/accessor_01.cc | 0 tests/{deal.II => dofs}/accessor_01.output | 0 tests/{deal.II => dofs}/accessor_02.cc | 0 tests/{deal.II => dofs}/accessor_02.output | 0 tests/{deal.II => dofs}/block_info.cc | 0 tests/{deal.II => dofs}/block_info.output | 0 tests/{deal.II => dofs}/block_info_02.cc | 0 tests/{deal.II => dofs}/block_info_02.output | 0 tests/{deal.II => dofs}/block_list.h | 0 tests/{deal.II => dofs}/block_list_01.cc | 0 tests/{deal.II => dofs}/block_list_01.output | 0 tests/{deal.II => dofs}/block_list_02.cc | 0 tests/{deal.II => dofs}/block_list_02.output | 0 tests/{deal.II => dofs}/block_list_03.cc | 0 tests/{deal.II => dofs}/block_list_03.output | 0 tests/{deal.II => dofs}/block_list_04.cc | 0 tests/{deal.II => dofs}/block_list_04.output | 0 .../{deal.II => dofs}/count_dofs_on_patch.cc | 0 .../count_dofs_on_patch.output | 0 tests/{deal.II => dofs}/mg_dof_handler.cc | 0 tests/{deal.II => dofs}/mg_dof_handler.output | 0 .../cell_similarity_crash_01.cc | 0 .../cell_similarity_crash_01.output | 0 .../cell_similarity_crash_02.cc | 0 .../cell_similarity_crash_02.output | 0 tests/{deal.II => fe}/curl_curl_01.cc | 0 tests/{deal.II => fe}/curl_curl_01.output | 0 tests/{deal.II => grid}/coarsening_02.cc | 0 tests/{deal.II => grid}/coarsening_02.output | 0 tests/{deal.II => grid}/coarsening_02_1d.cc | 0 .../{deal.II => grid}/coarsening_02_1d.output | 0 tests/{deal.II => grid}/coarsening_03.cc | 0 tests/{deal.II => grid}/coarsening_03.output | 0 tests/{deal.II => grid}/distort_random.cc | 0 tests/{deal.II => grid}/distort_random.output | 0 .../{deal.II => grid}/extent_in_direction.cc | 0 .../extent_in_direction.output | 0 .../{deal.II => grid}/face_orientations_3d.cc | 0 .../face_orientations_3d.output | 0 tests/{deal.II => grid}/filtered_iterator.cc | 0 .../filtered_iterator.output | 0 .../{deal.II => grid}/filtered_iterator_02.cc | 0 .../filtered_iterator_02.output | 0 .../{deal.II => grid}/filtered_iterator_03.cc | 0 .../filtered_iterator_03.output | 0 .../get_finest_common_cells_01.cc | 0 .../get_finest_common_cells_01.output | 0 .../get_finest_common_cells_02.cc | 0 .../get_finest_common_cells_02.output | 0 .../get_finest_common_cells_03.cc | 0 .../get_finest_common_cells_03.output | 0 .../get_finest_common_cells_04.cc | 0 .../get_finest_common_cells_04.output | 0 .../get_patch_around_cell_01.cc | 0 .../get_patch_around_cell_01.output | 0 .../have_same_coarse_mesh_01.cc | 0 .../have_same_coarse_mesh_01.output | 0 .../have_same_coarse_mesh_02.cc | 0 .../have_same_coarse_mesh_02.output | 0 .../have_same_coarse_mesh_03.cc | 0 .../have_same_coarse_mesh_03.output | 0 .../have_same_coarse_mesh_04.cc | 0 .../have_same_coarse_mesh_04.output | 0 .../intergrid_constraints.cc | 0 .../intergrid_constraints.output | 0 tests/{deal.II => grid}/intergrid_map.cc | 0 tests/{deal.II => grid}/intergrid_map.output | 0 tests/{deal.II => grid}/line_coarsening_3d.cc | 0 .../line_coarsening_3d.output | 0 .../maximal_cell_diameter.cc | 0 .../maximal_cell_diameter.output | 0 tests/{deal.II => grid}/measure_et_al.cc | 0 tests/{deal.II => grid}/measure_et_al.output | 0 tests/{deal.II => grid}/measure_et_al_02.cc | 0 .../{deal.II => grid}/measure_et_al_02.output | 0 .../measure_of_3d_face_01.cc | 0 .../measure_of_3d_face_01.output | 0 tests/{deal.II => grid}/mesh_smoothing_02.cc | 0 .../mesh_smoothing_02.output | 0 tests/{deal.II => grid}/mesh_smoothing_03.cc | 0 .../mesh_smoothing_03.output | 0 .../minimal_cell_diameter.cc | 0 .../minimal_cell_diameter.output | 0 .../normals_at_vertices_01.cc | 0 .../normals_at_vertices_01.output | 0 .../normals_at_vertices_02.cc | 0 .../normals_at_vertices_02.output | 0 tests/{deal.II => grid}/number_cache_01.cc | 0 .../{deal.II => grid}/number_cache_01.output | 0 tests/{deal.II => grid}/number_cache_02.cc | 0 .../{deal.II => grid}/number_cache_02.output | 0 tests/{deal.II => grid}/number_cache_03.cc | 0 .../{deal.II => grid}/number_cache_03.output | 0 tests/{deal.II => grid}/vertex_as_face_01.cc | 0 .../vertex_as_face_01.output | 0 tests/{deal.II => grid}/vertex_as_face_02.cc | 0 .../vertex_as_face_02.output | 0 tests/{deal.II => grid}/vertex_as_face_03.cc | 0 .../vertex_as_face_03.output | 0 tests/{deal.II => grid}/vertex_as_face_04.cc | 0 .../vertex_as_face_04.output | 0 tests/{deal.II => grid}/vertex_as_face_05.cc | 0 .../vertex_as_face_05.output | 0 tests/{deal.II => grid}/vertex_as_face_06.cc | 0 .../vertex_as_face_06.output | 0 tests/{deal.II => grid}/vertex_as_face_07.cc | 0 .../vertex_as_face_07.output | 0 tests/{deal.II => grid}/vertex_as_face_08.cc | 0 .../vertex_as_face_08.output | 0 tests/{deal.II => grid}/vertex_as_face_09.cc | 0 .../vertex_as_face_09.output | 0 tests/{deal.II => grid}/vertex_as_face_10.cc | 0 .../vertex_as_face_10.output | 0 tests/{deal.II => grid}/vertex_as_face_11.cc | 0 .../vertex_as_face_11.output | 0 tests/{deal.II => grid}/vertex_as_face_12.cc | 0 .../vertex_as_face_12.output | 0 tests/{deal.II => grid}/vertex_as_face_13.cc | 0 .../vertex_as_face_13.output | 0 .../block_matrices_04.cc} | 0 .../block_matrices_04.output} | 0 .../filtered_matrix_06.cc} | 0 .../filtered_matrix_06.output} | 0 .../matrices.cc => lac/matrices_02.cc} | 0 .../matrices_02.output} | 0 .../assemble_matrix_parallel_01.cc | 0 .../assemble_matrix_parallel_01.output | 0 .../assemble_matrix_parallel_02.cc | 0 .../assemble_matrix_parallel_02.output | 0 .../assemble_matrix_parallel_03.cc | 0 .../assemble_matrix_parallel_03.output | 0 .../assemble_matrix_parallel_04.cc | 0 .../assemble_matrix_parallel_04.output | 0 tests/{deal.II => numerics}/boundaries.cc | 0 tests/{deal.II => numerics}/boundaries.output | 0 tests/{deal.II => numerics}/data_out.cc | 0 tests/{deal.II => numerics}/data_out.output | 0 .../data_out_postprocessor_01.cc | 0 .../data_out_postprocessor_01.output | 0 .../data_out_postprocessor_scalar_01.cc | 0 .../data_out_postprocessor_scalar_01.output | 0 .../data_out_postprocessor_vector_01.cc | 0 .../data_out_postprocessor_vector_01.output | 0 .../derivative_approximation.cc | 0 .../derivative_approximation.output | 0 .../derivative_approximation_2.cc | 0 .../derivative_approximation_2.output | 0 tests/{deal.II => numerics}/derivatives.cc | 0 .../{deal.II => numerics}/derivatives.output | 0 .../{deal.II => numerics}/error_estimator.cc | 0 .../error_estimator.output | 0 .../error_estimator_01.cc | 0 .../error_estimator_01.output | 0 .../error_estimator_02.cc | 0 .../error_estimator_02.output | 0 tests/{deal.II => numerics}/kelly_crash_02.cc | 0 .../kelly_crash_02.output | 0 .../tangential_flux_inhom_01.cc | 0 .../tangential_flux_inhom_01.output | 0 .../vectors_boundary_rhs_01.cc | 0 .../vectors_boundary_rhs_01.output | 0 .../vectors_boundary_rhs_02.cc | 0 .../vectors_boundary_rhs_02.output | 0 .../vectors_boundary_rhs_03.cc | 0 .../vectors_boundary_rhs_03.output | 0 171 files changed, 160 deletions(-) rename tests/{deal.II => base}/memory_consumption_01.cc (100%) rename tests/{deal.II => base}/memory_consumption_01.with_64bit_indices=off.output (100%) rename tests/{deal.II => base}/memory_consumption_01.with_64bit_indices=on.output (100%) delete mode 100644 tests/deal.II/helmholtz.h delete mode 100644 tests/deal.II/helmholtz1.th delete mode 100644 tests/deal.II/helmholtz1mg.th rename tests/{deal.II => dofs}/accessor_01.cc (100%) rename tests/{deal.II => dofs}/accessor_01.output (100%) rename tests/{deal.II => dofs}/accessor_02.cc (100%) rename tests/{deal.II => dofs}/accessor_02.output (100%) rename tests/{deal.II => dofs}/block_info.cc (100%) rename tests/{deal.II => dofs}/block_info.output (100%) rename tests/{deal.II => dofs}/block_info_02.cc (100%) rename tests/{deal.II => dofs}/block_info_02.output (100%) rename tests/{deal.II => dofs}/block_list.h (100%) rename tests/{deal.II => dofs}/block_list_01.cc (100%) rename tests/{deal.II => dofs}/block_list_01.output (100%) rename tests/{deal.II => dofs}/block_list_02.cc (100%) rename tests/{deal.II => dofs}/block_list_02.output (100%) rename tests/{deal.II => dofs}/block_list_03.cc (100%) rename tests/{deal.II => dofs}/block_list_03.output (100%) rename tests/{deal.II => dofs}/block_list_04.cc (100%) rename tests/{deal.II => dofs}/block_list_04.output (100%) rename tests/{deal.II => dofs}/count_dofs_on_patch.cc (100%) rename tests/{deal.II => dofs}/count_dofs_on_patch.output (100%) rename tests/{deal.II => dofs}/mg_dof_handler.cc (100%) rename tests/{deal.II => dofs}/mg_dof_handler.output (100%) rename tests/{deal.II => fe}/cell_similarity_crash_01.cc (100%) rename tests/{deal.II => fe}/cell_similarity_crash_01.output (100%) rename tests/{deal.II => fe}/cell_similarity_crash_02.cc (100%) rename tests/{deal.II => fe}/cell_similarity_crash_02.output (100%) rename tests/{deal.II => fe}/curl_curl_01.cc (100%) rename tests/{deal.II => fe}/curl_curl_01.output (100%) rename tests/{deal.II => grid}/coarsening_02.cc (100%) rename tests/{deal.II => grid}/coarsening_02.output (100%) rename tests/{deal.II => grid}/coarsening_02_1d.cc (100%) rename tests/{deal.II => grid}/coarsening_02_1d.output (100%) rename tests/{deal.II => grid}/coarsening_03.cc (100%) rename tests/{deal.II => grid}/coarsening_03.output (100%) rename tests/{deal.II => grid}/distort_random.cc (100%) rename tests/{deal.II => grid}/distort_random.output (100%) rename tests/{deal.II => grid}/extent_in_direction.cc (100%) rename tests/{deal.II => grid}/extent_in_direction.output (100%) rename tests/{deal.II => grid}/face_orientations_3d.cc (100%) rename tests/{deal.II => grid}/face_orientations_3d.output (100%) rename tests/{deal.II => grid}/filtered_iterator.cc (100%) rename tests/{deal.II => grid}/filtered_iterator.output (100%) rename tests/{deal.II => grid}/filtered_iterator_02.cc (100%) rename tests/{deal.II => grid}/filtered_iterator_02.output (100%) rename tests/{deal.II => grid}/filtered_iterator_03.cc (100%) rename tests/{deal.II => grid}/filtered_iterator_03.output (100%) rename tests/{deal.II => grid}/get_finest_common_cells_01.cc (100%) rename tests/{deal.II => grid}/get_finest_common_cells_01.output (100%) rename tests/{deal.II => grid}/get_finest_common_cells_02.cc (100%) rename tests/{deal.II => grid}/get_finest_common_cells_02.output (100%) rename tests/{deal.II => grid}/get_finest_common_cells_03.cc (100%) rename tests/{deal.II => grid}/get_finest_common_cells_03.output (100%) rename tests/{deal.II => grid}/get_finest_common_cells_04.cc (100%) rename tests/{deal.II => grid}/get_finest_common_cells_04.output (100%) rename tests/{deal.II => grid}/get_patch_around_cell_01.cc (100%) rename tests/{deal.II => grid}/get_patch_around_cell_01.output (100%) rename tests/{deal.II => grid}/have_same_coarse_mesh_01.cc (100%) rename tests/{deal.II => grid}/have_same_coarse_mesh_01.output (100%) rename tests/{deal.II => grid}/have_same_coarse_mesh_02.cc (100%) rename tests/{deal.II => grid}/have_same_coarse_mesh_02.output (100%) rename tests/{deal.II => grid}/have_same_coarse_mesh_03.cc (100%) rename tests/{deal.II => grid}/have_same_coarse_mesh_03.output (100%) rename tests/{deal.II => grid}/have_same_coarse_mesh_04.cc (100%) rename tests/{deal.II => grid}/have_same_coarse_mesh_04.output (100%) rename tests/{deal.II => grid}/intergrid_constraints.cc (100%) rename tests/{deal.II => grid}/intergrid_constraints.output (100%) rename tests/{deal.II => grid}/intergrid_map.cc (100%) rename tests/{deal.II => grid}/intergrid_map.output (100%) rename tests/{deal.II => grid}/line_coarsening_3d.cc (100%) rename tests/{deal.II => grid}/line_coarsening_3d.output (100%) rename tests/{deal.II => grid}/maximal_cell_diameter.cc (100%) rename tests/{deal.II => grid}/maximal_cell_diameter.output (100%) rename tests/{deal.II => grid}/measure_et_al.cc (100%) rename tests/{deal.II => grid}/measure_et_al.output (100%) rename tests/{deal.II => grid}/measure_et_al_02.cc (100%) rename tests/{deal.II => grid}/measure_et_al_02.output (100%) rename tests/{deal.II => grid}/measure_of_3d_face_01.cc (100%) rename tests/{deal.II => grid}/measure_of_3d_face_01.output (100%) rename tests/{deal.II => grid}/mesh_smoothing_02.cc (100%) rename tests/{deal.II => grid}/mesh_smoothing_02.output (100%) rename tests/{deal.II => grid}/mesh_smoothing_03.cc (100%) rename tests/{deal.II => grid}/mesh_smoothing_03.output (100%) rename tests/{deal.II => grid}/minimal_cell_diameter.cc (100%) rename tests/{deal.II => grid}/minimal_cell_diameter.output (100%) rename tests/{deal.II => grid}/normals_at_vertices_01.cc (100%) rename tests/{deal.II => grid}/normals_at_vertices_01.output (100%) rename tests/{deal.II => grid}/normals_at_vertices_02.cc (100%) rename tests/{deal.II => grid}/normals_at_vertices_02.output (100%) rename tests/{deal.II => grid}/number_cache_01.cc (100%) rename tests/{deal.II => grid}/number_cache_01.output (100%) rename tests/{deal.II => grid}/number_cache_02.cc (100%) rename tests/{deal.II => grid}/number_cache_02.output (100%) rename tests/{deal.II => grid}/number_cache_03.cc (100%) rename tests/{deal.II => grid}/number_cache_03.output (100%) rename tests/{deal.II => grid}/vertex_as_face_01.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_01.output (100%) rename tests/{deal.II => grid}/vertex_as_face_02.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_02.output (100%) rename tests/{deal.II => grid}/vertex_as_face_03.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_03.output (100%) rename tests/{deal.II => grid}/vertex_as_face_04.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_04.output (100%) rename tests/{deal.II => grid}/vertex_as_face_05.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_05.output (100%) rename tests/{deal.II => grid}/vertex_as_face_06.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_06.output (100%) rename tests/{deal.II => grid}/vertex_as_face_07.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_07.output (100%) rename tests/{deal.II => grid}/vertex_as_face_08.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_08.output (100%) rename tests/{deal.II => grid}/vertex_as_face_09.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_09.output (100%) rename tests/{deal.II => grid}/vertex_as_face_10.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_10.output (100%) rename tests/{deal.II => grid}/vertex_as_face_11.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_11.output (100%) rename tests/{deal.II => grid}/vertex_as_face_12.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_12.output (100%) rename tests/{deal.II => grid}/vertex_as_face_13.cc (100%) rename tests/{deal.II => grid}/vertex_as_face_13.output (100%) rename tests/{deal.II/block_matrices.cc => lac/block_matrices_04.cc} (100%) rename tests/{deal.II/block_matrices.output => lac/block_matrices_04.output} (100%) rename tests/{deal.II/filtered_matrix.cc => lac/filtered_matrix_06.cc} (100%) rename tests/{deal.II/filtered_matrix.output => lac/filtered_matrix_06.output} (100%) rename tests/{deal.II/matrices.cc => lac/matrices_02.cc} (100%) rename tests/{deal.II/matrices.output => lac/matrices_02.output} (100%) rename tests/{deal.II => numerics}/assemble_matrix_parallel_01.cc (100%) rename tests/{deal.II => numerics}/assemble_matrix_parallel_01.output (100%) rename tests/{deal.II => numerics}/assemble_matrix_parallel_02.cc (100%) rename tests/{deal.II => numerics}/assemble_matrix_parallel_02.output (100%) rename tests/{deal.II => numerics}/assemble_matrix_parallel_03.cc (100%) rename tests/{deal.II => numerics}/assemble_matrix_parallel_03.output (100%) rename tests/{deal.II => numerics}/assemble_matrix_parallel_04.cc (100%) rename tests/{deal.II => numerics}/assemble_matrix_parallel_04.output (100%) rename tests/{deal.II => numerics}/boundaries.cc (100%) rename tests/{deal.II => numerics}/boundaries.output (100%) rename tests/{deal.II => numerics}/data_out.cc (100%) rename tests/{deal.II => numerics}/data_out.output (100%) rename tests/{deal.II => numerics}/data_out_postprocessor_01.cc (100%) rename tests/{deal.II => numerics}/data_out_postprocessor_01.output (100%) rename tests/{deal.II => numerics}/data_out_postprocessor_scalar_01.cc (100%) rename tests/{deal.II => numerics}/data_out_postprocessor_scalar_01.output (100%) rename tests/{deal.II => numerics}/data_out_postprocessor_vector_01.cc (100%) rename tests/{deal.II => numerics}/data_out_postprocessor_vector_01.output (100%) rename tests/{deal.II => numerics}/derivative_approximation.cc (100%) rename tests/{deal.II => numerics}/derivative_approximation.output (100%) rename tests/{deal.II => numerics}/derivative_approximation_2.cc (100%) rename tests/{deal.II => numerics}/derivative_approximation_2.output (100%) rename tests/{deal.II => numerics}/derivatives.cc (100%) rename tests/{deal.II => numerics}/derivatives.output (100%) rename tests/{deal.II => numerics}/error_estimator.cc (100%) rename tests/{deal.II => numerics}/error_estimator.output (100%) rename tests/{deal.II => numerics}/error_estimator_01.cc (100%) rename tests/{deal.II => numerics}/error_estimator_01.output (100%) rename tests/{deal.II => numerics}/error_estimator_02.cc (100%) rename tests/{deal.II => numerics}/error_estimator_02.output (100%) rename tests/{deal.II => numerics}/kelly_crash_02.cc (100%) rename tests/{deal.II => numerics}/kelly_crash_02.output (100%) rename tests/{deal.II => numerics}/tangential_flux_inhom_01.cc (100%) rename tests/{deal.II => numerics}/tangential_flux_inhom_01.output (100%) rename tests/{deal.II => numerics}/vectors_boundary_rhs_01.cc (100%) rename tests/{deal.II => numerics}/vectors_boundary_rhs_01.output (100%) rename tests/{deal.II => numerics}/vectors_boundary_rhs_02.cc (100%) rename tests/{deal.II => numerics}/vectors_boundary_rhs_02.output (100%) rename tests/{deal.II => numerics}/vectors_boundary_rhs_03.cc (100%) rename tests/{deal.II => numerics}/vectors_boundary_rhs_03.output (100%) diff --git a/tests/deal.II/memory_consumption_01.cc b/tests/base/memory_consumption_01.cc similarity index 100% rename from tests/deal.II/memory_consumption_01.cc rename to tests/base/memory_consumption_01.cc diff --git a/tests/deal.II/memory_consumption_01.with_64bit_indices=off.output b/tests/base/memory_consumption_01.with_64bit_indices=off.output similarity index 100% rename from tests/deal.II/memory_consumption_01.with_64bit_indices=off.output rename to tests/base/memory_consumption_01.with_64bit_indices=off.output diff --git a/tests/deal.II/memory_consumption_01.with_64bit_indices=on.output b/tests/base/memory_consumption_01.with_64bit_indices=on.output similarity index 100% rename from tests/deal.II/memory_consumption_01.with_64bit_indices=on.output rename to tests/base/memory_consumption_01.with_64bit_indices=on.output diff --git a/tests/deal.II/helmholtz.h b/tests/deal.II/helmholtz.h deleted file mode 100644 index d35cf8b080..0000000000 --- a/tests/deal.II/helmholtz.h +++ /dev/null @@ -1,26 +0,0 @@ - -#include "../tests.h" -#include - -/** - * Generator for system matrix and right hand side. - */ -class Helmholtz -{ -public: - template - void build_all(MATRIX &A, - VECTOR &f, - const DoFHandler &dof, - const Quadrature &quadrature, - const Function &rhs); - - template - void build_mgmatrix(MGLevelObject &A, - const DoFHandler &dof, - const Quadrature &quadrature); -}; - -#include "helmholtz1.th" -#include "helmholtz1mg.th" - diff --git a/tests/deal.II/helmholtz1.th b/tests/deal.II/helmholtz1.th deleted file mode 100644 index 75dc5f8878..0000000000 --- a/tests/deal.II/helmholtz1.th +++ /dev/null @@ -1,69 +0,0 @@ - -template -void -Helmholtz::build_all(MATRIX& A, - VECTOR& f, - const DoFHandler& dof, - const Quadrature& quadrature, - const Function& rhs) -{ - const unsigned int fe_dofs = dof.get_fe().n_dofs_per_cell(); - FEValues fe(dof.get_fe(), quadrature, - UpdateFlags(update_values | update_gradients | - update_JxW_values | update_q_points)); - - Vector elvec(fe_dofs); - FullMatrix elmat(fe_dofs); - std::vector indices(fe_dofs); - -//////////////////////////////////////////////////// -// Loop for computing the element matrices -//////////////////////////////////////////////////// - - for (DoFHandler::active_cell_iterator c = dof.begin_active() - ; c != dof.end() ; ++c) - { - elmat.clear(); - elvec.clear(); - fe.reinit(c); - - // Quadrature loop - - for (unsigned k=0;k loc = fe.quadrature_point(k); - - // Test function loop - for (unsigned i=0;i dv = fe.shape_grad(i,k); - double v = fe.shape_value(i,k); - - elvec(i) += fe.JxW(k) * - rhs.value(loc) * v; - - //Trial function loop - for (unsigned j=0;j du = fe.shape_grad(j,k); - double u = fe.shape_value(j,k); - - elmat(i,j) += fe.JxW(k) * - (0.1*u*v + du*dv) - ; - } - } - } - - c->get_dof_indices(indices); - for (unsigned i=0;i -void -Helmholtz::build_mgmatrix(MGLevelObject& A, - const DoFHandler& dof, - const Quadrature& quadrature) -{ - const unsigned int fe_dofs = dof.get_fe().n_dofs_per_cell(); - FEValues fe(dof.get_fe(), quadrature, - UpdateFlags(update_values | update_gradients | - update_JxW_values | update_q_points)); - - Vector elvec(fe_dofs); - FullMatrix elmat(fe_dofs); - std::vector indices(fe_dofs); - -//////////////////////////////////////////////////// -// Loop for computing the element matrices -//////////////////////////////////////////////////// - - DoFHandler::cell_iterator dc = dof.DoFHandler::begin(); - for (DoFHandler::cell_iterator c = dof.begin() - ; c != dof.end() ; ++c, ++dc) - { - elmat.clear(); - elvec.clear(); - fe.reinit(dc); - - // Quadrature loop - - for (unsigned k=0;k loc = fe.quadrature_point(k); - - // Test function loop - for (unsigned i=0;i dv = fe.shape_grad(i,k); - double v = fe.shape_value(i,k); - - //Trial function loop - for (unsigned j=0;j du = fe.shape_grad(j,k); - double u = fe.shape_value(j,k); - - elmat(i,j) += fe.JxW(k) * - (0.1*u*v + du*dv) - ; - } - } - } - - const unsigned int level = c->level(); - - c->get_mg_dof_indices(indices); - for (unsigned i=0;i