Areas of improvement
- `PDE` supports excluded points
- Add `PointSetBC` for training using exact solution.
- Add `Hypersphere.boundary_normal()`
- Add `Triangle.inside()`
API changes
- `PDE` doesn't support numpy.array `X` as the third parameter. To define pde with extra coefficients, use `PDE(..., auxiliary_var_function=...)`, see an example at Lorenz_inverse_forced_Colab.ipynb.