simplification in homogenize/matvec.py