anchor.diff:  		    (mem-aref jcol :int i) col)
   	      (incf i)))
   	  (assert (= i nele_hess)))
  -	))
  +	;; Return the values. Symmetric, fill lower left
  +	;; Fill objective portion
  +	(progn
  +	  (setf (mem-aref values :double 0) (* obj_factor 2 x3) ; 0,0
  +		(mem-aref values :double 1) (* obj_factor x3) ; 1,0
  +		(mem-aref values :double 2) 0 ; 1,1
  +		
  +		(mem-aref values :double 3) (* obj_factor x3) ; 2,0
  +		(mem-aref values :double 4) 0 ; 2,1
  +		(mem-aref values :double 5) 0 ; 2,2
  +		
  +		(mem-aref values :double 6) (* obj_factor (+ (* 2 x0) x1 x2)) ; 3,0
  +		(mem-aref values :double 7) (* obj_factor x0) ; 3,1
  +		(mem-aref values :double 8) (* obj_factor x0) ; 3,2
  +		(mem-aref values :double 9) 0) ; 3,3
  +	  ;; Add portion for the 1st constraint
  +	  (incf (mem-aref values :double 1) (* lmb0 x2 x3)) ; 1,0
  +	  
  +	  (incf (mem-aref values :double 3) (* lmb0 x1 x3)) ; 2,0
  +	  (incf (mem-aref values :double 4) (* lmb0 x0 x3)) ; 2,1
  +	  
  +	  (incf (mem-aref values :double 6) (* lmb0 x1 x2)) ; 3,0
  +	  (incf (mem-aref values :double 7) (* lmb0 x0 x2)) ; 3,1
  +	  (incf (mem-aref values :double 8) (* lmb0 x0 x1)) ; 3,2
  +	  
  +	  ;; Add portion for the 2nd constraint
  +	  (incf (mem-aref values :double 0) (* lmb1 2)) ; 0,0
  +	  (incf (mem-aref values :double 2) (* lmb1 2)) ; 1,1
  +	  (incf (mem-aref values :double 5) (* lmb1 2)) ; 2,2
  +	  (incf (mem-aref values :double 9) (* lmb1 2))))) ; 3,3
     1)
   
   (defun test ()
  
positive.code: 	 (x_u :double n)
  	 (g_l :double m)
  	 (g_u :pointer m)
  	 (status :int)
  	 (x :pointer n)
  	 (mult_x_l :double n)
  	 (mult_x_u :double n)
  	 (obj :double))
        )))
  
negative.code:   (let ((x0 (mem-aref x :double 0))
  	(x1 (mem-aref x :double 1))
  	(x2 (mem-aref x :double 2))
  	(x3 (mem-aref x :double 3))
  	(lmb0 (mem-aref lmb :double 0))
  	(lmb1 (mem-aref lmb :double 1)))
      (if (null-pointer-p values) ;; Return Hessian structure
  	(let ((i 0))
  	  (dotimes (row 4)
  