updated formulas in polarity function so that:

delta = (1 + c(x_{i}^A - x_{i}^D))(-1) / sqrt(3)
(1+ 0.8 (0 - 1))(-1)/sqrt(3)
-.2 /sqrt(3) = -0.115

issue #141