checkconsistency()

void stfinv::STFBaseEngine::checkconsistency ( ) const

private

Check consistency of data members.

Check consistency of time series data passed to the STFBaseEngine with respect to

1. number of samples
2. sampling interval

Definition at line 159 of file stfinvbase.cc.

References CWaveformHeader::dt, dt(), Mpairs, Mstf, Mtriples, CWaveformHeader::n, CTripleHeader::sampling, stfinv::Waveform::sampling, stfinv::Waveform::series, and STFINV_assert.

Referenced by initialize().

  {
    const unsigned int& n=Mstf.sampling.n;
    const double& dt=Mstf.sampling.dt;
    STFINV_assert(Mstf.sampling.n==Mstf.series.size(),
                  "inconsistent number of samples");
    const double tolerance=1.e-4;

    {
      stfinv::Tvectoroftriples::const_iterator I=Mtriples.begin();
      while (I!=Mtriples.end())
      {
        CTripleHeader header=I->header;
        STFINV_assert(header.sampling.n==n,
                      "inconsistent number of samples");
        STFINV_assert(header.sampling.n==I->data.size(),
                      "inconsistent number of samples");
        STFINV_assert(header.sampling.n==I->synthetics.size(),
                      "inconsistent number of samples");
        STFINV_assert(header.sampling.n==I->convolvedsynthetics.size(),
                      "inconsistent number of samples");
        STFINV_assert(std::abs(1.-header.sampling.dt/dt)<tolerance,
                      "inconsistent values of sampling interval");
        ++I;
      } // while (I!=Mtriples.end())
    }

    {
      stfinv::Tvectorofpairs::const_iterator I=Mpairs.begin();
      while (I!=Mpairs.end())
      {
        CWaveformHeader sampling=I->sampling;
        STFINV_assert(sampling.n==n,
                      "inconsistent number of samples");
        STFINV_assert(sampling.n==I->synthetics.size(),
                      "inconsistent number of samples");
        STFINV_assert(sampling.n==I->convolvedsynthetics.size(),
                      "inconsistent number of samples");
        STFINV_assert(std::abs(1.-sampling.dt/dt)<tolerance,
                      "inconsistent values of sampling interval");
        ++I;
      } // while (I!=Mpairs.end())
    }
  }

Here is the call graph for this function:

Here is the caller graph for this function: