AcaPostProcessingPABS< nDim > Class Template Reference

ACA post processing class for absoulte pressure calculation. More...

#include <acapostprocessing.hpp>

Inheritance diagram for AcaPostProcessingPABS< nDim >:

Collaboration diagram for AcaPostProcessingPABS< 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
std::vector< MFloat >	m_pabs

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 AcaPostProcessingPABS< nDim >

Definition at line 68 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 AcaPostProcessingPABS< nDim >::calc ( const MInt  observerId, const MFloat *const  data, const MFloat *const  dataComplex = nullptr  )

overridevirtual

Implements AcaPostProcessing.

Definition at line 87 of file acapostprocessing.hpp.

                                                                        {
  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];
    // product of p * p_complex_conjugate = p_re_sq + p_im_sq = [abs(p)]^2
    const MFloat p_magnSq = POW2(real) + POW2(imag);
    // Absolute pressure
    m_pabs[offset + t] = 2 * sqrt(p_magnSq);
  }
}

finish()

template<MInt nDim>

void AcaPostProcessingPABS< nDim >::finish

overridevirtual

Implements AcaPostProcessing.

Definition at line 101 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 dimSizes[] = {m_noGlobalObservers, m_noSamples};
  file.defineArray(PIO_FLOAT, "coordinates", 2, &dimSizesCoordinates[0]);
  file.defineArray(PIO_FLOAT, "pressABS", 2, &dimSizes[0]);
 
  const MBool isRoot = (m_rank == 0);
  file.setOffset((isRoot) ? m_noGlobalObservers : 0, 0, 2);
  file.writeArray(m_coords, "coordinates");
  ScratchSpace<MFloat> pressABS(m_noObservers, m_noSamples, AT_, "pressABS");
  file.setOffset(m_noObservers, m_offsetObserver, 2);
  for(MInt obs = 0; obs < m_noObservers; obs++) {
    const MInt offset = obs * m_noSamples;
    for(MInt t = 0; t < m_noSamples; t++) {
      pressABS(obs, t) = m_pabs[offset + t];
    }
  }
  file.writeArray(pressABS.data(), "pressABS");
}

init_()

template<MInt nDim>

void AcaPostProcessingPABS< nDim >::init_

overridevirtual

Implements AcaPostProcessing.

Definition at line 81 of file acapostprocessing.hpp.

                                        {
  m_fileNameIoParallel = m_outPath + "pressABS" + ParallelIo::fileExt();
  m_pabs.resize(m_noObservers * m_noSamples, 0.0);
}

Member Data Documentation

m_fileNameIoParallel

template<MInt nDim>

MString AcaPostProcessingPABS< nDim >::m_fileNameIoParallel

private

Definition at line 70 of file acapostprocessing.hpp.

m_pabs

template<MInt nDim>

std::vector<MFloat> AcaPostProcessingPABS< nDim >::m_pabs

private

Definition at line 71 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