cauchy

[JanuszBedkowski]

which one os correct?

inline double cauchy(double delta, double b){
return 1.0 / (M_PI * b *( 1.0 + ((delta)/b) * ((delta)/b) ) );
}

inline double get_cauchy_w(double r, double c){
return 1/(1 + (r/c)*(r/c));
}

[mwlasiuk]

which one os correct?

inline double cauchy(double delta, double b){ return 1.0 / (M_PI * b *( 1.0 + ((delta)/b) * ((delta)/b) ) ); }

inline double get_cauchy_w(double r, double c){ return 1/(1 + (r/c)*(r/c)); }

This is normalized Cauchy with peak at 1.0 / (M_PI * b) and integrral == 1:

inline double cauchy(double delta, double b){
return 1.0 / (M_PI * b *( 1.0 + ((delta)/b) * ((delta)/b) ) );
}

This is unnormalized (weighted) Cauchy where peak is equal to 1 but integral not equal to 1:

inline double get_cauchy_w(double r, double c){
return 1/(1 + (r/c)*(r/c));
}