All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Defines
Classes | Public Member Functions
QGpCoreWave::ModalStorage Class Reference

Storage for modal computation. More...

#include <ModalStorage.h>

Inheritance diagram for QGpCoreWave::ModalStorage:
QGpCoreWave::Dispersion QGpCoreWave::Ellipticity

List of all members.

Classes

class  Iterator

Public Member Functions

void abs ()
Curve< Point2D > curve (int iMode) const
void exp10 ()
ModalRefine & lastRefine ()
void log10 ()
 ModalStorage ()
 ModalStorage (int nModes, const QVector< double > *x)
RealValue * mode (int index)
const RealValue * mode (int index) const
int modeCount () const
ModalStorage & operator= (const ModalStorage &o)
ModalRefine & refineAdd (double x)
RefineIterator refineBegin () const
void refineClear ()
int refineCount () const
RefineIterator refineEnd () const
void refineSort ()
void set (int nModes, const QVector< double > *x)
void toStream (QTextStream &s, int modeIndex=-1) const
void writeReport (QDataStream &s) const
const QVector< double > * x () const
double x (int ix) const
int xCount () const
 ~ModalStorage ()

Detailed Description

Storage for modal computation.

Storage for modal computation such as dispersion curve, ellipticity or spac

Constructor & Destructor Documentation

QGpCoreWave::ModalStorage::ModalStorage ( )

Init a empty storage. Use set() before using this object.

References TRACE.

Referenced by QGpCoreWave::Dispersion::deltaK().

{
  TRACE;
  _nModes=0;
  _nX=0;
  _values=0;
  _x=0;
}
QGpCoreWave::ModalStorage::ModalStorage ( int  nModes,
const QVector< double > *  x 
) [inline]

Init storage with nModes and with the list of x o (not owned by this object) 

{set(nModes, x);}
QGpCoreWave::ModalStorage::~ModalStorage ( ) [inline]

Delete stored values but not the list of x provided in constructor's arguments 

{delete [] _values;}

Member Function Documentation

void QGpCoreWave::ModalStorage::abs ( )

Transforms all values into their absolute value

References QGpCoreTools::RealValue::abs().

Referenced by QGpCoreWave::DispersionFactory::calculate(), EllipticityReader::parse(), and QGpCoreWave::EllipticityFactory::peakMisfit().

{
  int nValues=_nModes * _nX;
  for(int i=0; i<nValues;i++) {
    _values[i].abs();
  }
  for(RefineList::iterator it=_refines.begin(); it!=_refines.end(); it++) {
    for(int i=0; i<_nModes; i++) {
      it->setValue(i, fabs(it->value(i).value()));
    }
  }
}
Curve< Point2D > QGpCoreWave::ModalStorage::curve ( int  iMode) const

Return a curve for mode iMode, mixing refines and original samples, removing invalid values

References QGpCoreTools::Curve< pointType >::append(), QGpCoreWave::ModalStorage::Iterator::atEnd(), QGpCoreTools::Value< numberType >::isValid(), and QGpCoreTools::Value< numberType >::value().

{
  double xFactor=0.5/M_PI;
  Curve<Point2D> points;
  for(Iterator it(this, iMode); !it.atEnd(); ++it) {
    const RealValue & val=it.value();
    if(val.isValid())
      points.append(Point2D( it.x() * xFactor, val.value()) );
  }
  return points;
}
void QGpCoreWave::ModalStorage::exp10 ( )

Transforms all values into their exp10

References QGpCoreTools::RealValue::exp10().

{
  int nValues=_nModes * _nX;
  for(int i=0; i<nValues;i++) {
    _values[i].exp10();
  }
  for(RefineList::iterator it=_refines.begin(); it!=_refines.end(); it++) {
    for(int i=0; i<_nModes; i++) {
      it->setValue(i, ::pow(10.0, it->value(i).value()));
    }
  }
}
ModalRefine& QGpCoreWave::ModalStorage::lastRefine ( ) [inline] 
{return _refines.last();}
void QGpCoreWave::ModalStorage::log10 ( )

Transforms all values into their log10.

References QGpCoreTools::RealValue::log10().

{
  int nValues=_nModes * _nX;
  for(int i=0; i<nValues;i++) {
    _values[i].log10();
  }
  for(RefineList::iterator it=_refines.begin(); it!=_refines.end(); it++) {
    for(int i=0; i<_nModes; i++) {
      it->setValue(i, ::log10(it->value(i).value()));
    }
  }
}
RealValue* QGpCoreWave::ModalStorage::mode ( int  index) [inline]

Referenced by QGpCoreWave::AutocorrFactory::calculateHorizontal(), QGpCoreWave::AutocorrFactory::calculateVertical(), QGpCoreWave::Dispersion::deltaK(), dispersion_curve_love_(), dispersion_curve_rayleigh_(), QGpCoreWave::ModalStorage::Iterator::Iterator(), QGpCoreWave::EllipticityFactory::mode(), QGpCoreWave::DispersionFactory::mode(), QGpCoreWave::AutocorrFactory::mode(), QGpCoreWave::Ellipticity::peakMisfit(), QGpCoreWave::Ellipticity::peaks(), Spac3CForward::radial(), QGpCoreWave::Dispersion::refine(), QGpCoreWave::Dispersion::setGroupSlowness(), Spac3CForward::setModel(), Spac3CForward::transverse(), and Spac3CForward::vertical().

{return _values+index*_nX;}
const RealValue* QGpCoreWave::ModalStorage::mode ( int  index) const [inline] 
{return _values+index*_nX;}
int QGpCoreWave::ModalStorage::modeCount ( ) const [inline]

Referenced by QGpCoreWave::AutocorrFactory::calculateHorizontal(), QGpCoreWave::AutocorrFactory::calculateVertical(), QGpCoreWave::Dispersion::deltaK(), QGpCoreWave::Dispersion::setGroupSlowness(), and QGpCoreWave::DispersionFactory::validate().

{return _nModes;}
ModalStorage & QGpCoreWave::ModalStorage::operator= ( const ModalStorage &  o)

Copy operator, only list of x remains shared.

References TRACE.

{
  TRACE;
  _nModes=o._nModes;
  _nX=o._nX;
  _x=o._x;
  int nValues=_nModes*_nX;
  delete [] _values;
  _values=new RealValue[nValues];
  memcpy (_values, o._values, nValues*sizeof(RealValue));
  return *this;
}
ModalRefine& QGpCoreWave::ModalStorage::refineAdd ( double  x) [inline]

References QGpCoreWave::ModalRefine::setValid().

Referenced by QGpCoreWave::Dispersion::refine().

                                   {
    _refines.append(ModalRefine(_nModes, x));
    _refines.last().setValid(false);
    return _refines.last();
  }
RefineIterator QGpCoreWave::ModalStorage::refineBegin ( ) const [inline]

Referenced by QGpCoreWave::ModalStorage::Iterator::Iterator().

{return _refines.begin();}
void QGpCoreWave::ModalStorage::refineClear ( ) [inline] 
{_refines.clear();}
int QGpCoreWave::ModalStorage::refineCount ( ) const [inline] 
{return _refines.count();}
RefineIterator QGpCoreWave::ModalStorage::refineEnd ( ) const [inline]

Referenced by QGpCoreWave::ModalStorage::Iterator::Iterator().

{return _refines.end();}
void QGpCoreWave::ModalStorage::refineSort ( ) [inline]

Referenced by QGpCoreWave::Ellipticity::peakMisfit().

{qSort(_refines);}
void QGpCoreWave::ModalStorage::set ( int  nModes,
const QVector< double > *  x 
)

Init storage with nModes and with the list of x o (not owned by this object)

References TRACE, and x().

{
  TRACE;
  _nModes=nModes;
  _nX=x->count();
  _values=new RealValue [ _nModes * _nX ];
  _x=x;
}
void QGpCoreWave::ModalStorage::toStream ( QTextStream &  s,
int  modeIndex = -1 
) const

Output to stream s, mixing refines and original samples, removing invalid values. If modeIndex is negative, all modes are ouput.

References QGpCoreWave::ModalStorage::Iterator::atEnd(), QGpCoreTools::Value< numberType >::isValid(), TRACE, and QGpCoreTools::Value< numberType >::value().

Referenced by EllipticityReader::parse(), DispersionReader::parse(), and SpacReader::parse().

{
  TRACE;
  static const QString fmt("%1 %2\n");
  double xFactor=0.5/M_PI;
  for (int im=0; im<_nModes; im++) {
    if(modeIndex<0 || im==modeIndex) {
      s << "# Mode " << im << "\n";
      for (Iterator it(this, im); !it.atEnd(); ++it) {
        const RealValue& val=it.value();
        if (val.isValid())
          s << fmt.arg(it.x()*xFactor, 0, 'g', 15).arg(val.value(), 0, 'g', 15);
      }
    }
  }
}
void QGpCoreWave::ModalStorage::writeReport ( QDataStream &  s) const

Save all to report (including frequency), mixing refines and original samples, removing invalid values

Begin with a table of the number of samples per mode

4 nm Modes For each mode: 4 number of samples for this mode mode, nx[im] nx[im]*2*8 values x,y

References QGpCoreWave::ModalStorage::Iterator::atEnd(), QGpCoreTools::Value< numberType >::isValid(), TRACE, and QGpCoreTools::Value< numberType >::value().

{
  TRACE;
  double xFactor=0.5/M_PI;
  s << _nModes;
  for(int im=0; im < _nModes;im++ ) {
    // room for sample per mode table
    int xCountReal=0;
    qint64 nSamplesPos=s.device() ->pos();
    s << xCountReal;
    for(Iterator it(this, im); !it.atEnd(); ++it) {
      const RealValue & val=it.value();
      if(val.isValid()) {
        s << it.x() * xFactor << val.value();
        xCountReal++;
      }
    }
    // Write the real number of x at the beginning
    qint64 curPos=s.device() ->pos();
    s.device() ->seek(nSamplesPos);
    s << xCountReal;
    s.device() ->seek(curPos);
  }
}
const QVector<double>* QGpCoreWave::ModalStorage::x ( ) const [inline]

Referenced by QGpCoreWave::Dispersion::deltaK(), QGpCoreWave::ModalStorage::Iterator::Iterator(), Spac3CForward::omega(), QGpCoreWave::Ellipticity::peakMisfit(), QGpCoreWave::Dispersion::refine(), set(), and QGpCoreWave::Dispersion::setGroupSlowness().

{return _x;}
double QGpCoreWave::ModalStorage::x ( int  ix) const [inline] 
{return _x->at(ix);}
int QGpCoreWave::ModalStorage::xCount ( ) const [inline]

Referenced by QGpCoreWave::Dispersion::deltaK(), DinverDCCore::TargetList::dispersionSampleCount(), QGpCoreWave::ModalStorage::Iterator::operator++(), QGpCoreWave::Ellipticity::peakMisfit(), QGpCoreWave::Ellipticity::peaks(), QGpCoreWave::Dispersion::refine(), and QGpCoreWave::Dispersion::setGroupSlowness().

{return _nX;}
The documentation for this class was generated from the following files:
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Defines

Modified 2013-10-22

Valid CSS! Valid XHTML 1.0!