# Triangular finite element mesh for the domain [0,1]^2
# for solving the equation -m*Laplace(u) + c*u = f 
# with homegeneous Dirichlet boundary conditions.
# The functions m,c,f are taken to be piecewise constant on elements.
#
#NODES:
#globalNodeNumber  #xCoordinate  #yCoordinate  #onBoundary
0 0.0 0.0 1
1 0.0 1.0 1
2 1.0 0.0 1
3 1.0 1.0 1
4 0.5 0.33333333 0
5 0.5 0.66666667 0
6 0.0 0.5 1
7 1.0 0.5 1
#ELEMENTS:
#globalElementNumber  #node1  #node2  #node3  #mValue  #cValue  #fValue
0 1 3 5 1.0 1.0 1.0
1 1 5 6 1.0 2.0 2.0
2 3 5 7 1.0 2.0 2.0
3 4 5 6 1.0 0.0 0.0
4 4 5 7 1.0 0.0 0.0
5 0 4 6 1.0 2.0 2.0
6 2 4 7 1.0 2.0 2.0
7 0 2 4 1.0 1.0 1.0
#