Hi,
I am trying to write a numeric code following Sloan's (Article A fast algorithm for constructing Delaunay triangulations in the plane
) article. He mentions that it is only necessary a set of N points (x,y) distributed in a plane.I elaborated the main program
program delau
!implicit real*8 (a-h,o-z)
parameter (n = 100)
INTEGER I,LIST(n+3),V(3,n+3),E(3,n+3),NUMTRI,BIN(n+3),NUMPTS
real X(n+3),Y(n+3)
rr = rand(78)
NUMPTS = n
do i =1,n
x(i) =rand()
y(i) =rand()
enddo
call deltri(numpts,n,x,y,list,bin,v,e,numtri)
end
the subroutines were copied and pasted from Sloan's article. But my code only triangulates 2 vertices and the EDG function displays
***ERROR IN FUNCTION EDG***
***ELEMENTS NOT ADJACENT***
Where am I going wrong??
best regards
Wanderson Wanzeller You appear to be attempting to implement the Delaunay triangulation procedure from Sloan's essay in Fortran 77, but are encountering difficulties. It is difficult to pinpoint the source of the problem without further information about the precise subroutines you are running and the error message you are receiving. However, here are a few things to think about that can help you fix the problem:
- Ascertain that the input data (the points x, y) are true and valid. Check for duplicate points and that the points are not collinear.
- Check that the subroutines are appropriately called and that the input and output variables are accurately supplied.
- Before calling deltri, ensure that the array dimensions are accurate and that the arrays are initialized with proper values.
- Check that the subroutines are properly built and linked.
- Examine the subroutines for any typographical mistakes or inconsistencies.
Compare your implementation to a reference implementation of the algorithm to determine if there are any disparities.
It's worth noting that the paper was published in 1986, and since then, computer science has produced many new techniques and libraries for Delaunay triangulation. In current Fortran, libraries such as Qhull, CGAL, and scipy.spatial are commonly used for this job.
I hope these solutions assist you in troubleshooting your issue. Please let me know if the situation persists, and I will do my best to assist you further.
Jay Patel the proposal is good but I don't have a license for MATLAB.
- Badges
- Science topic