ExpMesh

Generates exponential mesh.

Interface

܀ Interface

INTERFACE ExpMesh
  MODULE PURE FUNCTION ExpMesh(rmin, rmax, a, N) RESULT(Ans)
    REAL(REAL32), INTENT(IN) :: rmin
    !! left end of 1D domain
    REAL(REAL32), INTENT(IN) :: rmax
    !! right end of 1D domain
    REAL(REAL32), INTENT(IN) :: a
    !! Ratio of largest to smallest element, a should be positive
    !! a = 1, then we get uniform mesh
    INTEGER(I4B), INTENT(IN) :: N
    !! Number of elements present in mesh
    REAL(REAL32) :: ans(N + 1)
    !! Number of nodes in mesh
  END FUNCTION ExpMesh
END INTERFACE ExpMesh

️܀ See example

program main
use easifemBase
implicit none

call Display(expmesh(rmin=1.0, rmax=100.0, a=10.0, N=10), "ans: ", full=.true. )
end program main

See results

results

 ans:  
-------
  1.000
  3.422
  6.551
 10.592
 15.810
 22.551
 31.256
 42.500
 57.021
 75.777
100.000

↢