[gnuastro-commits] master 2c6c7e0 076/113: Imported recent work, no conflicts

[Mohammad Akhlaghi]

branch: master
commit 2c6c7e0dcda1a23ca039516cf85ec81703e8f8db
Merge: 5dac3c3 7cf87f9
Author: Mohammad Akhlaghi <mohammad@akhlaghi.org>
Commit: Mohammad Akhlaghi <mohammad@akhlaghi.org>

    Imported recent work, no conflicts
    
    There weren't any conflicts in this merge.
---
 NEWS                       |  5 +++++
 bin/mkcatalog/args.h       | 10 ++++-----
 bin/mkcatalog/columns.c    |  4 ++--
 bin/mkcatalog/main.h       |  2 +-
 bin/mkcatalog/mkcatalog.c  | 21 ++++++++++++-----
 bin/mkcatalog/ui.c         | 47 +++++++++++++++++++-------------------
 bin/mkcatalog/ui.h         |  4 ++--
 bin/mkcatalog/upperlimit.c | 30 +++++++++++++------------
 doc/gnuastro.texi          | 56 +++++++++++++++++++++++-----------------------
 lib/gnuastro/statistics.h  | 11 +--------
 lib/statistics.c           |  1 +
 11 files changed, 100 insertions(+), 91 deletions(-)

diff --git a/NEWS b/NEWS
index ef6d63d..5a359f9 100644
--- a/NEWS
+++ b/NEWS
@@ -12,6 +12,10 @@ GNU Astronomy Utilities NEWS                          -*- 
outline -*-
 
 ** Changed features
 
+  MakeCatalog:
+    - `--checkuplim': new name for `--checkupperlimit'.
+    - `--brightnessnoriver': new name for `--noriverbrightness'.
+
   Library:
     - gal_txt_write: new `colinfoinstdout' argument.
     - gal_table_write: new `colinfoinstdout' argument.
@@ -21,6 +25,7 @@ GNU Astronomy Utilities NEWS                          -*- 
outline -*-
   bug #54057: Building failure due to not finding gsl_interp_steffen.
   bug #54063: Match tests in make check fail randomly.
   bug #54186: MakeCatalog's --checkupperlimit not keeping output's name.
+  bug #54188: MakeCatalog's Upperlimit not being sigma-clipped properly.
 
 
 
diff --git a/bin/mkcatalog/args.h b/bin/mkcatalog/args.h
index feb5954..c06b3d9 100644
--- a/bin/mkcatalog/args.h
+++ b/bin/mkcatalog/args.h
@@ -334,13 +334,13 @@ struct argp_option program_options[] =
       GAL_OPTIONS_NOT_SET
     },
     {
-      "checkupperlimit",
-      UI_KEY_CHECKUPPERLIMIT,
+      "checkuplim",
+      UI_KEY_CHECKUPLIM,
       "INT[,INT]",
       0,
       "Check random distribution for one label.",
       UI_GROUP_UPPERLIMIT,
-      &p->checkupperlimit,
+      &p->checkuplim,
       GAL_TYPE_STRING,
       GAL_OPTIONS_RANGE_GT_0,
       GAL_OPTIONS_NOT_MANDATORY,
@@ -855,8 +855,8 @@ struct argp_option program_options[] =
       ui_column_codes_ll
     },
     {
-      "noriverbrightness",
-      UI_KEY_NORIVERBRIGHTNESS,
+      "brightnessnoriver",
+      UI_KEY_BRIGHTNESSNORIVER,
       0,
       0,
       "Sky (not river) subtracted clump brightness.",
diff --git a/bin/mkcatalog/columns.c b/bin/mkcatalog/columns.c
index 7fcf0b2..56f7881 100644
--- a/bin/mkcatalog/columns.c
+++ b/bin/mkcatalog/columns.c
@@ -922,7 +922,7 @@ columns_define_alloc(struct mkcatalogparams *p)
           oiflag[ OCOL_C_SUM ] = 1;
           break;
 
-        case UI_KEY_NORIVERBRIGHTNESS:
+        case UI_KEY_BRIGHTNESSNORIVER:
           name           = "NO_RIVER_BRIGHTNESS";
           unit           = MKCATALOG_NO_UNIT;
           ocomment       = NULL;
@@ -2035,7 +2035,7 @@ columns_fill(struct mkcatalog_passparams *pp)
                                         : NAN );
             break;
 
-          case UI_KEY_NORIVERBRIGHTNESS:
+          case UI_KEY_BRIGHTNESSNORIVER:
             ((float *)colarr)[cind] = ( ci[ CCOL_NUM ]>0.0f
                                         ? ci[ CCOL_SUM ] : NAN );
             break;
diff --git a/bin/mkcatalog/main.h b/bin/mkcatalog/main.h
index 3f6b0d8..dbf54e3 100644
--- a/bin/mkcatalog/main.h
+++ b/bin/mkcatalog/main.h
@@ -181,7 +181,7 @@ struct mkcatalogparams
   uint8_t             envseed;  /* Use the environment for random seed. */
   double       upsigmaclip[2];  /* Sigma clip to measure upper limit.   */
   float              upnsigma;  /* Multiple of sigma to define up-lim.  */
-  int32_t  checkupperlimit[2];  /* Object & clump ID to check dist.     */
+  int32_t       checkuplim[2];  /* Object & clump ID to check dist.     */
 
   /* Internal. */
   time_t              rawtime;  /* Starting time of the program.        */
diff --git a/bin/mkcatalog/mkcatalog.c b/bin/mkcatalog/mkcatalog.c
index 97ef3bc..5e3a8f1 100644
--- a/bin/mkcatalog/mkcatalog.c
+++ b/bin/mkcatalog/mkcatalog.c
@@ -125,18 +125,29 @@ mkcatalog_single_object(void *in_prm)
 
   /* If we have upper-limit mode, then allocate the container to keep the
      values to calculate the standard deviation. */
-  pp.up_vals = ( p->upperlimit
-                 ? gal_data_alloc(NULL, GAL_TYPE_FLOAT32, 1, &p->upnum,
+  if(p->upperlimit)
+    {
+      /* Allocate the space to keep the upper-limit values. */
+      pp.up_vals = gal_data_alloc(NULL, GAL_TYPE_FLOAT32, 1, &p->upnum,
                                   NULL, 0, p->cp.minmapsize, NULL, NULL,
-                                  NULL)
-                 : NULL );
+                                  NULL);
+
+      /* Set the blank checked flag to 1. By definition, this dataset won't
+         have any blank values. Also `flag' is initialized to `0'. So we
+         just have to set the checked flag (`GAL_DATA_FLAG_BLANK_CH') to
+         one to inform later steps that there are no blank values. */
+      pp.up_vals->flag |= GAL_DATA_FLAG_BLANK_CH;
+    }
+  else
+    pp.up_vals=NULL;
+
 
   /* Fill the desired columns for all the objects given to this thread. */
   for(i=0; tprm->indexs[i]!=GAL_BLANK_SIZE_T; ++i)
     {
       /* For easy reading. Note that the object IDs start from one while
          the array positions start from 0. */
-      pp.ci=NULL;
+      pp.ci     = NULL;
       pp.object = tprm->indexs[i] + 1;
       pp.tile   = &p->tiles[ tprm->indexs[i] ];
 
diff --git a/bin/mkcatalog/ui.c b/bin/mkcatalog/ui.c
index b7b9311..be48f60 100644
--- a/bin/mkcatalog/ui.c
+++ b/bin/mkcatalog/ui.c
@@ -121,15 +121,15 @@ ui_initialize_options(struct mkcatalogparams *p,
   cp->coptions           = gal_commonopts_options;
 
   /* Specific to this program. */
-  p->medstd              = NAN;
-  p->sfmagnsigma         = NAN;
-  p->sfmagarea           = NAN;
-  p->upnsigma            = NAN;
-  p->zeropoint           = NAN;
-  p->upsigmaclip[0]      = NAN;
-  p->upsigmaclip[1]      = NAN;
-  p->checkupperlimit[0]  = GAL_BLANK_INT32;
-  p->checkupperlimit[1]  = GAL_BLANK_INT32;
+  p->medstd         = NAN;
+  p->sfmagnsigma    = NAN;
+  p->sfmagarea      = NAN;
+  p->upnsigma       = NAN;
+  p->zeropoint      = NAN;
+  p->upsigmaclip[0] = NAN;
+  p->upsigmaclip[1] = NAN;
+  p->checkuplim[0]  = GAL_BLANK_INT32;
+  p->checkuplim[1]  = GAL_BLANK_INT32;
 
   /* Modify common options. */
   for(i=0; !gal_options_is_last(&cp->coptions[i]); ++i)
@@ -264,14 +264,13 @@ ui_check_upperlimit(struct argp_option *option, char *arg,
   /* Write. */
   if(lineno==-1)
     {
-      if(p->checkupperlimit[1]==GAL_BLANK_INT32)
+      if(p->checkuplim[1]==GAL_BLANK_INT32)
         {
-          if( asprintf(&str, "%d", p->checkupperlimit[0])<0 )
+          if( asprintf(&str, "%d", p->checkuplim[0])<0 )
             error(EXIT_FAILURE, 0, "%s: asprintf allocation", __func__);
         }
       else
-        if( asprintf(&str, "%d,%d", p->checkupperlimit[0],
-                     p->checkupperlimit[1])<0 )
+        if( asprintf(&str, "%d,%d", p->checkuplim[0], p->checkuplim[1])<0 )
           error(EXIT_FAILURE, 0, "%s: asprintf allocation", __func__);
       return str;
     }
@@ -322,8 +321,8 @@ ui_check_upperlimit(struct argp_option *option, char *arg,
         }
 
       /* Write the values in. */
-      p->checkupperlimit[0] = d[0];
-      p->checkupperlimit[1] = raw->size==2 ? d[1] : GAL_BLANK_INT32;
+      p->checkuplim[0] = d[0];
+      p->checkuplim[1] = raw->size==2 ? d[1] : GAL_BLANK_INT32;
 
       /* For no un-used variable warning. This function doesn't need the
          pointer.*/
@@ -365,14 +364,14 @@ ui_read_check_only_options(struct mkcatalogparams *p)
   /* If an upper-limit check table is requested with a specific clump, but
      no clump catalog has been requested, then abort and inform the
      user. */
-  if( p->checkupperlimit[1]!=GAL_BLANK_INT32 && p->clumpscat==0 )
+  if( p->checkuplim[1]!=GAL_BLANK_INT32 && p->clumpscat==0 )
     error(EXIT_FAILURE, 0, "no clumps catalog is requested, hence "
-          "`--checkupperlimit' is only available for objects (one value "
+          "`--checkuplim' is only available for objects (one value "
           "must be given to it).\n\n"
           "To ask for a clumps catalog, please append `--clumpscat' to the "
           "command calling MakeCatalog.\n\n"
           "If you want the upperlimit check table for an object, only give "
-          "one value (the object's label) to `--checkupperlimit'.");
+          "one value (the object's label) to `--checkuplim'.");
 
   /* See if `--skyin' is a filename or a value. When the string is ONLY a
      number (and nothing else), `tailptr' will point to the end of the
@@ -998,11 +997,11 @@ ui_preparations_read_inputs(struct mkcatalogparams *p)
     {
       /* If an upperlimit check was requested, make sure the object number
          is not larger than the maximum number of labels. */
-      if(p->checkupperlimit[0] != GAL_BLANK_INT32
-         && p->checkupperlimit[0] > p->numobjects)
+      if(p->checkuplim[0] != GAL_BLANK_INT32
+         && p->checkuplim[0] > p->numobjects)
         error(EXIT_FAILURE, 0, "%d (object identifier for the "
-              "`--checkupperlimit' option) is larger than the number of "
-              "objects in the input labels (%zu)", p->checkupperlimit[0],
+              "`--checkuplim' option) is larger than the number of "
+              "objects in the input labels (%zu)", p->checkuplim[0],
               p->numobjects);
 
       /* Read the mask file if it was given. */
@@ -1210,7 +1209,7 @@ ui_preparations_outnames(struct mkcatalogparams *p)
     }
 
   /* If an upperlimit check image is requsted, then set its filename. */
-  if(p->checkupperlimit)
+  if(p->checkuplim)
     {
       /* See if the directory should be respected. */
       p->cp.keepinputdir = p->cp.output ? 1 : p->cp.keepinputdir;
@@ -1640,7 +1639,7 @@ ui_free_report(struct mkcatalogparams *p, struct timeval 
*t1)
             "input is less than double their length. If the input is taken "
             "from a larger dataset, this issue can be solved by using a "
             "larger part of it. You can also run MakeCatalog with "
-            "`--checkupperlimit' to see the distribution for a special "
+            "`--checkuplim' to see the distribution for a special "
             "object or clump as a table and personally inspect its "
             "reliability. \n\n");
 
diff --git a/bin/mkcatalog/ui.h b/bin/mkcatalog/ui.h
index b711271..16e63a3 100644
--- a/bin/mkcatalog/ui.h
+++ b/bin/mkcatalog/ui.h
@@ -97,7 +97,7 @@ enum option_keys_enum
   UI_KEY_UPRANGE,
   UI_KEY_UPSIGMACLIP,
   UI_KEY_UPNSIGMA,
-  UI_KEY_CHECKUPPERLIMIT,
+  UI_KEY_CHECKUPLIM,
   UI_KEY_NOCLUMPSORT,
 
   UI_KEY_OBJID,                         /* Catalog columns. */
@@ -128,7 +128,7 @@ enum option_keys_enum
   UI_KEY_CLUMPSGEOW3,
   UI_KEY_BRIGHTNESSERR,
   UI_KEY_CLUMPSBRIGHTNESS,
-  UI_KEY_NORIVERBRIGHTNESS,
+  UI_KEY_BRIGHTNESSNORIVER,
   UI_KEY_MEAN,
   UI_KEY_MEDIAN,
   UI_KEY_CLUMPSMAGNITUDE,
diff --git a/bin/mkcatalog/upperlimit.c b/bin/mkcatalog/upperlimit.c
index a338a10..32dc016 100644
--- a/bin/mkcatalog/upperlimit.c
+++ b/bin/mkcatalog/upperlimit.c
@@ -418,12 +418,12 @@ upperlimit_write_check(struct mkcatalogparams *p, 
gal_list_sizet_t *check_x,
 
 
   /* Write exactly what object/clump this table is for. */
-  if( p->checkupperlimit[1]!=GAL_BLANK_INT32 )
-    if( asprintf(&tmp2, ", Clump %d", p->checkupperlimit[1]) <0 )
+  if( p->checkuplim[1]!=GAL_BLANK_INT32 )
+    if( asprintf(&tmp2, ", Clump %d", p->checkuplim[1]) <0 )
       error(EXIT_FAILURE, 0, "%s: asprintf allocation", __func__);
   if( asprintf(&tmp, "Upperlimit distribution for Object %d%s",
-               p->checkupperlimit[0],
-               ( p->checkupperlimit[1]==GAL_BLANK_INT32
+               p->checkuplim[0],
+               ( p->checkuplim[1]==GAL_BLANK_INT32
                  ? "" : tmp2) ) <0 )
     error(EXIT_FAILURE, 0, "%s: asprintf allocation", __func__);
   gal_list_str_add(&comments, tmp, 0);
@@ -494,7 +494,8 @@ upperlimit_measure(struct mkcatalog_passparams *pp, int32_t 
clumplab,
             case UI_KEY_UPPERLIMITSIGMA:
             case UI_KEY_UPPERLIMITONESIGMA:
 
-              /* We only need to do this once. */
+              /* We only need to do this once, but the columns can be
+                 requested in any order. */
               if(sigclip==NULL)
                 {
                   /* Calculate the sigma-clipped standard deviation. Since
@@ -504,7 +505,7 @@ upperlimit_measure(struct mkcatalog_passparams *pp, int32_t 
clumplab,
                   init_size=pp->up_vals->size;
                   sigclip=gal_statistics_sigma_clip(pp->up_vals,
                                                     p->upsigmaclip[0],
-                                                    p->upsigmaclip[1], 1, 1);
+                                                    p->upsigmaclip[1],1,1);
                   pp->up_vals->size=pp->up_vals->dsize[0]=init_size;
                   scarr=sigclip->array;
 
@@ -586,16 +587,16 @@ upperlimit_one_tile(struct mkcatalog_passparams *pp, 
gal_data_t *tile,
                                       "rcoord");
 
   /* See if a check table must be created for this distribution. */
-  if( p->checkupperlimit[0]==pp->object )
+  if( p->checkuplim[0]==pp->object )
     {
       /* We are on a clump */
       if( clumplab )
         {
-          if( p->checkupperlimit[1]==clumplab )
+          if( p->checkuplim[1]==clumplab )
             writecheck=1;
         }
       else
-        if( p->checkupperlimit[1]==GAL_BLANK_INT32 )
+        if( p->checkuplim[1]==GAL_BLANK_INT32 )
           writecheck=1;
     }
 
@@ -603,6 +604,7 @@ upperlimit_one_tile(struct mkcatalog_passparams *pp, 
gal_data_t *tile,
   /* Initializations. */
   tarray=tile->array;
   gsl_rng_set(pp->rng, seed);
+  pp->up_vals->flag &= ~GAL_DATA_FLAG_SORT_CH;
 
 
   /* Set the range of random values for this tile. */
@@ -774,13 +776,13 @@ upperlimit_calculate(struct mkcatalog_passparams *pp)
       /* If an upper-limit check image is requested, then make sure that
          the clump label is not more than the number of clumps in this
          object. */
-      if( p->checkupperlimit[0] == pp->object
-          && p->checkupperlimit[1] != GAL_BLANK_INT32
-          && p->checkupperlimit[1] > pp->clumpsinobj )
+      if( p->checkuplim[0] == pp->object
+          && p->checkuplim[1] != GAL_BLANK_INT32
+          && p->checkuplim[1] > pp->clumpsinobj )
         error(EXIT_FAILURE, 0, "object %d has %zu clumps, but an upperlimit "
-              "check table (using the `--checkupperlimit' option) has been "
+              "check table (using the `--checkuplim' option) has been "
               "requested for clump %d", pp->object, pp->clumpsinobj,
-              p->checkupperlimit[1]);
+              p->checkuplim[1]);
 
       /* Make tiles covering the clumps. */
       clumptiles=upperlimit_make_clump_tiles(pp);
diff --git a/doc/gnuastro.texi b/doc/gnuastro.texi
index 10b7e93..33e6412 100644
--- a/doc/gnuastro.texi
+++ b/doc/gnuastro.texi
@@ -736,6 +736,8 @@ history, before it has been fully documented here first. As 
discussed in
 @section Quick start
 
 @cindex Test
+@cindex Gzip
+@cindex Lzip
 @cindex Check
 @cindex Build
 @cindex Compile
@@ -749,18 +751,14 @@ and robust archival features, an 
@url{http://www.nongnu.org/lzip/lzip.html,
 Lzip} compressed tarball is also available at
 @url{http://ftp.gnu.org/gnu/gnuastro/gnuastro-latest.tar.lz,
 @file{gnuastro-latest.tar.lz}}, see @ref{Release tarball} for more details
-on the tarball release.
+on the tarball release@footnote{The Gzip library and program are commonly
+available on most systems. However, Gnuastro recommends Lzip as described
+above and the beta-releases are also only distributed in @file{tar.lz}. You
+can download and install Lzip's source (in @file{.tar.gz} format) from its
+webpage and follow the same process as below: Lzip has no dependencies, so
+simply decompress, then run @command{./configure}, @command{make},
+@command{sudo make install}.}.
 
-@cindex Gzip
-@cindex Lzip
-@cindex GNU Tar
-The Gzip library and program are commonly available on most
-systems. However, Gnuastro recommends Lzip as described above and the
-beta-releases are also only distributed in @file{tar.lz}. You can download
-and install Lzip's source (in @file{.tar.gz} format) from its webpage and
-follow the same process as below: Lzip has no dependencies, so simply
-decompress, then run @command{./configure}, @command{make}, @command{sudo
-make install}.
 
 Let's assume the downloaded tarball is in the @file{TOPGNUASTRO}
 directory. The first two commands below can be used to decompress the
@@ -786,9 +784,10 @@ $ tar xf gnuastro-latest.tar
 @end example
 
 Gnuastro has three mandatory dependencies and some optional dependencies
-for extra functionality, see @ref{Dependencies} for the full list. Links to
-the dependency sources and instructions on installing each are fully
-described there. When the mandatory dependencies are ready, you can
+for extra functionality, see @ref{Dependencies} for the full list. In
+@ref{Dependencies from package managers} we have prepared the command to
+easily install Gnuastro's dependencies using the package manager of some
+operating systems. When the mandatory dependencies are ready, you can
 configure, compile, check and install Gnuastro on your system with the
 following commands.
 
@@ -4992,12 +4991,11 @@ results and understanding (of what is going on, and 
thus how you interpret
 your scientific results) will also correspondingly improve.
 @end enumerate
 
-Based on your operating system and its package manager, you can use any of
-the following commands to install the mandatory and optional
-dependencies. If your operating system or package manager are not included,
-please send us the respective command, so we add it here. Gnuastro itself
-if also already packaged in some package managers (for example Debian or
-Homebrew).
+Based on your package manager, you can use any of the following commands to
+install the mandatory and optional dependencies. If your package manager
+isn't included in the list below, please send us the respective command, so
+we add it. Gnuastro itself if also already packaged in some package
+managers (for example Debian or Homebrew).
 
 As discussed above, we recommend installing the @emph{mandatory}
 dependencies manually from source. Therefore, in each command below, first
@@ -5007,7 +5005,7 @@ mandatory dependencies with your package manager, just 
ignore the empty
 line.
 
 @table @asis
-@item Debian-based OSs (@command{apt-get})
+@item @command{apt-get} (Debian-based OSs: Ubuntu, Linux Mint and etc)
 @cindex Debian
 @cindex Ubuntu
 @cindex @command{apt-get}
@@ -5025,7 +5023,7 @@ $ sudo apt-get install ghostscript libtool-bin  
libjpeg-dev \
                        libgsl0-dev libcfitsio-dev wcslib-dev
 @end example
 
-@item macOS (@command{brew})
+@item @command{brew} (macOS)
 @cindex macOS
 @cindex Homebrew
 @cindex MacPorts
@@ -5051,7 +5049,7 @@ $ brew tap brewsci/science
 $ brew install wcslib
 @end example
 
-@item Arch Linux (@command{pacman})
+@item @command{pacman} (Arch Linux)
 @cindex Arch Linux
 @cindex @command{pacman}
 Arch Linux is a smaller GNU/Linux distribution. As discussed in
@@ -5060,7 +5058,7 @@ KISS principle ("keep it simple, stupid") as the general 
guideline, and
 focuses on elegance, code correctness, minimalism and simplicity, and
 expects the user to be willing to make some effort to understand the
 system's operation''. Arch Linux uses ``Package manager'' (Pacman) to
-manage its packages.
+manage its packages/components.
 @example
 $ sudo pacman -S ghostscript libtool libjpeg libtiff libgit2 \
                                                              \
@@ -17665,7 +17663,7 @@ The multiple of the final (@mymath{\sigma}-clipped) 
standard deviation (or
 @mymath{\sigma}) used to measure the upper-limit brightness or
 magnitude.
 
-@item --checkupperlimit=INT[,INT]
+@item --checkuplim=INT[,INT]
 Print a table of positions and measured values for all the full random
 distribution used for one particular object or clump. If only one integer
 is given to this option, it is interpreted to be an object's label. If two
@@ -17874,7 +17872,7 @@ usable pixels (blank or below the threshold) are 
present over the clump or
 object, the stored value will be NaN, because zero (note that zero is
 meaningful).
 
-@item --noriverbrightness
+@item --brightnessnoriver
 [Clumps] The Sky (not river) subtracted clump brightness. By definition,
 for the clumps, the average brightness of the rivers surrounding it are
 subtracted from it for a first order accounting for contamination by
@@ -26875,9 +26873,11 @@ is larger than @code{1} (one), it must be an integer 
and will be
 interpreted as the number clips to do. If it is less than @code{1} (one),
 it is interpreted as the tolerance level to stop the iteration.
 
-The output dataset has the following elements with a
-@code{GAL_TYPE_FLOAT32} type:
+The returned dataset (let's call it @code{out}) contains a four-element
+array with type @code{GAL_TYPE_FLOAT32}. The final number of clips is
+stored in the @code{out->status}.
 @example
+float *array=out->array;
 array[0]: Number of points used.
 array[1]: Median.
 array[2]: Mean.
diff --git a/lib/gnuastro/statistics.h b/lib/gnuastro/statistics.h
index 92b7868..ad96d01 100644
--- a/lib/gnuastro/statistics.h
+++ b/lib/gnuastro/statistics.h
@@ -167,7 +167,7 @@ gal_statistics_cfp(gal_data_t *data, gal_data_t *bins, int 
normalize);
 
 
 /****************************************************************
- *****************        Sigma clip         ********************
+ *****************         Outliers          ********************
  ****************************************************************/
 gal_data_t *
 gal_statistics_sigma_clip(gal_data_t *input, float multip, float param,
@@ -177,15 +177,6 @@ gal_statistics_sigma_clip(gal_data_t *input, float multip, 
float param,
 
 
 
-/****************************************************************/
-/*************         Identify outliers         ****************/
-/****************************************************************/
-void
-gal_statistics_remove_outliers_flat_cdf(float *sorted, size_t *outsize);
-
-
-
-
 __END_C_DECLS    /* From C++ preparations */
 
 #endif           /* __GAL_STATISTICS_H__ */
diff --git a/lib/statistics.c b/lib/statistics.c
index a29303e..7d6e2f4 100644
--- a/lib/statistics.c
+++ b/lib/statistics.c
@@ -2068,6 +2068,7 @@ gal_statistics_sigma_clip(gal_data_t *input, float 
multip, float param,
 
 
 
+
 #if 0
 /* Using the cumulative distribution function this function will
    remove outliers from a dataset. */