default:
Assert (false, ExcNotImplemented());
return "";
- };
+ }
}
for (unsigned int i=0; i<vertices.size(); ++i)
if (vertex_used[i])
- {
out << '\t' << vertices[i] << '\n';
- };
// write cells or faces
const bool write_cells = dx_flags.write_cells;
for (unsigned int d=spacedim+1; d<=3; ++d)
out << " 0"; // fill with zeroes
out << std::endl;
- };
+ }
// Write cells preamble
out << "$ENDNOD" << std::endl
++vertex)
out << cell->vertex_index(GeometryInfo<dim>::ucd_to_deal[vertex])+1 << ' ';
out << std::endl;
- };
+ }
// write faces and lines with
// non-zero boundary indicator
<< "# For a description of the UCD format see the AVS Developer's guide."
<< '\n'
<< "#" << '\n';
- };
+ }
// start with ucd data
out << n_vertices << ' '
for (unsigned int d=spacedim+1; d<=3; ++d)
out << " 0"; // fill with zeroes
out << '\n';
- };
+ }
// write cells. Enumerate cells
// consecutively, starting with 1
break;
default:
Assert (false, ExcNotImplemented());
- };
+ }
// it follows a list of the
// vertices of each cell. in 1d
++vertex)
out << cell->vertex_index(GeometryInfo<dim>::ucd_to_deal[vertex])+1 << ' ';
out << '\n';
- };
+ }
// write faces and lines with
// non-zero boundary indicator
break;
default:
Assert (false, ExcNotImplemented ());
- };
+ }
out << "\n";
// (iii) write vertex ordering
break;
default:
Assert (false, ExcNotImplemented ());
- };
+ }
// (iv) write a list of vertices of cells
out << "\n#"
break;
default:
Assert (false, ExcNotImplemented());
- };
+ }
out << static_cast<unsigned int>(face->boundary_indicator())
<< ' '
<< static_cast<unsigned int>(face->boundary_indicator())
out << '\n';
++index;
- };
+ }
}
break;
default:
Assert (false, ExcNotImplemented());
- };
+ }
// note: vertex numbers are 1-base
for (unsigned int vertex=0; vertex<GeometryInfo<dim>::vertices_per_face; ++vertex)
out << face->vertex_index(GeometryInfo<dim-1>::ucd_to_deal[vertex])+1 << ' ';
out << '\n';
++index;
- };
+ }
}
{
Assert(false, ExcInternalError());
break;
- };
+ }
case 2:
{
line_list.push_back (LineEntry(p0, p1,
face->user_flag_set(),
cell->level()));
- };
- };
- };
+ }
+ }
+ }
break;
- };
+ }
case 3:
{
y_max = std::max (y_max, line->second(1));
max_level = std::max (max_level, line->level);
- };
+ }
// scale in x-direction such that
// in the output 0 <= x <= 300.
"/MCshow { currentpoint stroke m\n"
"exch dup MFwidth -2 div 3 -1 roll R MFshow } def\n")
<< '\n';
- };
+ }
out << "%%EndProlog" << '\n'
<< '\n';
// set fine lines
out << eps_flags_base.line_width << " setlinewidth" << '\n';
- };
+ }
// now write the lines
const Point<2> offset(x_min, y_min);
out << ")] "
<< "] -6 MCshow"
<< '\n';
- };
- };
+ }
+ }
// and the vertex numbers
if ((dim == 2) && (eps_flags_2.write_vertex_numbers == true))
<< ")] "
<< "] -6 MCshow"
<< '\n';
- };
- };
+ }
+ }
out << "showpage" << '\n';