Default: SurfaceTesselationRegularBoundary_cRadius = -1.0 (i.e. no surface tesselation is executed)
The surface tesselation is performed in the following way. First, each boundary points obtains a circular disc
with the radius SurfaceTesselationRegularBoundary_cRadius*SmoothingLength. Then the discs cut each other, such that, in the oprimal case,
we obtain a set of non-overlapping tesselation cells.