LobattoStiffnessMatrix

Evaluates the Stiffness matrix for Lobatto polynomial.

Interface

܀ Interface

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

️܀ See example

This example shows the usage of LobattoStiffnessMatrix 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	1	0	0
0	0	0	1	0
0	0	0	0	1

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

end program main

↢