reflectivitymodel.h

Go to the documentation of this file.

#pragma once
 
/*+
________________________________________________________________________
 
 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        Kristofer
 Date:          Jan 2011
 RCS:           $Id$
________________________________________________________________________
 
-*/
 
#include "algomod.h"
#include "odcomplex.h"
#include "sets.h"
 
 
mClass(Algo) ReflectivitySpike
{
public:
                        ReflectivitySpike()
                            : reflectivity_( mUdf(float), mUdf(float) )
                            , time_( mUdf(float) )
                            , correctedtime_( mUdf(float) )
                            , depth_( mUdf(float) )
                        {}
 
    inline bool         isDefined() const;
    inline float        time(bool isnmo=true) const;
 
    inline bool         operator==(const ReflectivitySpike& s) const;
    inline bool         operator!=(const ReflectivitySpike& s) const;
 
    float_complex       reflectivity_;
    float               time_;
    float               correctedtime_; 
    float               depth_;
};
 
 
typedef TypeSet<ReflectivitySpike> ReflectivityModel;
 
 
//Implementations
 
inline bool ReflectivitySpike::operator==(const ReflectivitySpike& s) const
{
    return mIsEqualWithUdf( reflectivity_.real(),s.reflectivity_.real(),1e-5) &&
           mIsEqualWithUdf( reflectivity_.imag(),s.reflectivity_.imag(),1e-5) &&
           mIsEqualWithUdf( time_, s.time_, 1e-5 ) &&
           mIsEqualWithUdf( correctedtime_, s.correctedtime_, 1e-5 ) &&
           mIsEqualWithUdf( depth_, s.depth_, 1e-5 );
}
 
inline bool ReflectivitySpike::operator!=(const ReflectivitySpike& s) const
{ return !(*this==s); }
 
 
inline bool ReflectivitySpike::isDefined() const
{
    return !mIsUdf(reflectivity_) && !mIsUdf(time_) &&
           !mIsUdf(correctedtime_) && !mIsUdf(depth_);
}
 
 
inline float ReflectivitySpike::time( bool isnmo ) const
{
    return isnmo ? correctedtime_ : time_;
}