escript  Revision_
FluxLimiter.h
Go to the documentation of this file.
1 
2 /*****************************************************************************
3 *
4 * Copyright (c) 2003-2018 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 by Centre for Geoscience Computing (GeoComp)
14 *
15 *****************************************************************************/
16 
17 #ifndef __PASO_FLUXLIMITER_H__
18 #define __PASO_FLUXLIMITER_H__
19 
20 #include "Transport.h"
21 
22 namespace paso {
23 
24 
26 struct FCT_FluxLimiter
27 {
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);
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:43
escript::DataTypes::real_t
double real_t
type of all real-valued scalars in escript
Definition: DataTypes.h:73
FluxLimiter.h
paso::FCT_FluxLimiter::addLimitedFluxes_Complete
void addLimitedFluxes_Complete(double *b)
Definition: FluxLimiter.cpp:218
paso::SystemMatrix
this class holds a (distributed) stiffness matrix
Definition: SystemMatrix.h:58
paso::FCT_FluxLimiter::setU_tilde
void setU_tilde(const double *Mu_tilde)
Definition: FluxLimiter.cpp:73
paso::Coupler
Definition: Coupler.h:52
paso::FCT_FluxLimiter::u_tilde_coupler
Coupler_ptr< real_t > u_tilde_coupler
Definition: FluxLimiter.h:76
paso::SystemMatrix_ptr
boost::shared_ptr< SystemMatrix > SystemMatrix_ptr
Definition: SystemMatrix.h:51
paso::SystemMatrixPattern_ptr
boost::shared_ptr< SystemMatrixPattern > SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:50
escript::DataTypes::real_t_max
real_t real_t_max()
Returns the maximum finite value for the real_t type.
Definition: DataTypes.h:111
paso::const_SystemMatrixPattern_ptr
boost::shared_ptr< const SystemMatrixPattern > const_SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:52
paso::FCT_FluxLimiter::getTotalNumRows
dim_t getTotalNumRows() const
Definition: FluxLimiter.h:54
paso::FCT_FluxLimiter::mpi_info
escript::JMPI mpi_info
Definition: FluxLimiter.h:69
escript::DataTypes::dim_t
index_t dim_t
Definition: DataTypes.h:87
escript::JMPI
boost::shared_ptr< JMPI_ > JMPI
Definition: EsysMPI.h:70
paso::Coupler_ptr
boost::shared_ptr< Coupler< T > > Coupler_ptr
Definition: Coupler.h:53
paso::FCT_FluxLimiter::R
double * R
Definition: FluxLimiter.h:73
paso::FCT_FluxLimiter::borrowed_lumped_mass_matrix
double * borrowed_lumped_mass_matrix
Definition: FluxLimiter.h:77
paso::FCT_FluxLimiter::u_tilde
double * u_tilde
Definition: FluxLimiter.h:71
paso::FCT_FluxLimiter::R_coupler
Coupler_ptr< real_t > R_coupler
Definition: FluxLimiter.h:75
paso::FCT_FluxLimiter::antidiffusive_fluxes
SystemMatrix_ptr antidiffusive_fluxes
Definition: FluxLimiter.h:68
paso::const_SystemMatrix_ptr
boost::shared_ptr< const SystemMatrix > const_SystemMatrix_ptr
Definition: SystemMatrix.h:53
paso::const_TransportProblem_ptr
boost::shared_ptr< const TransportProblem > const_TransportProblem_ptr
Definition: Transport.h:32
Transport.h
paso::FCT_FluxLimiter::addLimitedFluxes_Start
void addLimitedFluxes_Start()
Definition: FluxLimiter.cpp:148
escript::DataTypes::index_t
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:82
paso::FCT_FluxLimiter::~FCT_FluxLimiter
~FCT_FluxLimiter()
Definition: FluxLimiter.cpp:64
PASO_DLL_API
#define PASO_DLL_API
Definition: Paso.h:58
paso::FCT_FluxLimiter::getFluxPattern
SystemMatrixPattern_ptr getFluxPattern() const
Definition: FluxLimiter.h:59
paso
Definition: BiCGStab.cpp:25
paso::FCT_FluxLimiter::dt
double dt
Definition: FluxLimiter.h:70
paso::FCT_FluxLimiter::MQ
double * MQ
Definition: FluxLimiter.h:72