advection.py

Go to the documentation of this file.

# This file is part of the ExaHyPE2 project. For conditions of distribution and
# use, please see the copyright notice at www.peano-framework.org
from .equation import Equation
 
 
class Advection(Equation):
    def __init__(self, dimensions, adv_speed=1.0):
        self.dimensionsdimensionsdimensions = dimensions
        self.num_unknownsnum_unknownsnum_unknowns = 2 if dimensions == 2 else 3
        self.num_auxiliary_variablesnum_auxiliary_variablesnum_auxiliary_variables = 0
        self.adv_speedadv_speed = adv_speed
        self.is_linearis_linearis_linear = True
 
    def eigenvalues(self):
        return (
            """
  constexpr double v = """
            + str(self.adv_speedadv_speed)
            + """;
  return v;
"""
        )
 
    def flux(self):
        return (
            """
    F[0] = 0.0;
    F[1] = 0.0;
#if DIMENSIONS == 3
    F[2] = 0.0;
#endif
 
    F[normal] = """
            + str(self.adv_speedadv_speed)
            + """ * Q[normal];
"""
        )