| 17 | |
| 18 | template <typename Tx> |
| 19 | INLINE auto NodalShapePowerL2(Tx lam, int i, int order) |
| 20 | { |
| 21 | Tx prod = Tx(1.0); |
| 22 | for (int j = 0; j < i; j++) |
| 23 | prod *= ((order+1)*lam - j-0.5) / (i-j); |
| 24 | return prod; |
| 25 | } |
| 26 | |
| 27 | |
| 28 | template <> template<typename Tx, typename TFA> |