changeset 3656:96679fb690a4

[project @ 2000-04-04 06:16:22 by jwe]
author jwe
date Tue, 04 Apr 2000 06:16:23 +0000
parents 6324a4994278
children a908150a3a32
files scripts/ChangeLog scripts/statistics/distributions/beta_inv.m scripts/statistics/distributions/f_inv.m src/ChangeLog src/dynamic-ld.cc
diffstat 5 files changed, 14 insertions(+), 9 deletions(-) [+]
line wrap: on
line diff
--- a/scripts/ChangeLog
+++ b/scripts/ChangeLog
@@ -1,3 +1,8 @@
+2000-04-04  John Smith <john@arrows.demon.co.uk>
+
+	* statistics/distributions/beta_inv.m: Provide better(?) initial
+	guess for iteration.
+
 2000-03-31   Paul Kienzle  <pkienzle@kienzle.powernet.co.uk>
 
 	* image/image.m: Allow image (A) or image (x, y, A).
--- a/scripts/statistics/distributions/beta_inv.m
+++ b/scripts/statistics/distributions/beta_inv.m
@@ -57,7 +57,7 @@
     a = a (k);
     b = b (k);
     x = x (k);
-    y = a ./ b;
+    y = a ./ (a + b);
     l = find (y < eps);
     if (any (l))
       y(l) = sqrt (eps) * ones (length (l), 1);
--- a/scripts/statistics/distributions/f_inv.m
+++ b/scripts/statistics/distributions/f_inv.m
@@ -54,11 +54,7 @@
 
   k = find ((x > 0) & (x < 1) & (m > 0) & (n > 0));
   if (any (k))
-    fprintf (stderr, "n1: %f\n", n(k));
-    fprintf (stderr, "n2: %f\n", m(k));
-    t = beta_inv (1 - x(k), n(k) / 2, m(k) / 2)
-    fprintf (stderr, "qbeta: %f\n", t);
-    inv(k) = ((1 ./ t - 1)
+    inv(k) = ((1 ./ beta_inv (1 - x(k), n(k) / 2, m(k) / 2) - 1)
 	      .* n(k) ./ m(k));
   endif
 
--- a/src/ChangeLog
+++ b/src/ChangeLog
@@ -1,3 +1,7 @@
+2000-04-04  John W. Eaton  <jwe@bevo.che.wisc.edu>
+
+	* dynamic-ld.cc (octave_dynamic_loader::do_load): Undo previous change.
+
 2000-04-03  John W. Eaton  <jwe@bevo.che.wisc.edu>
 
 	* dynamic-ld.cc (octave_dynamic_loader::do_load): Also fail with
--- a/src/dynamic-ld.cc
+++ b/src/dynamic-ld.cc
@@ -267,11 +267,11 @@
 	= X_CAST (octave_dld_fcn_installer, function);
 
       retval = f (oct_file);
+
+      if (! retval)
+	::error ("failed to install dld function `%s'", fcn_name.c_str ());
     }
   
-  if (! retval)
-    ::error ("failed to install dld function `%s'", fcn_name.c_str ());
-
   unwind_protect::run_frame ("octave_dynamic_loader::do_load");
 
   return retval;