Skip to content
Snippets Groups Projects
Commit 15e13c02 authored by Robert K's avatar Robert K
Browse files

Make euler example work with FV space.

parent 8dbfe3ef
No related branches found
No related tags found
1 merge request!48[feature][FVOperator] A specialized FV operator for finite volume.
......@@ -63,7 +63,8 @@ namespace Fem
typedef AlgorithmConfigurator< GridSelectorGridType,
Galerkin::Enum::dg,
Adaptivity::Enum::yes,
DiscreteFunctionSpaces::Enum::orthonormal,
DiscreteFunctionSpaces::Enum::finitevolume,
//DiscreteFunctionSpaces::Enum::orthonormal,
//DiscreteFunctionSpaces::Enum::lagrange,
//DiscreteFunctionSpaces::Enum::gausslobatto,
//DiscreteFunctionSpaces::Enum::gausslegendre,
......@@ -125,8 +126,7 @@ namespace Fem
};
private:
typedef typename AC::template DefaultOpTraits< ModelType, FunctionSpaceType, polOrd >
OpTraits;
typedef typename AC::template DefaultOpTraits< ModelType, FunctionSpaceType, polOrd > OpTraits;
typedef DGAdaptationIndicatorOperator< OpTraits > IndicatorType;
typedef Estimator< DiscreteFunctionType, typename ModelType::ProblemType > GradientIndicatorType ;
public:
......
......@@ -7,7 +7,7 @@ set( GRIDTYPE ALUGRID_CUBE )
#set( GRIDTYPE SPGRID_COUNT_FLOPS )
#set( GRIDTYPE SPGRID )
set( GRIDDIM 2 )
set( POLORDER 2 )
set( POLORDER 0 )
add_definitions( "-D${GRIDTYPE}" )
add_definitions( "-DGRIDDIM=${GRIDDIM}" )
......
fem.verboserank: 0
# verbosity level
fem.verbositylevel: 3
# PROBLEM SELECTION
#------------------
# possible: "sod" , "withman", "withmansmooth", "smooth1d" , "ffs" , "diffractio
......@@ -22,8 +25,8 @@ femdg.limiter.indicator: jump
#femdg.nonblockingcomm: false
femdg.stepper.printcount: 1000
fem.io.outputformat: none
fem.io.savestep: 10.0
fem.io.outputformat: vtk-cell
fem.io.savestep: 0.01
fem.ode.odesolver: EX
fem.ode.order: 3
......
# toggle verbosity of on process
# (-1: off, p: process p is verbose, in general 0 is a good choice)
fem.verboserank: 0
fem.verbositylevel: 1
fem.verbositylevel: 3
# OMP THREADS
#------------
......
......@@ -81,7 +81,9 @@ namespace Fem
//! Lagrange space with GaussLobatto interpolation points
gausslobatto,
//! Lagrange space with GaussLegendre interpolation points
gausslegendre
gausslegendre,
//! FiniteVolume space
finitevolume
};
}
......
......@@ -46,6 +46,7 @@
#include <dune/fem-dg/assemble/primalmatrix.hh>
#include <dune/fem/space/discontinuousgalerkin.hh>
#include <dune/fem/space/lagrange.hh>
#include <dune/fem/space/finitevolume.hh>
#include <dune/fem-dg/operator/dg/operatortraits.hh>
......@@ -525,6 +526,12 @@ namespace Fem
};
#endif
template< class FunctionSpaceImp, class GridPartImp, int polOrder>
struct DiscreteFunctionSpaceSelector< FunctionSpaceImp, GridPartImp, polOrder, DiscreteFunctionSpaces::Enum::finitevolume, Galerkin::Enum::dg >
{
typedef FiniteVolumeSpace< FunctionSpaceImp, GridPartImp, polOrder > type;
};
template< class ModelImp, class DiscreteFunctionSpaceImp,
DiffusionFlux::Enum diffFluxId, Formulation::Enum form >
struct DiffusionFluxSelector;
......
......@@ -30,7 +30,7 @@
// otherwise optional values from MIN_POLORD to MAX_POLORD are provided.
// This will increase compile time and program size and should only be chosen
// for special purpose
#if POLORDER
#if POLORDER >= 0
#define MIN_POLORD POLORDER
#define MAX_POLORD POLORDER
#else
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment