added ability to work with NAs as input
fixed problem with nrOps = 1
first version of a GEVP mass determination for symmetric matrices