diff options
Diffstat (limited to 'noao/astutil/pdm/pdmundelete.x')
| -rw-r--r-- | noao/astutil/pdm/pdmundelete.x | 124 |
1 files changed, 124 insertions, 0 deletions
diff --git a/noao/astutil/pdm/pdmundelete.x b/noao/astutil/pdm/pdmundelete.x new file mode 100644 index 00000000..2f5f0e31 --- /dev/null +++ b/noao/astutil/pdm/pdmundelete.x @@ -0,0 +1,124 @@ +include <mach.h> +include <ctype.h> +include <error.h> +include <gset.h> +include "pdm.h" + +define MSIZE 2.0 # Mark size + +# PDM_UNDELETE -- Undelete the nearest deleted point from the plot +# and set it's inuse array entry to one (in-use). + +int procedure pdm_undelete (pdmp, cx, cy, ptype) + +pointer pdmp # pointer to PDM data structure +double cx, cy # device coordinates of point to undelete +int ptype # plot type flag + +pointer gp +real x, y +int npt, i, j, index +real x0, y0, r2, r2min + +begin + gp = PDM_GP(pdmp) + npt = 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 not in-use. + j = 0 + r2min = MAX_REAL + do i = 1, npt { + if (PDM_INUSE(pdmp,i) == 1) + 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 + } + } + + # Unmark the deleted point and reset the weight. + if (j != 0) { + call gscur (gp, real(PDM_X(pdmp,j)), real(PDM_DY(pdmp,j))) + call gseti (gp, G_PMLTYPE, GL_CLEAR) + call gmark (gp, real(PDM_X(pdmp,j)), real(PDM_DY(pdmp,j)), + GM_CROSS, MSIZE, MSIZE) + call gseti (gp, G_PMLTYPE, GL_SOLID) + if (PDM_NPT(pdmp) <= PDM_PLUSPOINT(pdmp)) + call gpmark (gp, real(PDM_X(pdmp,j)), real(PDM_DY(pdmp,j)), + 1, GM_PLUS, MSIZE, MSIZE) + else + call gpmark (gp, real(PDM_X(pdmp,j)), real(PDM_DY(pdmp,j)), + 1, GM_POINT, 1.0, 1.0) + PDM_INUSE(pdmp,j) = 1 + 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 not in-use. + j = 0 + r2min = MAX_REAL + do i = 1, npt { + index = PDM_SORT(pdmp,i) + if (PDM_INUSE(pdmp,index) == 1) + 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 + } + } + + # Unmark the deleted point and reset the weight. + if (j != 0) { + call gscur (gp, real(PDM_XPH(pdmp,j))+1.0, + real(PDM_YPH(pdmp,j))) + call gseti (gp, G_PMLTYPE, GL_CLEAR) + call gmark (gp, real(PDM_XPH(pdmp,j))+1.0, + real(PDM_YPH(pdmp,j)), GM_CROSS, MSIZE, MSIZE) + call gseti (gp, G_PMLTYPE, GL_SOLID) + if (PDM_NPT(pdmp) <= PDM_PLUSPOINT(pdmp)) + call gpmark (gp, real(PDM_XPH(pdmp,j))+1.0, + real(PDM_YPH(pdmp,j)), 1, GM_PLUS, MSIZE, MSIZE) + else + call gpmark (gp, real(PDM_XPH(pdmp,j))+1.0, + real(PDM_YPH(pdmp,j)), 1, GM_POINT, 1.0, 1.0) + + call gscur (gp, real(PDM_XPH(pdmp,j)), real(PDM_YPH(pdmp,j))) + call gseti (gp, G_PMLTYPE, GL_CLEAR) + call gmark (gp, real(PDM_XPH(pdmp,j)), real(PDM_YPH(pdmp,j)), + GM_CROSS, MSIZE, MSIZE) + call gseti (gp, G_PMLTYPE, GL_SOLID) + if (PDM_NPT(pdmp) <= PDM_PLUSPOINT(pdmp)) + call gpmark (gp, real(PDM_XPH(pdmp,j)), + real(PDM_YPH(pdmp,j)), 1, GM_PLUS, MSIZE, MSIZE) + else + call gpmark (gp, real(PDM_XPH(pdmp,j)), + real(PDM_YPH(pdmp,j)), 1, GM_POINT, 1.0, 1.0) + + # Calculate which point this corresponds to. + index = PDM_SORT(pdmp,j) + PDM_INUSE(pdmp,index) = 1 + call gflush (gp) + } + + return (index) + } else + return (0) +end |