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