splines Module Source file:subroutines

Dependencies

Subroutines
  • splder(n, y, g, dy)
    Splder fits a cubic spline to y and returns the first derivatives at the grid points in dy. Dy is equivalent to a 4th-order Pade difference formula for dy/di.
    • integer intent(in) :: n
    • real(dl) intent(in) :: y(n)
    • real(dl) intent(in) :: g(n)
    • real(dl) intent(out) :: dy(n)
  • spline_def(n, x, y, d2)
    Low-level initialize spline arrays with default boundary conditions
    • integer intent(in) :: n
    • real(sp_acc) intent(in) :: x(n)
    • real(sp_acc) intent(in) :: y(n)
    • real(sp_acc) intent(out) :: d2(n)
  • spline_deriv(n, x, y, y2, y1)
    Get derivative y1 given array of x, y and y''
    • INTEGER intent(in) :: n
    • real(dl) intent(in) :: x(n)
    • real(dl) intent(in) :: y(n)
    • real(dl) intent(in) :: y2(n)
    • real(dl) intent(out) :: y1(n)
  • spline_integrate(n, x, y, y2, yint)
    Cumulative integral of cubic spline
    • integer intent(in) :: n
    • real(dl) intent(in) :: x(n)
    • real(dl) intent(in) :: y(n)
    • real(dl) intent(in) :: y2(n)
    • real(dl) intent(out) :: yint(n)
  • splini(n, g)
    Splini must be called before splder to initialize array g in common.
    • integer intent(in) :: n
    • real(dl) intent(out) :: g(n)
  • splint(n, y, z)
    Splint integrates a cubic spline, providing the output value z = integral from 1 to n of s(i)di, where s(i) is the spline fit to y(i).
    • integer intent(in) :: n
    • real(dl) intent(in) :: y(n)
    • real(dl) intent(out) :: z