| 7 | namespace seissol::functions { |
| 8 | |
| 9 | uint64_t rangeProduct(uint64_t from, uint64_t to) { |
| 10 | uint64_t product = 1; |
| 11 | for (; from <= to; ++from) { |
| 12 | product *= from; |
| 13 | } |
| 14 | return product; |
| 15 | } |
| 16 | |
| 17 | double JacobiP(unsigned n, unsigned a, unsigned b, double x) { |
| 18 | if (n == 0) { |