| 588 | } |
| 589 | |
| 590 | double EEDFTwoTermApproximation::norm(const Eigen::VectorXd& f, const Eigen::VectorXd& grid) |
| 591 | { |
| 592 | string m_quadratureMethod = "simpson"; |
| 593 | Eigen::VectorXd p(f.size()); |
| 594 | for (int i = 0; i < f.size(); i++) { |
| 595 | p[i] = f(i) * pow(grid[i], 0.5); |
| 596 | } |
| 597 | return numericalQuadrature(m_quadratureMethod, p, grid); |
| 598 | } |
| 599 | |
| 600 | } |