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.
|