GeometryNodes1D Subroutine

  1. SELF_Mesh_3D_t.f90
  2. SELF_Mesh_3D_t
  3. GeometryNodes1D

private subroutine GeometryNodes1D(nGeo, quadrature, nodes)

Returns the nGeo+1 one-dimensional geometry node positions, in [-1,1], of a mesh whose node coordinates are sampled on quadrature points. These are the points GenerateFromMesh interpolates the node coordinates from (see Lagrange_t%Init), so a coordinate sampled here is reproduced exactly by the geometry interpolant.

Arguments

TypeIntentOptionalAttributesName
integer, intent(in) :: nGeo
integer, intent(in) :: quadrature
real(kind=prec), intent(out) :: nodes(1:nGeo+1)

Calls

proc~~geometrynodes1d~~CallsGraph proc~geometrynodes1d GeometryNodes1D proc~legendrequadrature LegendreQuadrature proc~geometrynodes1d->proc~legendrequadrature proc~chebyshevquadrature ChebyshevQuadrature proc~geometrynodes1d->proc~chebyshevquadrature proc~uniformpoints UniformPoints proc~geometrynodes1d->proc~uniformpoints proc~legendregausslobatto LegendreGaussLobatto proc~legendrequadrature->proc~legendregausslobatto proc~legendregauss LegendreGauss proc~legendrequadrature->proc~legendregauss proc~chebyshevgausslobatto ChebyshevGaussLobatto proc~chebyshevquadrature->proc~chebyshevgausslobatto proc~chebyshevgauss ChebyshevGauss proc~chebyshevquadrature->proc~chebyshevgauss proc~legendreqandl LegendreQandL proc~legendregausslobatto->proc~legendreqandl proc~legendrepolynomial LegendrePolynomial proc~legendregauss->proc~legendrepolynomial

Called by

proc~~geometrynodes1d~~CalledByGraph proc~geometrynodes1d GeometryNodes1D proc~extrudedmesh_mesh3d_t ExtrudedMesh_Mesh3D_t proc~extrudedmesh_mesh3d_t->proc~geometrynodes1d

Contents

Source Code

GeometryNodes1D

Source Code

  subroutine GeometryNodes1D(nGeo,quadrature,nodes)
  !! Returns the nGeo+1 one-dimensional geometry node positions, in [-1,1], of a mesh
  !! whose node coordinates are sampled on `quadrature` points. These are the points
  !! GenerateFromMesh interpolates the node coordinates from (see Lagrange_t%Init), so
  !! a coordinate sampled here is reproduced exactly by the geometry interpolant.
    implicit none
    integer,intent(in) :: nGeo
    integer,intent(in) :: quadrature
    real(prec),intent(out) :: nodes(1:nGeo+1)
    ! Local
    real(prec) :: weights(1:nGeo+1)

    if(quadrature == GAUSS .or. quadrature == GAUSS_LOBATTO) then
      call LegendreQuadrature(nGeo,nodes,weights,quadrature)
    elseif(quadrature == CHEBYSHEV_GAUSS .or. quadrature == CHEBYSHEV_GAUSS_LOBATTO) then
      call ChebyshevQuadrature(nGeo,nodes,weights,quadrature)
    elseif(quadrature == UNIFORM) then
      nodes = UniformPoints(-1.0_prec,1.0_prec,0,nGeo)
    else
      error stop "GeometryNodes1D : unrecognized mesh quadrature type"
    endif

  endsubroutine GeometryNodes1D