interp.h

浏览该文件的文档.

#ifndef LMMC_INTERP_H
#define LMMC_INTERP_H
 
#include "lmmc/config.h"
#include "lmmc/status.h"
#include <stddef.h>
 
#ifdef __cplusplus
extern "C" {
#endif

lmmc_status_t lmmc_interp_linear(
    const lmmc_real_t* xs,
    const lmmc_real_t* ys,
    size_t n,
    lmmc_real_t query_x,
    lmmc_real_t* out_y
);

typedef struct lmmc_interp_cspline_t lmmc_interp_cspline_t;

typedef enum {
    LMMC_SPLINE_NATURAL = 0,   
    LMMC_SPLINE_CLAMPED,       
    LMMC_SPLINE_NOT_A_KNOT,    
    LMMC_SPLINE_PERIODIC       
} lmmc_spline_bc_t;

lmmc_status_t lmmc_interp_cspline_create(
    const lmmc_real_t* xs,
    const lmmc_real_t* ys,
    size_t n,
    lmmc_interp_cspline_t** out_spline
);

lmmc_status_t lmmc_interp_cspline_create_ex(
    const lmmc_real_t* xs,
    const lmmc_real_t* ys,
    size_t n,
    lmmc_spline_bc_t bc,
    lmmc_real_t deriv_left,
    lmmc_real_t deriv_right,
    lmmc_interp_cspline_t** out_spline
);

lmmc_status_t lmmc_interp_cspline_eval(
    const lmmc_interp_cspline_t* spline,
    lmmc_real_t query_x,
    lmmc_real_t* out_y
);

void lmmc_interp_cspline_destroy(lmmc_interp_cspline_t* spline);

typedef struct lmmc_interp_pchip_t lmmc_interp_pchip_t;

lmmc_status_t lmmc_interp_pchip_create(
    const lmmc_real_t* xs,
    const lmmc_real_t* ys,
    size_t n,
    lmmc_interp_pchip_t** out
);

lmmc_status_t lmmc_interp_pchip_eval(
    const lmmc_interp_pchip_t* p,
    lmmc_real_t x,
    lmmc_real_t* out_y
);

void lmmc_interp_pchip_destroy(lmmc_interp_pchip_t* p);

typedef struct lmmc_interp_akima_t lmmc_interp_akima_t;

lmmc_status_t lmmc_interp_akima_create(
    const lmmc_real_t* xs,
    const lmmc_real_t* ys,
    size_t n,
    lmmc_interp_akima_t** out
);

lmmc_status_t lmmc_interp_akima_eval(
    const lmmc_interp_akima_t* a,
    lmmc_real_t x,
    lmmc_real_t* out_y
);

void lmmc_interp_akima_destroy(lmmc_interp_akima_t* a);

typedef struct lmmc_interp_lagrange_t lmmc_interp_lagrange_t;

lmmc_status_t lmmc_interp_lagrange_create(
    const lmmc_real_t* xs,
    const lmmc_real_t* ys,
    size_t n,
    lmmc_interp_lagrange_t** out_lagrange
);

lmmc_status_t lmmc_interp_lagrange_eval(
    const lmmc_interp_lagrange_t* lagrange,
    lmmc_real_t query_x,
    lmmc_real_t* out_y
);

void lmmc_interp_lagrange_destroy(lmmc_interp_lagrange_t* lagrange);

lmmc_status_t lmmc_interp_bilinear(
    const lmmc_real_t* xs, size_t nx,
    const lmmc_real_t* ys, size_t ny,
    const lmmc_real_t* zs,
    lmmc_real_t qx, lmmc_real_t qy,
    lmmc_real_t* out_z
);

lmmc_status_t lmmc_interp_bicubic(
    const lmmc_real_t* xs, size_t nx,
    const lmmc_real_t* ys, size_t ny,
    const lmmc_real_t* zs,
    lmmc_real_t qx, lmmc_real_t qy,
    lmmc_real_t* out_z
);
 
#ifdef __cplusplus
}
#endif
 
#endif