UnscaledLobattoStiffnessMatrix

Evaluates the Stiffness matrix for UnscaledLobatto polynomial.

Interface

܀ Interface

INTERFACE
  MODULE PURE FUNCTION UnscaledLobattoStiffnessMatrix(n) RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: n
    REAL(DFP) :: ans(n + 1, n + 1)
  END FUNCTION UnscaledLobattoStiffnessMatrix
END INTERFACE

️܀ See example

This example shows the usage of UnscaledLobattoStiffnessMatrix method.

This routine evaluates stiffness matrix for lobatto polynomials.

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

  n = 0
  call callme

0.5

  n = 1
  call callme

0.5	-0.5
-0.5	0.5

  n = 4
  call callme

0.5	-0.5	0	0	0
-0.5	0.5	0	0	0
0	0	0.66667	0	0
0	0	0	0.4	0
0	0	0	0	0.28571

  contains
  subroutine callme
    ans= UnscaledLobattoStiffnessMatrix( n=n )
    astr = MdEncode(ans)
    CALL Display( "UnscaledLobatto mass matrix = " )
    CALL Display(astr, "" )
    call blanklines(nol=2)
  end subroutine callme

end program main

↢