Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/glentner/cubicsplinelib
C++ facility for natural cubic spline interpolation in 1D, 2D, and 3D.
https://github.com/glentner/cubicsplinelib
Last synced: 6 days ago
JSON representation
C++ facility for natural cubic spline interpolation in 1D, 2D, and 3D.
- Host: GitHub
- URL: https://github.com/glentner/cubicsplinelib
- Owner: glentner
- License: gpl-2.0
- Created: 2015-01-22T02:36:27.000Z (almost 10 years ago)
- Default Branch: master
- Last Pushed: 2015-02-10T21:42:27.000Z (almost 10 years ago)
- Last Synced: 2024-05-01T15:26:43.331Z (7 months ago)
- Language: C++
- Size: 637 KB
- Stars: 33
- Watchers: 4
- Forks: 10
- Open Issues: 1
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
Cubic Spline Library
====================
This small C++ library contains my working code for performing
natural cubic spline interpolation in one, two, and three
dimensions.
**dependencies:** OpenMP (omp.h)
The CubicSpline class is fully functional. A test program is
available in *CubicSpineTest.cc*, the results of which I
have plotted in *CubicSplineTest.png* using Python.
BiCubicSpline and TriCubicSpline are still being developed.
Documentation and usage will be provided in the future.
##Building
You should be able to just run 'make'. This compiles and links
the programs. The user will either need to make this available
on their standard include path or link to this directory.
I recommend something like:
```
export CSLDIR=/path/to/CubicSplineLib
export CSLINC="-I$CSLDIR -L$CSLDIR -lCubicSpline -fopenmp"
```
After which, programs can be compiled similar to:
```bash
g++ your_program.cc $CSLINC
```
##Short Example Usage
```C++
// create spline object
CubicSpline raw_data( x, y, 2 );
// resample 'x' data
std::vector new_y = raw_data.interpolate( new_x );
```
The user can construct the spline object with float,
double, or long double types. The third argument in
the constructor is optional and has a default value
of one. This says what the maximum allowed threads
should be for OpenMP. There is a fourth optional argument
that takes either true or false. By default this
implementation uses a quicksort algorithm to ensure
the horizontal axis is in ascending order before
building the spline polynomials. This behavior can
be supressed by specifying false. Also, in the
interpolation function, an optional second and third
argument taking either true or false and a numeric
type the same as the template respectively allow
for returning fixed values outside the domain of the
original data.
e.g.,
```C++
raw_data.interpolate( new_x, true, 0.0 );
```
would return a std::vector containing 0.0 for all positions
outside the original domain of 'x'.