39 int main(
int argc,
char *argv[])
41 const char *new_mesh_file =
"mobius-strip.mesh";
47 bool visualization =
true;
50 args.
AddOption(&new_mesh_file,
"-m",
"--mesh-out-file",
51 "Output Mesh file to write.");
52 args.
AddOption(&nx,
"-nx",
"--num-elements-x",
53 "Number of elements in x-direction.");
54 args.
AddOption(&ny,
"-ny",
"--num-elements-y",
55 "Number of elements in y-direction.");
56 args.
AddOption(&order,
"-o",
"--mesh-order",
57 "Order (polynomial degree) of the mesh elements.");
58 args.
AddOption(&close_strip,
"-c",
"--close-strip",
59 "How to close the strip: 0 - open, 1 - closed, 2 - twisted.");
60 args.
AddOption(&dg_mesh,
"-dm",
"--discont-mesh",
"-cm",
"--cont-mesh",
61 "Use discontinuous or continuous space for the mesh nodes.");
63 "Number of twists of the strip.");
64 args.
AddOption(&visualization,
"-vis",
"--visualization",
"-no-vis",
66 "Enable or disable GLVis visualization.");
79 mesh =
new Mesh(nx, ny, el_type, 1, 2*M_PI, 2.0);
86 for (
int i = 0; i < v2v.Size(); i++)
91 for (
int j = 0; j <= ny; j++)
93 int v_old = nx + j * (nx + 1);
94 int v_new = ((close_strip == 1) ? j : (ny - j)) * (nx + 1);
98 for (
int i = 0; i < mesh->
GetNE(); i++)
103 for (
int j = 0; j < nv; j++)
109 for (
int i = 0; i < mesh->
GetNBE(); i++)
114 for (
int j = 0; j < nv; j++)
131 for (
int i = 0; i < nodes.
Size(); i++)
133 if (std::abs(nodes(i)) < 1e-12)
139 ofstream ofs(new_mesh_file);
146 char vishost[] =
"localhost";
149 sol_sock.precision(8);
150 sol_sock <<
"mesh\n" << *mesh << flush;
158 double a = 1.0 + 0.5 * (x[1] - 1.0) * cos(
num_twists * x[0] );
161 p[0] = a * cos( x[0] );
162 p[1] = a * sin( x[0] );
163 p[2] = 0.5 * (x[1] - 1.0) * sin(
num_twists * x[0] );
Class for grid function - Vector with associated FE space.
void mobius_trans(const Vector &x, Vector &p)
virtual void GetVertices(Array< int > &v) const =0
Returns element's vertices.
void SetSize(int s)
Resizes the vector if the new size is different.
int GetNBE() const
Returns number of boundary elements.
int Size() const
Returns the size of the vector.
void Transform(void(*f)(const Vector &, Vector &))
int GetNE() const
Returns number of elements.
void RemoveInternalBoundaries()
void SetCurvature(int order, bool discont=false, int space_dim=-1, int ordering=1)
Type
Constants for the classes derived from Element.
const Element * GetElement(int i) const
void PrintUsage(std::ostream &out) const
virtual void Print(std::ostream &out=std::cout) const
Print the mesh to the given stream using the default MFEM mesh format.
int main(int argc, char *argv[])
void RemoveUnusedVertices()
Remove unused vertices and rebuild mesh connectivity.
void AddOption(bool *var, const char *enable_short_name, const char *enable_long_name, const char *disable_short_name, const char *disable_long_name, const char *description, bool required=false)
void PrintOptions(std::ostream &out) const
virtual int GetNVertices() const =0
void GetNodes(Vector &node_coord) const
Abstract data type element.
const Element * GetBdrElement(int i) const