make mvn plain normal instead of f-normal to prepare for derivatives
add operator== to Scalar