# Executable: SolveXcts # Version: 2022.02.17 # Branch with additional optimizations for Fig. 18: # https://github.com/nilsvu/spectre/tree/elliptic_solver_paper_optimizations # The commented-out options in this file are only available on that branch. Background: Kerr: Mass: 1. Spin: [0., 0., 0.] Center: [0., 0., 0.] InitialGuess: Flatness DomainCreator: Shell: InnerRadius: 2. OuterRadius: 10. RadialPartitioning: [] RadialDistribution: [Logarithmic] # Placeholders for mesh refinement: # - L: h-refinement level # - P: Polynomial degree InitialRefinement: {{ L }} InitialGridPoints: [{{ P + 1 }}, {{ P + 1 }}] UseEquiangularMap: True EquatorialCompression: None WhichWedges: All TimeDependence: None BoundaryConditions: InnerBoundary: AnalyticSolution: ConformalFactor: Neumann LapseTimesConformalFactor: Dirichlet ShiftExcess: Dirichlet OuterBoundary: AnalyticSolution: ConformalFactor: Dirichlet LapseTimesConformalFactor: Dirichlet ShiftExcess: Dirichlet Discretization: DiscontinuousGalerkin: PenaltyParameter: 1. Massive: True Observers: VolumeFileName: "KerrSchildVolume" ReductionFileName: "KerrSchildReductions" NonlinearSolver: NewtonRaphson: ConvergenceCriteria: MaxIterations: 20 RelativeResidual: 0. AbsoluteResidual: 1.e-10 SufficientDecrease: 1.e-4 MaxGlobalizationSteps: 40 DampingFactor: 1. Verbosity: Verbose LinearSolver: Gmres: ConvergenceCriteria: MaxIterations: 300 RelativeResidual: 1.e-4 AbsoluteResidual: 9.e-11 Verbosity: Quiet Multigrid: # Set to 0 to disable preconditioning altogether (Fig. 15, dotted line) Iterations: 1 # Set to 1 to disable multigrid (Fig. 15, dashed line) MaxLevels: Auto PreSmoothing: True # Set to False when MaxLevels is 1 to run the Schwarz smoother only once # (Fig. 15, dashed line) PostSmoothingAtBottom: True Verbosity: Silent OutputVolumeData: False # This optimization option is irrelevant, since the grid is p-uniform # ElementAllocation: # WeightByNumPoints: True SchwarzSmoother: MaxOverlap: 2 Iterations: 3 Verbosity: Silent SubdomainSolver: Gmres: ConvergenceCriteria: # Set MaxIteration to >20 for Fig. 17 MaxIterations: 3 RelativeResidual: 1.e-4 AbsoluteResidual: 1.e-10 Verbosity: Silent Restart: None Preconditioner: MinusLaplacian: Solver: ExplicitInverse: # Enable optimizations for Fig. 18 # FillFactor: 1 # Verbosity: Silent # Set to 'Dirichlet' to run with only Dirichlet-Laplacians (Fig. 17, # black dots and shaded region) BoundaryConditions: Auto ObservePerCoreReductions: False SkipResets: True EventsAndTriggers: ? EveryNIterations: N: 1 Offset: 0 : - ObserveErrorNorms: SubfileName: ErrorNorms