GaussianVariogram Class Reference

#include <GaussianVariogram.h>

Inheritance diagram for GaussianVariogram:

Public Member Functions
	GaussianVariogram ()
	GaussianVariogram (Options *options)
	~GaussianVariogram ()
Object *	copy ()
void	DeepEcho ()
void	Echo ()
int	Id ()
void	Marshall (char **pmarshalled_data, int *pmarshalled_data_size, int marshall_direction)
int	ObjectEnum ()
double	Covariance (double deltax, double deltay)
double	SemiVariogram (double deltax, double deltay)
Public Member Functions inherited from Variogram
virtual	~Variogram ()
Public Member Functions inherited from Object
virtual	~Object ()

Data Fields
double	nugget
double	sill
double	range

Detailed Description

Definition at line 12 of file GaussianVariogram.h.

Constructor & Destructor Documentation

GaussianVariogram() [1/2]

GaussianVariogram::GaussianVariogram

(

)

Definition at line 15 of file GaussianVariogram.cpp.

                                    {/*{{{*/
   this->nugget = 0.2;
   this->sill   = 1;
   this->range  = SQRT3;
   return;
}

GaussianVariogram() [2/2]

GaussianVariogram::GaussianVariogram

(

Options *

options

)

Definition at line 22 of file GaussianVariogram.cpp.

                                                    {/*{{{*/
 
   /*Defaults*/
   this->nugget = 0.2;
   this->sill   = 1;
   this->range  = SQRT3;
 
   /*Overwrite from options*/
   if(options->GetOption("nugget")) options->Get(&this->nugget,"nugget");
   if(options->GetOption("sill"))   options->Get(&this->sill,"sill");
   if(options->GetOption("range"))  options->Get(&this->range,"range");
 
   /*Checks*/
   if(nugget==sill) _error_("nugget and sill cannot be equal (constant semivariogram not allowed)");
}

~GaussianVariogram()

GaussianVariogram::~GaussianVariogram

(

)

Definition at line 38 of file GaussianVariogram.cpp.

                                     {/*{{{*/
   return;
}

Member Function Documentation

copy()

Object * GaussianVariogram::copy ( void  )

virtual

Implements Object.

Definition at line 44 of file GaussianVariogram.cpp.

                                   {/*{{{*/
      return new GaussianVariogram(*this);
}

DeepEcho()

void GaussianVariogram::DeepEcho ( void  )

inlinevirtual

Implements Object.

Definition at line 26 of file GaussianVariogram.h.

{_error_("Not implemented yet");};

Echo()

void GaussianVariogram::Echo ( void  )

virtual

Implements Object.

Definition at line 48 of file GaussianVariogram.cpp.

                                {/*{{{*/
   _printf_("GaussianVariogram\n");
   _printf_("   nugget: " << this->nugget << "\n");
   _printf_("   sill  : " << this->sill << "\n");
   _printf_("   range : " << this->range << "\n");
}

Id()

int GaussianVariogram::Id ( void  )

inlinevirtual

Implements Object.

Definition at line 28 of file GaussianVariogram.h.

{_error_("Not implemented yet");}; 

Marshall()

void GaussianVariogram::Marshall ( char **  pmarshalled_data, int *  pmarshalled_data_size, int  marshall_direction  )

inlinevirtual

Implements Object.

Definition at line 29 of file GaussianVariogram.h.

{ _error_("not implemented yet!"); };

ObjectEnum()

int GaussianVariogram::ObjectEnum ( void  )

inlinevirtual

Implements Object.

Definition at line 30 of file GaussianVariogram.h.

{_error_("Not implemented yet");};

Covariance()

double GaussianVariogram::Covariance ( double  deltax, double  deltay  )

virtual

Implements Variogram.

Definition at line 57 of file GaussianVariogram.cpp.

                                                               {/*{{{*/
   /*The covariance can be deduced from the variogram from the following
    * relationship:
    *    2 gamma = C(x,x) + C(y,y) -2 C(x,y)
    * so
    *    C(h) = sill - gamma                                            */
   double h2,a,cova;
 
   /*Calculate length square*/
   h2=deltax*deltax + deltay*deltay;
 
   /*If h is too small, return sill*/
   if(h2<0.0000001) return sill;
 
   /*compute covariance*/
   a     = 1./3.;
   cova = (sill-nugget)*exp(-h2/(a*range*range));
 
   return cova;
}

SemiVariogram()

double GaussianVariogram::SemiVariogram ( double  deltax, double  deltay  )

virtual

Implements Variogram.

Definition at line 78 of file GaussianVariogram.cpp.

                                                                  {/*{{{*/
   /*http://en.wikipedia.org/wiki/Variogram*/
   double h2,a,gamma;
 
   /*Calculate length square*/
   h2=deltax*deltax + deltay*deltay;
 
   /*return semi-variogram*/
   a     = 1./3.;
   gamma = (sill-nugget)*(1.-exp(-h2/(a*range*range))) + nugget;
 
   //if(h2>1000*1000) _printf_("gamma = " << gamma << " h= " << sqrt(h2) << "\n");
   _printf_("h = " << sqrt(h2) << " gamma = " << gamma << "\n");
   return gamma;
}

Field Documentation

nugget

double GaussianVariogram::nugget

Definition at line 15 of file GaussianVariogram.h.

sill

double GaussianVariogram::sill

Definition at line 16 of file GaussianVariogram.h.

range

double GaussianVariogram::range

Definition at line 17 of file GaussianVariogram.h.

---

The documentation for this class was generated from the following files:

• src/c/classes/kriging/GaussianVariogram.h
• src/c/classes/kriging/GaussianVariogram.cpp