Executable: EvolveGhSingleBlackHole

---

Parallelization:
  ElementDistribution: NumGridPointsAndGridSpacing


InitialData: &InitialData
  GeneralizedHarmonic(KerrSchild):
    Mass: 1.0
    Spin: [0.0, 0.0, 0.0] 
    Center: [0.0, 0.0, 0.0]

DomainCreator:
  Sphere:
    InnerRadius: &InnerRadius 1.8 
    OuterRadius: 1000.0
    Interior:
      ExciseWithBoundaryCondition:
        DemandOutgoingCharSpeeds:
    InitialRefinement:
      Shell0: [2, 2, 4] 
      Shell1: [1, 1, 5] 
    InitialGridPoints: 10
    UseEquiangularMap: True
    EquatorialCompression: None
    RadialPartitioning: [50.0]
    RadialDistribution: [Logarithmic, Linear]
    WhichWedges: All
    TimeDependentMaps: 
      InitialTime: 0.0
      ShapeMap:
        LMax: &LMax 16
        InitialValues:
          Mass: 1.0
          Spin: [0.0, 0.0, 0.0] 
        SizeInitialValues: [0.0, 0.0, 0.0]
      RotationMap: None
      ExpansionMap:
        InitialValues: [1.0, 0.0, 0.0]
        InitialValuesOuterBoundary: [1.0, -1.0e-6, 0.0]
        DecayTimescaleExpansion: 1.0
        DecayTimescaleExpansionOuterBoundary: 1.0
      TranslationMap: None
    OuterBoundaryCondition:
      ConstraintPreservingBjorhus:
        Type: ConstraintPreservingPhysical

Evolution:
  InitialTime: 0.0
  InitialTimeStep: 0.0002
  InitialSlabSize: 0.1
  StepChoosers:
    - Increase:
        Factor: 2
    - ElementSizeCfl:
        SafetyFactor: 0.95
    - ErrorControl:
        AbsoluteTolerance: 1e-8
        RelativeTolerance: 1e-6
        MaxFactor: 2
        MinFactor: 0.25
        SafetyFactor: 0.95
  TimeStepper:
    AdamsMoultonPcMonotonic:
      Order: 4

EvolutionSystem:
  GeneralizedHarmonic:
    GaugeCondition: Harmonic


    DampingFunctionGamma0:
      TimeDependentTripleGaussian:
        Constant: 0.001
        Gaussian1:
          Amplitude: 7.0
          Width: 2.5                 
          Center: [0.0, 0.0, 0.0]
        Gaussian2:
          Amplitude: 0.1             
          Width: 100.0                 
          Center: [0.0, 0.0, 0.0]
        Gaussian3:
          Amplitude: 0.0
          Width: 100.0
          Center: [0.0, 0.0, 0.0]
    DampingFunctionGamma1:
      Constant:
        Value: -1.0
    DampingFunctionGamma2:
      TimeDependentTripleGaussian:
        Constant: 0.001
        Gaussian1:
          Amplitude: 7.0
          Width: 2.5                 
          Center: [0.0, 0.0, 0.0]
        Gaussian2:
          Amplitude: 0.1             
          Width: 100.0                 
          Center: [0.0, 0.0, 0.0]
        Gaussian3:
          Amplitude: 0.0
          Width: 100.0
          Center: [0.0, 0.0, 0.0]

Filtering:
  ExpFilter0:
    Alpha: 64.0
    HalfPower: 210
    Enable: true
    BlocksToFilter: All

SpatialDiscretization:
  BoundaryCorrection:
    UpwindPenalty:
  DiscontinuousGalerkin:
    Formulation: StrongInertial
    Quadrature: GaussLobatto

Observers:
  VolumeFileName: "BlackHoleVolume"
  ReductionFileName: "BlackHoleReductions"
  SurfaceFileName: "BlackHoleSurfaces"

Amr:
  Criteria:
  Policies:
    EnforceTwoToOneBalanceInNormalDirection: true
    Isotropy: Anisotropic
    Limits:
      RefinementLevel: Auto
      NumGridPoints: Auto
  Verbosity: Quiet

PhaseChangeAndTriggers:
  - Trigger:
      Slabs:
       EvenlySpaced:
         Interval: 100
         Offset: 0
    PhaseChanges:
      - CheckpointAndExitAfterWallclock:
          WallclockHours: 10000

EventsAndTriggers:
  - Trigger:
      Slabs:
        EvenlySpaced:
          Interval: 10
          Offset: 0
    Events:
      - ObserveTimeStep:
          SubfileName: TimeSteps
          PrintTimeToTerminal: True
          ObservePerCore: False
      - ObserveNorms:
          SubfileName: Norms
          TensorsToObserve:
          - Name: Lapse
            NormType: L2Norm
            Components: Individual
          - Name: PointwiseL2Norm(GaugeConstraint)
            NormType: L2Norm
            Components: Sum
          - Name: PointwiseL2Norm(ThreeIndexConstraint)
            NormType: L2Norm
            Components: Sum
  - Trigger:
      Slabs:
        EvenlySpaced:
          Interval: 10
          Offset: 0
    Events:
      - ErrorIfDataTooBig:
          Threshold: 100
          VariablesToCheck: [SpacetimeMetric]
  - Trigger:
      Slabs:
        EvenlySpaced:
          Interval: 1
          Offset: 0
    Events:
      - ApparentHorizon
      - ExcisionBoundary
  - Trigger:
      TimeCompares:
        Comparison: GreaterThanOrEqualTo
        Value: 2000.0
    Events:
      - Completion

EventsAndDenseTriggers:
  - Trigger:
      Times:
        EvenlySpaced:
          Interval: 0.1
          Offset: 0.0
    Events:
      - BondiSachsInterpolation

Interpolator:
  DumpVolumeDataOnFailure: false

ApparentHorizons:
  ApparentHorizon: &Ah
    InitialGuess:
      LMax: *LMax 
      Radius: 2.2
      Center: [0., 0., 0.]
    FastFlow:
      Flow: Fast
      Alpha: 1.0
      Beta: 0.5
      AbsTol: 1e-12
      TruncationTol: 1e-2
      DivergenceTol: 1.2
      DivergenceIter: 5
      MaxIts: 100
    Verbosity: Verbose
  ControlSystemSingleAh: *Ah
  ControlSystemCharSpeedAh: *Ah

InterpolationTargets:
  BondiSachsInterpolation:
    LMax: *LMax
    Radius: [100, 150, 200]
    Center: [0, 0, 0]
    AngularOrdering: Cce
  ExcisionBoundary: &ExBdry
    LMax: *LMax
    Center: [0., 0., 0.]
    Radius: *InnerRadius
    AngularOrdering: "Strahlkorper"
  ControlSystemCharSpeedExcision: *ExBdry

Cce:
  BondiSachsOutputFilePrefix: "BondiSachs"

ControlSystems:
  WriteDataToDisk: true
  MeasurementsPerUpdate: 4
  Verbosity: Silent
  Shape:
    IsActive: false 
    Averager:
      AverageTimescaleFraction: 0.25
      Average0thDeriv: false
    Controller:
      UpdateFraction: 0.2 
    TimescaleTuner:
      InitialTimescales: 0.2
      MinTimescale: 1.0e-2
      MaxTimescale: 10.0
      IncreaseThreshold: 2.5e-4
      DecreaseThreshold: 1.0e-3
      IncreaseFactor: 1.01
      DecreaseFactor: 0.98
    ControlError:
  Size:
    IsActive: true
    Averager:
      AverageTimescaleFraction: 0.25
      Average0thDeriv: true
    Controller:
      UpdateFraction: 0.1
    TimescaleTuner:
      InitialTimescales: 0.2
      MinTimescale: 1.0e-4
      MaxTimescale: 20.0
      IncreaseThreshold: 2.5e-4
      IncreaseFactor: 1.01
    ControlError:
      MaxNumTimesForZeroCrossingPredictor: 4
      SmoothAvgTimescaleFraction: 0.25
      DeltaRDriftOutwardOptions: None
      InitialState: DeltaR
      SmootherTuner:
        InitialTimescales: [0.2]
        MinTimescale: 1.0e-4
        MaxTimescale: 20.0
        IncreaseThreshold: 2.5e-4
        DecreaseThreshold: 1.0e-3
        IncreaseFactor: 1.01
        DecreaseFactor: 0.9

ResourceInfo:
  AvoidGlobalProc0: false
  Singletons: Auto