Procedures

Coefficient alpha from the recursive formulation of Legendre polynomials, for the Legendre mode 'mode'.

Prodcut of alpha(numerator) * beta(denominator) / alpha(denominator) as needed by the Clenshaw algorithm in ply_split_legendre_matrix.

Quotient of two alpha values.

Return the penalty parameter for the IP discretizations of higher order equations.

append a value to the dynamic array and return its position.

Append conservative variables for Euler equations.

appending a value to the dynamic array

appending a sorted list of values to the dynamic array

Append conservative variables for Heat equations.

Append conservative variables for Heat equations.

Append conservative variables for Heat equations.

Interface for fluxes of acoustic equations.

calculate flux of pure acoustic equation directly on the face-vector

Function for physical flux of the acoustic equation F, 1D? Since it is 1d, there need to be passed the correct background velocity (u0 for F - flux in x direction, v0 for G - flux in y direction, w0 for H - flux in z direction)

Interface for fluxes of acoustic equations.

calculate flux of pure acoustic equation directly on the face-vector

Function for physical flux of the acoustic equation F, 1D? Since it is 1d, there need to be passed the correct background velocity (u0 for F - flux in x direction, v0 for G - flux in y direction, w0 for H - flux in z direction)

Setup timers to assess the runtime of various parts of Ateles.

Setup timers to assess the runtime of various parts of Ateles.

Allocate temporary arrays.

\brief append the variables for acoustic simulations the acoustic equation only has 'primitive' variables or different speaking, the equation describe the pertubation in primitive variables

append the variables for acoustic simulations

Append conservative variables for Euler equations.

append / set data of derived quantities

Append primitive variables for euler equation

Append conservative variables for Euler equations.

append / set data of derived quantities

Append primitive variables for euler equation

Append conservative variables for Euler equations.

Append / set methods and data to compute derived quantities to the variable system.

Append primitive variables for euler equation

append / set data of derived quantities

append the variables for LinearEuler simulations

append / set data of derived quantities

append the variables for LinearEuler simulations

append the variables for electrodynamic simulations

append the variables for electrodynamic simulations

append the variables for electrodynamic simulations that include

append the variables for diffusion simulations

Adds the configured material variables to the variable system.

Append / set methods and data to compute derived quantities to the variable system.

Assign reference to spacetime functions to the affected elements. This means: The position of the variable in the variable system, which reflects the spacetime function is determined and assigned to the element_list.

Define material properties for the faces of all fluid elements and inherit the face material property to all ghost elements on the finer level.

A most basic flux function which uses the average of both sides for 1D Euler.

Prominently let the user now, what he actually is running right now.

Define the method to set the solverData_evalElem routine for stfuns.

Bubble sorting of array of real numbers of size n

Routine to check if the physicle values of the state are physically meaningful or not for the Euler and Linear Euler equation, checking the cfl for a fixed timestep

Routine to check if the physicle values of the state are physically meaningful or not for the Navier-Stokes equation, checking the cfl for a fixed timestep

Damp the modal coefficients of the state vector by a given spectral viscosity method.

Routine to check if the physical values of a state are physically meaningful or not.

Apply conservative limitation of denisty and energy to ensure positivity for density and pressure.

Apply conservative limitation of denisty and energy to ensure positivity for density and pressure.

Covolume filtering for 3D equation.

Covolume filtering for 1D equation.

Covolume filtering for 2D equation.

Project two co-volume elements onto single a single element.

Project two co-volume elements onto single a single element.

Project two co-volume elements onto single a single element.

Recursive routine to project the state from primal mesh to covolume mesh.

Recursive routine to project the state from primal mesh to covolume mesh.

Recursive routine to project the state from primal mesh to covolume mesh.

This routine creates hard coded variables which are required by ateles.

Create separate compute list for constant-constant, constant-variable (or vice versa) and variable-variable material parameter compute faces on this rank.

Create separate compute list for constant-constant, constant-variable (or vice versa) and variable-variable material parameter compute faces on this rank.

Subroutine to create element index list for constant and non-constant material parameters.

Deallocates the array for storing the sourceData for the currentLevel

Subroutine do define a specific stencil for a certain scheme.

This routine evaluates the lambda2 criterion. The input is the nodal value

This routine evaluates the q_criterion. The input is the nodal value of

Dump material description to lua file

Performance results are written to a file for statistical review The file-format is simple can be evaluated with gnuplot

Dump weights to a file.

Initialization of the Acoustic equation.

Reading boundary conditions for the acoustic equations.

Initialization of the linearized Euler equations.

Reading boundary conditions for the advection equations in 1D.

Convert conservative to primitive variables.

Convert conservative to primitive variables.

Convert primitive varibales to conservative variables.

Convert primitive varibales to conservative variables.

Convert conservative to primitive variables.

Convert conservative to primitive variables.

Convert conservative to primitive variables including the gradients. (instate(npnts, nScalars+nScalars)

Convert conservative to primitive variables (including temperature instead of pressure).

Convert conservative to conservative variables (including velocity instead of momentum).

Convert primitive varibales to conservative variables.

Convert primitive varibales to conservative variables.

Convert primitive varibales (including temperature instead of pressure) to conservative variables.

Convert conservative varibales (including velocity instead of temperature) to conservative variables.

Convert conservative to primitive variables.

Convert conservative to primitive variables.

Convert conservative to primitive variables including the gradients.

Convert conservative to primitive variables (including temperature instead of pressure).

Convert conservative to conservative variables (including velocity instead of momentum).

Solve the equation system with just the penalization terms to find an implicit update for the IMEX timestepping procedure.

Initialization of the Euler equation.

Reading boundary conditions for the euler equations.

Convert primitive varibales to conservative variables.

Convert primitive varibales to conservative variables.

Convert primitive varibales to conservative variables including their gradients.

Convert primitive varibales (including temperature instead of pressure) to conservative variables.

Convert conservative varibales (including velocity instead of temperature) to conservative variables.

Initialization of the Navier-Stokes equations.

Initialization of the Heat equation.

Reading boundary conditions for the Heat equation.

Initialization of the linearized Euler equations.

Reading boundary conditions for the LinearEuler equations.

Convert conservative to primitive variables.

Convert primitive varibales to conservative variables.

Convert conservative to primitive variables.

Solve the equation system with just the penalization terms to find an implicit update for the IMEX timestepping procedure.

Initialization of the Maxwell equations.

Subroutine to load boundary conditions for Maxwell equations from a Lua configuration file.

Convert primitive varibales to conservative variables.

Subroutine to load boundary conditions for Maxwell equations with divergence correction from a Lua configuration file. For the correction in E-field dirichlet bc and for correction in B-field neumann bc are defined.

Convert conservative to primitive variables.

Convert primitive varibales to conservative variables.

Initialization of the Navier-Stokes equations.

VK! density for viscous terms VKbc_state_gradient(1,1) = bc_state(1) VKkvp%key = trim(bc_state_gradient(1,1)%state_name) VKcall append( me = bc_varDict_gradient, val = kvp )

Convert conservative to primitive variables.

Convert primitive varibales to conservative variables.

Convert conservative to primitive variables.

VK allocate(bc_normal_vec_gradient(2)) VK allocate(bc_trafo_gradient(2))

Convert primitive varibales to conservative variables.

Routine which updates the background since it is a temporal function and

Evaluate the state on the edge using the exact riemann solver

Compute the exact solution of the Riemann problem for the Euler equation.

Set the general constants for the exact riemann solver.

Evaluate the solution to the 1D Riemann problem for a given sample point s.

evaluate "currentDensity" source

evaluate "currentDensity" source

Evaluates the material properties for all elements contained in the computeElems variable of the material_desc datatype.

Subroutine to evaluate the material properties on the nodes of the compute faces (left and right element's trace).

Lift face data to volume data.

Lift face data to volume data.

Lift face data to volume data.

Create source elements list for given source variable

Routine to obtain a modal representation for a variable, which is only available in nodal form, like space-time functions.

This routine prepares the data for variable derivation or operators. It gathers all input variables from the variable system, oversamples and projects them into nodal space, calls the function with the actual calculation and projects the results back into modal space. As these projections are common to all elementwise variable accesses, this generic routine does all necessary operations in a generic way.

Routine to obtain a modal representation for a variable, which is only available in nodal form, like space-time functions.

This function returns local modular variable atl_elemTimers to apesmate

This function returns local modular variable atl_elemTimers to apesmate

Routine to get a pointer to a new instance of atl_varSys_data_type to be used as method data for a variable in the variable system.

Subroutine to count the number of elements per level.

Get the number of (virtual) boundary elements for each level and each direction.

Return the penalty parameter for the IP discretizations of higher order equations.

Get all the surface points for a specific boundary.

This function returns local modular variable atl_timerHandles to apesmate

Subrountine which gather all calls to get the timestep for the current iteration

Getting the linearization indicator for Euler equations from the config.

This routine implements the getElement interface for material variables.

This routine implements the getPoint interface for material variables.

Determines maximum propagation speed, i.e. the speed of light depends only on material parameters.

Determines maximum propagation speed for Maxwell equation with divergence cleaning (hyperbolic), i.e. the speed of light depends only on material parameters.

boundary timer should only be measured for boundary elements VK setBnd = tem_getTimerVal( timerHandle = atl_timerHandles%setBnd) VK readBC = tem_getTimerVal( timerHandle = atl_timerHandles%readBC)

Routine to load the timestepping scheme from the config.

Godunov flux for the Euler equation.

Godunov flux for the 1D Euler equation.

Godunov flux for the 2D Euler equation.

Calculate the HLL flux given the left an right state.

Calculate the 1D HLL flux given the left an right state.

Calculate the 2D HLL flux given the left an right state.

Subroutine to load configuration file options for the conservative positivity preserving limiter.

Subroutine to load configuration file options for the covolume filter method.

Subroutine to load configuration file options for the spectral viscosity method.

Subroutine to load configuration file options for the spectral viscosity method.

Init source terms for flow simulations.

\brief init the variables for acoustic equation

Init source terms for flow simulations.

init the variables for acoustic equation

Init the variable system for simulations of advection equation.

Subroutine to create the levelwise list of boundaries.

Initialize the parallel module to make it usable in ATELES.

Method to initialize a cube mesh by tree and boundary definitions obtained by treelm.

Initialize the cubic elements.

Creates boundary informations for the faces.

Initialize the container module.

This subroutine reads the equation system to solve from the configuration.

This routine initializes possible source variables and returns the filled up list of the poss_srcVars

Init the variable system for Euler (inviscid) flow simulations.

This routine initializes possible source variables and returns the filled up list of the poss_srcVars

Init the variable system for Euler (inviscid) flow simulations.

Adds the properties of the expected source terms to the list of possible variables to extract these expected variables later on from the configuration file.

Init source terms for flow simulations. This routine initializes possible source variables and returns the filled up list of the poss_srcVars

Init the variable system for Euler (inviscid) flow simulations.

Routine to initialize explicit runge kutta scheme for timestepping.

Initialize explicit euler scheme for timestepping.

Routine to initialize explicit runge kutta scheme for timestepping.

Routine to initialize explicit runge kutta scheme for timestepping.

Routine to initialize explicit runge kutta taylor scheme for timestepping.

Creates boundary informations for the faces.

Initializes the face data by a given set of faces (asymmetrically).

Initializes the face data by a given set of faces.

Initializes the face data by a given set of faces.

initializes the face buffers for communication.

Init source terms for flow simulations. This routine initializes possible source variables and returns the filled up list of the poss_srcVars

Routine to init the timestepping scheme.

Init the variable system for simulations of Heat equation.

Init the variable system for simulations of Heat equation.

Init the variable system for simulations of Heat equation.

Routine to initialize IMEX Runge-Kutta scheme for timestepping.

init routine for the kerneldata type.

init the variables for LinearEuler equation

init the variables for LinearEuler equation

Init source terms for flow simulations. This routine initializes possible source variables and returns the filled up list of the poss_srcVars

Read the configuration for the material paramters for Maxwell equations from configuration files and init the material parameter datatype.

init source terms for electrodynamic simulations.

init the variables for maxwell equation (2D, TE-mode formulation).

Adds the properties of the expected source terms to the list of possible variables to extract these expected variables later on from the configuration file.

init source terms for electrodynamic simulations.

init the variables for maxwell equation

Adds the properties of the expected source terms to the list of possible variables to extract these expected variables later on from the configuration file.

init source terms for electrodynamic simulations with divergence

init the variables for maxwell equation

Initiate the MODG kernel for cubic elements on all levels.

Initiate the MODG kernel for cubic elements on all levels.

Initiate the MODG kernel for cubic elements on all levels.

init the variables for nernst-planck equation

initialize the parallel module to make it usable in ATELES

Routine to init container for penalization data.

Read the info for the physical checks from the configuration file

Read from preCICE the requested points

Init the variable system for filtered NavierStokes equation.

subroutine to intialize a scheme as specified by a given lua script file.

Initialize the space basis, this subroutine has to be called before using the module variable space_basis#basis.

Initialize the statedata.

Routine to initialize the complete Ateles solver.

This routine is used to initialize an array in an OpenMP PARALLEL region. Usually this is done using a WORKSHARE directive, but due to a bug in Intel 15 we cannot make use of WORKSHARE.

Initializes writing restart files, if activated.

Subroutine to initialize the timestep information for the first iteration

Recursive interpolation of element states among the levels of the mesh.

Find the maximal deviation for the polynomials representing the state in each element.

Find the maximal gradient estimation for the polynomials representing the state in each element.

Lax-Friedrich flux (in fully conservative variables) for the Acoustic equation

Lax-Friedrich flux (in fully conservative variables) for the Euler equation

Lax-Friedrich flux (in fully conservative variables) for the 1D Euler equation

Lax-Friedrich flux (in fully conservative variables) for the 2D Euler equation

Lax-Friedrich flux (in fully conservative variables) for the

Lax-Friedrich flux (in fully conservative variables) for the

Lax-Friedrich flux (in fully conservative variables) for the Euler equation

Load the definition of a Legendre polynomial variable from a Lua script.

Function for physical flux of the LinearEuler equation F, 1D? Since it is 1d, there need to be passed the correct background velocity (u0 for F - flux in x direction, v0 for G - flux in y direction, w0 for H - flux in z direction)

Function for physical flux of the LinearEuler equation F, 1D? Since it is 1d, there need to be passed the correct background velocity (u0 for F - flux in x direction, v0 for G - flux in y direction, w0 for H - flux in z direction)

Load the configuration for acoustic equations from the Lua script.

subroutine to initialize an equation of type advection equation as defined in the configuration file

subroutine to intialize BBM

Get the boundary configuration.

Load the boundary condition for state variables.

Subroutine to initialize an equation of type euler equation as defined in the configuration file

subroutine to initialize an equation of type filtered-navier-stokes equation (turbulence modelling) as defined in the configuration file

subroutine to initialize an equation of type heat equation as defined in the configuration file

subroutine to load the initial conditions from a lua configuration file.

subroutine to initialize an equation of type linear euler equation as defined in the configuration file

Subroutine to initialize an equation of type navier stokes equation as defined in the configuration file

subroutine to intialize Nernst-Planck equation with constant

Routine to initialize the global parameters, sets the solver module variable.

Interface for fluxes of pure Maxwell equations.

Interface for fluxes of pure Maxwell equations.

Interface for fluxes of Maxwell equations with hyperbolic divergence cleaning.

Calculate the numerical flux for the advection equation and MODG scheme

Calculate the physical flux for the MODG scheme and advection equation.

Subroutine to create the modal representation for a ceratin boundary face.

\brief Interpolate modal face representation from coarse to next finer faces (level difference between coarser and finer faces has to be 1).

Lift the value on the face to a positive value if necessary.

Calculate the numerical flux for Euler equation and MODG scheme

Calculate the physical flux for the MODG scheme and Euler equation with constant characteristic (mask function) in the element.

Calculate the physical flux for the MODG scheme and Euler equation with variable characteristic (mask function) in the element.

Interpolate modal face representation from next finer faces to coarse level (level

Calculate the numerical flux for the Heat equation and MODG scheme

Calculate the physical flux for the MODG scheme and Heat equation.

Applies the inverse of the mass matrix for a 3D scheme.

Calculate the physical flux for the MODG scheme and Linearized euler equation.

Projects modal representation of each cell to its faces, i.e. this subroutine creates a modal representation on the faces.

Projects modal representation of each cell to its faces, i.e. this subroutine creates a modal representation on the faces.

Subroutine to project modal representations of physical flux, numerical flux and source terms onto test functions.

Subroutine to project modal representations of physical flux, numerical flux and source terms onto test functions.

Subroutine to set face values to impose boundary conditions.

Project modal representation of an element to one of its faces for Q space.

Calculate the numerical flux for acoustic equation and MODG scheme

Calculate the physical flux for the MODG scheme and Acoustic equation.

Subroutine to create the modal representation for a ceratin boundary face.

Interpolate modal face representation from coarse to next finer faces (level

Lift the mean on the face to a positive value if necessary.

Calculate the numerical flux for Euler equation and MODG scheme

Numerical flux calculation for Euler equation across the faces in a single spatial direction (with constant penalization parameters).

Numerical flux calculation for Euler equation across the faces in a single spatial direction (with non-constant penalization parameters).

Calculate the physical flux for the MODG scheme and Euler equation.

Calculate the numerical flux for Navier-Stokes equation and MODG scheme

Calculate the physical flux for the MODG scheme and Navier-Stokes equation (with constant penalizations).

Calculate the physical flux for the MODG scheme and Navier-Stokes equation (with non-constant penalizations).

Interpolate modal face representation from next finer faces to coarse level (level

Calculate the numerical flux for Heat equation and MODG scheme

Calculate the physical flux for the MODG scheme and Heat equation.

Applies the inverse of the mass matrix for a 2D scheme.

Calculate the numerical flux for LinearEuler equation and MODG scheme

Calculate the physical flux for the MODG scheme and Linearized euler equation.

Calculate the physical flux for the MODG scheme and Linearized euler equation.

Projects modal representation of each cell to its faces, i.e. this subroutine creates a modal representation on the faces.

Calculate the numerical flux for Navier-Stokes equation and MODG scheme

Calculate the physical flux for the MODG scheme and Navier-Stokes equation (with constant penalizations).

Calculate the physical flux for the MODG scheme and Navier-Stokes equation (with non-constant penalizations).

Subroutine to project modal representations of numerical flux and source terms onto test functions.

Subroutine to project modal representations of physical flux, numerical flux and source terms onto test functions.

Projection of the source terms (in modal representation) to the test functions.

Subroutine to set face values to impose boundary conditions at a certain point of the domain. The subroutine is operating levelwise.

Project modal representation of gradients of an element to one of its faces for Q space.

Project modal representation of an element to one of its faces for P space.

Project modal representation of an element to one of its faces for Q space.

Calculate the numerical flux for acoustic equation and MODG scheme

Calculate the physical flux for the MODG scheme and Acoustic equation.

Subroutine to create the modal representation for a ceratin boundary face.

Project coarse parent element to its 8 finer child elements by a simple L2 projection.

Project coarse parent element to its 8 finer child elements by a simple L2 projection.

Project coarse parent element to its 8 finer child elements by a simple L2 projection.

Interpolate modal face representation from coarse to next finer faces (level

Lift the mean on the face to a positive value if necessary.

Calculate the numerical flux for Euler equation and MODG scheme

Numerical flux calculation for Euler equation across the faces in a single spatial direction.

Numerical flux calculation for Euler equation across the faces in a single spatial direction.

Calculate the physical flux for the MODG scheme and Euler equation.

Calculate the physical flux for the MODG scheme and Euler equation.

Calculate the numerical flux for Navier-Stokes equation and MODG scheme

Calculate the physical flux for the MODG scheme and Navier-Stokes equation (with constant penalizations).

Calculate the physical flux for the MODG scheme and Navier-Stokes equation (with non-constant penalizations).

Project data from 8 smaller elements to its parent element in terms of L2 projections.

Project data from 8 smaller elements to its parent element in terms of L2 projections.

Project data from 8 smaller elements to its parent element in terms of L2 projections.

Interpolate modal face representation from next finer faces to coarse level (level

Calculate the numerical flux for Heat equation and MODG scheme

Numerical flux calculation for Heat equation across the faces in a single spatial direction.

Calculate the physical flux for the MODG scheme and Heat equation.

Applies the inverse of the mass matrix for a 3D scheme.

Subroutine to test the various routines of this module.

Calculate the numerical flux for LinearEuler equation and MODG scheme

Calculate the physical flux for the MODG scheme and Linearized euler equation.

Calculate the physical flux for the MODG scheme and Linearized euler equation.

Calculate the numerical flux for Maxwell equation and MODG scheme

Calculate the physical flux for the MODG scheme and Maxwell equation. used for both space P and Q

Calculate the physical flux for the MODG scheme and Maxwell equation. used for both space P and Q

Calculate the numerical flux for Maxwell equation and MODG scheme

Calculate the numerical flux for Maxwell equation with hyperbolic divergence cleaning and MODG scheme

Calculate the physical flux for the MODG scheme and Maxwell equation with hyperbolic divergenc cleaning.

Calculate the physical flux for the MODG scheme and Maxwell equation with hyperbolic divergenc cleaning.

Projects modal representation of each cell to its faces, i.e. this subroutine creates a modal representation on the faces.

Calculate the numerical flux for Navier-Stokes equation and MODG scheme

Calculate the physical flux for the MODG scheme and Navier-Stokes equation (with constant penalizations).

Calculate the physical flux for the MODG scheme and Navier-Stokes equation (with non-constant penalizations).

Subroutine to project modal representations of physical flux, numerical flux and source terms onto test functions.

Subroutine to project modal representations of physical flux, numerical flux and source terms onto test functions.

Projection of the source terms (in modal representation) to the test functions.

Projection of the source terms (in modal representation) to the test functions.

Projection of the source terms (in modal representation) to the test functions.

Applies a scaled transposed inverse of the mass matrix for a 3D scheme. The result is a transformation of the polynomial basis from the ansatz- to the test-polynomials. This is useful for a local predictor

Applies a scaled transposed inverse of the mass matrix for a 3D scheme. The result is a transformation of the polynomial basis from the ansatz- to the test-polynomials. This is useful for a local predictor

Projection of the physical flux for the local predictor: This function expects the physical flux transformed in the basis of test-polynomials and projects it onto the ansatz-functions. Thus, the result is directly given in the basis of ansatz-polynomials (no invMassMatrix-call needed afterwards)

Projection of the physical flux for the local predictor: This function expects the physical flux transformed in the basis of test-polynomials and projects it onto the ansatz-functions. Thus, the result is directly given in the basis of ansatz-polynomials (no invMassMatrix-call needed afterwards)

Subroutine to set face values to impose boundary conditions at a certain point of the domain. The subroutine is operating levelwise.

Numerical flux calculation for stab-viscous part of the Navier-Stokes equation across the faces in a single spatial direction.

Numerical flux calculation for viscous part of the Navier-Stokes equation across the faces in a single spatial direction.

Project modal representation of gradients of an element to one of its faces for Q space.

Project modal representation of an element to one of its faces for P space.

Project modal representation of an element to one of its faces for Q space.

For a 2D testcase just the edges are of importance, in 3D also triangles have to be provided loop over all faces and edges. Set an offset, to get from one face to the other For a 3D testcase in addition to the edges triangles have to be created

This routine takes in a variable and differentiates it in a modal way

to do

This routine takes in a variable and differentiates it in a modal way

This routine takes in a variable and differentiates it in a modal way

This routine takes in a variable and differentiates it in a modal way

This routine takes points coordinates, pass them tp the input variables the opertaion depends and return indices

Routine to check if the physical values of a state are physically meaningful or not for the acoustic equation.

Routine to check if the physical values of a state are physically meaningful or not for the acoustic 2d equation.

Routine to check if the physical values of a state are physically meaningful or not for the Euler equation.

Routine to check if the physical values of a state are physically meaningful or not for the Euler equation.

Routine to check if the physical values of a state are physically meaningful or not for the Euler equation.

Routine to check if the physical values of a state are physically meaningful or not for the linear euler equation.

Routine to check if the physical values of a state are physically meaningful

Routine to check if the physical values of a state are physically meaningful or not for the Maxwell equation.

Routine to check if the physical values of a state are

Physical flux calculation along x direction for Euler equation.

Physical flux calculation along x direction for Euler equation.

Physical flux calculation along x direction for Euler equation.

Physical flux calculation along x direction for Euler equation.

Function for physical flux of the Maxwell equations in terms of D and B.

Function for physical flux of the Maxwell equations in terms of D and B.

Function for physical flux of the Maxwell equations in terms of D and B.

Physical flux calculation along x direction for the

Apply pointwise limitation of denisty and energy to ensure positivity for density and pressure.

Apply pointwise limitation of denisty and energy to ensure positivity for density and pressure.

Subroutine to execute the preprocessing for the MODG kernels. Currently this includes: Convert external source terms to modal representation.

Subroutine to execute the preprocessing for the MODG kernels. Currently this includes: Convert external source terms to modal representation.

interface for preprocessing the data for the kernel

call varSys%method%val(fun%inpuit_varPos(1))%get_point( & & varSys = varSys, & & point = point, & & time = time, & & tree = tree, & & nPnts = nPnts, & & res = density ) call varSys%method%val(fun%input_varPos(2))%get_point( & & varSys = varSys, & & point = point, & & time = time, & & tree = tree, & & nPnts = nPnts, & & res = momentum ) call varSys%method%val(fun%input_varPos(3))%get_point( & & varSys = varSys, & & point = point, & & time = time, & & tree = tree, & & nPnts = nPnts, & & res = energy )

Project two elements onto single co-volume element.

Project two elements onto single co-volume element.

Project two elements onto single co-volume element.

loop over linked list of spacetimefunction set the current to head loop over all shapes of that spacetimefunction check if there is a precice kind get the number of unquie points ( xyz should be same number)

Read the serialized restart file into the state vectors

Reassigns the spacetime function pointers based on the material variable position and the spacetime function position.

Transform reference points to physical points in the element of the tree identified by the provided elempos.

Subroutine to move points defined on the reference element [-1,+1] to the physical element coordinates.

dump the equation variables into the lua file

Routine to set boundary values for the covolume filter.

This routine sets elementTimers passed by apesmate

This routine sets elementTimers passed by apesmate

This routine sets atl_timerHandles passed by apesmate

NE ! Dump and reset the timer after each Iteration NE call atl_dumpTimers( general = params%general, & NE & nElems = tree%global%nElems, & NE & nDofs = element_container%cubes & NE & %scheme_list(tree%global%minLevel) & NE & %nDofs, & NE & nVars = equation%varSys%nScalars ) NE call atl_resetTimers()

This routine converts primitive variables in source terms to convervative

Damp )the modal coefficients of the state vector by a given spectral viscosity method.

Damp )the modal coefficients of the state vector by a given spectral viscosity method.

Damp )the modal coefficients of the state vector by a given spectral viscosity method.

Subroutine to apply the stabilization procedure to the state vector.

Routine to store position of user variable defined state and state_gradient boundary variable in bc(iBC)%state(iVar)%varPos and bc(iBC)%state_gradient

Get the subresolution data for all elements for a given color and in the requested format.

Subroutine to load subresolution information for a given tree.

Read the configuration for the material paramters for Maxwell equations from configuration files and init the material parameter datatype.

subroutine to calculate the RHS of the PDE from the sum of all

To obtain values of a given variable, it is necessary to state the treeID and time at which the variable should be evaluated. The interface is nDofs values to cover the all degrees of freedoms in the element. Of course the variable system itself also needs to be passed in, to allow the computation of other derived quantities as needed. The method description itself is passed in automatically, and has not to be provided explicitly.

Interface description for a variable access method (single point).

Routine for gettint the actual value for a given array of indices. The indices belong to the grwarray of points storing levelwise in Pointdata%pntLvl(iLevel). Hence this routines takes the indeices as input, can refer to the pointData and evaluate the variable and returns the values

Method to load user defined variables from the configuration.

This routine takes points coordinates, stores them in the method_data and return indices where points are located in the growing array of points or values ( sometimes we do not need to store the points ) It is need to setup points for every variable. Points will be provided by boundaries or sources depends on what uses the variable. This points do not change with time . This indices will be stored in corresponding boundary or source to evaluate a value on that point later using tem_varSys_proc_getValOfIndex.

Physical flux calculation along x direction for Euler equation.

Physical flux calculation along x direction for Euler equation.

Physical flux calculation along x direction for Euler equation.

Project elemental state to a particular face.

Project elemental state to a particular face.

Project elemental state to a particular face.

ceate equidistant points for write to precice and use the number of these points for the setup_indices routine. If its set by the user, the Nearest-Pojection can be used for these points

setup_indices for boundary variables are read_var in precice space-time function we are not able to use 'idx' created here as input for get_valOfIndex for write_vars.Therefore we have to do setup_indices for write_vars once during initialization and then store the idx in precice_coupling type. This is done for the non-equidistant points. We also set the vertices and the edges as well as the triangles for the Interpolation between the domains.

loop over linked list of spacetimefunction set the current to head loop over all shapes of that spacetimefunction check if there is a precice kind

loop over linked list of spacetimefunction set the current to head check if there is a precice kind loop over all shapes of that spacetimefunction

Writes a restart file for the current point in time.

Writes a restart file, if necessary.

Write solver specific info to scratch file

We only want to have one component, thus we don't need to loop over them

We only want to have one component, thus we don't need to loop over them

We only want to have one component, thus we don't need to loop over them

We only want to have one component, thus we don't need to loop over them

We only want to have one component, thus we don't need to loop over them

We only want to have one component, thus we don't need to loop over them

Coefficient beta from the recursive formulation of Legendre polynomials, for the Legendre mode 'mode'.

Calculate barycentric coordinates of the tree ids given in treeids.

Function to find a single global time step for all levels and processes.

Calculate time step based on a given CFL condition for a cube in a acoustic simulation.

Calculate time step based on a given CFL condition for a cube in a acoustic simulation.

Calculate time step based on a given CFL number of cubes in a advection simulation.

Calculate time step based on a given CFL number of cubes in a electrodynamic simulation.

Calculate time step based on a given CFL condition for a cube in a flow simulation.

Calculate time step based on a given CFL condition for a cube in a flow simulation.

Calculate time step based on a given CFL condition for a cube in a flow simulation.

Calculate time step based on a given CFL condition for a cube in a flow simulation.

Calculate time step based on a given CFL condition for a cube in a linear Euler 2d simulation.

Calculate time step based on a given CFL condition for a cube in a linear Euler simulation.

Calculate time step based on a given CFL number of cubes in a Nernst-Planck simulation.

Calculate time step based on a given CFL condition for a cube in a viscous flow simulation.

Calculate time step based on a given CFL condition for a cube in a viscous flow simulation.

Calculate the timestep for a whole part of a cubic mesh by a CFL condition.

Compute the physical flux in x direction.

Compute the physical flux in x direction.

Compute the physical flux in x direction.

Compute the physical flux in x direction.

Compute the physical flux in x direction. For other directions a properly defined variable permutation can be used. This routine covers only constant material parameters.

Compute the physical flux in x direction. For other directions a properly defined variable permutation can be used. This routine covers non-constant material parameters.

compute the right hand side of your discrete equation.

Computes the right hand side for cubical elements and MODG scheme.

Computes the right hand side for cubical elements and 1D MODG scheme.

Computes the right hand side for cubical elements and 2D MODG scheme.

precompute the scalar products of the anstaz and test function

Copy the FPT header information.

\brief subroutine to create a single global timestep.

destroy the dynamic array

destruction of a dynamic array

Interface definition for elementwise timestepping routine.

Elemental operation for timestepping IMEX-RK.

Elemental operation for timestepping of order 4.

Elemental operation for timestepping of order 4.

Elemental operation for timestepping of type Runge Kutta Taylor

Elemental operation for timestepping of order 2.

Interface definition for elementwise timestepping routine.

Elemental operation for timestepping of order 4.

Elemental operation for timestepping of order 4.

Elemental operation for timestepping of order 4.

Elemental operation for timestepping of type Runge Kutta Taylor

Elemental operation for timestepping of order 2.

empty the array, reset nvals to be 0

empty all contents of the array without changing the size or status of any array

Subroutine to load boundary conditions for 2D Maxwell equations from a Lua configuration file.

evaluate "currentDensity" source

evaluate "currentDensity" source

evaluate "charge" source

evaluate "currentDensity" source

evaluate "pml" source (uniaxial PML)

increase the size of the container for the array.

expanding the dynamic array

current state

Get sampled data.

Get sampled data.

Subroutine calculates a substep of corrector, this is the same as a usual substep of the RKDG

Subroutine calculates a substep of the IMEX Runge-Kutta timestepping scheme.

Subroutine calculates the final update step of the Runge-Kutta method. It is performing levelwise.

and the code from here should move into its hlp module.

initialize the dynamic array

initialization of a dynamic array

Initialize the kernel states for all parts of the mesh.

routine to init the timestepping scheme.

Deactivate adaptive inviscous computations.

Estimate the impact of viscous terms in 2D.

Estimate the impact of viscous terms in 3D.

This function provides the test for equality of two projections.

This function provides the test for equality of two projections.

This function provides the test for equality of two nodes descriptions.

This function provides the test for equality of two projections.

This function provides the test for equality of the header for two projections.

This function provides the test for equality of two projections.

This function provides a > comparison of two projections.

This function provides a > comparison of two projections.

This function provides a > comparison of nodes descriptions.

This function provides a > comparison of two projections.

This function provides a > comparison of the header of two projections.

This function provides a > comparison of two projections.

This function provides a >= comparison of two projections.

This function provides a >= comparison of two projections.

This function provides a >= comparison of two nodes descriptions.

This function provides a >= comparison of two projections.

This function provides a >= comparison of the header of two projections.

This function provides a >= comparison of two projections.

This function provides a < comparison of two projections.

This function provides a < comparison of two projections.

This function provides a < comparison of two nodes descriptions.

This function provides a < comparison of two projections.

This function provides a < comparison of the header of two projections.

This function provides a < comparison of two projections.

This function provides a <= comparison of two projections.

This function provides a <= comparison of two projections.

This function provides a <= comparison of two nodes descriptions.

This function provides a <= comparison of two projections.

This function provides a <= comparison of the header of two projections.

This function provides a <= comparison of two projections.

This function provides the test for unequality of two projections.

This function provides the test for unequality of two projections.

This function provides the test for unequality of two nodes descriptions.

This function provides the test for unequality of two projections.

This function provides the test for unequality of the header of two projections.

This function provides the test for unequality of two projections.

An indicator that completely deactivates linearization.

An indicator to decide whether linearization of fluxes is tolerable based on the density.

An indicator to decide whether linearization of fluxes is tolerable based on the energy in 2D.

An indicator to decide whether linearization of fluxes is tolerable based on the energy in 2D.

An indicator to decide whether linearization of fluxes is tolerable based on the energy in 3D.

An indicator to decide whether linearization of fluxes is tolerable based on the error estimate.

An indicator to decide whether linearization of fluxes is tolerable based on the error estimate.

An indicator to decide whether linearization of fluxes is tolerable based on the error estimate.

Subroutine calculates a substep of the local predictor

Subroutine to calculate the flux for pure Maxwell equations without any divergence cleaning on the reference cubic face.

Subroutine to calculate the flux for pure Maxwell equations without

calculate flux of pure maxwell equation directly on the

calculate flux of pure maxwell equation directly on the face-vector

calculate flux of pure maxwell equation directly on the

calculate flux of pure maxwell equation directly on the face-vector

Subroutine to calculate the flux for pure Maxwell equations with

calculate flux of maxwell equation with hyperbolic divergence

calculate flux of maxwell equation with hyperbolic divergence

subroutine for timestepping with explicit euler.

Subroutine for timestepping with IMEX Runge-Kutta.

Subroutine for timestepping with explicit runge kutta of order 4.

Subroutine for timestepping with explicit runge kutta of order 4. It contains four Euler forwarding substeps: 0.1 temp = u0 1.1 u1 = rk4_substep(temp) = rhs(temp) 1.2 temp = u0 + dt/2 * u1 2.1 u2 = rk4_substep(temp) = rhs(temp) 2.2 temp = u0 + dt/2 * u2 3.1 u3 = rk4_substep(temp) = rhs(temp) 3.2 temp = u0 + dt * u3 4.1 u4 = rk4_substep(temp) = rhs(temp) 4.2 u0 = u0 + dt/6 * ( u1 + 2*(u2+u3) + u4 ) After each substep, the results are stabilized

Subroutine for timestepping with explicit runge kutta taylor

Subroutine for timestepping with explicit runge kutta of order 2.

Numerical flux calculation for the advection equation across the faces in a single spatial direction.

Function to extrapolate face values for a given boundary condition in physical or modal space.

Function to mirror a modal representation around a given boundary condition in modal space.

Function to mirror pointvalues for a given boundary conditions.

Numerical flux calculation for Euler equation across the faces in a single spatial direction.

Numerical flux calculation for the Heat equation across the faces in a single spatial direction.

Applies the inverse of the mass matrix for a 3D scheme.

Set boundary values in a modal way

Set boundary values in a nodal way

Projection of the numerical flux in x direction onto the differentiated

Projection of the numerical flux in x direction onto the testfunctions.

Projection of the penalization terms (in modal representation) to the test functions.

Projection of the physical flux in x direction onto the testfunctions.

Projection of the source terms (in modal representation) to the test functions.

Subroutine to project modal representations of physical flux, numerical flux and source terms onto test functions.

Projects derivative of modal representation of each cell to its faces, i.e. this subroutine creates a modal representation on the faces. Project modal representation of an element to one of its faces for Q space.

Function to extrapolate face values for a given boundary condition in physical or modal space.

Function to mirror pointvalues for a given boundary conditions.

This subroutine computes the physical fluxes for various equation system

Applies the inverse of the mass matrix for a 2D scheme in P_space.

Applies the inverse of the mass matrix for a 2D scheme in Q_space.

Set boundary values in a modal way

Set boundary values in a nodal way

Projection of the numerical flux in x direction onto the testfunctions.

Projection of the numerical flux in x direction onto the testfunctions for Q_space.

Projection of the numerical flux in y direction onto the testfunctions.

Projection of the numerical flux in y direction onto the testfunctions for Q_space.

Projection of the penalization terms (in modal representation) to the test functions.

Projection of the physical flux in x direction onto the testfunctions.

Projection of the physical flux in x direction onto the testfunctions.

Projection of the physical flux in y direction onto the testfunctions.

Projection of the physical flux in y direction onto the testfunctions.

Projection of the source terms (in modal representation) to the test functions.

Projection of the source terms (in modal representation) to the test functions.

Projection of the numerical flux in x direction onto the testfunctions.

Projection of the numerical flux in y direction onto the testfunctions.

Numerical flux calculation for Rans 2D equation across the faces in a single spatial direction (with constant penalization parameters).

Numerical flux calculation for Rans equation across the faces in a single spatial direction (with non-constant penalization parameters).

Project modal representation of a semi-refined face (i.e. on a face

Project modal representation on a face to a semi-refined face (i.e. on a face

Numerical flux calculation for stab-viscous part of the Navier-Stokes equation across the faces in a single spatial direction.

Numerical flux calculation for stab-viscous part of the RANS equation across the faces in a single spatial direction.

Numerical flux calculation for viscous part of the Navier-Stokes equation across the faces in a single spatial direction.

Numerical flux calculation for viscous part of the RANS equation across the faces in a single spatial direction.

Function to extrapolate face values for a given boundary condition in physical or modal space.

Function to mirror pointvalues for a given boundary conditions.

TODO NA - Move this routine to the atl_modg_kernel_module

Applies the inverse of the mass matrix for a 3D scheme.

Applies the inverse of the mass matrix for a 3D scheme.

Applies the inverse of the mass matrix for a 3D scheme.

Set boundary values in a modal way

Calculate the projection of the physical flux for the MODG scheme and Nernst-Planck equation.

Numerical flux calculation for Nernst-Planck equation across the faces in X direction.

Numerical flux calculation for Nernst-Planck equation across the faces in Y direction.

Numerical flux calculation for Nernst-Planck equation across the faces in Z direction.

Set boundary values in a nodal way

X direction for 5 scalars projection of the physical flux onto the testfunctions, with unrolled loops

X direction for 6 scalars projection of the physical flux onto the testfunctions, with unrolled loops

Y direction for 5 scalars projection of the physical flux onto the testfunctions, with unrolled loops

Y direction for 6 scalars projection of the physical flux onto the testfunctions, with unrolled loops

Z direction for 5 scalars projection of the physical flux onto the testfunctions, with unrolled loops

Z direction for 6 scalars projection of the physical flux onto the testfunctions, with unrolled loops

Projection of the numerical flux onto the differentiated testfunctions.

Projection of the numerical flux onto the testfunctions for P_space.

Projection of the numerical flux onto the testfunctions for Q_space.

Projection of the penalization terms (in modal representation) to the test functions.

Projection of the physical flux onto the testfunctions, with unrolled loops => fewer loop-overhead/instructions, but more "random" memory accesses MZ: perhaps this version is faster for low order (or always, depending on the machine?)

Projection of the physical flux onto the testfunctions, with unrolled loops

Projection of the numerical flux in x direction onto the testfunctions.

Projection of the numerical flux in y direction onto the testfunctions.

Projection of the numerical flux in y direction onto the testfunctions.

Subroutine to semi-coarsen an element with modal polynomial representation to its semi-parent.

Subroutine to semi-coarsen an element with modal polynomial representation to its semi-parent.

Subroutine to semi-coarsen an element with modal polynomial representation to its semi-parent.

Project modal representation of a semi-refined face (i.e. on a face

Subroutine to semi-refine an element with modal polynomial representation into its semi-children.

Subroutine to semi-refine an element with modal polynomial representation into its semi-children.

Subroutine to semi-refine an element with modal polynomial representation into its semi-children.

Project modal representation on a face to a semi-refined face (i.e. on a face

Project modal representation of an element to its two faces in X.

Project modal representation of an element to its two faces in Y.

Project modal representation of an element to its two faces in Y.

Subroutine to calculate the numerical flux for the first equation (u)

Subroutine to calculate the flux for Nernst-Planck equations.

Subroutine to calculate the numerical flux for the first equation (u)

Subroutine to calculate the flux for Nernst-Planck equations.

Compute one-sided iteration

Initial setup for the iterative computation of the solution to the Riemann problem.

Outputs the result of appending a variable. If appending was not successful, the program is aborted.

Coefficients from the recursive formulation of legendre polynomials. L_n = alpha * x * L_n-1 + beta * L_n-2

Prodcut of alpha(numerator) * beta(denominator) / alpha(denominator)

Quotient of two alpha values.

Coefficients from the recursive formulation of legendre polynomials. L_n = alpha * x * L_n-1 + beta * L_n-2

do IDeg1 = 1, mPd+1 do IDeg2 = 1, mPd=1 !! iDeg2 = mod(iDeg-1,mpd+1)+1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Compute the derivative in X direction for 3D Legendre polynomial.

Compute the derivative in Y direction for 3D Legendre polynomial.

Compute the derivative in Y direction for 3D Legendre polynomial.

Subroutine to change the polynomial space (Q or P) of an atl_statedata_type from Q-space to P-space and vice versa.

Copy a single element state into a larger array and pad it with zeroes.

Copy a single 1D element state into a larger array and pad it with zeroes.

Copy a single 2D element state into a larger array and pad it with zeroes.

Copy a single element state into a larger array and pad it with zeroes.

Truncating an oversampled polynomial representation back to the original representation.

Truncating an oversampled 1D polynomial representation back to the original representation.

Truncating an oversampled 2D polynomial representation back to the original representation.

Truncating an oversampled polynomial representation back to the original representation.

Subroutine to convert linearized dof index to ansatz function number for Q-Polynomials.

Evaluate three-dimensional tensor product Legendre polynomials (not-normalized) at a given set of coordinates.

Returns the value of the non-normalized Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the non-normalized Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the non-normalized differentiated Legendre polynomial at the leftboundary of the reference element, i.e. at -1.

Returns the value of the gradient of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the gradient of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.Vectorized version.

Returns the value of the derivaitve of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the derivaitve of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.Vectorized version.

Returns the value of the non-normalized differentiated Legendre polynomial at the right boundary of the reference element, i.e. at +1.

Returns the value of the gradient of dual Legendre polynomial at the right boundary of the reference element, i.e. at +1.

Returns the value of the gradient of dual Legendre polynomial at the right boundary of the reference element, i.e. at +1. Vectorized version.