From c45fb96de45971f7c02dc943d58d596746a3a023 Mon Sep 17 00:00:00 2001 From: bangerth Date: Tue, 18 Aug 2009 21:24:31 +0000 Subject: [PATCH] Remove assertions. It is true that they are satisfied all the time but just not on very coarse meshes where sometimes we have quadrature points inside the inner ring because we only use a Q1 mapping. git-svn-id: https://svn.dealii.org/trunk@19307 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-32/step-32.cc | 13 +++++-------- 1 file changed, 5 insertions(+), 8 deletions(-) diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc index e757edb002..f1f46a3a61 100644 --- a/deal.II/examples/step-32/step-32.cc +++ b/deal.II/examples/step-32/step-32.cc @@ -138,9 +138,6 @@ namespace EquationData const double rho = (r-R0)/h; - Assert (rho >= -.001, ExcInternalError()); - Assert (rho <= 1.001, ExcInternalError()); - return T1+(T0-T1)*((1-rho)*(1-rho)); } @@ -1306,7 +1303,7 @@ compute_viscosity (const std::vector &old_temperature, // side is cheap anyway so we won't even // notice that this part is not parallized // by threads. - // + // // Regarding the implementation of // inhomogeneous Dirichlet boundary // conditions: Since we use the temperature @@ -2462,7 +2459,7 @@ local_assemble_temperature_rhs (const std::pair global_T_range, time_step * nu * ext_grad_T * scratch.grad_phi_T[i] + - time_step * + time_step * scratch.gamma_values[q] * scratch.phi_T[i]) * scratch.temperature_fe_values.JxW(q); @@ -2470,15 +2467,15 @@ local_assemble_temperature_rhs (const std::pair global_T_range, if (temperature_constraints.is_inhomogeneously_constrained(data.local_dof_indices[i])) { for (unsigned int j=0; j