Inverse matrix diagonalization for photonic band calculation