A lot of changes.
 - the output did work
 - routines for elementwise (triangles) and pointwise representation
 of the calculation domain
 - two ways of output