../../lib/numru/derivative.rb

# module NumRu::Derivative in derivative.rb

## todo

- decide argument of b_expand_linear_ext is Symbol or Numeric.
  - now is Numeric.
  - it denpends the treatment of dRuby.
- support other boundary conditions.

## Index

- module NumRu::Derivative
  - threepoint_O2nd_deriv
    - First derivative (2nd Order difference use three point.)
  - cderiv
    - First derivative (center difference use two point.)
  - b_expand_linear_ext
    - return array extended boundaries with linear extention.
  - cdiff
    - return difference. (center difference)

# module NumRu::Derivative

Module functions of Derivative Operater for NArray.

`threepoint_O2nd_deriv(z, x, dim, bc=LINEAR_EXT)`

  Derivate `z` respect to `dim` th dimension with 2nd Order difference. return an NArray which result of the difference z divided difference `x` (in other wards,

      (s**2*z_{i+1} + (t**2 - s**2)*f_{i} - t**2*f_{i-1}) / (s*t*(s + t)):
      now s represents (x_{i} - x_{i-1}) ,t represents (x_{i+1} - x_{i})
      and _{i} represents the suffix of {i} th element in the ((<dim>)) th
      dimension of array. ).

  ARGUMENTS

  - z (NArray): a NArray which you want to derivative.
  - x (NArray): a NArray represents the dimension which derivative respect to. z.rank must be 1.
  - dim (Numeric): a Numeric represents the dimention which derivative respect to. you can give number count backward (dim<0), but z.rank ¡Üdim must be > 0.
  - bc (Numeric) : a Numeric which represent boundary condition. now only LINEAR_EXT(=1) supported. LINEAR_EXT load b_expand_linear_ext which extend boundary with lenear value.

  RETURN VALUE

  - O2nd_deriv_data (NArray): (s**2*z_{i+1} + (t**2 - s**2)*f_{i} - t**2*f_{i-1}) / (s*t*(s + t))
`cderiv(z, x, dim, bc=LINEAR_EXT)`

  Derivate `z` respect to `dim` th dimension with center difference. return an NArray which result of the difference z divided difference `x` ( in other wards, (z_{i+1} - z_{i-1}) / (x_{i+1} - x_{i-1}): now _{i} represents the suffix of {i} th element in the dim th dimension of array. ).

  ARGUMENTS

  - z (NArray): a NArray which you want to derivative.
  - x (NArray): a NArray represents the dimension which derivative respect to. z.rank must be 1.
  - dim (Numeric): a Numeric represents the dimention which derivative respect to. you can give number count backward (dim<0), but z.rank ¡Üdim must be > 0.
  - bc (Numeric) : a Numeric which represent boundary condition. now only LINEAR_EXT(=1) supported. LINEAR_EXT load b_expand_linear_ext which extend boundary with lenear value.

  RETURN VALUE

  - cderiv_data (NArray): (z_{i+1} - z_{i-1}) / (x_{i+1} - x_{i-1})
`b_expand_linear_ext(z, dim)`

  expand boundary with linear value. extend array with 1 grid at each boundary with dim th dimension, and assign th value which diffrential value between a grid short of boundary and boundary grid in original array. (on other wards, 2*z_{0}-z_{1} or 2*z_{n-1}-z_{n-2}: now _{i} represents the

      suffix of {i} th element in the ((<dim>)) th dimension of array. ).

  ARGUMENTS

  - z (NArray): a NArray which you want to expand boundary.
  - dim (Numeric): a Numeric represents the dimention which derivative respect to. you can give number count backward (dim<0), but z.rank ¡Üdim must be > 0.

  RETURN VALUE

  - expand_data (NArray):
`cdiff(x, dim)`

  Diffrence operater. return an NArray which a difference x ( in other wards, (x_{i+1} - x_{i-1}): now _{i} represents the suffix of {i} th element in the dim th dimension of array. ).

  ARGUMENTS

  - x (NArray): a NArray which you want to get difference.
  - dim (Numeric): a Numeric representing the dimention which derivative respect to. you can give number count backward (dim<0), but z.rank ¡Üdim must be > 0.

  RETURN VALUE

  - cdiff_data (NArray): (x_{i+1} - x_{i-1})