#include <postprocessingfvlpt.h>

Inheritance diagram for PostProcessingFvLPT< nDim, SysEqn >:

Collaboration diagram for PostProcessingFvLPT< nDim, SysEqn >:

Public Types
using	lptSolver = LPT< nDim >

using	fv = FvCartesianSolverXD< nDim, SysEqn >

using	BaseFv = PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >

using	BaseLpt = PostProcessing< nDim, PostProcessingLPT< nDim > >

using	ppFv = PostProcessingFv< nDim, SysEqn >

using	ppLpt = PostProcessingLPT< nDim >

Public Types inherited from PostProcessingLPT< nDim >
using	SolverType = LPT< nDim >

using	Base = PostProcessing< nDim, PostProcessingLPT< nDim > >

Public Types inherited from PostProcessing< nDim, PostProcessingLPT< nDim > >
using	PPGrid = CartesianGrid< nDim >

using	Cell = typename PPGrid::Cell

using	GridProxy = typename maia::grid::Proxy< nDim >

using	Data = PostData< nDim >

Public Member Functions
	PostProcessingFvLPT (MInt postprocessingId_, PostData< nDim > data, fv ppFvSolver_, lptSolver *ppLptSolver_)

virtual	~PostProcessingFvLPT ()

Solver *	mSolver () const override

void	initPostProcessing () override

void	postprocessPreSolve () override

void	postprocessPostSolve () override

void	postprocessInSolve (const MBool finalTimeStep) override

void	postprocessSolution () override

void	initSprayData () override

Public Member Functions inherited from PostProcessingLPT< nDim >
	PostProcessingLPT (MInt postprocessingId_, PostData< nDim > data, SolverType ppSolver_)

virtual	~PostProcessingLPT ()

SolverType &	solver () const

void	initSprayData () override
	init arrays More...

Public Member Functions inherited from PostProcessing< nDim, PostProcessingLPT< nDim > >
	PostProcessing (MInt postprocessingId_, PostData< nDim > *data)

	~PostProcessing ()
	Destructor for the massive paralle postprocessing solver. More...

void	initPostProcessing () override
	Reads all required properties in and prepares for postprocessing. More...

void	postprocessPreSolve () override

void	postprocessPostSolve () override

void	postprocessInSolve (MBool finalTimeStep) override

void	postprocessSolution () override

Solver *	mSolver () const override

Data &	postData () const

Public Member Functions inherited from PostProcessingInterface
	PostProcessingInterface (const MInt postprocessingId_)

virtual	~PostProcessingInterface ()=default

	PostProcessingInterface (const PostProcessingInterface &)=delete

PostProcessingInterface &	operator= (const PostProcessingInterface &)=delete

virtual void	initPostProcessing ()=0

virtual void	postprocessPreSolve ()=0

virtual void	postprocessPostSolve ()=0

virtual void	postprocessInSolve (MBool finalTimeStep)=0

virtual void	postprocessSolution ()=0

virtual Solver *	mSolver () const =0

MInt	a_postprocessingId () const

Static Public Attributes
static MString	s_ppType

Protected Member Functions
void	computeSprayData () override
	base function which is called every time-Step in postprocessInSolve More...

void	writeSprayData () override

void	initLPTSolutionFile () override

void	writeLPTSolutionFile () override

Protected Member Functions inherited from PostProcessingLPT< nDim >
void	initLPTSolutionFile () override
	init function for LPT particle solution file More...

void	writeLPTSolutionFile () override
	base function which is called every time-Step in postprocessInSolve More...

void	initParticleStatistics () override
	init function for LPT particle statistics file More...

void	computeParticleStatistics () override
	compute average quantites of particle phase average particle Reynolds number ReP the root mean square particle velocity vrms the total kinetic energy of the particles More...

void	initParticleLog ()
	initialize particle.log More...

void	writeParticleLog ()
	write particle Log file, containing particle position and velocity at timestep t More...

void	particleMass ()
	computes conervation variables (mass/momentum/energy!) More...

void	parcelStatistics ()
	computes parcel statistics More...

void	particlePenetration ()
	calculate vertical and horizontal penetration from given coordinate More...

MInt	getInjectionData ()

void	advanceDataStep ()

void	resetDataStep ()

Protected Member Functions inherited from PostProcessing< nDim, PostProcessingLPT< nDim > >
MBool	isActive () const

void	transferSensorWeights ()

void	initReduceToLevelAvg ()
	Initializes properties for grid level reduction. More...

void	initTimeStepProperties ()
	Initializes timestep properties for postprocessing. More...

void	initTimeStepPropertiesSlice ()

void	initTimeStepPropertiesLine ()

virtual void	initAveragingProperties ()
	Initialize properties relevant for temporal averaging. More...

void	initProbePoint ()
	Initializes properties for point probing. More...

void	initProbeLine ()
	Initializes properties and memory for line probing. More...

void	initProbeSlice ()
	Initializes properties and memory for slice probing. More...

void	initSpatialAveraging ()

void	initProbeArbitraryLine ()
	initializes properties and data for arbitrary line probing More...

void	initProbeArbitrarySlice ()
	initializes properties and data for arbitrary slice probing More...

void	initProbeLinePeriodic ()

virtual void	probeLinePeriodic ()

virtual void	probeLinePeriodicPost ()

virtual void	initSprayData ()

virtual void	computeSprayData ()

virtual void	writeSprayData ()

virtual void	initLPTSolutionFile ()

virtual void	writeLPTSolutionFile ()

virtual void	initParticleStatistics ()

virtual void	computeParticleStatistics ()

virtual void	initPLIsoTurbulenceStatistics ()

virtual void	computePLIsoTurbulenceStatistics ()

virtual void	initAverageVariables ()
	allocates memory for averageSolutions() and averageSolutionsInSolve() More...

void	pp_saveCoarseSolution ()

void	averageSolutions ()
	Averages solutions. More...

virtual void	averageSolutionsInSolve ()

void	computeAndSaveDivergence ()

void	computeAndSaveMean ()

void	initAverageSolutionsSlice ()
	Initializes properties for slice averaging. More...

void	averageSolutionsSlice ()

void	probeLocation ()

void	findClosestProbePointsGridCells ()

virtual void	initPointSamplingData ()

virtual void	savePointSamplingData ()

virtual void	initSurfaceSamplingData ()

virtual void	saveSurfaceSamplingData ()

virtual void	initVolumeSamplingData ()

virtual void	saveVolumeSamplingData ()

virtual void	initIsoTurbulenceStatistics ()

virtual void	computeIsoTurbulenceStatistics ()

void	initWritePointData ()

void	writePointData ()

void	spatialAveraging ()

void	spatialAveragingPost ()

void	createCellToMap1D (std::map< MFloat, MInt, coord_comp_1d_ > &coordinates, std::map< MFloat, MFloat, coord_comp_1d_ > &cell_to_map)

void	createCellToMap2D (std::map< std::pair< MFloat, MFloat >, MInt, coord_comp_2d_ > &coordinates, std::map< std::pair< MFloat, MFloat >, std::pair< MFloat, MFloat >, coord_comp_2d_ > &cell_to_map)

void	probeLine ()

void	probeLinePre ()

void	probeLinePost ()

void	probeSlice ()

void	probeSliceIn ()

void	probeSlicePre ()

void	probeSlicePost ()

void	collectMinLvlCells ()

void	findContainingCell (const MFloat *coord, MInt &cellId)

void	probeArbitraryLine ()

void	probeArbitraryLinePost ()

void	probeArbitrarySlice ()

void	probeArbitrarySlicePost ()

virtual void	initMovingAverage ()
	Initializes properties and allocates memory for moving averaging. More...

void	movingAverage ()

void	movingAveragePost ()

void	initCorrelation ()

void	initPeriodicSliceAverage ()
	Initializes the periodic averaging on a slice. More...

void	periodicSliceAverage ()

MInt	findNearestGridCell (const MFloat *coord)

void	movePointsToGrid (MFloat *in_points, MInt in_noPoints, MInt in_moveType)

void	getSampleVariables (MInt cellId, const MFloat *&vars, MBool mode)

void	getSampleVariables (MInt const cellId, std::vector< MFloat > &vars)

void	calcSamplingVar (const MInt cellId, const MInt sampleVarId, MFloat *const vars)

void	saveSliceAiaFileFormat (const MInt step, const MInt noVars, MFloatScratchSpace &vars, const MInt sliceId)

MFloat	calcDivergence (const MInt cellIdSolver)

virtual void	calcVorticity (const MFloatTensor &deriv, MFloat vorticity[nDim *2 - 3])

virtual void	getVorticity (MFloat *const vorticity)

virtual void	getVorticityT (MFloat *const vorticity)

virtual void	getPrimitiveVariables (MInt, MFloat , MFloat , MInt)

virtual void	computeAcousticSourceTermQe (MFloatScratchSpace &, MFloatScratchSpace &, MFloatScratchSpace &, MFloatScratchSpace &)

void	getSampleVarsDerivatives (const MInt cellId, const MFloat *&vars)

MBool	getSampleVarsDerivatives (const MInt, std::vector< MFloat > &)

virtual MFloat &	vorticityAtCell (const MInt cellId, const MInt dir)

virtual MFloat	getBoundaryHeatFlux (const MInt cellId) const

virtual void	getPrimitiveVarNames (MString *names) const

void	neededMeanVarsForSourceTerm (const MInt sourceTerm, std::vector< MInt > &meanVars) const

MBool	isMeanFile ()

Private Member Functions
lptSolver &	lpt () const

fv &	fvSolver () const

void	updateData ()
	updates all post-processing data from postprocessing-fv or postprocessing-lpt More...

void	advanceDataStep ()

void	resetDataStep ()

Private Member Functions inherited from PostProcessingFv< nDim, SysEqn >
	PostProcessingFv (MInt postprocessingId_, PostData< nDim > data, SolverType ppSolver_)

virtual	~PostProcessingFv ()

void	initPostProcessing () override

SolverType &	solver () const

void	probeLinePeriodicPost () override

void	probeLinePeriodic () override

void	initMovingAverage () override
	Initializes properties and allocates memory for moving averaging. More...

void	initAveragingProperties () override
	Initialize properties relevant for temporal averaging. More...

void	initPointSamplingData () override

void	savePointSamplingData () override

void	initSurfaceSamplingData () override

void	saveSurfaceSamplingData () override

void	initVolumeSamplingData () override

void	saveVolumeSamplingData () override

void	initSprayData () override

void	calcVorticity (const MFloatTensor &deriv, MFloat vorticity[nDim *2 - 3]) override

void	getVorticity (MFloat *const vorticity) override

void	getVorticityT (MFloat *const vorticity) override

void	getSampleVarsDerivatives (MInt cellId, const MFloat *&vars)

MBool	getSampleVarsDerivatives (const MInt cellId, std::vector< MFloat > &vars)

MFloat &	vorticityAtCell (const MInt cellId, const MInt dir) override

void	getPrimitiveVariables (MInt cellId, MFloat Xp, MFloat vars, MInt order) override

void	computeAcousticSourceTermQe (MFloatScratchSpace &QeI, MFloatScratchSpace &QeIII, MFloatScratchSpace &cSquared, MFloatScratchSpace &drhodt) override

void	getPrimitiveVarNames (MString *names) const override

MFloat	getBoundaryHeatFlux (const MInt cellId) const override

void	vapourPenetration (MFloat spawnCoord[nDim])

void	vapourMass (const MInt)

void	advanceDataStep ()

void	resetDataStep ()

Private Member Functions inherited from PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >
	PostProcessing (MInt postprocessingId_, PostData< nDim > *data)

	~PostProcessing ()
	Destructor for the massive paralle postprocessing solver. More...

void	initPostProcessing () override
	Reads all required properties in and prepares for postprocessing. More...

void	postprocessPreSolve () override

void	postprocessPostSolve () override

void	postprocessInSolve (MBool finalTimeStep) override

void	postprocessSolution () override

Solver *	mSolver () const override

Data &	postData () const

MBool	isActive () const

void	transferSensorWeights ()

void	initReduceToLevelAvg ()
	Initializes properties for grid level reduction. More...

void	initTimeStepProperties ()
	Initializes timestep properties for postprocessing. More...

void	initTimeStepPropertiesSlice ()

void	initTimeStepPropertiesLine ()

virtual void	initAveragingProperties ()
	Initialize properties relevant for temporal averaging. More...

void	initProbePoint ()
	Initializes properties for point probing. More...

void	initProbeLine ()
	Initializes properties and memory for line probing. More...

void	initProbeSlice ()
	Initializes properties and memory for slice probing. More...

void	initSpatialAveraging ()

void	initProbeArbitraryLine ()
	initializes properties and data for arbitrary line probing More...

void	initProbeArbitrarySlice ()
	initializes properties and data for arbitrary slice probing More...

void	initProbeLinePeriodic ()

virtual void	probeLinePeriodic ()

virtual void	probeLinePeriodicPost ()

virtual void	initSprayData ()

virtual void	computeSprayData ()

virtual void	writeSprayData ()

virtual void	initLPTSolutionFile ()

virtual void	writeLPTSolutionFile ()

virtual void	initParticleStatistics ()

virtual void	computeParticleStatistics ()

virtual void	initPLIsoTurbulenceStatistics ()

virtual void	computePLIsoTurbulenceStatistics ()

virtual void	initAverageVariables ()
	allocates memory for averageSolutions() and averageSolutionsInSolve() More...

void	pp_saveCoarseSolution ()

void	averageSolutions ()
	Averages solutions. More...

virtual void	averageSolutionsInSolve ()

void	computeAndSaveDivergence ()

void	computeAndSaveMean ()

void	initAverageSolutionsSlice ()
	Initializes properties for slice averaging. More...

void	averageSolutionsSlice ()

void	probeLocation ()

void	findClosestProbePointsGridCells ()

virtual void	initPointSamplingData ()

virtual void	savePointSamplingData ()

virtual void	initSurfaceSamplingData ()

virtual void	saveSurfaceSamplingData ()

virtual void	initVolumeSamplingData ()

virtual void	saveVolumeSamplingData ()

virtual void	initIsoTurbulenceStatistics ()

virtual void	computeIsoTurbulenceStatistics ()

void	initWritePointData ()

void	writePointData ()

void	spatialAveraging ()

void	spatialAveragingPost ()

void	createCellToMap1D (std::map< MFloat, MInt, coord_comp_1d_ > &coordinates, std::map< MFloat, MFloat, coord_comp_1d_ > &cell_to_map)

void	createCellToMap2D (std::map< std::pair< MFloat, MFloat >, MInt, coord_comp_2d_ > &coordinates, std::map< std::pair< MFloat, MFloat >, std::pair< MFloat, MFloat >, coord_comp_2d_ > &cell_to_map)

void	probeLine ()

void	probeLinePre ()

void	probeLinePost ()

void	probeSlice ()

void	probeSliceIn ()

void	probeSlicePre ()

void	probeSlicePost ()

void	collectMinLvlCells ()

void	findContainingCell (const MFloat *coord, MInt &cellId)

void	probeArbitraryLine ()

void	probeArbitraryLinePost ()

void	probeArbitrarySlice ()

void	probeArbitrarySlicePost ()

virtual void	initMovingAverage ()
	Initializes properties and allocates memory for moving averaging. More...

void	movingAverage ()

void	movingAveragePost ()

void	initCorrelation ()

void	initPeriodicSliceAverage ()
	Initializes the periodic averaging on a slice. More...

void	periodicSliceAverage ()

MInt	findNearestGridCell (const MFloat *coord)

void	movePointsToGrid (MFloat *in_points, MInt in_noPoints, MInt in_moveType)

void	getSampleVariables (MInt cellId, const MFloat *&vars, MBool mode)

void	getSampleVariables (MInt const cellId, std::vector< MFloat > &vars)

void	calcSamplingVar (const MInt cellId, const MInt sampleVarId, MFloat *const vars)

void	saveSliceAiaFileFormat (const MInt step, const MInt noVars, MFloatScratchSpace &vars, const MInt sliceId)

MFloat	calcDivergence (const MInt cellIdSolver)

virtual void	calcVorticity (const MFloatTensor &deriv, MFloat vorticity[nDim *2 - 3])

virtual void	getVorticity (MFloat *const vorticity)

virtual void	getVorticityT (MFloat *const vorticity)

virtual void	getPrimitiveVariables (MInt, MFloat , MFloat , MInt)

virtual void	computeAcousticSourceTermQe (MFloatScratchSpace &, MFloatScratchSpace &, MFloatScratchSpace &, MFloatScratchSpace &)

void	getSampleVarsDerivatives (const MInt cellId, const MFloat *&vars)

MBool	getSampleVarsDerivatives (const MInt, std::vector< MFloat > &)

virtual MFloat &	vorticityAtCell (const MInt cellId, const MInt dir)

virtual MFloat	getBoundaryHeatFlux (const MInt cellId) const

virtual void	getPrimitiveVarNames (MString *names) const

void	neededMeanVarsForSourceTerm (const MInt sourceTerm, std::vector< MInt > &meanVars) const

MBool	isMeanFile ()

Private Member Functions inherited from PostProcessingInterface
	PostProcessingInterface (const MInt postprocessingId_)

virtual	~PostProcessingInterface ()=default

	PostProcessingInterface (const PostProcessingInterface &)=delete

PostProcessingInterface &	operator= (const PostProcessingInterface &)=delete

virtual void	initPostProcessing ()=0

virtual void	postprocessPreSolve ()=0

virtual void	postprocessPostSolve ()=0

virtual void	postprocessInSolve (MBool finalTimeStep)=0

virtual void	postprocessSolution ()=0

virtual Solver *	mSolver () const =0

MInt	a_postprocessingId () const

Private Attributes
lptSolver *	m_ppSolverLpt

fv *	m_ppSolverFv

MInt	m_fvLPTSpeciesId = 0

MInt	m_sprayDataStep = 0

MFloat **	m_cInjData

Private Attributes inherited from PostProcessingFv< nDim, SysEqn >
std::unique_ptr< PointData< nDim, SolverType > >	m_pointData

std::unique_ptr< SurfaceData< nDim, SolverType > >	m_surfaceData

std::unique_ptr< VolumeData< nDim, SolverType > >	m_volumeData

Private Attributes inherited from PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >
MInt	m_noLocalVars

MInt	m_restartTimeStep

MBool	m_statisticCombustionAnalysis

MBool	m_acousticAnalysis

MInt	m_noPostprocessingOps

MString *	m_postprocessingOps

tvecpost	m_postprocessingSolution

std::vector< tvecpost >	m_postprocessingMethods

std::vector< std::vector< MString > >	m_postprocessingMethodsDesc

Data *	m_postData

MInt	m_noVariables

MBool	m_square

MBool	m_skewness

MBool	m_kurtosis

MBool	m_computeAndSaveMean

MBool	m_twoPass

MBool	m_useKahan

MBool	m_correlation

std::vector< MInt > *	m_cell2globalIndex

MInt *	m_noPeriodicSliceCells

MFloat **	m_sliceAverage

MInt	m_globalnoSlicePositions

std::multimap< MFloat, MFloat >	m_sliceGlobalPositions

MBool	m_averageVorticity

MBool	m_averageSpeedOfSound

MFloat	m_gamma

MInt	m_movingAverageStartTimestep

MInt	m_movingAverageStopTimestep

MInt	m_movingAverageInterval

MInt	m_movingAverageDataPoints

MInt	m_movingAverageCounter

MFloat **	m_movAvgVariables

MInt	m_movAvgNoVariables

MString *	m_movAvgVarNames

MInt	m_spatialDirection1

MInt	m_spatialDirection2

MFloat	m_spatialWeightSinglePoint

MFloat	m_spatialLvlWeight [20]

MInt *	m_spatialDispls

MInt *	m_spatialVarsDispls

MInt *	m_spatialRecvcnts

MInt *	m_spatialVarsRecvcnts

MInt	m_spatialCoordSum

std::map< MFloat, MInt, coord_comp_1d_ >	m_spatialLineCoordinates

std::map< MFloat, MFloat, coord_comp_1d_ >	m_spatialLineCellToMap

std::map< MFloat, MInt, coord_comp_1d_ >	m_spatialGlobalLineCoordinates

std::map< MFloat, MFloat, coord_comp_1d_ >	m_spatialGlobalLineCellToMap

MInt	m_spatialLineNoCells

MFloat *	m_spatialLineAllVars

MFloat *	m_spatialLineAllCoord

std::map< std::pair< MFloat, MFloat >, MInt, coord_comp_2d_ >	m_spatialPlaneCoordinates

std::map< std::pair< MFloat, MFloat >, std::pair< MFloat, MFloat >, coord_comp_2d_ >	m_spatialPlaneCellToMap

std::map< std::pair< MFloat, MFloat >, MInt, coord_comp_2d_ >	m_spatialGlobalPlaneCoordinates

std::map< std::pair< MFloat, MFloat >, std::pair< MFloat, MFloat >, coord_comp_2d_ >	m_spatialGlobalPlaneCellToMap

MInt	m_spatialPlaneNoCells

MFloat *	m_spatialPlaneAllVars

MFloat *	m_spatialPlaneAllCoord

MInt	m_noProbeLines

MInt *	m_probeLineDirection

MInt *	m_probeLinePeriodic

MFloat **	m_probeLineCoordinates

MInt **	m_probeLineIds

MInt *	m_noProbeLineIds

MFloat **	m_probeLinePositions

MInt *	m_globalNoProbeLineIds

MInt *	m_probeLineOffsets

MInt **	m_noGlobalProbeLineIds

MInt	m_probeLineInterval

MInt	m_probeLineStartTimestep

MInt	m_probeLineStopTimestep

MInt *	m_correlationDirection

MInt *	m_correlationVariableIndex

MFloat **	m_correlationCoordinates

MInt	m_noCorrelationLines

MInt *	m_noCorrelationIds

MInt **	m_correlationIds

MInt **	m_globalCorrelationIds

MFloat **	m_globalCorrelationPositions

MInt *	m_globalNoCorrelationIds

MFloat **	m_correlationPositions

MInt **	m_correlationIndexMapping

MFloat **	m_correlationExchangeVar

MFloat **	m_correlationExchangeVarMean

MFloat **	m_correlationExchangeVarRMS

MFloat **	m_globalCorrelationExchangeVar

MFloat **	m_globalCorrelationExchangeVarMean

MFloat **	m_globalCorrelationExchangeVarRMS

MInt *	m_probeLineAverageDirection

MFloat **	m_probeLineAverageCoordinates

MInt **	m_probeLineAverageCoordinatesSign

MInt **	m_noProbeLineAverageIds

MFloat ***	m_globalProbeLineAverageVars

MInt	m_noProbeLineAverageSteps

MInt **	m_globalNoProbeLineAverageIds

MInt ***	m_probeLineAverageIds

MFloat ***	m_probeLineAveragePositions

MInt	m_noProbeSlices

MInt *	m_probeSliceDir

MFloat *	m_probeSliceCoordinate

MInt **	m_probeSliceIds

MInt *	m_noProbeSliceIds

MFloat **	m_probeSlicePositions

MInt *	m_globalNoProbeSliceIds

MInt *	m_probeSliceOffsets

MInt **	m_noGlobalProbeSliceIds

MString *	m_probeSliceGridNames

MString *	m_sliceAxis

MFloat *	m_sliceIntercept

MBool	m_sliceAiaFileFormat

MBool	m_optimizedSliceIo

MInt *	m_noProbeSliceNoHilbertIds

MInt *	m_noProbeSliceMaxNoHilbertIds

MInt **	m_noProbeSliceHilbertInfo

MInt *	m_noProbeSliceNoContHilbertIds

MInt *	m_noProbeSliceMaxNoContHilbertIds

MInt **	m_noProbeSliceContHilbertInfo

MInt	m_probeSliceInterval

MInt	m_probeSliceStartTimestep

MInt	m_probeSliceStopTimestep

std::vector< MInt >	m_sliceVarIds
	List of slice variables. More...

std::vector< MInt >	m_noSliceVars
	Number of variables for each slice variable. More...

std::vector< std::vector< MString > >	m_sliceVarNames
	List of variable names for each slice variable. More...

MInt	m_noMinLvlIds

MInt *	m_minLvlIds

MInt	m_movePointsToGrid

MBool	m_spatialAveraging

MFloat *	m_arbLinePoints

MInt	m_noArbLines

MInt *	m_noArbLineIds

MFloat *	m_arbLinePointsDistribution

MInt *	m_arbLineHasNewDistribution

MInt *	m_arbLineFull

MInt *	m_moveLinePointsToGrid

MInt *	m_globalNoArbLineIds

MInt **	m_arbLineIds

MInt *	m_arbLineOffsets

MFloat **	m_arbLineCoordinates

MFloat *	m_arbSlicePoints

MInt	m_noArbSlices

MInt *	m_noArbSlicePoints

MFloat *	m_arbSlicePointsDistribution

MInt *	m_globalNoArbSlicePoints

MInt **	m_arbSliceIds

MInt *	m_arbSliceOffsets

MFloat **	m_arbSliceCoordinates

MString	m_postprocessFileName

MFloat **	m_localVars

MBool	m_averageRestart

MInt	m_averageInterval

MInt	m_averageStartTimestep

MInt	m_averageStopTimestep

MInt	m_noAveragedVorticities

MInt	m_noSpeedOfSoundVars

MInt	m_noSourceVars

MBool	m_needVorticity

MString	m_ReStressesAverageFileName

MFloat **	m_probeCoordinates

MString	m_probePath

MInt *	m_probeCellIds

std::ofstream *	m_probeFileStreams

MInt	m_noProbePoints

MInt	m_probeInterval

MInt	m_probeWriteInterval

MString	m_pdFileName

MInt	m_pdStartTimestep

MInt	m_pdStopTimestep

MInt	m_pdRestartInterval

std::vector< MFloat >	m_pdPoints

std::vector< MInt >	m_pdCells

std::vector< MFloat >	m_pdVars

MInt	m_pdNoPoints

MInt	m_sprayComputeInterval

MInt	m_sprayWriteInterval

MInt	m_sprayDataSize

MInt	m_sprayDataStep

MFloat **	m_particleCV

MFloat **	m_particlePen

MFloat **	m_sprayStat

MFloat **	m_sprayTimes

MFloat **	m_injectionData

MFloat **	m_vapourCV

MFloat **	m_vapourPen

MInt	m_LPTSolutionInterval

MBool	m_forceOutput

std::vector< MInt >	m_activeSourceTerms

std::array< MString, ST::totalNoSourceTerms >	s_sourceTermNames

std::set< MInt >	m_activeMeanVars

Private Attributes inherited from PostProcessingInterface
MBool	m_finalTimeStep = false

MInt	m_postprocessingId

Friends
template<MInt nDim_, class ppType >
class	PostProcessing

template<MInt nDim_>
class	LPT

Additional Inherited Members
Public Attributes inherited from PostProcessing< nDim, PostProcessingLPT< nDim > >
MInt	m_noLocalVars

MInt	m_restartTimeStep

MBool	m_statisticCombustionAnalysis

MBool	m_acousticAnalysis

Public Attributes inherited from PostProcessingInterface
MBool	m_finalTimeStep = false

Protected Types inherited from PostProcessing< nDim, PostProcessingLPT< nDim > >
typedef void(PostProcessing::*	tpost) ()

typedef std::vector< tpost >	tvecpost

Protected Attributes inherited from PostProcessingLPT< nDim >
SolverType *	m_ppSolver

MFloat **	m_cInjData = nullptr

MFloat	m_AvgRe_p = F0

MFloat	m_VpRMS = F0

MFloat	m_EkinP = F0

MBool	m_writeParLog = false

MInt	m_parLogInterval = 50

MBool	m_parLogApp = false

Protected Attributes inherited from PostProcessing< nDim, PostProcessingLPT< nDim > >
MInt	m_noPostprocessingOps

MString *	m_postprocessingOps

tvecpost	m_postprocessingSolution

std::vector< tvecpost >	m_postprocessingMethods

std::vector< std::vector< MString > >	m_postprocessingMethodsDesc

Data *	m_postData

MInt	m_noVariables

MBool	m_square

MBool	m_skewness

MBool	m_kurtosis

MBool	m_computeAndSaveMean

MBool	m_twoPass

MBool	m_useKahan

MBool	m_correlation

std::vector< MInt > *	m_cell2globalIndex

MInt *	m_noPeriodicSliceCells

MFloat **	m_sliceAverage

MInt	m_globalnoSlicePositions

std::multimap< MFloat, MFloat >	m_sliceGlobalPositions

MBool	m_averageVorticity

MBool	m_averageSpeedOfSound

MFloat	m_gamma

MInt	m_movingAverageStartTimestep

MInt	m_movingAverageStopTimestep

MInt	m_movingAverageInterval

MInt	m_movingAverageDataPoints

MInt	m_movingAverageCounter

MFloat **	m_movAvgVariables

MInt	m_movAvgNoVariables

MString *	m_movAvgVarNames

MInt	m_spatialDirection1

MInt	m_spatialDirection2

MFloat	m_spatialWeightSinglePoint

MFloat	m_spatialLvlWeight [20]

MInt *	m_spatialDispls

MInt *	m_spatialVarsDispls

MInt *	m_spatialRecvcnts

MInt *	m_spatialVarsRecvcnts

MInt	m_spatialCoordSum

std::map< MFloat, MInt, coord_comp_1d_ >	m_spatialLineCoordinates

std::map< MFloat, MFloat, coord_comp_1d_ >	m_spatialLineCellToMap

std::map< MFloat, MInt, coord_comp_1d_ >	m_spatialGlobalLineCoordinates

std::map< MFloat, MFloat, coord_comp_1d_ >	m_spatialGlobalLineCellToMap

MInt	m_spatialLineNoCells

MFloat *	m_spatialLineAllVars

MFloat *	m_spatialLineAllCoord

std::map< std::pair< MFloat, MFloat >, MInt, coord_comp_2d_ >	m_spatialPlaneCoordinates

std::map< std::pair< MFloat, MFloat >, std::pair< MFloat, MFloat >, coord_comp_2d_ >	m_spatialPlaneCellToMap

std::map< std::pair< MFloat, MFloat >, MInt, coord_comp_2d_ >	m_spatialGlobalPlaneCoordinates

std::map< std::pair< MFloat, MFloat >, std::pair< MFloat, MFloat >, coord_comp_2d_ >	m_spatialGlobalPlaneCellToMap

MInt	m_spatialPlaneNoCells

MFloat *	m_spatialPlaneAllVars

MFloat *	m_spatialPlaneAllCoord

MInt	m_noProbeLines

MInt *	m_probeLineDirection

MInt *	m_probeLinePeriodic

MFloat **	m_probeLineCoordinates

MInt **	m_probeLineIds

MInt *	m_noProbeLineIds

MFloat **	m_probeLinePositions

MInt *	m_globalNoProbeLineIds

MInt *	m_probeLineOffsets

MInt **	m_noGlobalProbeLineIds

MInt	m_probeLineInterval

MInt	m_probeLineStartTimestep

MInt	m_probeLineStopTimestep

MInt *	m_correlationDirection

MInt *	m_correlationVariableIndex

MFloat **	m_correlationCoordinates

MInt	m_noCorrelationLines

MInt *	m_noCorrelationIds

MInt **	m_correlationIds

MInt **	m_globalCorrelationIds

MFloat **	m_globalCorrelationPositions

MInt *	m_globalNoCorrelationIds

MFloat **	m_correlationPositions

MInt **	m_correlationIndexMapping

MFloat **	m_correlationExchangeVar

MFloat **	m_correlationExchangeVarMean

MFloat **	m_correlationExchangeVarRMS

MFloat **	m_globalCorrelationExchangeVar

MFloat **	m_globalCorrelationExchangeVarMean

MFloat **	m_globalCorrelationExchangeVarRMS

MInt *	m_probeLineAverageDirection

MFloat **	m_probeLineAverageCoordinates

MInt **	m_probeLineAverageCoordinatesSign

MInt **	m_noProbeLineAverageIds

MFloat ***	m_globalProbeLineAverageVars

MInt	m_noProbeLineAverageSteps

MInt **	m_globalNoProbeLineAverageIds

MInt ***	m_probeLineAverageIds

MFloat ***	m_probeLineAveragePositions

MInt	m_noProbeSlices

MInt *	m_probeSliceDir

MFloat *	m_probeSliceCoordinate

MInt **	m_probeSliceIds

MInt *	m_noProbeSliceIds

MFloat **	m_probeSlicePositions

MInt *	m_globalNoProbeSliceIds

MInt *	m_probeSliceOffsets

MInt **	m_noGlobalProbeSliceIds

MString *	m_probeSliceGridNames

MString *	m_sliceAxis

MFloat *	m_sliceIntercept

MBool	m_sliceAiaFileFormat

MBool	m_optimizedSliceIo

MInt *	m_noProbeSliceNoHilbertIds

MInt *	m_noProbeSliceMaxNoHilbertIds

MInt **	m_noProbeSliceHilbertInfo

MInt *	m_noProbeSliceNoContHilbertIds

MInt *	m_noProbeSliceMaxNoContHilbertIds

MInt **	m_noProbeSliceContHilbertInfo

MInt	m_probeSliceInterval

MInt	m_probeSliceStartTimestep

MInt	m_probeSliceStopTimestep

std::vector< MInt >	m_sliceVarIds
	List of slice variables. More...

std::vector< MInt >	m_noSliceVars
	Number of variables for each slice variable. More...

std::vector< std::vector< MString > >	m_sliceVarNames
	List of variable names for each slice variable. More...

MInt	m_noMinLvlIds

MInt *	m_minLvlIds

MInt	m_movePointsToGrid

MBool	m_spatialAveraging

MFloat *	m_arbLinePoints

MInt	m_noArbLines

MInt *	m_noArbLineIds

MFloat *	m_arbLinePointsDistribution

MInt *	m_arbLineHasNewDistribution

MInt *	m_arbLineFull

MInt *	m_moveLinePointsToGrid

MInt *	m_globalNoArbLineIds

MInt **	m_arbLineIds

MInt *	m_arbLineOffsets

MFloat **	m_arbLineCoordinates

MFloat *	m_arbSlicePoints

MInt	m_noArbSlices

MInt *	m_noArbSlicePoints

MFloat *	m_arbSlicePointsDistribution

MInt *	m_globalNoArbSlicePoints

MInt **	m_arbSliceIds

MInt *	m_arbSliceOffsets

MFloat **	m_arbSliceCoordinates

MString	m_postprocessFileName

MFloat **	m_localVars

MBool	m_averageRestart

MInt	m_averageInterval

MInt	m_averageStartTimestep

MInt	m_averageStopTimestep

MInt	m_noAveragedVorticities

MInt	m_noSpeedOfSoundVars

MInt	m_noSourceVars

MBool	m_needVorticity

MString	m_ReStressesAverageFileName

MFloat **	m_probeCoordinates

MString	m_probePath

MInt *	m_probeCellIds

std::ofstream *	m_probeFileStreams

MInt	m_noProbePoints

MInt	m_probeInterval

MInt	m_probeWriteInterval

MString	m_pdFileName

MInt	m_pdStartTimestep

MInt	m_pdStopTimestep

MInt	m_pdRestartInterval

std::vector< MFloat >	m_pdPoints

std::vector< MInt >	m_pdCells

std::vector< MFloat >	m_pdVars

MInt	m_pdNoPoints

MInt	m_sprayComputeInterval

MInt	m_sprayWriteInterval

MInt	m_sprayDataSize

MInt	m_sprayDataStep

MFloat **	m_particleCV

MFloat **	m_particlePen

MFloat **	m_sprayStat

MFloat **	m_sprayTimes

MFloat **	m_injectionData

MFloat **	m_vapourCV

MFloat **	m_vapourPen

MInt	m_LPTSolutionInterval

MBool	m_forceOutput

std::vector< MInt >	m_activeSourceTerms

std::array< MString, ST::totalNoSourceTerms >	s_sourceTermNames

std::set< MInt >	m_activeMeanVars

Protected Attributes inherited from PostProcessingInterface
MInt	m_postprocessingId

Private Types inherited from PostProcessingFv< nDim, SysEqn >
using	SolverType = FvCartesianSolverXD< nDim, SysEqn >

using	Base = PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >

Private Types inherited from PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >
using	PPGrid = CartesianGrid< nDim >

using	Cell = typename PPGrid::Cell

using	GridProxy = typename maia::grid::Proxy< nDim >

using	Data = PostData< nDim >

typedef void(PostProcessing::*	tpost) ()

typedef std::vector< tpost >	tvecpost

Detailed Description

template<MInt nDim, class SysEqn>
class PostProcessingFvLPT< nDim, SysEqn >

Definition at line 26 of file postprocessingfvlpt.h.

Member Typedef Documentation

◆ BaseFv

template<MInt nDim, class SysEqn >

using PostProcessingFvLPT< nDim, SysEqn >::BaseFv = PostProcessing<nDim, PostProcessingFv<nDim, SysEqn> >

Definition at line 36 of file postprocessingfvlpt.h.

◆ BaseLpt

template<MInt nDim, class SysEqn >

using PostProcessingFvLPT< nDim, SysEqn >::BaseLpt = PostProcessing<nDim, PostProcessingLPT<nDim> >

Definition at line 37 of file postprocessingfvlpt.h.

◆ fv

template<MInt nDim, class SysEqn >

using PostProcessingFvLPT< nDim, SysEqn >::fv = FvCartesianSolverXD<nDim, SysEqn>

Definition at line 35 of file postprocessingfvlpt.h.

◆ lptSolver

template<MInt nDim, class SysEqn >

using PostProcessingFvLPT< nDim, SysEqn >::lptSolver = LPT<nDim>

Definition at line 34 of file postprocessingfvlpt.h.

◆ ppFv

template<MInt nDim, class SysEqn >

using PostProcessingFvLPT< nDim, SysEqn >::ppFv = PostProcessingFv<nDim, SysEqn>

Definition at line 38 of file postprocessingfvlpt.h.

◆ ppLpt

template<MInt nDim, class SysEqn >

using PostProcessingFvLPT< nDim, SysEqn >::ppLpt = PostProcessingLPT<nDim>

Definition at line 39 of file postprocessingfvlpt.h.

Constructor & Destructor Documentation

◆ PostProcessingFvLPT()

template<MInt nDim, class SysEqn >

PostProcessingFvLPT< nDim, SysEqn >::PostProcessingFvLPT	(	MInt	postprocessingId_,
		PostData< nDim > *	data,
		fv *	ppFvSolver_,
		lptSolver *	ppLptSolver_
	)

Definition at line 16 of file postprocessingfvlpt.cpp.

  : PostProcessingInterface(postprocessingId_),
    PostProcessingLPT<nDim>(postprocessingId_, data, ppLptSolver_),
    PostProcessingFv<nDim, SysEqn>(postprocessingId_, data, ppFvSolver_) {
  m_ppSolverLpt = ppLptSolver_;
  m_ppSolverFv = ppFvSolver_;
 
  m_fvLPTSpeciesId = 0;
  m_fvLPTSpeciesId = Context::getSolverProperty<MInt>("fvLptSpeciesId", lpt().solverId(), AT_, &m_fvLPTSpeciesId);
}

◆ ~PostProcessingFvLPT()

template<MInt nDim, class SysEqn >

virtual PostProcessingFvLPT< nDim, SysEqn >::~PostProcessingFvLPT ( )

inlinevirtual

Definition at line 44 of file postprocessingfvlpt.h.

44{};

Member Function Documentation

◆ advanceDataStep()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::advanceDataStep ( )

inlineprivate

Definition at line 98 of file postprocessingfvlpt.h.

                         {
    ppFv::advanceDataStep();
    ppLpt::advanceDataStep();
    m_sprayDataStep++;
  };

◆ computeSprayData()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::computeSprayData

overrideprotectedvirtual

Author: Tim Wegmann

Date: Jan. 2021

Reimplemented from PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >.

Definition at line 37 of file postprocessingfvlpt.cpp.

                                                         {
  TRACE();
 
 
  if(globalTimeStep % this->m_sprayComputeInterval == 0) {
    updateData();
    advanceDataStep();
  }
 
 
  if(globalTimeStep % this->m_sprayWriteInterval == 0) {
    writeSprayData();
    resetDataStep();
  }
}

◆ fvSolver()

template<MInt nDim, class SysEqn >

fv & PostProcessingFvLPT< nDim, SysEqn >::fvSolver ( ) const

inlineprivate

Definition at line 74 of file postprocessingfvlpt.h.

74{ return *m_ppSolverFv; }

◆ initLPTSolutionFile()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::initLPTSolutionFile ( )

inlineoverrideprotectedvirtual

Reimplemented from PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >.

Definition at line 66 of file postprocessingfvlpt.h.

66{ ppLpt::initLPTSolutionFile(); }

PostProcessingLPT::initLPTSolutionFile

void initLPTSolutionFile() override

init function for LPT particle solution file

Definition: postprocessinglpt.cpp:459

◆ initPostProcessing()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::initPostProcessing ( )

inlineoverridevirtual

Reimplemented from PostProcessingFv< nDim, SysEqn >.

Definition at line 51 of file postprocessingfvlpt.h.

51{ BaseFv::initPostProcessing(); };

PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >::initPostProcessing

void initPostProcessing() override

Reads all required properties in and prepares for postprocessing.

◆ initSprayData()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::initSprayData ( )

inlineoverridevirtual

Reimplemented from PostProcessingFv< nDim, SysEqn >.

Definition at line 58 of file postprocessingfvlpt.h.

                                {
    ppFv::initSprayData();
    ppLpt::initSprayData();
  }

◆ lpt()

template<MInt nDim, class SysEqn >

lptSolver & PostProcessingFvLPT< nDim, SysEqn >::lpt ( ) const

inlineprivate

Definition at line 73 of file postprocessingfvlpt.h.

73{ return *m_ppSolverLpt; }

◆ mSolver()

template<MInt nDim, class SysEqn >

Solver * PostProcessingFvLPT< nDim, SysEqn >::mSolver ( ) const

inlineoverridevirtual

Implements PostProcessingInterface.

Definition at line 46 of file postprocessingfvlpt.h.

46{ return static_cast<Solver*>(m_ppSolverFv); };

Solver

Parent class of all solvers This class is the base for all solvers. I.e. all solver class (e....

Definition: solver.h:29

◆ postprocessInSolve()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::postprocessInSolve ( const MBool finalTimeStep )

inlineoverridevirtual

Implements PostProcessingInterface.

Definition at line 54 of file postprocessingfvlpt.h.

54{ BaseFv::postprocessInSolve(finalTimeStep); };

PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >::postprocessInSolve

void postprocessInSolve(MBool finalTimeStep) override

Definition: postprocessing.cpp:185

◆ postprocessPostSolve()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::postprocessPostSolve ( )

inlineoverridevirtual

Implements PostProcessingInterface.

Definition at line 53 of file postprocessingfvlpt.h.

53{ BaseFv::postprocessPostSolve(); };

PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >::postprocessPostSolve

void postprocessPostSolve() override

Definition: postprocessing.cpp:206

◆ postprocessPreSolve()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::postprocessPreSolve ( )

inlineoverridevirtual

Implements PostProcessingInterface.

Definition at line 52 of file postprocessingfvlpt.h.

52{ BaseFv::postprocessPreSolve(); };

PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >::postprocessPreSolve

void postprocessPreSolve() override

Definition: postprocessing.cpp:175

◆ postprocessSolution()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::postprocessSolution ( )

inlineoverridevirtual

Implements PostProcessingInterface.

Definition at line 55 of file postprocessingfvlpt.h.

55{ BaseFv::postprocessSolution(); };

PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >::postprocessSolution

void postprocessSolution() override

Definition: postprocessing.cpp:196

◆ resetDataStep()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::resetDataStep ( )

inlineprivate

Definition at line 103 of file postprocessingfvlpt.h.

                       {
    ppFv::resetDataStep();
    ppLpt::resetDataStep();
    m_sprayDataStep = 0;
  };

◆ updateData()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::updateData

private

Author: Tim Wegmann

Date: Jan. 2021

Definition at line 60 of file postprocessingfvlpt.cpp.

                                                   {
  TRACE();
 
 
  m_sprayTimes[m_sprayDataStep][0] = globalTimeStep;
  m_sprayTimes[m_sprayDataStep][1] = fvSolver().physicalTime();
  m_sprayTimes[m_sprayDataStep][2] = fvSolver().crankAngle(fvSolver().physicalTime(), 0);
 
  fvSolver().stopLoadTimer(AT_);
  if(lpt().isActive()) {
    lpt().startLoadTimer(AT_);
    this->particleMass();
    this->particlePenetration();
    this->parcelStatistics();
    lpt().stopLoadTimer(AT_);
  }
 
  fvSolver().startLoadTimer(AT_);
 
  if(fvSolver().isActive()) {
    this->vapourMass(m_fvLPTSpeciesId);
    this->vapourPenetration(lpt().m_spawnCoord);
  }
}

◆ writeLPTSolutionFile()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::writeLPTSolutionFile ( )

inlineoverrideprotectedvirtual

Reimplemented from PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >.

Definition at line 67 of file postprocessingfvlpt.h.

67{ ppLpt::writeLPTSolutionFile(); };

PostProcessingLPT::writeLPTSolutionFile

void writeLPTSolutionFile() override

base function which is called every time-Step in postprocessInSolve

Definition: postprocessinglpt.cpp:475

◆ writeSprayData()

template<MInt nDim, class SysEqn >

void PostProcessingFvLPT< nDim, SysEqn >::writeSprayData

overrideprotectedvirtual

Reimplemented from PostProcessing< nDim, PostProcessingFv< nDim, SysEqn > >.

Definition at line 86 of file postprocessingfvlpt.cpp.

                                                       {
  TRACE();
 
  static MBool firstTime = true;
 
  if(this->m_finalTimeStep && fvSolver().domainId() == 0) {
    cerr << "Writing final postProcessing Data" << endl;
  }
 
  if(firstTime && fvSolver().domainId() == 0) {
    cerr0 << "Writing first postProcessing Data" << endl;
  }
 
  fvSolver().stopLoadTimer(AT_);
 
  const MBool timerEnabled = fvSolver().dlbTimersEnabled();
 
  if(timerEnabled) {
    maia::dlb::g_dlbTimerController.disableAllDlbTimers();
  }
 
  const MBool hasInjection = lpt().m_sprayModel != nullptr ? true : false;
  MInt injDataSize = this->getInjectionData();
  MPI_Allreduce(MPI_IN_PLACE, &injDataSize, 1, MPI_INT, MPI_MAX, fvSolver().mpiComm(), AT_, "INPLACE", "injDataSize");
  const MInt dataSize = (this->m_finalTimeStep || firstTime) ? m_sprayDataStep : m_sprayDataSize;
  const MInt vapourPenSize = nDim == 3 ? 16 : 12;
 
  // rank with injector writes data
  const MInt injectorCellId = lpt().isActive() ? lpt().injectorCellId() : -1;
 
  if(hasInjection && injDataSize != dataSize) {
    cerr0 << "Spray-PP data differs: " << injDataSize << " inj-Size " << dataSize << " sprayData-Size" << endl;
  }
 
  if(dataSize > 0) {
    // exchange data:
    MPI_Allreduce(MPI_IN_PLACE, &m_particleCV[0][0], m_sprayDataSize * (2 + nDim), MPI_DOUBLE, MPI_SUM,
                  fvSolver().mpiComm(), AT_, "INPLACE", "m_particleCV");
 
    MPI_Allreduce(MPI_IN_PLACE, &m_sprayStat[0][0], m_sprayDataSize * 5, MPI_DOUBLE, MPI_SUM, fvSolver().mpiComm(), AT_,
                  "INPLACE", "m_sprayStat");
 
    MPI_Allreduce(MPI_IN_PLACE, &m_particlePen[0][0], m_sprayDataSize * 19, MPI_DOUBLE, MPI_MAX, fvSolver().mpiComm(),
                  AT_, "INPLACE", "m_particlePen");
 
    MPI_Allreduce(MPI_IN_PLACE, &m_vapourCV[0][0], m_sprayDataSize * (5 + 2 * nDim), MPI_DOUBLE, MPI_SUM,
                  fvSolver().mpiComm(), AT_, "INPLACE", "m_vapourCV");
 
    MPI_Allreduce(MPI_IN_PLACE, &m_vapourPen[0][0], m_sprayDataSize * vapourPenSize, MPI_DOUBLE, MPI_MAX,
                  fvSolver().mpiComm(), AT_, "INPLACE", "m_vapourPen");
  }
 
  if(injectorCellId >= 0) {
 
    ofstream ofl;
 
    // backup a few older files
    if(firstTime) {
      // file with conservation variables
      MString fileNameCheck = lpt().outputDir() + "sprayCons";
      if(fileExists(fileNameCheck)) {
        const MString fileName = lpt().outputDir() + "sprayCons_" + to_string(lpt().m_restartTimeStep);
        if(!fileExists(fileName)) copyFile(fileNameCheck, fileName);
      }
      if(!lpt().m_restartFile) {
        ofl.open(fileNameCheck, ios_base::out | ios_base::trunc);
        ofl << "#0:ts 1:time 2:CAD 3:pMass 4:pEnergy 5-7:pMomentum 8:vMass 9:aMass 10:gMass 11:gEnergy 12:gMomX "
               "13:gMomY 14:gMomZ";
        ofl << endl;
        ofl.close();
      }
 
      // file with liquid penetration lengths
      fileNameCheck = lpt().outputDir() + "liquidPen";
      if(fileExists(fileNameCheck)) {
        const MString fileName = lpt().outputDir() + "liquidPen_" + to_string(lpt().m_restartTimeStep);
        if(!fileExists(fileName)) copyFile(fileNameCheck, fileName);
      }
      if(!lpt().m_restartFile) {
        ofl.open(fileNameCheck, ios_base::out | ios_base::trunc);
        ofl << "#0:ts 1:time 2:CAD 3-5:liquidPen"
               "6-8:liquidWidth";
        ofl << endl;
        ofl.close();
      }
 
      // file with vapour penetration lengths
      fileNameCheck = lpt().outputDir() + "vapourPen";
      if(fileExists(fileNameCheck)) {
        const MString fileName = lpt().outputDir() + "vapourPen_" + to_string(lpt().m_restartTimeStep);
        if(!fileExists(fileName)) copyFile(fileNameCheck, fileName);
      }
      if(!lpt().m_restartFile) {
        ofl.open(fileNameCheck, ios_base::out | ios_base::trunc);
        ofl << "#0:ts 1:time 2:CAD 3-5:v01Pen 6-8:v1Pen 9-11:v5Pen";
        ofl << endl;
        ofl.close();
      }
 
      // file with spray statistics
      fileNameCheck = lpt().outputDir() + "sprayStats";
      if(fileExists(fileNameCheck)) {
        const MString fileName = lpt().outputDir() + "sprayStats_" + to_string(lpt().m_restartTimeStep);
        if(!fileExists(fileName)) copyFile(fileNameCheck, fileName);
      }
 
      if(!lpt().m_restartFile) {
        ofl.open(fileNameCheck, ios_base::out | ios_base::trunc);
        ofl << "#0:ts 1:time 2:CAD 3:noParcles 4:noDroplets 5:SMS 6:avgDiameter";
        ofl << endl;
        ofl.close();
      }
 
      // file with injection statistics
      if(hasInjection) {
        fileNameCheck = lpt().outputDir() + "injStats";
        if(fileExists(fileNameCheck)) {
          const MString fileName = lpt().outputDir() + "injStats_" + to_string(lpt().m_restartTimeStep);
          if(!fileExists(fileName)) copyFile(fileNameCheck, fileName);
        }
        if(!lpt().m_restartFile) {
          ofl.open(fileNameCheck, ios_base::out | ios_base::trunc);
          ofl << "#0:ts 1:timeSinceSOI 2:injProgress 3:rampFactor 4:noInjPart 5:noInjDrop 6:injDiameter 7:injMass "
                 "8:injXMom 9:injYMom 10:injZMom 11:injEnergy";
          ofl << endl;
          ofl.close();
        }
 
        // file with sum of conservative variables from injection
        fileNameCheck = lpt().outputDir() + "injSums";
        if(fileExists(fileNameCheck)) {
          const MString fileName = lpt().outputDir() + "injSums_" + to_string(lpt().m_restartTimeStep);
          if(!fileExists(fileName)) copyFile(fileNameCheck, fileName);
        }
        if(!lpt().m_restartFile) {
          ofl.open(fileNameCheck, ios_base::out | ios_base::trunc);
          ofl << "#0:ts 1:injMass 2-4:injMomentum 5:injEnergy 6:evapMass";
          ofl << endl;
          ofl.close();
        }
      }
    }
 
 
    MString fileName = lpt().outputDir() + "sprayCons";
    ofl.open(fileName, ios_base::out | ios_base::app);
 
    if(ofl.is_open() && ofl.good()) {
      for(MInt i = 0; i < dataSize; i++) {
        ofl << m_sprayTimes[i][0] << setprecision(7) << " " << m_sprayTimes[i][1] << " " << m_sprayTimes[i][2] << " "
            << setprecision(12);
        for(MInt v = 0; v < 2 + nDim; v++) {
          ofl << m_particleCV[i][v] << " ";
        }
        for(MInt v = 0; v < 5 + 2 * nDim; v++) {
          ofl << m_vapourCV[i][v] << " ";
        }
        ofl << endl;
      }
    }
    ofl.close();
 
    fileName = lpt().outputDir() + "liquidPen";
    ofl.open(fileName, ios_base::out | ios_base::app);
    if(ofl.is_open() && ofl.good()) {
      for(MInt i = 0; i < dataSize; i++) {
        ofl << m_sprayTimes[i][0] << setprecision(7) << " " << m_sprayTimes[i][1] << " " << setprecision(7)
            << m_sprayTimes[i][2] << " " << setprecision(8);
        for(MInt j = 0; j < 19; j++) {
          ofl << m_particlePen[i][j] << " ";
        }
        ofl << endl;
      }
    }
    ofl.close();
 
    fileName = lpt().outputDir() + "vapourPen";
    ofl.open(fileName, ios_base::out | ios_base::app);
    if(ofl.is_open() && ofl.good()) {
      for(MInt i = 0; i < dataSize; i++) {
        ofl << m_sprayTimes[i][0] << setprecision(7) << " " << m_sprayTimes[i][1] << " " << setprecision(7)
            << m_sprayTimes[i][2] << " " << setprecision(8);
        for(MInt j = 0; j < vapourPenSize; j++) {
          ofl << m_vapourPen[i][j] << " ";
        }
        ofl << endl;
      }
    }
    ofl.close();
 
    fileName = lpt().outputDir() + "sprayStats";
    ofl.open(fileName, ios_base::out | ios_base::app);
    if(ofl.is_open() && ofl.good()) {
      for(MInt i = 0; i < dataSize; i++) {
        ofl << m_sprayTimes[i][0] << setprecision(7) << " " << m_sprayTimes[i][1] << " " << m_sprayTimes[i][2] << " "
            << setprecision(12);
        for(MInt j = 0; j < 5; j++) {
          ofl << m_sprayStat[i][j] << " ";
        }
        ofl << endl;
      }
    }
    ofl.close();
 
    if(hasInjection) {
      fileName = lpt().outputDir() + "injStats";
      ofl.open(fileName, ios_base::out | ios_base::app);
      if(ofl.is_open() && ofl.good()) {
        for(MInt i = 0; i < injDataSize; i++) {
          ofl << setprecision(12);
          for(MInt j = 0; j < 15; j++) {
            ofl << m_injectionData[i][j] << " ";
          }
          ofl << endl;
        }
      }
      ofl.close();
 
      fileName = lpt().outputDir() + "injSums";
      ofl.open(fileName, ios_base::out | ios_base::app);
      if(ofl.is_open() && ofl.good()) {
        for(MInt i = 0; i < injDataSize; i++) {
          ofl << m_injectionData[i][0] << " " << setprecision(12);
          for(MInt j = 0; j < nDim + 5; j++) {
            ofl << m_cInjData[i][j] << " ";
          }
          ofl << endl;
        }
      }
      ofl.close();
    }
  }
  firstTime = false;
 
  // reset data
  for(MInt i = 0; i < m_sprayDataSize; i++) {
    for(MInt j = 0; j < 5; j++) {
      m_sprayStat[i][j] = 0;
    }
    for(MInt j = 0; j < 2 + nDim; j++) {
      m_particleCV[i][j] = 0;
    }
    for(MInt j = 0; j < 19; j++) {
      m_particlePen[i][j] = 0;
    }
    if(hasInjection) {
      for(MInt j = 0; j < 15; j++) {
        m_injectionData[i][j] = 0;
      }
    }
  }
 
  if(timerEnabled) {
    maia::dlb::g_dlbTimerController.enableAllDlbTimers();
  }
  fvSolver().startLoadTimer(AT_);
}

Friends And Related Function Documentation

◆ LPT

template<MInt nDim, class SysEqn >

template<MInt nDim_>

friend class LPT

friend

Definition at line 31 of file postprocessingfvlpt.h.

◆ PostProcessing

template<MInt nDim, class SysEqn >

template<MInt nDim_, class ppType >

friend class PostProcessing

friend

Definition at line 29 of file postprocessingfvlpt.h.

Member Data Documentation

◆ m_cInjData

template<MInt nDim, class SysEqn >

MFloat** PostProcessingLPT< nDim >::m_cInjData

private

Definition at line 71 of file postprocessinglpt.h.

◆ m_fvLPTSpeciesId

template<MInt nDim, class SysEqn >

MInt PostProcessingFvLPT< nDim, SysEqn >::m_fvLPTSpeciesId = 0

private

Definition at line 76 of file postprocessingfvlpt.h.

◆ m_ppSolverFv

template<MInt nDim, class SysEqn >

fv* PostProcessingFvLPT< nDim, SysEqn >::m_ppSolverFv

private

Definition at line 71 of file postprocessingfvlpt.h.

◆ m_ppSolverLpt

template<MInt nDim, class SysEqn >

lptSolver* PostProcessingFvLPT< nDim, SysEqn >::m_ppSolverLpt

private

Definition at line 70 of file postprocessingfvlpt.h.

◆ m_sprayDataStep

template<MInt nDim, class SysEqn >

MInt PostProcessingFvLPT< nDim, SysEqn >::m_sprayDataStep = 0

private

Definition at line 96 of file postprocessingfvlpt.h.

◆ s_ppType

template<MInt nDim, class SysEqn >

MString PostProcessingFvLPT< nDim, SysEqn >::s_ppType

static

Definition at line 48 of file postprocessingfvlpt.h.

The documentation for this class was generated from the following files:

/home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/POST/postprocessingfvlpt.h
/home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/POST/postprocessingfvlpt.cpp

Public Types

Public Member Functions

Static Public Attributes

Protected Member Functions

Private Member Functions

Private Attributes

Friends

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ BaseFv

◆ BaseLpt

◆ fv

◆ lptSolver

◆ ppFv

◆ ppLpt

Constructor & Destructor Documentation

◆ PostProcessingFvLPT()

◆ ~PostProcessingFvLPT()

Member Function Documentation

◆ advanceDataStep()

◆ computeSprayData()

◆ fvSolver()

◆ initLPTSolutionFile()

◆ initPostProcessing()

◆ initSprayData()

◆ lpt()

◆ mSolver()

◆ postprocessInSolve()

◆ postprocessPostSolve()

◆ postprocessPreSolve()

◆ postprocessSolution()

◆ resetDataStep()

◆ updateData()

◆ writeLPTSolutionFile()

◆ writeSprayData()

Friends And Related Function Documentation

◆ LPT

◆ PostProcessing

Member Data Documentation

◆ m_cInjData

◆ m_fvLPTSpeciesId

◆ m_ppSolverFv

◆ m_ppSolverLpt

◆ m_sprayDataStep

◆ s_ppType