Function compute_gpt_center

Function Documentation

void compute_gpt_center(double alpha0, const double *r0, double alpha1, const double *r1, double gamma_inv, double *gpt_center)

The sum of two Gaussian centres (via the gaussian product theorem).

\[ \frac{\alpha_{1}R_{0}+\alpha_{2}R_{1}}{\alpha_{1}+\alpha_{2}} \]

See eqn 2.3 and preceeding line in Taketa et al. (1966)

Parameters
  • alpha0: Exponent on centre 0
  • r0: Centre 0
  • alpha1: Exponent on centre 1
  • r1: Centre 1
  • gamma_inv: \( \frac{1}{\alpha_{1}+\alpha_{2}} \)
  • gpt_center: Output variable. The resultant gaussian centre.