Reference documentation for deal.II version 8.1.0
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
step-38.h
1 
1073  triangulation.set_boundary (0); /* ** */
1074  GridTools::transform (&warp<spacedim>, triangulation); /* ** */
1075  std::ofstream x("x"), y("y");
1076  GridOut().write_gnuplot (volume_mesh, x);
1077  GridOut().write_gnuplot (triangulation, y);
1078  }
1079 
1080  std::cout << "Surface mesh has " << triangulation.n_active_cells()
1081  << " cells."
1082  << std::endl;
1083 
1084  ...
1085 @endcode
1086 
1087 Note that the only essential addition has been the three lines marked with
1088 asterisks. It is worth pointing out one other thing here, though: because we
1089 detach the manifold description from the surface mesh, whenever we use a
1090 mapping object in the rest of the program, it has no curves boundary
1091 description to go on any more. Rather, it will have to use the implicit,
1092 StraightBoundary class that is used on all parts of the boundary not
1093 explicitly assigned a different manifold object. Consequently, whether we use
1094 MappingQ(2), MappingQ(15) or MappingQ1, each cell of our mesh will be mapped
1095 using a bilinear approximation.
1096 
1097 All these drawbacks aside, the resulting pictures are still pretty. The only
1098 other differences to what's in @ref step_38 "step-38" is that we changed the right hand side
1099 to @f$f(\mathbf x)=\sin x_3@f$ and the boundary values (through the
1100 <code>Solution</code> class) to @f$u(\mathbf x)|_{\partial\Omega}=\cos x_3@f$. Of
1101 course, we now non longer know the exact solution, so the computation of the
1102 error at the end of <code>LaplaceBeltrami::run</code> will yield a meaningless
1103 number.
1104 <a name="PlainProg"></a>
1105 <h1> The plain program</h1>
1106 @include "step-38.cc"
1107  */
void write_gnuplot(const Triangulation< dim, spacedim > &tria, std::ostream &out, const Mapping< dim, spacedim > *mapping=0) const
VectorizedArray< Number > sin(const ::VectorizedArray< Number > &x)
void transform(const Transformation &transformation, Triangulation< dim, spacedim > &triangulation)
VectorizedArray< Number > cos(const ::VectorizedArray< Number > &x)