|
1 | | -// Copyright 2022 -2024 Thijs Janzen |
| 1 | +// Copyright 2022 -2025 Thijs Janzen |
2 | 2 | // This program is free software: you can redistribute it and/or modify |
3 | 3 | // it under the terms of the GNU General Public License as published by |
4 | 4 | // the Free Software Foundation, either version 3 of the License, or |
@@ -428,13 +428,14 @@ namespace correction { |
428 | 428 | doublecorrect_pda(double Ic,size_t num_tips) {// colless corrections |
429 | 429 | double denom =powf(num_tips,1.5f); |
430 | 430 | return1.0 * Ic / denom; |
431 | | -} |
| 431 | +} |
432 | 432 |
|
433 | | -doublecorrect_yule(double Ic,size_t num_tips) { |
| 433 | +doublecorrect_yule(double Ic,size_t num_tips) { |
434 | 434 | staticconstdouble g =0.577215664901532; |
435 | 435 | auto output = (Ic - |
436 | 436 | num_tips *log(num_tips) - |
437 | 437 | num_tips * (g -1 -log(2))) / num_tips; |
438 | 438 | return output; |
439 | | - } |
| 439 | +} |
| 440 | + |
440 | 441 | }// namespace correction |