Current State:
 Working program:
 adaptively polygonizes a parametric curve to minimize open angle between two consecutive line segments, currently uses random numbers so that no two approximations are alike.
 using a galerkin finite element method, computes general relativistic initial data for axisymmetric spacetime using the above curve as a generator for a marginally trapped surface of revolution.
 uses bisection method to find point on the actual parametric curve that is closest to a given boundary point on the approximating polygon
 Parameters:
 theta0 = opening angle criteria, the closer this is to pi or 180 degrees, the more accurate the approximating polygon will be.
 hmax = "maximum size" of a given finite element, problems occur if this is taken too small
 grade = maximum ratio of the areas of two adjacent elements
 rmax = radius at which the discretization truncates spacetime.
 Interface:
 program is called either from the command line (usage statement is included) or from a perl script which reads a config file and reads the above parameters.
 Specification of curve:
 hard coded in source code and specified through the make file.
