ply_split_element_1D_test Subroutine

private subroutine ply_split_element_1D_test(nModes, success)

Testing the 1D splitting.

We test all combinations, even though a higher number of modes in the children is probably not that relevant, it should still be possible.

Arguments

TypeIntentOptionalAttributesName
integer, intent(in) :: nModes

Number of modes in the (1D) polynomials to use in the check.

logical, intent(out) :: success

Indication whether the tests were completed successfully.


Calls

proc~~ply_split_element_1d_test~~CallsGraph proc~ply_split_element_1d_test ply_split_element_1D_test proc~ply_split_element_init ply_split_element_init proc~ply_split_element_1d_test->proc~ply_split_element_init proc~ply_legendre_1d ply_legendre_1D proc~ply_split_element_1d_test->proc~ply_legendre_1d proc~ply_split_element_singled ply_split_element_singleD proc~ply_split_element_1d_test->proc~ply_split_element_singled proc~ply_split_legendre_matrix ply_split_legendre_matrix proc~ply_split_element_init->proc~ply_split_legendre_matrix proc~alpha_beta alpha_beta proc~ply_split_legendre_matrix->proc~alpha_beta proc~beta beta proc~ply_split_legendre_matrix->proc~beta proc~alpha alpha proc~ply_split_legendre_matrix->proc~alpha proc~alpha_frac alpha_frac proc~ply_split_legendre_matrix->proc~alpha_frac

Called by

proc~~ply_split_element_1d_test~~CalledByGraph proc~ply_split_element_1d_test ply_split_element_1D_test proc~ply_split_element_test ply_split_element_test proc~ply_split_element_test->proc~ply_split_element_1d_test

Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: parentModes
integer, private :: childmodes
integer, private :: iPoint
integer, private :: iElem
real(kind=rk), private :: xi(nModes)
real(kind=rk), private :: x_left(nModes)
real(kind=rk), private :: x_right(nModes)
real(kind=rk), private :: legchild(nModes,nModes)
real(kind=rk), private :: legleft(nModes,nModes)
real(kind=rk), private :: legright(nModes,nModes)
real(kind=rk), private :: rootval(nModes,2)
real(kind=rk), private :: childval
real(kind=rk), private, allocatable:: rootelem(:,:)
real(kind=rk), private, allocatable:: childelem(:,:)
real(kind=rk), private :: tolerance