grid_Inv_Grow - fabiankindermann/ce-fortran GitHub Wiki

function grid_Inv_Grow(x, left, right, growth, n)

Inverts the formula

for gridpoints with growing distance on the interval , such that

If is equal to a particular gridpoint created by the subroutine grid_Cons_Grow, then the return value of this function is the integer index . If lies in between two gridpoints values, the function will return a value in between the two integer indices of the gridpoints according to the above formula.

• real*8 :: x    or    real*8 :: x(:)
  The gridpoint of which to calculate the index, corresponds to . This can be either a scalar or a one-dimensional array containing many different gridpoints.
• real*8 :: left
  The left interval endpoint, corresponds to .
• real*8 :: right
  The right interval endpoint, corresponds to .
• real*8 :: growth
  The growth rate of gridpoints, corresponds to .
• integer :: n
  The number of gridpoints to use, corresponds to .

• real*8 :: grid_Inv_Grow    or    real*8 :: grid_Inv_Grow(:)
  The index of the gridpoint i. This is a scalar if x is a scalar and a one-dimensional array if x is an array.

• For further reading refer to:
  • Süli, E. & Mayers, D. F. (2003). An Introduction to Numerical Analysis. Cambridge: Cambridge University Press.
  • Powell, M.J.D. (1996). Approximation Theory and Methods. Cambridge: Cambridge University Press.