/*********************************************************************
tile -- work with tesselations over a host dataset.
This is part of GNU Astronomy Utilities (Gnuastro) package.

Original author:
     Mohammad Akhlaghi <mohammad@akhlaghi.org>
Contributing author(s):
Copyright (C) 2017-2021, Free Software Foundation, Inc.

Gnuastro is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation, either version 3 of the License, or (at your
option) any later version.

Gnuastro is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
General Public License for more details.

You should have received a copy of the GNU General Public License
along with Gnuastro. If not, see <http://www.gnu.org/licenses/>.
**********************************************************************/
#ifndef __GAL_INTERPOLATE_H__
#define __GAL_INTERPOLATE_H__

/* Include other headers if necessary here. Note that other header files
   must be included before the C++ preparations below */
#include <gnuastro/data.h>
#include <gnuastro/tile.h>
#include <gsl/gsl_spline.h>

/* C++ Preparations */
#undef __BEGIN_C_DECLS
#undef __END_C_DECLS
#ifdef __cplusplus
# define __BEGIN_C_DECLS extern "C" {
# define __END_C_DECLS }
#else
# define __BEGIN_C_DECLS                /* empty */
# define __END_C_DECLS                  /* empty */
#endif
/* End of C++ preparations */

/* Actual header contants (the above were for the Pre-processor). */
__BEGIN_C_DECLS  /* From C++ preparations */




/* Metrics to use for nearest-neighbor  */
enum gal_interpolate_neighbors_metric
{
 GAL_INTERPOLATE_NEIGHBORS_METRIC_INVALID,

 GAL_INTERPOLATE_NEIGHBORS_METRIC_RADIAL,
 GAL_INTERPOLATE_NEIGHBORS_METRIC_MANHATTAN,
};



/* Types of interpolation for nearest-neighbor. */
enum gal_interpolate_neighbors_func
{
 GAL_INTERPOLATE_NEIGHBORS_FUNC_INVALID,

 GAL_INTERPOLATE_NEIGHBORS_FUNC_MIN,
 GAL_INTERPOLATE_NEIGHBORS_FUNC_MAX,
 GAL_INTERPOLATE_NEIGHBORS_FUNC_MEDIAN,
};



/* Types of interpolation. */
enum gal_interpolate_1D_types
{
 GAL_INTERPOLATE_1D_INVALID,

 GAL_INTERPOLATE_1D_LINEAR,
 GAL_INTERPOLATE_1D_POLYNOMIAL,
 GAL_INTERPOLATE_1D_CSPLINE,
 GAL_INTERPOLATE_1D_CSPLINE_PERIODIC,
 GAL_INTERPOLATE_1D_AKIMA,
 GAL_INTERPOLATE_1D_AKIMA_PERIODIC,
 GAL_INTERPOLATE_1D_STEFFEN,
};



gal_data_t *
gal_interpolate_neighbors(gal_data_t *input,
                          struct gal_tile_two_layer_params *tl,
                          uint8_t metric, size_t numneighbors,
                          size_t numthreads, int onlyblank,
                          int aslinkedlist, int function);

gsl_spline *
gal_interpolate_1d_make_gsl_spline(gal_data_t *X, gal_data_t *Y, int type_1d);

void
gal_interpolate_1d_blank(gal_data_t *in, int type_1d);


__END_C_DECLS    /* From C++ preparations */

#endif
