Skip to content
Snippets Groups Projects

Exchange troubled cell

Merged Robert K requested to merge feature/exchange-troubled-cell into master
Compare and Show latest version
18 files
+ 333
167
Compare changes
  • Side-by-side
  • Inline
Files
18
@@ -98,6 +98,7 @@ namespace Fem
typedef Estimator< DestinationType, typename ModelType::ProblemType > GradientIndicatorType ;
typedef AdaptIndicator< IndicatorType, GradientIndicatorType > AdaptIndicatorType;
typedef typename FullOperatorType :: TroubledCellIndicatorType TroubledCellIndicatorType;
// solver selection, available fem, istl, petsc, ...
typedef typename MatrixFreeSolverSelector< solverId, symmetric > :: template LinearInverseOperatorType< DiscreteFunctionSpaceType, DiscreteFunctionSpaceType > LinearSolverType ;
@@ -183,6 +184,13 @@ namespace Fem
}
}
void setTroubledCellIndicator(TroubledCellIndicatorType indicator)
{
fullOperator_.setTroubledCellIndicator(indicator);
explOperator_.setTroubledCellIndicator(indicator);
}
/** \copydoc SpaceOperatorInterface::setTime */
void setTime( const double time )
{
@@ -191,6 +199,12 @@ namespace Fem
double timeStepEstimate() const { return fullOperator_.timeStepEstimate(); }
//! return number of interior elements visited by the operator
inline size_t gridSizeInterior() const
{
return std::max( fullOperator_.numberOfElements(), explOperator_.numberOfElements() );
}
//// End Methods from SpaceOperatorInterface /////
protected:
@@ -222,7 +236,6 @@ namespace Fem
mutable double fixedTimeStep_ ;
mutable bool initialized_;
};
} // end namespace Fem
Loading