Operator precedence bug in calculation of R13

octue / es-flow

Atmospheric and marine flow characterisation tools.

Other

3 stars 0 forks source link

Operator precedence bug in calculation of R13 #39

Closed thclark closed 5 years ago

thclark commented 5 years ago

In reynolds_stress_13 there is an operator precedence bug computing the N term.

Once validated against MATLAB implementation, correct to n = get_coles_wake(1.0, pi_coles) + pi_coles * (wc_plus - wc_minus) / (2.0 * d_pi); because without hte final set of brackets, this presently gives the wrong answer.

thclark commented 5 years ago

Closing as validated code correct.