Optimization loop over mesh subsets nodes

Substitute find_if with sorted vector of nodes.
Again this trades space vs. time.

This trades a copy of nodes pointers vs. O(#Elements * #Nodes/2) complexity.

Calculate surface flux; small fixes

Fix a bug in the initialization of the secondary variable

Solve Axially Symmetric Problems