MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

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
==========================================================


  • Prev by Date: Problem in Solving Double Integral for PDF transformation
  • Next by Date: Re: is there a better way to do constraint logic programming in Mathematica?
  • Previous by thread: making a fast projection of a crytsallographic unit cell
  • Next by thread: Re: [Mathematica 6] Bug in Labeled, Center align option