SHORELINE
Triangulate the Domain where 0 <= F(X,Y)

SHORELINE, a MATLAB library which tries to identify and triangulate the 2D domain over which some function F(X,Y) is nonnegative.

It does this by starting with a set of sample points, where the function is evaluated. Triangles are dropped from the triangulation if all 3 vertices have negative functions. Triangles with 3 positive vertices are assumed (to within the requested resolution) to be part of the region. Therefore, a more accurate answer can be sought by refining those triangles whose vertices include both positive and negative values.

This results in a refined "shoreline" around the positive triangles. The function is evaluated at the new nodes, and again the negative triangles are dropped, the positive ones kept, and the intermediate triangles can be further refined, if desired.

It is also possible to estimate the integral of some function g(x,y) over the region implicitly defined by 0 <= f(x,y).

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Languages:

SHORELINE is available in a MATLAB version.

Related Data and Programs:

CONTOUR_GRADIENT, a MATLAB program which displays contours and gradient vectors for a function f(x,y).

EDGE, a MATLAB library which defines some test functions in 1D, 2D and 3D for the detection of edges.

HYPERSPHERE_SURFACE, a MATLAB library which illustrates a procedure for estimating the location of a hypersurface defined by a characteristic function or a signed function.

LEVELS, a MATLAB library which makes a contour plot, choosing the contour levels using random sampling.

MATLAB_SURF, a MATLAB library which demonstrates the MATLAB surf() function for displaying a 3D surface of the form Z=F(X,Y).

shoreline_test

TOMS626, a FORTRAN77 library which creates a contour plot of data associated with a triangulated set of points;
this is ACM TOMS algorithm 626.

TRIANGULATION_ORDER3_CONTOUR, a MATLAB program which makes contour plot of scattered data, or of data defined on an order 3 triangulation.

TWOD_PLOTC, a MATLAB program which makes a color contour plot of scalar data defined on a 2D finite element mesh, by Jeff Borggaard.

Source Code:

• circular_normal_pdf.m, evaluates the Circular Normal PDF.
• cull_negative_triangles.m, removes from the triangulation those triangles for which all vertices have a negative function value.
• i4mat_write.m, writes an I4MAT to a file.
• product_grid.m, samples a rectangle with a product grid.
• r8mat_write.m, writes an R8MAT to a file.
• refine_compute.m, computes the refinement of a triangulation.
• refine_size.m, determines the "size" of the refinement of a triangulation.
• shoreline.m, traces out a triangulated region for which 0 <= F(X,Y).
• timestamp.m, prints the current YMDHMS date as a time stamp.
• triangulation_area.m, returns the area of a triangulation.
• triangulation_quad.m, estimates the integral of a function over a triangulated region.