NormalsRegularizationBorders

This program sets up two regularization matrices for two different regional areas. For a given set of points defined by grid it is evaluated, whether each point (corresponding to an unknown parameter of a respective parameterization by space localizing basis functions) is inside or outside a certain area given by border. Each regularization matrix is a diagonal matrix, one of them features a one if the point is inside, and a zero if the point lies outside the area. The other matrix features a zero if the point is inside, and a one if the point lies outside the area This results in two regularization matrices with \[ \M R_1+\M R_2=\M I. \]The two matrices are provided as vectors of the diagonal in the output files outputfileOutside and outputfileInside. The regularization matrices are then used by normalEquation:regularization. As an example, the two different areas could be oceanic regions on the one hand and continental areas on the other hand.

Name	Type	Annotation
outputfileInside	filename
outputfileOutside	filename
grid	grid	nodal point distribution of parameters, e.g harmonics splines
border	border	regularization areas, e.g land and ocean
R	double	reference radius for ellipsoidal coordinates for border
inverseFlattening	double	reference flattening for ellipsoidal coordinates for border, 0: spherical coordinates