AnisoSurfaceWaveObjective.h

Go to the documentation of this file.

#ifndef ANISOSURFACEWAVEOBJECTIVE_H_
#define ANISOSURFACEWAVEOBJECTIVE_H_

#include "PlottableObjective.h"
#include "ParkSurfaceWaveData.h"
#include "AnisoSurfaceWaveSynthetic.h"
#include <string>

namespace gplib
  {
    /** \addtogroup seistools Seismic data analysis and modeling */
    /* @{ */
//! This class calculates the misfit for anisotropic surface wave dispersion data
class AnisoSurfaceWaveObjective: public PlottableObjective
  {
private:
  ParkSurfaceWaveData MeasuredData;
  ParkSurfaceWaveData SynthData;
  AnisoSurfaceWaveSynthetic Synthetic;
  double errorlevel;
  double poisson;
  double backazimuth;
  double avelratio;
public:

  //! Write the current model to ascii file for calculations
  virtual void WriteModel(const std::string &filename)
    {
      Synthetic.WriteModel(filename);
    }
  //! Write the current model to ascii file for plotting
  virtual void WritePlot(const std::string &filename)
    {
      Synthetic.WritePlot(filename);
    }
  //! Write synthetic data as ascii file
  virtual void WriteData(const std::string &filename)
        {
           SynthData.WriteAscii(filename);
    }
  //! Provide read only access to the synthetic data
  const ParkSurfaceWaveData &GetSynthetic() const
    {
      return SynthData;
    }
  AnisoSurfaceWaveObjective(const AnisoSurfaceWaveObjective &Old);
  AnisoSurfaceWaveObjective& operator=(const AnisoSurfaceWaveObjective& source);
  AnisoSurfaceWaveObjective(const ParkSurfaceWaveData &Data, const double ba, const double avel);
  virtual ~AnisoSurfaceWaveObjective();
  //! We need clone and create for building an array of derived objects, see FAQ lite 20.8, the return type depends on the derived class
  virtual AnisoSurfaceWaveObjective *clone() const
    {
      return new AnisoSurfaceWaveObjective(*this);
    }
  //! Some operations cannot be done in parallel, these are done before
  virtual void PreParallel(const ttranscribed &member);
  //! Some operations cannot be done in parallel, these are done after, returns the misfit value
  virtual double PostParallel(const ttranscribed &member);
  //! The core performance calculation, has to be safe to be done in parallel
  virtual void SafeParallel(const ttranscribed &member);
  };
  /* @} */
  }
#endif /*ANISOSURFACEWAVEOBJECTIVE_H_*/