share/minpack/rtest_minpack.mac

   1 (kill(all), 0);
   2 0$
   3
   4 (load(minpack), float_approx_equal_tolerance : 1e-12, 0);
   5 0$
   6
   7 minpack_lsquares(powell(x1,x2,x3,x4), [x1,x2,x3,x4],
   8   [3,-1,0,1]), powell(x1,x2,x3,x4) :=
   9                    [x1+10*x2, sqrt(5)*(x3-x4), (x2-2*x3)^2,
  10                    sqrt(10)*(x1-x4)^2];
  11 [[1.652117596168394e-17, - 1.652117596168393e-18,
  12                   2.643388153869468e-18, 2.643388153869468e-18],
  13                   6.109327859207777e-34, 4]$
  14
  15 minpack_lsquares(powell(x1,x2,x3,x4), [x1,x2,x3,x4],
  16   [3,-1,0,1], jacobian = false), powell(x1,x2,x3,x4) :=
  17                    [x1+10*x2, sqrt(5)*(x3-x4), (x2-2*x3)^2,
  18                    sqrt(10)*(x1-x4)^2];
  19 [[5.060282149485331e-11, - 5.060282149491206e-12,
  20                   2.179447843547218e-11, 2.179447843547218e-11],
  21                   3.534491794847031e-21, 5]$
  22
  23 minpack_solve ( [x1 + x2 - 2, x1 - x2 - 2], [x1, x2], [0, 0]);
  24 [[2.0,-1.479114197289397e-31],2.091783358105259e-31,1]$
  25
  26 minpack_solve ( [x1 + x2 - 2, x1 - x2 - 2], [x1, x2], [0, 0], jacobian=true);
  27 [[2.0, - 1.479114197289397e-31], 2.091783358105259e-31, 1]$
  28
  29 /* An example from User Guide for MINPACK-1, Chapter 4
  30 http://www.mcs.anl.gov/~more/ANL8074b.pdf
  31 */
  32 block([eq, varlist, gu, x], local(eq, x),
  33   eq: makelist([], 9),
  34   eq[1]: (3-2*x[1])*x[1]            - 2*x[2]                         + 1,
  35   for i: 2 thru 8 do
  36   eq[i]:         -x[i-1] + (3-2*x[i])*x[i]            - 2*x[i+1]     + 1,
  37   eq[9]:                             -x[8]   + (3-2*x[9])*x[9]       + 1,
  38   varlist: makelist(x[i], i, 9),
  39   gu : makelist(0,       9),
  40   minpack_solve (eq, varlist, gu));
  41 [[-0.5706545124664266,-0.6816283412828047,-0.7017324513732917,
  42   -0.7042129396771661,-0.7013690482934034,-0.6918656444557993,
  43   -0.6657920125574163,-0.5960342005611687,-0.4164120628161226],
  44  2.306068468708184e-10,1]$
  45
  46 (reset (float_approx_equal_tolerance), 0);
  47 0$