LobattoEvalAll

Evaluate all Lobatto polynomials upto order n.

Interface 1

܀ Interface

INTERFACE
  MODULE PURE FUNCTION LobattoEvalAll(n, x) RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: n
    REAL(DFP), INTENT(IN) :: x
    REAL(DFP) :: ans(n + 1)
    !! Evaluate Lobatto polynomial of order = 0 to n (total n+1)
    !! at point x
  END FUNCTION LobattoEvalAll
END INTERFACE

️܀ See example

This example shows the usage of LobattoEvalAll method.

This routine evaluates Lobatto polynomial upto order n, for many points

program main
  use easifembase
  implicit none
  integer( i4b ) :: n
  real( dfp ), allocatable :: ans( :, : ), x( : )
  type(string) :: astr

  x = [-1.0, 0.0, 1.0]
  n = 0; call callme

P0
---
1
0.5
0

  n = 1; call callme

l0	l1
1	0
0.5	0.5
0	1

  x = [-1.0, 0.0, 1.0, 0.5]
  n = 5; call callme

l0	l1	l2	l3	l4	l5
1	0	0	0	0	0
0.5	0.5	-0.61237	-0	0.23385	0
0	1	0	0	0	0
0.25	0.75	-0.45928	-0.29646	-4.38475E-02	0.1243

  contains
  subroutine callme
    ans= LobattoEvalAll( n=n, x=x )
    astr = MdEncode( ans )
    call display( astr%chars(), "" )
  end subroutine callme

end program main

↢

Interface 2

܀ Interface

INTERFACE
  MODULE PURE FUNCTION LobattoEvalAll(n, x) RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: n
    REAL(DFP), INTENT(IN) :: x(:)
    REAL(DFP) :: ans(SIZE(x), n + 1)
    !! Evaluate Lobatto polynomial of order = 0 to n (total n+1)
    !! at point x
  END FUNCTION LobattoEvalAll
END INTERFACE

️܀ See example

This example shows the usage of LobattoEvalAll method.

This routine evaluates Lobatto polynomial upto order n, for many points.

program main
  use easifembase
  implicit none
  integer( i4b ) :: n
  real( dfp ), allocatable :: ans( : ), x
  type(string) :: astr

  x = -1.0_DFP
  n = 5; call callme

l0	l1	l2	l3	l4	l5
1	0	0	0	0	0

  x = 1.0_DFP
  n = 5; call callme

l0	l1	l2	l3	l4	l5
0	1	0	0	0	0

  x = 0.5_DFP
  n = 5; call callme

l0	l1	l2	l3	l4	l5
0.25	0.75	-0.45928	-0.29646	-4.38475E-02	0.1243

  contains
  subroutine callme
    ans= LobattoEvalAll( n=n, x=x )
    astr = MdEncode( ans )
    call display( astr%chars(), "" )
  end subroutine callme

end program main

↢