# Parameter definitions for thorn IDBrillData
# $Header$


shares:einstein

EXTENDS KEYWORD initial_data 
{
   "brilldata"   :: "Brill wave initial data"
} ""

EXTENDS BOOLEAN use_conformal 
{
}


private:

# Parameters for elliptic solve

KEYWORD brill_solver "Which elliptic solver to use"
{
  "sor"    ::  "Use SOR solver"
  "petsc"  ::  "Use PETSc solver"
  "bam"    ::  "Use bam solver"
} "sor"

KEYWORD brill_bound "Which boundary condition to use"
{
  "const" :: "constant boundary: set brill_const_v0"
  "robin" :: "Robin boundary: set brill_robin_falloff, brill_robin_inf"
} "const"

INT brill_robin_falloff "Fall-off of Robin BC"
{
  0: :: "any positive integer value"
} 1

REAL brill_const_v0 "Value of constant BC"
{
  : :: "anything goes"
} 1.0

REAL brill_robin_inf "Value at infinity of Robin BC"
{
  : :: "anything goes"
} 1.0

REAL brill_thresh "How far (absolute norm) to go"
{
  0.0: :: "Positive number please"
} 0.00001


# Brill wave parameters

BOOLEAN axisym "Axisymmetric initial data?"
{
} "yes"

INT brill_q  "Form of function q [0,1,2]"
{
  0:2 :: "Only cases 0,1,2 defined at the moment"
} 2

REAL brill_a "Brill wave: Amplitude"
{
  : :: "Anything"
} 0.0

INT brill_b "Brill wave: rho^b"
{
  : :: "Anything"
} 2

INT brill_c "Brill wave: (r^2 - r0^2)^(c/2)"
{  
  : :: "Anything"
} 2

REAL brill_rho0 "Brill wave: radius of torus in rho"
{
  : :: "Anything"
} 0.0

REAL brill_r0 "Brill wave: radius of torus in r"
{
  : :: "Anything" 
} 0.0

REAL brill_srho "Brill wave: sigma in rho"
{
  : :: "Anything"
} 1.0

REAL brill_sr "Brill wave: sigma in r"
{
  : :: "Anything"
} 1.0

REAL brill_sz "Brill wave: sigma in z"
{
  : :: "Anything"
} 1.0

# 3D Brill wave parameters

REAL brill_d "3D Brill wave:  d rho^m cos^2(n (phi + phi0))"
{
  : :: "Anything"
} 0.0

REAL brill_e "3D Brill wave:  d rho^m cos^2(n (phi + phi0))"
{
  : :: "Anything"
} 1.0

REAL brill_m "3D Brill wave:  d rho^m cos^2(n (phi + phi0))"
{
  : :: "Anything"
} 2.0

REAL brill_n "3D Brill wave:  d rho^m cos^2(n (phi + phi0))"
{
  : :: "Anything"
} 2.0

REAL brill_phi0 "3D Brill wave:  d rho^m cos^2(n (phi + phi0))"
{
  : :: "Anything"
} 0.0


# Additional parameters

REAL brill_rhofudge "delta rho for axis fudge"
{
  0: :: "Positive please"
} 0.00001

INT sor_maxit "Maximum number of iterations"
{
0:* :: "Positive"
} 100