Fixed bug in subSurface().

The actual bug was inside the LRSplineSurface constructor (called from
subSurface() only) which did not update the LRBSpline2D basis functions
with the mesh set in LRSplineSurface.
Added support for unit testing in lrsplines2D, with the bugfix as the first unit test.