read()

void datrw::ihpmostream::read ( std::istream &  is, const bool &  verbose = false  )

private

read one data file

Definition at line 101 of file ihpmostream.cc.

References DATRW_assert, datrw::hpmo::dt(), Mdatafile, Mindex, Mnblocks, Mnextblock, Mnextchannel, datrw::hpmo::MinuteBlock::Mtime, datrw::hpmo::nminutes, datrw::hpmo::nsamples, datrw::hpmo::qualityreports(), datrw::hpmo::readfile(), and datrw::idatstream::setfilefree().

Referenced by ihpmostream().

  {
    // simply read the thing
    datrw::hpmo::Tvecofblocks blocks(datrw::hpmo::readfile(is));
    if (verbose)
    {
      std::cout << "(ihpmostream::read): read "
        << blocks.size() << " blocks of data" << std::endl;
    }
    int nblocksread=blocks.size();
    DATRW_assert((nblocksread <= datrw::hpmo::nminutes),
                   "ERROR (ihpmostream::read): " 
                   "illegal number of minute blocks in file!" "\n"
                   "  Use ihpmostream only together with original HPMO data");
    // extract the stuff now
    if (verbose) { std::cout << "(ihpmostream::read): copy block #"; }
    libtime::TRelativeTime length(datrw::hpmo::dt()*datrw::hpmo::nsamples);
    for (int iblock=0; iblock<nblocksread; iblock++)
    {
      this->Mdatafile[iblock]=blocks[iblock];
      datrw::hpmo::MinuteBlock& block=this->Mdatafile[iblock];
      if (iblock > 0)
      {
        datrw::hpmo::MinuteBlock& prevblock=this->Mdatafile[iblock-1];
        if (prevblock.Mtime+length != block.Mtime)
        {
          Mindex[Mnblocks]=iblock;
          ++Mnblocks;
          if (verbose) { 
            std::cout << std::endl 
              << "(ihpmostream::read): found gap" << std::endl; 
            std::cout << "  between block # " << iblock-1 << " @ "
              << prevblock.Mtime.timestring() << std::endl;
            std::cout << "      and block # " << iblock << " @ "
              << block.Mtime.timestring() << std::endl;
            std::cout << "(ihpmostream::read): copy block #"; 
          }
        }
      }
      else
      {
        Mindex[Mnblocks]=0;
        ++Mnblocks;
      }
      if (verbose) { std::cout << " " << iblock; }
    }
    if (verbose) 
    { 
      std::cout << std::endl; 
      std::cout << "(ihpmostream::read): extracted "
        << Mnblocks << " contiguous data block(s)" << std::endl; 
    }
    // marks last block
    Mindex[Mnblocks]=nblocksread;
    // init reading pointers
    Mnextblock=0;
    Mnextchannel=1;
    // prepare file FREE block
    std::ostringstream oss;
    oss << nblocksread;
    sff::FREE filefree;
    filefree.append("ihpmostream read " + oss.str() + " blocks from file");
    oss.str("");
    oss << Mnblocks;
    filefree.append("found " + oss.str() + " contiguous sets of blocks");
    filefree.append(qualityreports(Mdatafile, nblocksread));
    this->setfilefree(filefree);
  } // void ihpmostream::read(std::istream& is, const bool& verbose)

Here is the call graph for this function:

Here is the caller graph for this function: