1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
|
include <mach.h>
include <ctype.h>
include <error.h>
include <gset.h>
include "pdm.h"
define MSIZE 2.0 # Mark size
# PDM_DELETE -- Delete the point from the plot and set it's inuse array
# entry to zero.
int procedure pdm_delete (pdmp, cx, cy, ptype)
pointer pdmp # pointer to PDM data structure
double cx, cy # device cursor coordinates
int ptype # plot type
pointer gp
real x, y
int npts, i, j, index
real x0, y0, r2, r2min
begin
gp = PDM_GP(pdmp)
npts = PDM_NPT(pdmp)
if (ptype == DATAPLOT) {
# Transform world cursor coordinates to NDC.
call gctran (gp, real(cx), real(cy), x, y, 1, 0)
# Search for nearest point in-use.
j = 0
r2min = MAX_REAL
do i = 1, npts {
if (PDM_INUSE(pdmp,i) == 0)
next
call gctran (gp, real(PDM_X(pdmp,i)), real(PDM_DY(pdmp,i)),
x0, y0, 1, 0)
r2 = (x0 - x) ** 2 + (y0 - y) ** 2
if (r2 < r2min) {
r2min = r2
j = i
}
}
# Mark the deleted point with a cross and set the weight to zero.
if (j != 0) {
call gscur (gp, real(PDM_X(pdmp,j)), real(PDM_DY(pdmp,j)))
call gmark (gp, real(PDM_X(pdmp,j)), real(PDM_DY(pdmp,j)),
GM_CROSS, MSIZE, MSIZE)
PDM_INUSE(pdmp,j) = 0
call gflush(gp)
}
return (j)
} else if (ptype == PHASEPLOT) {
# Transform world cursor coordinates to NDC.
call gctran (gp, real(cx), real(cy), x, y, 1, 0)
# Search for nearest point in-use.
j = 0
r2min = MAX_REAL
do i = 1, npts {
index = PDM_SORT(pdmp,i)
if (PDM_INUSE(pdmp,index) == 0)
next
call gctran (gp, real(PDM_XPH(pdmp,i)), real(PDM_YPH(pdmp,i)),
x0, y0, 1, 0)
r2 = (x0 - x) ** 2 + (y0 - y) ** 2
if (r2 < r2min) {
r2min = r2
j = i
}
}
# Mark the deleted point with a cross and set the weight to zero.
# We mark two points since all of the points are displayed twice
# on the phase plot.
if (j != 0) {
call gscur (gp, real(PDM_XPH(pdmp,j))+1.0,
real(PDM_YPH(pdmp,j)))
call gmark (gp, real(PDM_XPH(pdmp,j))+1.0,
real(PDM_YPH(pdmp,j)), GM_CROSS, MSIZE, MSIZE)
call gscur (gp, real(PDM_XPH(pdmp,j)), real(PDM_YPH(pdmp,j)))
call gmark (gp, real(PDM_XPH(pdmp,j)), real(PDM_YPH(pdmp,j)),
GM_CROSS, MSIZE, MSIZE)
# Calculate which point this corresponds to.
index = PDM_SORT(pdmp,j)
PDM_INUSE(pdmp,index) = 0
call gflush (gp)
}
return (index)
} else
return (0)
end
|
