From aef3c3ee32736c8c9d99dba7d9d4724916759739 Mon Sep 17 00:00:00 2001 From: wolf Date: Mon, 2 Jul 2001 16:11:15 +0000 Subject: [PATCH] Undo last check-in, as it was with the wrong message. git-svn-id: https://svn.dealii.org/trunk@4792 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_values.cc | 8 ++++---- deal.II/deal.II/source/fe/mapping_c1.cc | 16 ++++++++-------- deal.II/deal.II/source/fe/mapping_q1.cc | 4 ++-- 3 files changed, 14 insertions(+), 14 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_values.cc b/deal.II/deal.II/source/fe/fe_values.cc index 2014f91d77..bb2f61d398 100644 --- a/deal.II/deal.II/source/fe/fe_values.cc +++ b/deal.II/deal.II/source/fe/fe_values.cc @@ -141,7 +141,7 @@ void FEValuesBase::get_function_values (const InputVector &fe_fu present_cell->get_interpolated_dof_values(fe_function, dof_values); // initialize with zero - std::fill_n (values.begin(), n_quadrature_points, 0); + fill_n (values.begin(), n_quadrature_points, 0); // add up contributions of trial // functions @@ -265,7 +265,7 @@ get_function_grads (const InputVector &fe_function, present_cell->get_interpolated_dof_values(fe_function, dof_values); // initialize with zero - std::fill_n (gradients.begin(), n_quadrature_points, Tensor<1,dim>()); + fill_n (gradients.begin(), n_quadrature_points, Tensor<1,dim>()); // add up contributions of trial // functions @@ -306,7 +306,7 @@ get_function_grads (const InputVector // initialize with zero for (unsigned i=0;i()); + fill_n (gradients[i].begin(), gradients[i].size(), Tensor<1,dim>()); // add up contributions of trial // functions @@ -387,7 +387,7 @@ get_function_2nd_derivatives (const InputVector // initialize with zero for (unsigned i=0;i()); + fill_n (second_derivs[i].begin(), second_derivs[i].size(), Tensor<2,dim>()); // add up contributions of trial // functions diff --git a/deal.II/deal.II/source/fe/mapping_c1.cc b/deal.II/deal.II/source/fe/mapping_c1.cc index 797d0fe3c0..cc20ddb6f3 100644 --- a/deal.II/deal.II/source/fe/mapping_c1.cc +++ b/deal.II/deal.II/source/fe/mapping_c1.cc @@ -142,14 +142,14 @@ MappingC1<2>::add_line_support_points (const Triangulation<2>::cell_iterator &ce coordinate_axis /= h; const double alpha = std::atan2(coordinate_axis[1], coordinate_axis[0]); - const double b = ((face_vertex_normals[0][1] * std::sin(alpha) - +face_vertex_normals[0][0] * std::cos(alpha)) / - (face_vertex_normals[0][1] * std::cos(alpha) - -face_vertex_normals[0][0] * std::sin(alpha))), - c = ((face_vertex_normals[1][1] * std::sin(alpha) - +face_vertex_normals[1][0] * std::cos(alpha)) / - (face_vertex_normals[1][1] * std::cos(alpha) - -face_vertex_normals[1][0] * std::sin(alpha))); + const double b = ((face_vertex_normals[0][1] * sin(alpha) + +face_vertex_normals[0][0] * cos(alpha)) / + (face_vertex_normals[0][1] * cos(alpha) + -face_vertex_normals[0][0] * sin(alpha))), + c = ((face_vertex_normals[1][1] * sin(alpha) + +face_vertex_normals[1][0] * cos(alpha)) / + (face_vertex_normals[1][1] * cos(alpha) + -face_vertex_normals[1][0] * sin(alpha))); // next evaluate the so // determined cubic diff --git a/deal.II/deal.II/source/fe/mapping_q1.cc b/deal.II/deal.II/source/fe/mapping_q1.cc index 78873e1739..e992094f51 100644 --- a/deal.II/deal.II/source/fe/mapping_q1.cc +++ b/deal.II/deal.II/source/fe/mapping_q1.cc @@ -679,8 +679,8 @@ MappingQ1::compute_fill_face (const typename DoFHandler::cell_iterator | update_JxW_values)) for (unsigned int i=0;i