Main.py

# -*- coding: utf-8 -*-
"""
@author: Laura C. Kühle

Code-Style: E226, W503
Docstring-Style: D200, D400

"""
from DG_Approximation import DGScheme
import timeit

tic = timeit.default_timer()

alpha = 1
p = 2
cfl = 0.2
N = 256    # 40 elements work well for Condition 3
finalTime = 1
xL = -1
xR = 1

detector_config = {}
fold_len = 16
whisker_len = 3

detector_config['fold_len'] = fold_len
detector_config['whisker_len'] = whisker_len

init_config = {}
init_config['factor'] = 4
init_config['left_factor'] = 3

limiter_config = {}
limiter_config['mod_factor'] = 0
limiter_config['erase_factor'] = 1

quadrature_config = {}
quadrature_config['num_eval_points'] = 12  # 12

detector = 'Theoretical'
update_scheme = 'SSPRK3'
init_cond = 'Box'
limiter = 'ModifiedMinMod'
quadrature = 'Gauss'
dg_scheme = DGScheme(detector, detector_config=detector_config, init_cond=init_cond, init_config=init_config,
                     limiter=limiter, limiter_config=limiter_config, quadrature=quadrature,
                     quadrature_config=quadrature_config, update_scheme=update_scheme, wave_speed=alpha,
                     polynom_degree=p, cfl_number=cfl, num_grid_cells=N, final_time=finalTime,
                     show_plot=True, left_bound=xL, right_bound=xR)

# __, __, troubled_cells, __ =
dg_scheme.approximate()
dg_scheme.save_plots()

toc = timeit.default_timer()
print("Time:", toc-tic)
# =============================================================================
#
# print(troubled_cells)
#
# print(dg_scheme.update_scheme)
# print(dg_scheme.init_cond)
# print(dg_scheme.detector)
# print(dg_scheme.limiter)
# # print(dg_scheme.projection)
# =============================================================================

# if __name__ == '__main__':
#     pass