Ice Sheet System Model  4.18
Code documentation
Public Member Functions | Data Fields
SphericalVariogram Class Reference

#include <SphericalVariogram.h>

Inheritance diagram for SphericalVariogram:
Variogram Object

Public Member Functions

 SphericalVariogram ()
 
 SphericalVariogram (Options *options)
 
 ~SphericalVariogram ()
 
Objectcopy ()
 
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 11 of file SphericalVariogram.h.

Constructor & Destructor Documentation

◆ SphericalVariogram() [1/2]

SphericalVariogram::SphericalVariogram ( )

Definition at line 15 of file SphericalVariogram.cpp.

15  {/*{{{*/
16  this->nugget = 0.2;
17  this->sill = 1;
18  this->range = SQRT3;
19  return;
20 }

◆ SphericalVariogram() [2/2]

SphericalVariogram::SphericalVariogram ( Options options)

Definition at line 22 of file SphericalVariogram.cpp.

22  {/*{{{*/
23 
24  /*Defaults*/
25  this->nugget = 0.2;
26  this->sill = 1;
27  this->range = SQRT3;
28 
29  /*Overwrite from options*/
30  if(options->GetOption("nugget")) options->Get(&this->nugget,"nugget");
31  if(options->GetOption("sill")) options->Get(&this->sill,"sill");
32  if(options->GetOption("range")) options->Get(&this->range,"range");
33 
34  /*Checks*/
35  if(nugget==sill) _error_("nugget and sill cannot be equal (constant semivariogram not allowed)");
36 }

◆ ~SphericalVariogram()

SphericalVariogram::~SphericalVariogram ( )

Definition at line 38 of file SphericalVariogram.cpp.

38  {/*{{{*/
39  return;
40 }

Member Function Documentation

◆ copy()

Object * SphericalVariogram::copy ( void  )
virtual

Implements Object.

Definition at line 44 of file SphericalVariogram.cpp.

44  {/*{{{*/
45  return new SphericalVariogram(*this);
46 }

◆ DeepEcho()

void SphericalVariogram::DeepEcho ( void  )
inlinevirtual

Implements Object.

Definition at line 25 of file SphericalVariogram.h.

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

◆ Echo()

void SphericalVariogram::Echo ( void  )
virtual

Implements Object.

Definition at line 48 of file SphericalVariogram.cpp.

48  {/*{{{*/
49  _printf_("SphericalVariogram\n");
50  _printf_(" nugget: " << this->nugget << "\n");
51  _printf_(" sill : " << this->sill << "\n");
52  _printf_(" range : " << this->range << "\n");
53 }

◆ Id()

int SphericalVariogram::Id ( void  )
inlinevirtual

Implements Object.

Definition at line 27 of file SphericalVariogram.h.

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

◆ Marshall()

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

Implements Object.

Definition at line 28 of file SphericalVariogram.h.

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

◆ ObjectEnum()

int SphericalVariogram::ObjectEnum ( void  )
inlinevirtual

Implements Object.

Definition at line 29 of file SphericalVariogram.h.

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

◆ Covariance()

double SphericalVariogram::Covariance ( double  deltax,
double  deltay 
)
virtual

Implements Variogram.

Definition at line 57 of file SphericalVariogram.cpp.

57  {/*{{{*/
58  /*The covariance can be deduced from the variogram from the following
59  * relationship:
60  * 2 gamma = C(x,x) + C(y,y) -2 C(x,y)
61  * so
62  * C(h) = sill - gamma */
63  double h,cova;
64 
65  /*Calculate length square*/
66  h=sqrt(deltax*deltax + deltay*deltay);
67 
68  /*return covariance*/
69  if(h<=range)
70  cova = (sill-nugget)*(1 - (3*h)/(2*range) + pow(h,3)/(2*pow(range,3)) );
71  else
72  cova = 0.;
73 
74  return cova;
75 }

◆ SemiVariogram()

double SphericalVariogram::SemiVariogram ( double  deltax,
double  deltay 
)
virtual

Implements Variogram.

Definition at line 77 of file SphericalVariogram.cpp.

77  {/*{{{*/
78  /*http://en.wikipedia.org/wiki/Variogram*/
79  double h,gamma;
80 
81  /*Calculate length square*/
82  h=sqrt(deltax*deltax + deltay*deltay);
83 
84  /*return semi-variogram*/
85  if(h<=range)
86  gamma = (sill-nugget)*( (3*h)/(2*range) - pow(h,3)/(2*pow(range,3)) ) + nugget;
87  else
88  gamma = sill;
89 
90  return gamma;
91 }

Field Documentation

◆ nugget

double SphericalVariogram::nugget

Definition at line 14 of file SphericalVariogram.h.

◆ sill

double SphericalVariogram::sill

Definition at line 15 of file SphericalVariogram.h.

◆ range

double SphericalVariogram::range

Definition at line 16 of file SphericalVariogram.h.


The documentation for this class was generated from the following files:
SphericalVariogram::nugget
double nugget
Definition: SphericalVariogram.h:14
_printf_
#define _printf_(StreamArgs)
Definition: Print.h:22
SphericalVariogram::SphericalVariogram
SphericalVariogram()
Definition: SphericalVariogram.cpp:15
Options::Get
void Get(OptionType *pvalue, const char *name)
Definition: Options.h:21
SphericalVariogram::range
double range
Definition: SphericalVariogram.h:16
Options::GetOption
Option * GetOption(const char *name)
Definition: Options.cpp:67
_error_
#define _error_(StreamArgs)
Definition: exceptions.h:49
SQRT3
#define SQRT3
Definition: constants.h:10
SphericalVariogram::sill
double sill
Definition: SphericalVariogram.h:15