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
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
// initialize with zero
for (unsigned i=0;i<gradients.size();++i)
- std::fill_n (gradients[i].begin(), gradients[i].size(), Tensor<1,dim>());
+ fill_n (gradients[i].begin(), gradients[i].size(), Tensor<1,dim>());
// add up contributions of trial
// functions
// initialize with zero
for (unsigned i=0;i<second_derivs.size();++i)
- std::fill_n (second_derivs[i].begin(), second_derivs[i].size(), Tensor<2,dim>());
+ fill_n (second_derivs[i].begin(), second_derivs[i].size(), Tensor<2,dim>());
// add up contributions of trial
// functions
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
| update_JxW_values))
for (unsigned int i=0;i<boundary_forms.size();++i)
{
- double f = std::sqrt(contract(boundary_forms[i],
- boundary_forms[i]));
+ double f = sqrt(contract(boundary_forms[i],
+ boundary_forms[i]));
if (update_flags & update_JxW_values)
{
JxW_values[i] = f * weights[i];