add test for polarization calculation
fix issue #32 but probably also breaks some tests ...
this is TDD