comparison splines.cpp @ 5:5f3c32dc6e17

* Adjust comment syntax to permit Doxygen to generate HTML documentation; add Doxyfile

 //---------------------------------------------------------------------------
 
-
 #include "splines.h"
 
+/** \file splines.h */
+
 //---------------------------------------------------------------------------
-/*
+/**
   function tridiagonal: solves linear system A[N][N]x[N]=d[N] where A is tridiagonal
 
   In: tridiagonal matrix A[N][N] gives as three vectors - lower subdiagonal
         a[1:N-1], main diagonal b[0:N-1], upper subdiagonal c[0:N-2]
       vector d[N]

   d[N-1]=d[N-1]/b[N-1];
   for (int k=N-2; k>=0; k--)
     d[k]=(d[k]-c[k]*d[k+1])/b[k];
 }//tridiagonal
 
-/*
+/**
   function CubicSpline: computing piece-wise trinomial coefficients of a cubic spline
 
   In: x[N+1], y[N+1]: knots
       bordermode: boundary mode of cubic spline
         bordermode=0: natural: y''(x0)=y''(x_N)=0

       data[i]=d[n]+xx*(c[n]+xx*(b[n]+xx*a[n]));
     }
   }
 }//CubicSpline
 
-/*
+/**
   function CubicSpline: computing piece-wise trinomial coefficients of a cubic spline with uniformly placed knots
 
   In: y[N+1]: spline values at knots (0, h, 2h, ..., Nh)
       bordermode: boundary mode of cubic spline
         bordermode=0: natural: y''(x0)=y''(x_N)=0

       data[i]=d[n]+xx*(c[n]+xx*(b[n]+xx*a[n]));
     }
   }
 }//CubicSpline
 
-/*
+/**
   function Smooth_Interpolate: smoothly interpolate/extrapolate P-piece spline from P-1 values. The
   interpolation scheme is chosen according to P:
     P-1=1:  constant
     P-1=2:  linear function
     P-1=3:  quadratic function