ytBSpline.h

/*
  math/ytBSpline.{h,c} : B-spline routines.
  Copyright (C) 2018, Yoshinori Tamada <tamada A T ytlab.jp>
  All rights reserved.
 
  See LICENSE.txt for details of the licensing agreement.
*/
 
#ifndef __YTLIB_BSPLINE_H
#define __YTLIB_BSPLINE_H
 
#include <stdio.h>
 
void ytBSpline_designMatrix(const double * x, const int n, 
                            const double xl, const double xr,
                            double * B, int * index,
                            const int ndx, const int degree, double * work);
void ytBSpline_designMatrix2(const double * x, const int n, double * B, int * index,
                            const int ndx, const int degree,  double * work);
void ytBSpline_Dk(const int m, const int k, double * D);
void ytBSpline_K(int m, int k, double * K, double * work);
double ytBSpline_det_K_plus(const double * K, const int m);
void ytBSpline_syBtB(const double * B, const int * index, const int n,
                     const int ndx, const int degree, double * BtB);
void ytBSpline_syBtB2(const double * B, const int * index, const int n,
                      const int ndx, const int degree, double * BtB);
void ytBSpline_dsycopy(const double * X, double * Y, int m);
void ytBSpline_dsycopy20(const double * X, double * Y);
void ytBSpline_syMBt(const double * restrict A,
                     const double * restrict B, const int* index,
                     const int n, const int ndx, const int degree, 
                     double * restrict C);
void ytBSpline_syMBt_D3(const double * restrict A,
                        const double * restrict B, const int * index,
                        const int n, const int ndx,
                        double * restrict C);
void ytBSpline_syMBt_D3_2(const double * restrict A,
                          const double * restrict B, const int * index,
                          const int n, const int ndx,
                          double * restrict C);
void ytBSpline_syMBt_D3_3(double * restrict A,
                          double * restrict B, const int * index,
                          const int n, const int ndx,
                          double * restrict C);
void ytBSpline_dxpy_mask(int n, double * x, const double * y, const char * mask);
void ytBSpline_dxmy_mask(int n, double * x, const double * y, const char * mask);
void ytBSpline_Mv(int n, int m, const double * A, const double * b, double * c);
void ytBSpline_Mv20(int n, const double * A, const double * b, double * c);
void ytBSpline_Mv20_mask(int n, const double * A, const double * b, double * c, const char * mask);
void ytBSpline_Bx3(const double * restrict B, const int * index,
                   const double * x, const int n, double * restrict y);
void ytBSpline_BM3(const double * restrict B, const int * index,
           const double * X, const int n, const int ndx, const int L,
           double * restrict Y);
double ytBSpline_var_mask(int n, const double * T, const double * g, 
                          const char * mask);
double ytBSpline_gKgL(int M, const double * K, const double * g);
double ytBSpline_gKg20L(const double * K, const double * g);
void ytBSpline_printB(FILE * fp, const double * B, const int * index,
                      int n, int M, int deg);
double ytBSpline_gamma20(double x, double xl, double xr, double * gamma);
#endif /* __YTLIB_BSPLINE_H */