Peano
Loading...
Searching...
No Matches
523_higher-order-fd.py File Reference

Go to the source code of this file.

Namespaces

namespace  523_higher-order-fd
 
namespace  _higher
 -order-fd
 

Functions

 523_higher-order-fd.create_scheme (name, solver)
 The numerical scheme, grid construction and plotter of one run.
 
 523_higher-order-fd.build_and_run (name, directory, peano_installation, end_time)
 Generate, compile and run one setup within directory.
 
 523_higher-order-fd.probe (directory, point, output_file)
 Probe all snapshots at one point.
 
 523_higher-order-fd.profile (directory)
 Probe lapse and phi in the centre of each volume or grid cell along the line y=z=h/2, i.e.
 
 523_higher-order-fd.radius (x)
 
 523_higher-order-fd.time_steps (directory)
 Number of time steps and average wall clock time per time step.
 
 523_higher-order-fd.plot_over_time_and_profile (runs, names, file_name)
 Lapse next to the puncture over time (left) and along the x-axis at the last snapshot all runs share (right)
 
 523_higher-order-fd.plot_fourier_analysis ()
 Fourier analysis of the fourth-order central differences and of the Kreiss-Oliger dissipation.
 

Variables

 523_higher-order-fd.Here = os.path.dirname(os.path.abspath(__file__))
 
 523_higher-order-fd.PeanoRoot = os.path.abspath(os.path.join(Here, "..", "..", "..", ".."))
 
str 523_higher-order-fd.PatchFile = "solution-CCZ4.peano-patch-file"
 
float 523_higher-order-fd.Mass = 1.0
 
float 523_higher-order-fd.DomainOffset = -9.0
 
float 523_higher-order-fd.DomainSize = 18.0
 
int 523_higher-order-fd.CellsPerAxis = 9
 
int 523_higher-order-fd.PatchSize = 10
 
int 523_higher-order-fd.VolumesPerAxis = CellsPerAxis * PatchSize
 
float 523_higher-order-fd.H = DomainSize / VolumesPerAxis
 
int 523_higher-order-fd.Threads = 4
 
float 523_higher-order-fd.TimeInBetweenPlots = 0.5
 
 523_higher-order-fd.DissipationSubset = list(range(0, 23)) + [53, 54]
 Unknowns which the restricted Kreiss-Oliger dissipation damps: Q0-Q22 (G, K, theta, Z, lapse, shift, b) and Q53-Q54 (traceK, phi)
 
dict 523_higher-order-fd.Runs
 The runs: name -> (Kreiss-Oliger sigma or None, lopsided advection, quantities of the dissipation or None for all).
 
dict 523_higher-order-fd.Style
 Colours, markers and labels per run, so all pictures use the same ones.
 
 523_higher-order-fd.parser
 
 523_higher-order-fd.default
 
 523_higher-order-fd.help
 
 523_higher-order-fd.type
 
 523_higher-order-fd.float
 
 523_higher-order-fd.nargs
 
 523_higher-order-fd.args = parser.parse_args()
 
 523_higher-order-fd.work_directory
 
dict 523_higher-order-fd.runs = {}
 
 523_higher-order-fd.directory = os.path.join(os.path.abspath(work_directory), name)
 
 523_higher-order-fd.runtime
 
 523_higher-order-fd.success
 
 523_higher-order-fd.x
 
 523_higher-order-fd.t
 
 523_higher-order-fd.lapse
 
 523_higher-order-fd.phi
 
int 523_higher-order-fd.centre = VolumesPerAxis // 2
 
 523_higher-order-fd.statistics = time_steps(directory)
 
list 523_higher-order-fd.ko_runs = [name for name in Style if name.startswith("fd4-ko") and name in runs]
 
list 523_higher-order-fd.lopsided_runs