variogrammodels.h

Go to the documentation of this file.

#pragma once

/*+
________________________________________________________________________

 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        A. Huck & H.Huck
 Date:          27-09-2011
________________________________________________________________________


-*/

#include <math.h>


inline void getVariogramModel( const char* typestr, float nugget, float sill,
                              float range, int size, float* in, float* out )
{
    const FixedString type( typestr );
    if ( type=="exponential" )
    {
        for ( int idx=0; idx<size; idx++ )
            out[idx] = (sill-nugget)*(1-exp(-3*in[idx]/range))+nugget;
    }
    else if ( type=="gaussian" )
    {
        for ( int idx=0; idx<size; idx++ )
        {
            out[idx] = (sill-nugget)*
                           (1-exp(-3*((in[idx]*in[idx])/(range*range))))+
                           nugget;
        }
    }
    else if ( type=="spherical" )
    {
        for ( int idx=0; idx<size; idx++ )
        {
            if ( in[idx] < range )
            {
                out[idx] = (sill-nugget)*(((3*in[idx])/(2*range))-
                           ((in[idx]*in[idx]*in[idx])/(2*range*range*range)))
                           +nugget;
            }
            else
                out[idx] = sill;
        }
    }
}