setup_profiles.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2009 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/Base/Formats/PSRFITS/setup_profiles.h
 
#ifndef __Pulsar_setup_profiles_h
#define __Pulsar_setup_profiles_h
 
template<typename I, typename P>
void setup_profiles_dat (I subint, P& profiles)
{
  const unsigned npol = subint->get_npol ();
  const unsigned nchan = subint->get_nchan ();
 
  profiles.resize (npol * nchan);
 
  for (unsigned ipol=0; ipol<npol; ipol++)
    for (unsigned ichan=0; ichan<nchan; ichan++)
      profiles[ipol*nchan + ichan] = subint->get_Profile (ipol, ichan);
}
 
template<class E, typename I, typename P>
void setup_profiles (I subint, P& profiles)
{
  E* ext = subint->get_Profile(0,0)->template get<E>();
  if (!ext)
    throw Error (InvalidState, "setup_profiles<Extension>",
                 "first profile is missing required Extension");
 
  const unsigned naux = ext->get_size();
  const unsigned nchan = subint->get_nchan ();
 
  profiles.resize (naux * nchan);
 
  for (unsigned ichan=0; ichan<nchan; ichan++)
  {
    ext = subint->get_Profile(0,ichan)->template get<E>();
    if (!ext)
      throw Error (InvalidState, "setup_profiles<Extension>",
                   "profile[%u] is missing required Extension", ichan);
 
    for (unsigned iaux=0; iaux<naux; iaux++)
      profiles[iaux*nchan + ichan] = ext->get_Profile (iaux);
  }
}
 
#endif