Loading [MathJax]/extensions/TeX/newcommand.js
 Reference documentation for deal.II version 9.1.1
\newcommand{\dealcoloneq}{\mathrel{\vcenter{:}}=}
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Modules Pages
step-38.h
1 
1060  std::ofstream x("x"), y("y");
1061  GridOut().write_gnuplot(volume_mesh, x);
1062  GridOut().write_gnuplot(triangulation, y);
1063  }
1064 
1065  std::cout << "Surface mesh has " << triangulation.n_active_cells()
1066  << " cells."
1067  << std::endl;
1068  ...
1069 }
1070 @endcode
1071 
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
1079 MappingQ(2), MappingQ(15) or MappingQ1, each cell of our mesh will be mapped
1080 using a bilinear approximation.
1081 
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
1088 number.
1089 <a name="PlainProg"></a>
1090 <h1> The plain program</h1>
1091 @include "step-38.cc"
1092  */
void write_gnuplot(const Triangulation< dim, spacedim > &tria, std::ostream &out, const Mapping< dim, spacedim > *mapping=nullptr) const
Definition: grid_out.cc:4221
VectorizedArray< Number > sin(const ::VectorizedArray< Number > &x)
VectorizedArray< Number > cos(const ::VectorizedArray< Number > &x)