do generic defect calculation the way one expects it (that does work for non-symmetric matrices too)