* don't need them after calling this
* function, or if so store them.
*/
- virtual void distribute_dofs (const FiniteElement<dim> &, unsigned int offset=0);
+ virtual void distribute_dofs (const FiniteElement<dim> &,
+ const unsigned int offset = 0);
/**
* Actually do the renumbering based on
Tensor<1,dim> tangential;
tangential[(nindex-1)%dim] = 1.;
data.unit_tangentials[i+nfaces].resize(n_original_q_points);
- fill (data.unit_tangentials[i+nfaces].begin(),
- data.unit_tangentials[i+nfaces].end(),
- tangential);
+ std::fill (data.unit_tangentials[i+nfaces].begin(),
+ data.unit_tangentials[i+nfaces].end(),
+ tangential);
}
}
}
= (Point<dim>(1,0,0)
- ((Point<dim>(1,0,0) * view_direction) * view_direction)
- ((Point<dim>(1,0,0) * unit_vector1) * unit_vector1));
- const Point<dim> unit_vector2 = vector2 / sqrt(vector2.square());
+ const Point<dim> unit_vector2 = vector2 / std::sqrt(vector2.square());
for (; line!=endline; ++line)
line_list.push_back (LineEntry(Point<2>(line->vertex(0) * unit_vector2,
if (compute_radius_automatically)
{
const Point<dim> vertex_relative = p0 - center;
- r = sqrt(vertex_relative.square());
+ r = std::sqrt(vertex_relative.square());
}
else
r=radius;
Assert(std::fabs(v0.square()-r*r)<eps, ExcInternalError());
Assert(std::fabs(v1.square()-r*r)<eps, ExcInternalError());
- const double alpha=std::acos((v0*v1)/sqrt(v0.square()*v1.square()));
+ const double alpha=std::acos((v0*v1)/std::sqrt(v0.square()*v1.square()));
const double d_alpha=alpha/(n+1);
const Point<dim> pm=0.5*(v0+v1);
Assert(right_index<n, ExcInternalError());
Assert(left_index+1>=1, ExcInternalError());
- double d=h*tan(beta);
+ double d=h*std::tan(beta);
points[right_index]=pm+d/length*(v1-v0);
points[left_index]=pm-d/length*(v1-v0);
}
template <int dim>
void MGDoFHandler<dim>::distribute_dofs (const FiniteElement<dim> &fe,
- unsigned int offset)
+ const unsigned int offset)
{
// first distribute global dofs
DoFHandler<dim>::distribute_dofs (fe, offset);
double max_diag_entry = 0;
for (unsigned int i=0; i<dofs_per_cell; ++i)
if (std::fabs(cell_matrix(i,i)) > max_diag_entry)
- max_diag_entry = fabs(cell_matrix(i,i));
+ max_diag_entry = std::fabs(cell_matrix(i,i));
#endif
mutex.acquire ();