escript  Revision_
FluxLimiter.h
Go to the documentation of this file.
1 
2 /*****************************************************************************
3 *
4 * Copyright (c) 2003-2020 by The University of Queensland
5 * http://www.uq.edu.au
6 *
7 * Primary Business: Queensland, Australia
8 * Licensed under the Apache License, version 2.0
9 * http://www.apache.org/licenses/LICENSE-2.0
10 *
11 * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12 * Development 2012-2013 by School of Earth Sciences
13 * Development from 2014-2017 by Centre for Geoscience Computing (GeoComp)
14 * Development from 2019 by School of Earth and Environmental Sciences
15 **
16 *****************************************************************************/
17 
18 #ifndef __PASO_FLUXLIMITER_H__
19 #define __PASO_FLUXLIMITER_H__
20 
21 #include "Transport.h"
22 
23 namespace paso {
24 
25 
27 struct FCT_FluxLimiter
28 {
31 
32  inline dim_t getTotalNumRows() const
33  {
34  return antidiffusive_fluxes->getTotalNumRows();
35  }
36 
38  {
39  return antidiffusive_fluxes->pattern;
40  }
41 
42  void setU_tilde(const double* Mu_tilde);
44  void addLimitedFluxes_Complete(double* b);
45 
48  double dt;
49  double* u_tilde;
50  double* MQ; // (M_C* Q_min, M_C* Q_max)
51  double* R; // (R-, R+)
52  //Coupler_ptr MQ_coupler;
55  double* borrowed_lumped_mass_matrix; // borrowed reference
56 };
57 
58 } // namespace paso
59 
60 #endif // __PASO_FLUXLIMITER_H__
61 
paso::LARGE_POSITIVE_FLOAT
static const real_t LARGE_POSITIVE_FLOAT
Definition: FCT_Solver.cpp:44
paso::FCT_FluxLimiter::FCT_FluxLimiter
FCT_FluxLimiter(const_TransportProblem_ptr tp)
Definition: FluxLimiter.cpp:45
escript::DataTypes::real_t
double real_t
type of all real-valued scalars in escript
Definition: DataTypes.h:76
FluxLimiter.h
paso::FCT_FluxLimiter::addLimitedFluxes_Complete
void addLimitedFluxes_Complete(double *b)
Definition: FluxLimiter.cpp:220
paso::SystemMatrix
this class holds a (distributed) stiffness matrix
Definition: SystemMatrix.h:60
paso::FCT_FluxLimiter::setU_tilde
void setU_tilde(const double *Mu_tilde)
Definition: FluxLimiter.cpp:75
paso::FCT_FluxLimiter::u_tilde_coupler
Coupler_ptr< real_t > u_tilde_coupler
Definition: FluxLimiter.h:79
paso::SystemMatrix_ptr
boost::shared_ptr< SystemMatrix > SystemMatrix_ptr
Definition: SystemMatrix.h:53
paso::SystemMatrixPattern_ptr
boost::shared_ptr< SystemMatrixPattern > SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:52
escript::DataTypes::real_t_max
real_t real_t_max()
Returns the maximum finite value for the real_t type.
Definition: DataTypes.h:114
paso::const_SystemMatrixPattern_ptr
boost::shared_ptr< const SystemMatrixPattern > const_SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:54
paso::FCT_FluxLimiter::getTotalNumRows
dim_t getTotalNumRows() const
Definition: FluxLimiter.h:57
paso::FCT_FluxLimiter::mpi_info
escript::JMPI mpi_info
Definition: FluxLimiter.h:72
escript::DataTypes::dim_t
index_t dim_t
Definition: DataTypes.h:90
escript::JMPI
boost::shared_ptr< JMPI_ > JMPI
Definition: EsysMPI.h:71
paso::Coupler_ptr
boost::shared_ptr< Coupler< T > > Coupler_ptr
Definition: Coupler.h:55
paso::FCT_FluxLimiter::R
double * R
Definition: FluxLimiter.h:76
paso::FCT_FluxLimiter::borrowed_lumped_mass_matrix
double * borrowed_lumped_mass_matrix
Definition: FluxLimiter.h:80
paso::FCT_FluxLimiter::u_tilde
double * u_tilde
Definition: FluxLimiter.h:74
paso::FCT_FluxLimiter::R_coupler
Coupler_ptr< real_t > R_coupler
Definition: FluxLimiter.h:78
paso::FCT_FluxLimiter::antidiffusive_fluxes
SystemMatrix_ptr antidiffusive_fluxes
Definition: FluxLimiter.h:71
paso::const_SystemMatrix_ptr
boost::shared_ptr< const SystemMatrix > const_SystemMatrix_ptr
Definition: SystemMatrix.h:55
paso::const_TransportProblem_ptr
boost::shared_ptr< const TransportProblem > const_TransportProblem_ptr
Definition: Transport.h:33
Transport.h
paso::FCT_FluxLimiter::addLimitedFluxes_Start
void addLimitedFluxes_Start()
Definition: FluxLimiter.cpp:150
escript::DataTypes::index_t
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:85
paso::FCT_FluxLimiter::~FCT_FluxLimiter
~FCT_FluxLimiter()
Definition: FluxLimiter.cpp:66
PASO_DLL_API
#define PASO_DLL_API
Definition: Paso.h:59
paso::FCT_FluxLimiter::getFluxPattern
SystemMatrixPattern_ptr getFluxPattern() const
Definition: FluxLimiter.h:62
paso
Definition: BiCGStab.cpp:26
paso::FCT_FluxLimiter::dt
double dt
Definition: FluxLimiter.h:73
paso::FCT_FluxLimiter::MQ
double * MQ
Definition: FluxLimiter.h:75