[S] Summary: Delunay.triangle

Steve Friedman (friedman@gis.umn.edu)
Tue, 27 Oct 1998 12:00:41 -0600


This is a followup to the summary I posted a few days ago regarding
the delunay.triangle function. It appears that an unrecognized programming
error was at fault.

I appreciate everyone's assistance, Rolf Turner who wrote the original function
Bill Dunlap who found the coding error and provided the attached delunay.patch
Stephen Kaluzny, and everyone else who found this request for help curious enough
to respond.

Thanks,
Steve Friedman

Bill Dunlap wrote:

> On Mon, 26 Oct 1998, Steve Friedman wrote:
>
> > Bill,
> >
> > I should have saved the last email a bit to confirm that the changes
> > will help.
> >
> > I've edited the file delunay.triangle as follows:
> >
> > a) changed all occurrences of inside to rolfinside
> > b) tried the function with all references to dyn.load as is, as well as
> > by changing this to dyn.load.shared.
> >
> > The program continues to bail.
> > What am I doing wrong that you've managed to fix?
>
> Rolf may well have a fixed version now, but here is what I did to make
> it work with Splus3.4 on a IRIX Release 6.2 IP22 machine (uname -a says
> "IRIX hostname 6.2 03131015 IP22"). I did the same on a IRIX Release
> 5.3 IP22 machine and got identical results.
>
> (a) unshar delaunay.shar from statlib (dated 6 December 1996)
> and patch it with
> % patch < delaunay.patch
> I've attached delaunay.patch; if you don't have the patch program get
> it from some GNU server or apply the patches by hand, there are not too
> many of them. One patch to Makefile changes /usr/local/lib/Splus to
> Splus3.4 as the command to start up Splus. Use whatever your local
> command is there and for the rest of this note. I actually don't make
> use of his Makefile.
>
> (b) make the shared library deldir.so with
> % Splus3.4 SHLIB -o deldir.so *.r
> You might as well check that inside has really been renamed with
> % Splus3.4 NM deldir.so | grep inside
> deldir.so: 5ffe3700 T rolfinside_
>
> (c) load the Splus functions with
> % cat *.sf | Splus3.4
>
> (d) test it with
> % Splus3.4
> S-PLUS : Copyright (c) 1988, 1996 MathSoft, Inc.
> S : Copyright AT&T.
> Version 3.4 Release 1 for Silicon Graphics Iris, IRIX 5.3 : 1996
> Working data will be in .Data
> > sc <- deldir(state.center)
> > printer(height=30) # or motif() or postscript() ...
> > usa(states=F)
> > plot(sc, add=T)
> Line segments out of bounds X= -128.2442 Y= 40.4146
> Line segments out of bounds X= -79.9160 Y= 51.38756
> Line segments out of bounds X= -97.463 Y= 46.1092
> Line segments out of bounds X= -104.9118 Y= 51.38756
> Line segments out of bounds X= -70.7289 Y= 37.0564
> Line segments out of bounds X= -63.1531 Y= 39.9366
> Line segments out of bounds X= -70.8911 Y= 45.27073
> Line segments out of bounds X= -77.1438 Y= 51.38756
> Line segments out of bounds X= -122.7552 Y= 51.38756
> Line segments out of bounds X= -114.4167 Y= 49.1438
> Line segments out of bounds X= -87.6798 Y= 42.62956
> Line segments out of bounds X= -89.9974 Y= 51.38756
> > show() # this looks better on motif()
>
>
>
>
> . . . . . .
> ..................................... .... . . .
> .... ........................................ . . . ...
> . .... .... .......... .. ..................................
> . .. ........ .. .. ......... .............................
> .. . .................. ....... ........................
> . . ....... .... .. ...... ..... ............. ..........
> ... . ..... ....... .. . ..... .............. ........... ..
> .. ......... ............ ... .. ..... .................... .
> . ... ...... ...................... .... .............
> .. .. .... .......... ..... ............. ........ .
> . . .... ...... ............... . .......... ..... .
> .. ..... ................ . ................. .... .
> ... ......... ............ .................. .. .
> ................................ ...... ......... . . .
> . .................. ................... . . .
> . . . ......................... .... . .
> . . . .. ... .......... ... ..
> . . . ... . . .. . ....
> . . . . ... .
> ..
>
> If it crashes while running deldir() (or other function) perhaps there is some other
> naming conflict. You can check for them before loading ./deldir.so with
>
> temp <- tempfile()
> unix(paste("NM deldir.so >", temp))
> nm <- read.table(temp, header = F, row.names = NULL, col.names = c("file",
> "addr", "type", "symbol"), as.is = T)
> nm$symbol[nm$type != "U" & is.loaded(nm$symbol)]
>
> If the last command prints out any names aside from some boilerplate ones
> starting with "." or "_" then they will need to be renamed as well.
>
> ----------------------------------------------------------------------------
> Bill Dunlap 22461 Mt Vernon-Big Lake Rd
> Data Analysis Products Div. of MathSoft, Inc. Mount Vernon, WA 98274
> bill@statsci.com 360-428-8146
>
> "All statements in this message represent the opinions of the author and do
> not necessarily reflect MathSoft policy or position."

-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu. To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message: unsubscribe s-news