"PROJECTIONS METHOD" PEA
1. They use the tensor method to approximate the conditional expectation with orthogonal Chebyshev polynomials.
2. The coefficients of the approximating function are such that they minimize the distance between the approximating function and the numerically calculated conditional expectation at a set of grid points.
3. The grid points are Chebyshev nodes.
4. The numerical integration procedure used to calculate the conditional expectation is Hermite Gaussian Quadrature. In my experience it is easier to obtain an accurate solution fast with quadrature methods than with Monte Carlo methods. An example of a PEA algorithm that uses Monte Carlo methods can be found at Monte Carlo PEA
5. The "iterative" programs iterate on a projection procedure to find the coefficients of the approximating function.
6. The "equation-solver" programs use a nonlinear equation solver to find the
value of the coefficients at which the approximating function equal the
numerically calculated conditional expectation.