ply_l2p_header_module Module

Parameters for the plain L2 projection method to transform between Legendre modes and nodal representation.

This method utilizes the L2 projection from Legendre to Lagrange polynomials or the other way around. A numerical Gauss-Legendre Quadrature is used to compute the integral over the product of both functions. The Lagrange polynomials can be defined on any nodeset, see also ply_nodeset_module.

Available options for the nodes to project onto are:

  • 'gauss-legendre' these are the Gauss-Legendre integration points that are also used for the numerical integration (this is the default).
  • 'chebyshev' these are the nodes from the Chebyshev integration.

The set of nodes to use is configured by the nodes_kind option, and if nodes_kind = 'chebyshev' it is also possible to make use of Lobatto points to include the interval boundaries in the nodal representation. This is achieved by setting lobattoPoints = true, by default this is false.

The configuration table for a projection with L2P may, for example, look as follows:

  projection = {
    kind = 'l2p',
    factor = 1.5,
    nodes_kind = 'chebyshev',
    lobattoPoints = true
  }

The example illustrates the three possible settings for the L2P transformation method:

  • factor - Oversampling factor to avoid aliasing.
  • nodes_kind - Selection of set of nodes to use in the nodal representation.
  • lobattoPoints - Whether to include interval bounds, only available for Chebyshev nodes.

Uses

  • module~~ply_l2p_header_module~~UsesGraph module~ply_l2p_header_module ply_l2p_header_module module~aotus_module aotus_module module~ply_l2p_header_module->module~aotus_module module~aot_out_module aot_out_module module~ply_l2p_header_module->module~aot_out_module module~tem_logging_module tem_logging_module module~ply_l2p_header_module->module~tem_logging_module module~tem_aux_module tem_aux_module module~ply_l2p_header_module->module~tem_aux_module module~tem_float_module tem_float_module module~ply_l2p_header_module->module~tem_float_module module~ply_nodes_header_module ply_nodes_header_module module~ply_l2p_header_module->module~ply_nodes_header_module module~tem_tools_module tem_tools_module module~ply_l2p_header_module->module~tem_tools_module module~env_module env_module module~ply_l2p_header_module->module~env_module module~ply_nodes_header_module->module~env_module

Used by

  • module~~ply_l2p_header_module~~UsedByGraph module~ply_l2p_header_module ply_l2p_header_module module~ply_l2p_module ply_l2p_module module~ply_l2p_module->module~ply_l2p_header_module module~ply_prj_header_module ply_prj_header_module module~ply_prj_header_module->module~ply_l2p_header_module module~ply_poly_project_module ply_poly_project_module module~ply_poly_project_module->module~ply_l2p_module module~ply_poly_project_module->module~ply_prj_header_module module~ply_dynarray_project_module ply_dynarray_project_module module~ply_dynarray_project_module->module~ply_prj_header_module

Contents


Interfaces

public interface assignment(=)

public interface operator(==)

  • private pure function isEqual(left, right) result(equality)

    This function provides the test for equality of two projections.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    type(ply_l2p_header_type), intent(in) :: left

    projection to compare

    type(ply_l2p_header_type), intent(in) :: right

    projection to compare against

    Return Value logical

    is equal??

public interface operator(/=)

  • private pure function isUnequal(left, right) result(unequality)

    This function provides the test for unequality of two projections.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    type(ply_l2p_header_type), intent(in) :: left

    projection to compare

    type(ply_l2p_header_type), intent(in) :: right

    projection to compare against

    Return Value logical

    is unequal??

public interface operator(<)

  • private pure function isSmaller(left, right) result(small)

    This function provides a < comparison of two projections.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    type(ply_l2p_header_type), intent(in) :: left

    projection to compare

    type(ply_l2p_header_type), intent(in) :: right

    projection to compare against

    Return Value logical

    is smaller??

public interface operator(<=)

  • private pure function isSmallerOrEqual(left, right) result(small)

    This function provides a <= comparison of two projections.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    type(ply_l2p_header_type), intent(in) :: left

    projection to compare

    type(ply_l2p_header_type), intent(in) :: right

    projection to compare against

    Return Value logical

    is smaller??

public interface operator(>)

  • private pure function isGreater(left, right) result(great)

    This function provides a > comparison of two projections.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    type(ply_l2p_header_type), intent(in) :: left

    projection to compare

    type(ply_l2p_header_type), intent(in) :: right

    projection to compare against

    Return Value logical

    is greater??

public interface operator(>=)

  • private pure function isGreaterOrEqual(left, right) result(great)

    This function provides a >= comparison of two projections.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    type(ply_l2p_header_type), intent(in) :: left

    projection to compare

    type(ply_l2p_header_type), intent(in) :: right

    projection to compare against

    Return Value logical

    is greater??


Derived Types

type, public :: ply_l2p_header_type

l2p projection header type, consisting of the node header which give information about the type and number of points for the projection

Components

TypeVisibilityAttributesNameInitial
type(ply_nodes_header_type), private :: nodes_header
real(kind=rk), private :: factor

Functions

private pure function isEqual(left, right) result(equality)

This function provides the test for equality of two projections.

Read more…

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(in) :: left

projection to compare

type(ply_l2p_header_type), intent(in) :: right

projection to compare against

Return Value logical

is equal??

private pure function isUnequal(left, right) result(unequality)

This function provides the test for unequality of two projections.

Read more…

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(in) :: left

projection to compare

type(ply_l2p_header_type), intent(in) :: right

projection to compare against

Return Value logical

is unequal??

private pure function isSmaller(left, right) result(small)

This function provides a < comparison of two projections.

Read more…

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(in) :: left

projection to compare

type(ply_l2p_header_type), intent(in) :: right

projection to compare against

Return Value logical

is smaller??

private pure function isSmallerOrEqual(left, right) result(small)

This function provides a <= comparison of two projections.

Read more…

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(in) :: left

projection to compare

type(ply_l2p_header_type), intent(in) :: right

projection to compare against

Return Value logical

is smaller??

private pure function isGreater(left, right) result(great)

This function provides a > comparison of two projections.

Read more…

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(in) :: left

projection to compare

type(ply_l2p_header_type), intent(in) :: right

projection to compare against

Return Value logical

is greater??

private pure function isGreaterOrEqual(left, right) result(great)

This function provides a >= comparison of two projections.

Read more…

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(in) :: left

projection to compare

type(ply_l2p_header_type), intent(in) :: right

projection to compare against

Return Value logical

is greater??


Subroutines

public subroutine ply_l2p_header_load(me, conf, thandle)

Load settings to describe a projection method from a Lua table.

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(out) :: me
type(flu_State), intent(inout) :: conf
integer, intent(in) :: thandle

public subroutine ply_l2p_header_define(me, factor, nodes_kind, lobattoPoints)

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(out) :: me

L2P header to define.

integer, intent(in), optional :: factor

Oversampling factor to use in the projection, defaults to 1.

character(len=*), optional :: nodes_kind

Set of nodes to use in the nodal representation.

May be 'gauss-legendre' or 'chebyshev', defaults to 'gauss-legendre'

logical, intent(in), optional :: lobattoPoints

Wether to use Lobatto points (include interval bounds) when using the chebyshev nodes, defaults to .false..

public subroutine ply_l2p_header_out(me, conf)

Write L2P settings into a Lua table.

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(in) :: me
type(aot_out_type), intent(inout) :: conf

public subroutine ply_l2p_header_display(me)

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(in) :: me

private pure subroutine Copy_l2p_header(left, right)

Arguments

TypeIntentOptionalAttributesName
type(ply_l2p_header_type), intent(out) :: left

fpt to copy to

type(ply_l2p_header_type), intent(in) :: right

fpt to copy from