A numpy implementation of PMVS