Triangulation<dim> tr;
if (dim==2)
{
- GridGenerator::hyper_ball(tr);
+ GridGenerator::hyper_ball(tr, Point<dim>(), 1);
}
else
GridGenerator::hyper_cube(tr, -1./sqrt(dim),1./sqrt(dim));
deallog << "Interpolated boundary values" << std::endl;
std::map<unsigned int,double> interpolated_bv;
VectorTools::interpolate_boundary_values (mapping, dof, function_map,
- interpolated_bv);
+ interpolated_bv, std::vector<bool>());
write_map (interpolated_bv);
// project boundary values
Triangulation<dim> tr;
if (dim==2)
- GridGenerator::hyper_ball(tr);
+ GridGenerator::hyper_ball(tr, Point<dim>(), 1);
else
GridGenerator::hyper_cube(tr, -1,1);
tr.refine_global (1);
Triangulation<dim> tr;
if (dim==2)
- GridGenerator::hyper_ball(tr);
+ GridGenerator::hyper_ball(tr, Point<dim>(), 1);
else
GridGenerator::hyper_cube(tr, -1,1);
tr.refine_global (1);
Triangulation<dim> tr;
if (dim==2)
- GridGenerator::hyper_ball(tr);
+ GridGenerator::hyper_ball(tr, Point<dim>(), 1);
else
GridGenerator::hyper_cube(tr, -1,1);
tr.refine_global (1);
CosineFunction<dim> cosine;
if (dim==2)
- GridGenerator::hyper_ball(tr);
+ GridGenerator::hyper_ball(tr, Point<dim>(), 1);
else
GridGenerator::hyper_cube(tr, -1,1);
// interpolate boundary values
std::map<unsigned int,double> bv;
- VectorTools::interpolate_boundary_values (mapping, dof, 0, cosine, bv);
+ VectorTools::interpolate_boundary_values (mapping, dof, 0, cosine, bv, std::vector<bool>());
// the cosine has too many zero
// values on the boundary of the
// domain, so reset the elements to
{
Triangulation<dim> tr;
if (dim==2)
- GridGenerator::hyper_ball(tr);
+ GridGenerator::hyper_ball(tr, Point<dim>(), 1);
else
GridGenerator::hyper_cube(tr, -1,1);
tr.refine_global (1);
Triangulation<dim> tr;
if (dim==2)
{
- GridGenerator::hyper_ball(tr);
+ GridGenerator::hyper_ball(tr, Point<dim>(), 1);
}
else
GridGenerator::hyper_cube(tr, -1./sqrt(dim),1./sqrt(dim));