Reduce use of dynamically allocated vector: the algorithm only uses the first two points of the spacepoint vector to calculate a new theta, but was allocating a vector of dimension equal to the number of spacepoints.
Moved this function to anonymous namespace, as thetaFromSpacePoints.
@christos , relating to general task of reducing dynamic allocations.