GetNearestNode

Get node coordinate and global node number of nearest node of a give node in the domain.

This method uses Kdtree

Make sure you call InitiateKdtree method before calling it. If Kdtree is not built then this method will build it for you. Therefore, the first call can take longer time to execute.

info

In Interface 1 you can get a single nearest node. In the Interface 2 you can get nn nearest nodes.

Interface 1

Interface

INTERFACE
  MODULE SUBROUTINE obj_GetNearestNode1(obj, qv, x, globalNode)
    CLASS(AbstractDomain_), INTENT(INOUT) :: obj
    REAL(DFP), INTENT(IN) :: qv(:)
    !! Query vector
    REAL(DFP), INTENT(INOUT) :: x(:)
    !! node coord of nearest node
    INTEGER(I4B), INTENT(OUT) :: globalNode
    !! globalNode number
  END SUBROUTINE obj_GetNearestNode1
END INTERFACE

example

PROGRAM main
USE easifemBase
USE easifemClasses
USE FEDomain_Class

IMPLICIT NONE

TYPE(FEDomain_) :: obj
TYPE(HDF5File_) :: meshfile
CHARACTER(*), PARAMETER :: filename = &
                  "../../Mesh/examples/meshdata/small_tri3_mesh_two_region.h5"
INTEGER(I4B) :: globalNode
REAL(DFP) :: qv(3), x(3)

CALL meshfile%Initiate(filename, "READ")
CALL meshfile%OPEN()
CALL obj%Initiate(meshfile, '')

qv = [1.0, 0.0, 0.0]
CALL obj%GetNearestNode(qv=qv, x=x, globalNode=globalNode)
CALL Display(globalNode, "globalNode: ")
CALL Display(x, "x: ")

qv = [1.5, 0.3, 0.0]
CALL obj%GetNearestNode(qv=qv, x=x, globalNode=globalNode)
CALL Display(globalNode, "globalNode: ")
CALL Display(x, "x: ")

CALL meshfile%DEALLOCATE()
CALL obj%DEALLOCATE()

END PROGRAM main

↢ close

Interface 2

In this interface you can get nn nearest node.

Interface

 INTERFACE
  MODULE SUBROUTINE obj_GetNearestNode2(obj, qv, x, globalNode, nn)
    CLASS(AbstractDomain_), INTENT(INOUT) :: obj
    REAL(DFP), INTENT(IN) :: qv(:)
    !! Query vector
    REAL(DFP), INTENT(INOUT) :: x(:, :)
    !! node coord of nearest node
    !! the size(x, 2) should be atleast nn
    INTEGER(I4B), INTENT(INOUT) :: globalNode(:)
    !! globalNode number, size of globalNode should be atleast nn
    INTEGER(I4B), INTENT(IN) :: nn
    !! number of nearest points
  END SUBROUTINE obj_GetNearestNode2
END INTERFACE

example

PROGRAM main
USE easifemBase
USE easifemClasses
USE FEDomain_Class

IMPLICIT NONE

TYPE(FEDomain_) :: obj
TYPE(HDF5File_) :: meshfile
CHARACTER(*), PARAMETER :: filename = &
                  "../../Mesh/examples/meshdata/small_tri3_mesh_two_region.h5"
INTEGER(I4B) :: globalNode(5)
REAL(DFP) :: qv(3), x(3, 5)

CALL meshfile%Initiate(filename, "READ")
CALL meshfile%OPEN()
CALL obj%Initiate(meshfile, '')

qv = [1.0, 0.0, 0.0]
CALL obj%GetNearestNode(qv=qv, x=x, nn=5, globalNode=globalNode)
CALL Display(globalNode, "globalNode: ")
CALL Display(x, "x: ")

qv = [1.5, 0.3, 0.0]
CALL obj%GetNearestNode(qv=qv, x=x, nn=5, globalNode=globalNode)
CALL Display(globalNode, "globalNode: ")
CALL Display(x, "x: ")

CALL meshfile%DEALLOCATE()
CALL obj%DEALLOCATE()

END PROGRAM main

↢ close