From f9817637169ff56aba3bc6766e18767fbc009516 Mon Sep 17 00:00:00 2001 From: Luca Heltai Date: Tue, 31 Mar 2015 21:07:32 +0200 Subject: [PATCH] Fixed step-34. --- examples/step-34/step-34.cc | 30 ++++++++++++++++-------------- 1 file changed, 16 insertions(+), 14 deletions(-) diff --git a/examples/step-34/step-34.cc b/examples/step-34/step-34.cc index 09641afb9d..84e6c1f8c4 100644 --- a/examples/step-34/step-34.cc +++ b/examples/step-34/step-34.cc @@ -42,7 +42,7 @@ #include #include #include -#include +#include #include #include @@ -183,7 +183,7 @@ namespace Step34 // // Experimenting a little with the computation of the angles gives very // accurate results for simpler geometries. To verify this you can comment - // out, in the read_domain() method, the tria.set_boundary(1, boundary) + // out, in the read_domain() method, the tria.set_manifold(1, manifold) // line, and check the alpha that is generated by the program. By removing // this call, whenever the mesh is refined new nodes will be placed along // the straight lines that made up the coarse mesh, rather than be pulled @@ -481,21 +481,22 @@ namespace Step34 // normals to the mesh are external to both the circle in 2d or the sphere // in 3d. // - // The other detail that is required for appropriate refinement of the - // boundary element mesh, is an accurate description of the manifold that - // the mesh is approximating. We already saw this several times for the - // boundary of standard finite element meshes (for example in step-5 and - // step-6), and here the principle and usage is the same, except that the - // HyperBallBoundary class takes an additional template parameter that - // specifies the embedding space dimension. The function object still has to - // be static to live at least as long as the triangulation object to which - // it is attached. + // The other detail that is required for appropriate refinement of + // the boundary element mesh, is an accurate description of the + // manifold that the mesh is approximating. We already saw this + // several times for the boundary of standard finite element meshes + // (for example in step-5 and step-6), and here the principle and + // usage is the same, except that the SphericalManifold class takes + // an additional template parameter that specifies the embedding + // space dimension. The function object still has to be static to + // live at least as long as the triangulation object to which it is + // attached. template void BEMProblem::read_domain() { static const Point center = Point(); - static const HyperBallBoundary boundary(center,1.); + static const SphericalManifold manifold(center); std::ifstream in; switch (dim) @@ -515,8 +516,9 @@ namespace Step34 GridIn gi; gi.attach_triangulation (tria); gi.read_ucd (in); - - tria.set_boundary(1, boundary); + + tria.set_all_manifold_ids(1); + tria.set_manifold(1, manifold); } -- 2.39.5