debugging ongoing
Symmetries still alter results
todo:
work on flow eq implementation
impl product G x S
implement output hdf5 and plotting of it