SeisTools.h

Go to the documentation of this file.

#ifndef SEISTOOLS_H_
#define SEISTOOLS_H_

#include <algorithm>
#include <cassert>
#include <boost/bind.hpp>
#include <vector>
#include <gsl/gsl_math.h>
/** \addtogroup seistools Seismic data analysis and modeling */
/* @{ */

//! Calculate density from a given S-velocity, the formula is taken from Owen et al. JGR 89,7783-7795 and modified for vs
struct CalcDensity{
        double operator()(const double svel)
        {
                return (0.77 + 0.554 * svel);
        }
};

inline void Normalize(std::vector<double> &Trace)
{
        const double maxamp = 1./ *max_element(Trace.begin(),Trace.end());
        assert(gsl_finite(maxamp));
        std::transform(Trace.begin(),Trace.end(),Trace.begin(),boost::bind(std::multiplies<double>(),_1,maxamp));
}
/* @} */
#endif /*SEISTOOLS_H_*/

---