1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
|
// Import the data
electrondensity = Import("watermolecule");
// Partition the data
electrondensity = Partition(electrondensity);
// Compute the gradient
gradient = Gradient(electrondensity);
// Create a camera
camera = AutoCamera(electrondensity,"off-front");
// Make about 10 samples
starts = Sample(electrondensity, 10);
// Compute the divergence and curl
div,curl = DivCurl(gradient);
// Create streamlines using the samples as starting points
// Use the curl field to make the tubes twist
lines = Streamline(gradient, starts,curl=curl);
// Map the original field onto the lines
streamlinesmapped = Map(lines,electrondensity);
// Color the lines
colored = AutoColor(streamlinesmapped);
// Turn the colored lines into colored tubes, with a diameter of 0.4
// in world coordinates, and three sides
tubes = Tube(colored,diameter=0.4,ngon=3);
Display(tubes,camera);
|
