@amabdelrehim I don't really understand where the weighing by the inverse eigenvalue is taking place. It seems to me like initwork[i] should be multiplied by inv[H][i], no?