AcaPostProcessingSPL< nDim > Class Template Reference

ACA post processing class for sound pressure level calculation. More...

#include <acapostprocessing.hpp>

Inheritance diagram for AcaPostProcessingSPL< nDim >:

Collaboration diagram for AcaPostProcessingSPL< nDim >:

Public Member Functions
void	init_ () override
void	calc (const MInt observerId, const MFloat *const data, const MFloat *const dataComplex=nullptr) override
void	finish () override
	AcaPostProcessing (const MPI_Comm comm)
Public Member Functions inherited from AcaPostProcessing
	AcaPostProcessing (const MPI_Comm comm)
virtual	~AcaPostProcessing ()=default
virtual void	init_ ()=0
void	init (const MInt noObservers, const MInt noGlobalObservers, const MInt offsetObserver, const MInt noSamples, const MFloat *const coords, const MFloat *const frequencies, const MString outPath)
virtual void	calc (const MInt observerId, const MFloat *const data, const MFloat *const dataComplex=nullptr)=0
virtual void	finish ()=0

Private Attributes
MString	m_fileNameIoParallel
MFloat	m_pRefSq = -1.0
const MFloat	m_gamma = 1.4
std::vector< MFloat >	m_spl
	isentropic exponent More...

Additional Inherited Members
Protected Attributes inherited from AcaPostProcessing
MPI_Comm	m_mpiComm
MInt	m_rank
MInt	m_noObservers = -1
MInt	m_noGlobalObservers = -1
MInt	m_offsetObserver = 0
MInt	m_noSamples = -1
MString	m_outPath
const MFloat *	m_coords = nullptr
const MFloat *	m_frequencies = nullptr

Detailed Description

template<MInt nDim>
class AcaPostProcessingSPL< nDim >

Definition at line 208 of file acapostprocessing.hpp.

Member Function Documentation

AcaPostProcessing()

template<MInt nDim>

AcaPostProcessing::AcaPostProcessing ( const MPI_Comm  comm )

inline

Definition at line 25 of file acapostprocessing.h.

: m_mpiComm(comm) { MPI_Comm_rank(m_mpiComm, &m_rank); }

calc()

template<MInt nDim>

void AcaPostProcessingSPL< nDim >::calc ( const MInt  observerId, const MFloat *const  data, const MFloat *const  dataComplex = nullptr  )

overridevirtual

Implements AcaPostProcessing.

Definition at line 246 of file acapostprocessing.hpp.

                                                                       {
  // Calculation of sound pressure level (pressure at each frequency)
  const MInt offset = observerId * m_noSamples;
  for(MInt t = 1; t < m_noSamples; t++) { // Skip frequency 0
    const MFloat real = dataComplex[t * 2 + 0];
    const MFloat imag = dataComplex[t * 2 + 1];
    // p times p complex conjugate = p_re_sq + p_imag_sq
    const MFloat pSq = POW2(real) + POW2(imag);
    // Sound pressure for each observer point for each frequency
    // Ref.:  Mendez et al. (2013): https://doi.org/10.1260%2F1475-472X.12.1-2.1
    //        or Beranek and Ver (1992); Delfs, Basics of Aeroacoustics (2016); ...
    // Remind: 20 * log(a) = 10 * log(a^2), and pSq is only halved value because of FFT
    m_spl[offset + t] = 10 * log10(2 * pSq / m_pRefSq);
  }
}

finish()

template<MInt nDim>

void AcaPostProcessingSPL< nDim >::finish

overridevirtual

Implements AcaPostProcessing.

Definition at line 264 of file acapostprocessing.hpp.

                                        {
  using namespace maia::parallel_io;
  ParallelIo file(m_fileNameIoParallel, PIO_REPLACE, m_mpiComm);
  ParallelIo::size_type dimSizesCoordinates[] = {m_noGlobalObservers, nDim};
  ParallelIo::size_type dimSizesSpl[] = {m_noGlobalObservers, m_noSamples};
  file.defineArray(PIO_FLOAT, "coordinates", 2, &dimSizesCoordinates[0]);
  file.defineArray(PIO_FLOAT, "frequency", m_noSamples);
  file.defineArray(PIO_FLOAT, "SPL", 2, &dimSizesSpl[0]);
 
  const MBool isRoot = (m_rank == 0);
  file.setOffset((isRoot) ? m_noGlobalObservers : 0, 0, 2);
  file.writeArray(m_coords, "coordinates");
  file.setOffset((isRoot) ? m_noSamples : 0, 0);
  file.writeArray(m_frequencies, "frequency");
  file.setOffset(m_noObservers, m_offsetObserver, 2);
  file.writeArray(m_spl.data(), "SPL");
}

init_()

template<MInt nDim>

void AcaPostProcessingSPL< nDim >::init_

overridevirtual

Implements AcaPostProcessing.

Definition at line 226 of file acapostprocessing.hpp.

Member Data Documentation

m_fileNameIoParallel

template<MInt nDim>

MString AcaPostProcessingSPL< nDim >::m_fileNameIoParallel

private

Definition at line 210 of file acapostprocessing.hpp.

m_gamma

template<MInt nDim>

const MFloat AcaPostProcessingSPL< nDim >::m_gamma = 1.4

private

Definition at line 212 of file acapostprocessing.hpp.

m_pRefSq

template<MInt nDim>

MFloat AcaPostProcessingSPL< nDim >::m_pRefSq = -1.0

private

Definition at line 211 of file acapostprocessing.hpp.

m_spl

template<MInt nDim>

std::vector<MFloat> AcaPostProcessingSPL< nDim >::m_spl

private

Definition at line 213 of file acapostprocessing.hpp.

---

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

• /home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/ACA/acapostprocessing.hpp