# 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 1 1
1 0.333333 1 1
2 0.666667 1 1
3 1 1 1
4 0 0.666667 1
5 0.333333 0.666667 0
6 0.666667 0.666667 0
7 1 0.666667 1
8 0 0.333333 1
9 0.333333 0.333333 0
10 0.666667 0.333333 0
11 1 0.333333 1
12 0 0 1
13 0.333333 0 1
14 0.666667 0 1
15 1 0 1
#ELEMENTS:
#globalElementNumber  #node1  #node2  #node3  #mValue  #cValue  #fValue
0 0 4 5 1 0 1
1 0 1 5 1 0 1
2 1 5 6 1 0 1
3 1 2 6 1 0 1
4 2 6 7 1 0 1
5 2 3 7 1 0 1
6 4 8 9 1 0 1
7 4 5 9 1 0 1
8 5 9 10 1 0 1
9 5 6 10 1 0 1
10 6 10 11 1 0 1
11 6 7 11 1 0 1
12 8 12 13 1 0 1
13 8 9 13 1 0 1
14 9 13 14 1 0 1
15 9 10 14 1 0 1
16 10 14 15 1 0 1
17 10 11 15 1 0 1
#