| 339 | } |
| 340 | |
| 341 | inline void evaluateGradPolynomials(double* phis, double xi, double eta, int numPoly) { |
| 342 | assert(numPoly > 0); |
| 343 | unsigned idx = 0; |
| 344 | for (unsigned int d = 0; d <= static_cast<unsigned>(numPoly); ++d) { |
| 345 | for (unsigned int j = 0; j <= d; ++j) { |
| 346 | const auto grad = seissol::functions::gradTriDubinerP({d - j, j}, {xi, eta}); |
| 347 | for (const auto& g : grad) { |
| 348 | phis[idx++] = g; |
| 349 | } |
| 350 | } |
| 351 | } |
| 352 | } |
| 353 | } // namespace tri_dubiner |
| 354 | } // namespace basisFunction |
| 355 | } // namespace seissol |
no test coverage detected