Re: making a fast projection of a crytsallographic unit cell

*To*: mathgroup at smc.vnet.net*Subject*: [mg80891] Re: [mg80779] making a fast projection of a crytsallographic unit cell*From*: Curtis Osterhoudt <cfo at lanl.gov>*Date*: Wed, 5 Sep 2007 02:45:49 -0400 (EDT)*Organization*: LANL*References*: <200709010432.AAA26593@smc.twtelecom.net>*Reply-to*: cfo at lanl.gov

Huh. I can't get past the first Show[Graphics3D... statement, because Mathematica immediately segfaults. Anyone else have this problem? This is on a fast machine, with 2GB RAM, and In[1]:= $Version Out[1]= "6.0 for Linux x86 (32-bit) (June 19, 2007)" On Friday 31 August 2007 22:32:47 Roger Bagula wrote: > The problem is that I want to plot representations of crystallographic > cells fast in Mathematica. > > The face centerd cubic unit cell can be represented in 3d as a set of > points. > Hermann-Mauguin space groupo notation: Fm3m > Sodium Chloride and many metal crystals use this crystal cell. > (* points for cubic close pack/ face centered cubic*) > p = {{1, 1, 1}, {1, 0, 0}, {2, 1, 0}, {0, 0, 1}, {2, 2, 1}, {1, 2, > 2}, {0, 1, 2}, {0, 1, 0}, {1, 2, 0}, {0, 2, 1}, {2, 0, 1}, { > 2, 1, 2}, {1, 0, 2}}; > pieces = > Complement[ > Flatten[Table[{i, j, k}, {i, 0, 2}, {j, 0, 2}, {k, 0, 2}], > 2], p] > (* showing as points*) > Show[Graphics3D[Table[Point[pieces[[n]]], { > n, 1, Length[pieces]}]], Boxed -> False] > > (* creating the equivalent of a graphic primitive like a Cuboid[] for a > sphere*) > > Clear[sphere] > << RealTime3D` > (* RealTime3D allows you to drag the graphic around to see from > different view points*) > sphere[n_]:=ParametricPlot3D[{Sin[t]*Sin[p]/2-pieces[[n]][[ > 1]],Cos[t]*Sin[p]/2-pieces[[n]][[2]],Cos[p]/2-pieces[[n]][[3]]},{t,- > Pi,Pi},{p,0,Pi},Axes\[Rule]False,Boxed\[Rule]False] > > (* plotting them at the points them at the points*) > gr=Table[sphere[n],{n,1,Length[pieces]}]; > > Show[gr] > > The problem isn't doing it, but doing it faster. > There are 230 Federov space groups in 3d. > Roger Bagula -- ========================================================== Curtis Osterhoudt cfo at remove_this.lanl.and_this.gov PGP Key ID: 0x4DCA2A10 Please avoid sending me Word or PowerPoint attachments See http://www.gnu.org/philosophy/no-word-attachments.html ==========================================================

**References**:**making a fast projection of a crytsallographic unit cell***From:*Roger Bagula <rlbagula@sbcglobal.net>