Function periodic 

pub fn periodic(blocks: Blocks) -> Result<(DMatrix<Complex>, PeriodicResidual)>

Convert scattering blocks to the transfer matrix used for periodic bands.

Upstream: treams.SMatrices.periodic.