/*!
 * \file init2d.c
 *
 * \brief Initialization of interpolation library data structures
 *
 * \author H. Mitasova, I. Kosinovsky, D. Gerdes Fall 1993 (original authors)
 * \author modified by McCauley in August 1995
 * \author modified by Mitasova in August 1995
 * \author modified by Brown in June 1999 - added elatt & smatt
 *
 * \copyright
 * (C) 1993-1999 by Helena Mitasova and the GRASS Development Team
 *
 * \copyright
 * This program is free software under the
 * GNU General Public License (>=v2).
 * Read the file COPYING that comes with GRASS
 * for details.
 *
 */

#include <stdio.h>
#include <math.h>
#include <unistd.h>
#include <grass/gis.h>
#include <grass/interpf.h>


/*! Initializes parameters used by the library */
void IL_init_params_2d(struct interp_params *params,
                       FILE * inp,  /*!< input stream */
                       int elatt,  /*!< which fp att in sites file? 1 = first */
                       int smatt,  /*!< which fp att in sites file to use for 
                                         * smoothing? (if zero use sm) 1 = first */
                       double zm,  /*!< multiplier for z-values */
                       int k1,  /*!< min number of points per segment for interpolation */
                       int k2,  /*!< max number of points per segment */
                       char *msk,  /*!< name of mask */
                       int rows, int cols,  /*!< number of rows and columns */
                       DCELL * ar1, DCELL * ar2, DCELL * ar3, DCELL * ar4, DCELL * ar5,
                       DCELL * ar6,  /*!< arrays for interpolated values (ar1-ar6) */
                       double tension,  /*!< tension */
                       int k3,  /*!< max number of points for interpolation */
                       int sc1, int sc2, int sc3,  /*!< multipliers for interpolation values */
                       double sm,  /*!< smoothing */
                       char *f1, char *f2, char *f3, char *f4, char *f5,
                       char *f6,  /*!< output files (f1-f6) */
                       double dm,  /*!< min distance between points */
                       double x_or,  /*!< x of origin */
                       double y_or,  /*!< y of origin */
                       int der,  /*!< 1 if compute partial derivatives */
                       double tet,  /*!< anisotropy angle (0 is East, counter-clockwise) */
                       double scl,  /*!< anisotropy scaling factor */
                       FILE * t1, FILE * t2, FILE * t3, FILE * t4, FILE * t5,
                       FILE * t6,  /*!< temp files for writing interp. values (t1-t6) */
                       FILE * dev,  /*!< pointer to deviations file */
                       struct TimeStamp *ts,
                       int c,  /*!< cross validation */
                       const char *wheresql     /*!< SQL WHERE statement */
    )
{
    params->fdinp = inp;
    params->elatt = elatt;
    params->smatt = smatt;
    params->zmult = zm;
    params->kmin = k1;
    params->kmax = k2;
    params->maskmap = msk;
    params->nsizr = rows;
    params->nsizc = cols;
    params->az = ar1;
    params->adx = ar2;
    params->ady = ar3;
    params->adxx = ar4;
    params->adyy = ar5;
    params->adxy = ar6;
    params->fi = tension;
    params->KMAX2 = k3;
    params->scik1 = sc1;
    params->scik2 = sc2;
    params->scik3 = sc3;
    params->rsm = sm;
    params->elev = f1;
    params->slope = f2;
    params->aspect = f3;
    params->pcurv = f4;
    params->tcurv = f5;
    params->mcurv = f6;
    params->dmin = dm;
    params->x_orig = x_or;
    params->y_orig = y_or;
    params->deriv = der;
    params->theta = tet;
    params->scalex = scl;
    params->Tmp_fd_z = t1;
    params->Tmp_fd_dx = t2;
    params->Tmp_fd_dy = t3;
    params->Tmp_fd_xx = t4;
    params->Tmp_fd_yy = t5;
    params->Tmp_fd_xy = t6;
    params->fddevi = dev;
    params->ts = ts;
    params->cv = c;
    params->wheresql = wheresql;
}

/*! Initializes functions used by the library  */
void IL_init_func_2d(struct interp_params *params,
                     grid_calc_fn * grid_f,  /*!< calculates grid for given segment */
                     matrix_create_fn * matr_f,  /*!< creates matrix for a given segment */
                     check_points_fn * point_f,  /*!< checks interpolation function at points */
                     secpar_fn * secp_f,  /*!< calculates aspect, slope, curvature */
                     interp_fn * interp_f,  /*!< radial basis function */
                     interpder_fn * interpder_f,  /*!< derivatives of radial basis function */
                     wr_temp_fn * temp_f        /*!< writes temp files */
    )
{
    params->grid_calc = grid_f;
    params->matrix_create = matr_f;
    params->check_points = point_f;
    params->secpar = secp_f;
    params->interp = interp_f;
    params->interpder = interpder_f;
    params->wr_temp = temp_f;

}