clc; clear; close all;

Globals2D;

tic;
Problem;

N = 4; % degree of polynomial representing solution within an element (local solution representation)

StartUp2D; % computing local operators and mesh information

NonreflectingBCStartup;

[mapI, vmapI, mapO, vmapO] = BuildBCMaps();

NumericalFluxStartup;
preproccesing_time = toc;

tic;
[rho,u,v,p] = Timestepping_PML(T); % time-stepping for the ODE system
%[rho,u,v,p] = Timestepping_sponge(T); % time-stepping for the ODE system
computation_time = toc;