Unfortunately, solving this minimization problem is NP hard (Wagner and Wagner 1993;Von Luxburg 2007).
an approximate solution can be found by relaxing the constraints on the elements of H in constraint 4.