EquidistanceInPoint

This function returns the equidistance points in triangle.

All points are inside the triangle

Interface

܀ Interface

INTERFACE
  MODULE PURE FUNCTION EquidistanceInPoint_Quadrangle(order, xij) RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: order
  !! order
    REAL(DFP), OPTIONAL, INTENT(IN) :: xij(:, :)
  !! coordinates of point 1 and point 2 in $x_{iJ}$ format
  !! number of rows = nsd
  !! number of cols = 3
    REAL(DFP), ALLOCATABLE :: ans(:, :)
  !! returned coordinates in $x_{iJ}$ format
  END FUNCTION EquidistanceInPoint_Quadrangle
END INTERFACE

️܀ See example

program main
  use easifembase
  implicit none
  integer( i4b ) :: i1, i2, order
  real( dfp ), allocatable :: x(:,:)
  order=1
  x = EquidistanceInPoint_Quadrangle( order=order )
  call display( TRANSPOSE(x), "xij (order="//tostring(order)//")=" )
  call blanklines(nol=2)
  order=2
  x = EquidistanceInPoint_Quadrangle( order=order )
  call display( TRANSPOSE(x), "xij (order="//tostring(order)//")=" )
  call blanklines(nol=2)
  order=3
  x = EquidistanceInPoint_Quadrangle( order=order )
  call display( TRANSPOSE(x), "xij (order="//tostring(order)//")=" )
  call blanklines(nol=2)
end program main

See results

xij (order=1)=
--------------
 
 
 xij (order=2)= 
----------------
0.00000  0.00000
 
 
   xij (order=3)=   
--------------------
-0.333333  -0.333333
 0.333333  -0.333333
 0.333333   0.333333
-0.333333   0.333333

↢