From a3c5013be65ff20f4c8a91315e2e2be89c502e58 Mon Sep 17 00:00:00 2001 From: kronbichler Date: Sun, 20 Apr 2014 20:35:35 +0000 Subject: [PATCH] Next step in identifying the problem on the tester: Output the some of the values actually computed git-svn-id: https://svn.dealii.org/trunk@32796 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/matrix_free/step-48b.cc | 38 ++++++++++++++++++++++++++++++- tests/matrix_free/step-48b.output | 4 ++++ 2 files changed, 41 insertions(+), 1 deletion(-) diff --git a/tests/matrix_free/step-48b.cc b/tests/matrix_free/step-48b.cc index 4b28f2366b..1a2e35acd6 100644 --- a/tests/matrix_free/step-48b.cc +++ b/tests/matrix_free/step-48b.cc @@ -120,6 +120,7 @@ namespace Step48 { AssertDimension (src.size(), 2); FEEvaluationGL current (data), old (data); + deallog << "submitted values: "; for (unsigned int cell=cell_range.first; cell::n_array_elements; ++v) + deallog << current.begin_values()[0][v] << " "; current.integrate (true,true); current.distribute_local_to_global (dst); } + deallog << std::endl; + + deallog << "no sine term: "; + // now output similar terms again but without the sine term (and do not + // write anything back to the solution vector) + for (unsigned int cell=cell_range.first; cell current_value = current.get_value(q); + const VectorizedArray old_value = old.get_value(q); + + current.submit_value (2.*current_value - old_value, q); + current.submit_gradient (- delta_t_sqr * + current.get_gradient(q), q); + } + // output first value on quadrature point for all vector components + // (should be stable irrespective of vectorization width) + for (unsigned int v=0; v::n_array_elements; ++v) + deallog << current.begin_values()[0][v] << " "; + } + deallog << std::endl; } @@ -250,7 +286,7 @@ namespace Step48 QGaussLobatto<1> quadrature (fe_degree+1); typename MatrixFree::AdditionalData additional_data; additional_data.tasks_parallel_scheme = - MatrixFree::AdditionalData::partition_partition; + MatrixFree::AdditionalData::none; matrix_free_data.reinit (dof_handler, constraints, quadrature, additional_data); diff --git a/tests/matrix_free/step-48b.output b/tests/matrix_free/step-48b.output index 600934fc3d..20f0dcd6e1 100644 --- a/tests/matrix_free/step-48b.output +++ b/tests/matrix_free/step-48b.output @@ -4,7 +4,11 @@ DEAL:: Number of degrees of freedom: 289 DEAL:: Time step size: 0.2500, finest cell: 7.500 DEAL:: DEAL:: Time: -10.0 , solution norm: 7.11111164 +DEAL::submitted values: 0 0 0 0 0 0 0 0 0 0 0 0 0.569151584 0 0 0 +DEAL::no sine term: 0 0 0 0 0 0 0 0 0 0 0 0 0.562500000 0 0 0 DEAL:: Time: -9.75 , solution norm: 3.78861141 +DEAL::submitted values: 0 0 0 0 0 0 -0.00121093741 0 0 -0.00121093741 0 0 0.222882619 -0.00121093741 -0.00121093741 0 +DEAL::no sine term: 0 0 0 0 0 0 -0.00124999978 0 0 -0.00124999978 0 0 0.225803189 -0.00124999978 -0.00124999978 0 DEAL:: Time: -9.50 , solution norm: 0.869590104 DEAL:: DEAL:: Performed 3 time steps. -- 2.39.5