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 {
29  ~FCT_FluxLimiter();
30 
31  inline dim_t getTotalNumRows() const
32  {
33  return antidiffusive_fluxes->getTotalNumRows();
34  }
35 
37  {
38  return antidiffusive_fluxes->pattern;
39  }
40 
41  void setU_tilde(const double* Mu_tilde);
42  void addLimitedFluxes_Start();
43  void addLimitedFluxes_Complete(double* b);
44 
47  double dt;
48  double* u_tilde;
49  double* MQ; // (M_C* Q_min, M_C* Q_max)
50  double* R; // (R-, R+)
51  //Coupler_ptr MQ_coupler;
54  double* borrowed_lumped_mass_matrix; // borrowed reference
55 };
56 
57 } // namespace paso
58 
59 #endif // __PASO_FLUXLIMITER_H__
60 
paso::LARGE_POSITIVE_FLOAT
static const real_t LARGE_POSITIVE_FLOAT
Definition: FCT_Solver.cpp:43
paso::FCT_FluxLimiter::FCT_FluxLimiter
FCT_FluxLimiter(const_TransportProblem_ptr tp)
Definition: FluxLimiter.cpp:32
PASO_DLL_API
#define PASO_DLL_API
Definition: paso/src/system_dep.h:29
escript::DataTypes::real_t
double real_t
type of all real-valued scalars in escript
Definition: DataTypes.h:52
FluxLimiter.h
paso::FCT_FluxLimiter::addLimitedFluxes_Complete
void addLimitedFluxes_Complete(double *b)
Definition: FluxLimiter.cpp:207
paso::SystemMatrix
this class holds a (distributed) stiffness matrix
Definition: SystemMatrix.h:50
paso::FCT_FluxLimiter::setU_tilde
void setU_tilde(const double *Mu_tilde)
Definition: FluxLimiter.cpp:62
paso::Coupler
Definition: Coupler.h:100
paso::FCT_FluxLimiter::u_tilde_coupler
Coupler_ptr< real_t > u_tilde_coupler
Definition: FluxLimiter.h:53
paso::SystemMatrixPattern_ptr
boost::shared_ptr< SystemMatrixPattern > SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:40
paso::FCT_FluxLimiter
Definition: FluxLimiter.h:27
escript::DataTypes::real_t_max
real_t real_t_max()
Returns the maximum finite value for the real_t type.
Definition: DataTypes.h:90
paso::const_SystemMatrixPattern_ptr
boost::shared_ptr< const SystemMatrixPattern > const_SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:42
paso::FCT_FluxLimiter::getTotalNumRows
dim_t getTotalNumRows() const
Definition: FluxLimiter.h:31
paso::SystemMatrix_ptr
boost::shared_ptr< SystemMatrix< T > > SystemMatrix_ptr
Definition: SystemMatrix.h:42
paso::FCT_FluxLimiter::mpi_info
escript::JMPI mpi_info
Definition: FluxLimiter.h:46
escript::DataTypes::dim_t
index_t dim_t
Definition: DataTypes.h:66
escript::JMPI
boost::shared_ptr< JMPI_ > JMPI
Definition: EsysMPI.h:74
paso::Coupler_ptr
boost::shared_ptr< Coupler< T > > Coupler_ptr
Definition: Coupler.h:43
paso::FCT_FluxLimiter::R
double * R
Definition: FluxLimiter.h:50
paso::FCT_FluxLimiter::borrowed_lumped_mass_matrix
double * borrowed_lumped_mass_matrix
Definition: FluxLimiter.h:54
paso::FCT_FluxLimiter::u_tilde
double * u_tilde
Definition: FluxLimiter.h:48
paso::FCT_FluxLimiter::R_coupler
Coupler_ptr< real_t > R_coupler
Definition: FluxLimiter.h:52
paso::FCT_FluxLimiter::antidiffusive_fluxes
SystemMatrix_ptr< double > antidiffusive_fluxes
Definition: FluxLimiter.h:45
paso::const_TransportProblem_ptr
boost::shared_ptr< const TransportProblem > const_TransportProblem_ptr
Definition: Transport.h:34
Transport.h
paso::FCT_FluxLimiter::addLimitedFluxes_Start
void addLimitedFluxes_Start()
Definition: FluxLimiter.cpp:137
escript::DataTypes::index_t
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:61
paso::FCT_FluxLimiter::~FCT_FluxLimiter
~FCT_FluxLimiter()
Definition: FluxLimiter.cpp:53
paso::FCT_FluxLimiter::getFluxPattern
SystemMatrixPattern_ptr getFluxPattern() const
Definition: FluxLimiter.h:36
paso
Definition: BiCGStab.cpp:25
paso::FCT_FluxLimiter::dt
double dt
Definition: FluxLimiter.h:47
paso::const_SystemMatrix_ptr
boost::shared_ptr< const SystemMatrix< T > > const_SystemMatrix_ptr
Definition: SystemMatrix.h:43
paso::FCT_FluxLimiter::MQ
double * MQ
Definition: FluxLimiter.h:49