
ENHANCEMENTS

* calculate both s- and p-polarisation at once in recursive_fresnelcpp

CLEANUP

* tests using testthat

* simple demos

* check the derivation of fresnel coeffs for T, turn doc into vignette

* write up the matrix transfer method
