Re: [Q] ExtendGraphics ListSurfacePlot3D inconsistent error

*To*: mathgroup at smc.vnet.net*Subject*: [mg3313] Re: [Q] ExtendGraphics ListSurfacePlot3D inconsistent error*From*: ianc (Ian Collier)*Date*: Mon, 26 Feb 1996 02:54:45 -0500*Organization*: Wolfram Research, Inc.*Sender*: owner-wri-mathgroup at wolfram.com

In article <4gi5df$6n1 at dragonfly.wolfram.com>, adbened at alumnae.caltech.edu ( Alto D. Benedicto) wrote: > When two datasets sharing a common boundary were "Union"ed and > Flattened such that it now has the form > { {x1,y1,z1}, {x2,y2,z2}, {x3,y3,z3}, ... , {x99,y99,z99} } > ListSurfacePlot3D[theabovelist] gave the following error message: > LinkObject::linkd > LinkObject[delaunay, 2, 2] > is closed; the connection is dead. > > Is this error related to the illustrative test example below? > > In[1] := AppendTo[$Path, "/usr/local/TWJ_Packages"] > In[2] := <<ExtendGraphics`SurfaceGraphics3D` > In[3] := data1 = Table[{x/100, y/1000., 0.25 x}, {x,-1,2},{y,2,4}] > In[4] := data2 = Table[{x/100, y/1000., 0.25 x}, {x,2,2},{y,2,4}] > In[5] := dataunion = Union[data1, data2] > In[6] := data1flat = Flatten[data1,1] > In[7] := data2flat = Flatten[data2,1] > In[8] := dataunionflat = Flatten[dataunion,1] > > Now, ListSurfacePlot3D[data1] gives correct plot while > ListSurfacePlot3D[data1flat] gives wrong plot (opposite to our > intuition, since syntax for that of data1 is supposedly incorrect). > Curiously, ListSurfacePlot3D[data2flat] resulted in same error > message as my real data, though ListSurfacePlot3D[data2] has no > error message. > ListSurfacePlot3D[dataunion] is correct while > ListSurfacePlot3D[dataunionflat] is wrong. > > Even stranger,when the limit in data2 is changed to {x,2,4}, the > resultant ListSurfacePlot3D[dataunionflat] is now correct, inspite > of the fact that ListSurfacePlot3D[data1flat] is still wrong. > > The original data2 might result in an error because of {x,2,2}. > But what about the other errors? They don't even seem to be consistent. > And to reiterate, why did my real data give that error message when the > component datasets (also Flatten) making it give correct plots? > > THANK YOU very much in advance for any help suggested. > Here is a reply from Tom Wickham-Jones: There are some problems that I have found with the triangulation routine when it is applied to regular data. One way to avoid these problems is to slightly perturb the data. This can be done with the following function. PerturbData[ d_List] := Map[ PerturbData, d] PerturbData[ d_] := If[ NumberQ[ d], N[ d + d (Random[ ]-0.5)/1000000.], d] Then in this example the following will work. ListSurfacePlot3D[ PerturbData[ dataunionflat]] Sometime in the future I will updating these functions with this and a number of other improvements. Tom Wickham-Jones WRI ----------------------------------------------------------- Ian Collier Wolfram Research, Inc. ----------------------------------------------------------- tel:(217) 398-0700 fax:(217) 398-0747 ianc at wolfram.com Wolfram Research Home Page: http://www.wolfram.com/ ----------------------------------------------------------- ==== [MESSAGE SEPARATOR] ====