A subroutine to test the tem_polygon_material_value function in tem_polygon_material_test.
The only argument will be true, if all calculations result in the expected values.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| logical, | intent(out) | :: | success |
Indicator if all tests were computed correctly. |
subroutine tem_polygon_material_test_value(success) ! ---------------------------------------------------------------------- !> Indicator if all tests were computed correctly. logical, intent(out) :: success ! ---------------------------------------------------------------------- type(tem_polygon_material_type) :: polygon type(tem_polygon_vertex_type) :: me real(kind=rk) :: matval(1) ! ---------------------------------------------------------------------- ! Use simple polygon with 5 vertices to check the containment function. ! It has roughly this shape: ! ! x-------x ! \ | ! \ | ! x | ! / | ! / | ! x-------x ! ! With the 'x' indicating the 5 vertices. me%nVertices = 5 allocate(me%Vertex(me%nVertices,2)) me%Vertex(1,:) = [ 0.0_rk, 0.0_rk] me%Vertex(2,:) = [-1.0_rk, -1.0_rk] me%Vertex(3,:) = [ 1.0_rk, -1.0_rk] me%Vertex(4,:) = [ 1.0_rk, 1.0_rk] me%Vertex(5,:) = [-1.0_rk, 1.0_rk] ! Defining returned values for inside and outside the polygon: polygon%nComponents = 1 allocate(polygon%inval(polygon%nComponents)) allocate(polygon%outval(polygon%nComponents)) polygon%inval = 1.0_rk polygon%outval = 0.0_rk success = .true. ! (If any subsequent test fails, we set this to false.) write(*,*) 'Checking point containment for this 5 point polygon:' write(*,*) write(*,*) '(-1,1)---(1,1)' write(*,*) ' \ |' write(*,*) ' \ |' write(*,*) ' (0,0) |' write(*,*) ' / |' write(*,*) ' / |' write(*,*) '(-1,-1)--(1,-1)' write(*,*) ! Check a point inside the polygon: matval = tem_polygon_material_value( me = me, & & point = [0.5_rk, 0.5_rk], & & inval = polygon%inval, & & outval = polygon%outval, & & nComponents = polygon%nComponents ) if (matval(1) < 1.0_rk) success = .false. write(*,*) 'Point (0.5, 0.5) has value:', matval ! Check a point outside the polygon: matval = tem_polygon_material_value( me = me, & & point = [-1.0_rk, 0.0_rk], & & inval = polygon%inval, & & outval = polygon%outval, & & nComponents = polygon%nComponents ) if (matval(1) > epsilon(1.0_rk)) success = .false. write(*,*) 'Point (-1.0, 0.0) has value:', matval ! Check a point on a side of the polygon (should be considered inside): matval = tem_polygon_material_value( me = me, & & point = [-0.5_rk, -0.5_rk], & & inval = polygon%inval, & & outval = polygon%outval, & & nComponents = polygon%nComponents ) if (matval(1) < 1.0_rk) success = .false. write(*,*) 'Point (-0.5, -0.5) has value:', matval ! Check an outward pointing corner (should be outside): matval = tem_polygon_material_value( me = me, & & point = [1.0_rk, 1.0_rk], & & inval = polygon%inval, & & outval = polygon%outval, & & nComponents = polygon%nComponents ) if (matval(1) > epsilon(1.0_rk)) success = .false. write(*,*) 'Point (1.0, 1.0) has value:', matval ! Check an inward pointing corner (should be inside): matval = tem_polygon_material_value( me = me, & & point = [0.0_rk, 0.0_rk], & & inval = polygon%inval, & & outval = polygon%outval, & & nComponents = polygon%nComponents ) if (matval(1) < 1.0_rk) success = .false. write(*,*) 'Point (0.0, 0.0) has value:', matval end subroutine tem_polygon_material_test_value