>>804
∇f = 2 A' (A x - b)
H_f = 2 A' A