A method is disclosed. The method includes applying a human visual system (HVS) model to a Continuous Tone Image (CTI) and a initial Half Tone Image (HTI) to generate a perceived CTI and a perceived HTI and computing a change in pixel error for a first pixel by toggling the first pixel with all the possible output states and swapping the first pixel with all neighbor pixels.