Re: Re: Re: Conditinal Plots

*To*: mathgroup at smc.vnet.net*Subject*: [mg51148] Re: [mg51122] Re: [mg51037] Re: Conditinal Plots*From*: "David Park" <djmp at earthlink.net>*Date*: Wed, 6 Oct 2004 04:34:30 -0400 (EDT)*Sender*: owner-wri-mathgroup at wolfram.com

DrawGraphics from my web site has the ability to do this kind of thing. Here is an example. Needs["DrawGraphics`DrawingMaster`"] f[x_, y_] := Sin[3 Sqrt[x^2 + y^2]]/Sqrt[x^2 + y^2] grid = N[MakePolyGrid[{41, 41}, {{-Pi, -Pi}, {Pi, Pi}}]]; trimmedgrid = TrimPolygons[#1^2 + #2^2 & , {-1, Pi^2}][grid]; MakePolyGrid creates a 2D grid of rectangular polygons. TrimPolygons trimes the Polygons according to whether they exceed the limits (the second parameter) of a function (the first parameter). In the above case, the grid is trimmed to a circle. The following displays it. Draw2D[ {MediumAquamarine, trimmedgrid, trimmedgrid // PolygonOutline[Black]}, AspectRatio -> Automatic, Background -> Linen]; We can then plot the function in 3D by raising the trimmed grid to a 3D surface as with the following statement. Draw3DItems[ {SurfaceColor[MediumAquamarine], EdgeForm[ColorMix[MediumAquamarine, Black][0.4]], trimmedgrid // RaiseTo3D[f]}, NeutralLighting[0.3, 0.5, 0.1], PlotRange -> {-1, 3}, BoxRatios -> {1, 1, 1/2}, BoxStyle -> Gray, ViewPoint -> {1.737, -2.778, 0.845}, Background -> Linen, ImageSize -> 450]; NeutralLighting inserts a series of Lighting options that tone down the color saturation of the lights so surface colors will better display. I used EdgeForm to make the edges a darker shade of the surface color. David Park djmp at earthlink.net http://home.earthlink.net/~djmp/ From: Garry Helzer [mailto:gah at math.umd.edu] To: mathgroup at smc.vnet.net I missed the original posting by philippe so I don't know the precise question or philippe's email address. Clipping a graphic to one side of a plane is a basic computational geometry task but it is considered rather low level and Mathematica has no user level function that will do the clipping. With such a function one could plot the graphic over a sufficiently large domain and then clip by the vertical planes over the polygon edges. The lecture notes on Geometry for Computer Graphics at http://www.math.umd.edu/~gah/ contain an algorithm for such clipping which is implemented in the accompanying package Peano.m. However the implementation is for show rather than go?and slow. Further, the mathematics used is not widely known and some study of the lecture notes would be necessary before the function could be used with confidence (the trivector defining the plane is one of the arguments). The more general question of clipping a graphic at a surface is addressed in the article "Clipping Polygons" in The Mathematica Journal, volume 8, number 4. That article contains an illustration of a surface ( a monkey saddle) plotted over a rectangle, then clipped by the cylinder over the circle x^2+y^2=1. Unfortunately, when the issue is downloaded the package with the clipping function is not included. The package works under version 4 but needs modification to work under version 5 because of a syntax change to the function FindRoot. I intend to produce a version 5 form of the package but will not be able to do so for several months. On Oct 2, 2004, at 3:18 AM, philippe wrote: > If you really want to trim your plot with a cylinder you can try to > adapt the principle which is present in the "SurfaceColor" example in > the help browser. > To get the line drawn on the surface, first convert your plot to a > generic graphics (Graphics[....]), then type %//InputForm. > You will have the internal structure of the graphic. Then, with > replacement rules, you can keep only the part you are interested in. > You can try something like that. > > > > On 3 Sep 04 11:04:27 -0400 (EDT), matt wrote: >> Dear all, >> >> Í´m trying to produce a 3D plot with a 2D ploygon (within unit > circle) >> as its base. Equivalently, I want my 3D surface to sit inside (be >> bounded by) a cylinder with n-gon cross section. >> >> How can I trim off the edges? >> I´ve tried to use >> If[ (condition), Plot, Blah] but it doesn´t seem to work..... > > Garry Helzer gah at math.umd.edu