'Place H in the Hessian of all functions.'

def place_H(self, big_H, H, var_offsets):

    offset = 0
    for var in self.variables():
        var_size = (var.size[0] * var.size[1])
        var_H = H[offset:(offset + var_size), offset:(offset + var_size)]
        self.block_add(big_H, var_H, var_offsets[var.data], var_offsets[var.data], var_size, var_size)
        offset += var_size
