share/contrib/rtest_augmented_lagrangian.mac

   1 (kill(all), 'done);
   2 'done$
   3
   4 (load (augmented_lagrangian), float_approx_equal_tolerance : 1e-12, 'done);
   5 'done$
   6
   7 block([FOM: x^2 + 2*y^2, xx: [x, y], C: [x + y - 1], yy: [1, 1]],
   8   augmented_lagrangian_method(FOM, xx, C, yy, iprint=[-1,0]));
   9 [[x = 0.666659841080023,y = 0.333340272455448],%lambda = [-1.333337940892518]]$
  10
  11 block([FOM: x^2 + 2*y^2, xx: [x, y], grad : [2*x, 4*y], C: [x + y - 1], yy: [1, 1]],
  12   augmented_lagrangian_method ([FOM, grad], xx, C, yy,
  13     iprint = [-1, 0]));
  14 [[x = 0.6666598410800247,y = 0.3333402724554464],
  15 %lambda = [-1.333337940892525]]$
  16
  17 (reset(float_approx_equal_tolerance),0);
  18 0$