1056 triangulation.set_manifold (0);
1058 std::ofstream x(
"x"), y(
"y");
1063 std::cout <<
"Surface mesh has " << triangulation.n_active_cells()
1070 Note that the only essential addition has been the three lines marked with
1071 asterisks. It is worth pointing out one other thing here, though: because we
1072 detach the manifold description from the surface mesh, whenever we use a
1073 mapping
object in the rest of the program, it has no curves boundary
1074 description to go on any more. Rather, it will have to use the implicit,
1075 FlatManifold class that is used on all parts of the domain not
1076 explicitly assigned a different manifold object. Consequently, whether we use
1078 using a bilinear approximation.
1080 All these drawbacks aside, the resulting pictures are still pretty. The only
1081 other differences to what's in @ref step_38 "step-38" is that we changed the right hand side
1082 to @f$f(\mathbf x)=
\sin x_3@f$ and the boundary values (through the
1083 <code>Solution</code>
class) to @f$u(\mathbf x)|_{\partial\Omega}=
\cos x_3@f$. Of
1084 course, we now non longer know the exact solution, so the computation of the
1085 error at the end of <code>LaplaceBeltrami::run</code> will yield a meaningless
1087 <a name=
"PlainProg"></a>
1088 <h1> The plain program</h1>
1089 @include
"step-38.cc"
VectorizedArray< Number > sin(const ::VectorizedArray< Number > &x)
void write_gnuplot(const Triangulation< dim, spacedim > &tria, std::ostream &out, const Mapping< dim, spacedim > *mapping=nullptr) const
VectorizedArray< Number > cos(const ::VectorizedArray< Number > &x)