|
 |
"Cipped" away about half of 9000 cells. Increasing chance on a "grain"
along y axis, so the lower ones will be bigger.
---%<------%<------%<---
#version 3.7;
global_settings{assumed_gamma 1.0}
#default{finish{ambient 0.1 diffuse 0.9}}
#include "math.inc"
#include "transforms.inc"
#include "rand.inc"
// run: scene
// run: voro++pov.py -c "%i %q %w %p %s %t %l" -o -v -1 40 -1 40 -1 40
vorobox.txt
// run: scene
camera{
perspective angle 50
location <0,30,-100>
look_at <0,-5,0>
}
light_source{< -500,300,-500> rgb <0.3,0.3,0.6>}
light_source{<1500,3000,-30> rgb 1}
#macro VoroPolygons(Normals, FaceIndex, Vertices)
/*==
Creates a single voronoi cell from data generated by Voro++ by
constructing
a polygon per face.
Normals (arr, vec) : an array of normals
FaceIndex (arr, int) : an array with indices that points to vertices
Vertices (arr, vec) : an array with vertices of the cell faces
*/
#local Length = dimension_size(FaceIndex,1);
#local Voro = union{
#for(Face,0,Length-1,1)
#local VertexIndices = FaceIndex[Face];
#local Normal = Normals[Face];
#local Distance = vdot(Normal,Vertices[VertexIndices[0]]);
#local Transform = transform{translate Normal*Distance};//*scale
#local VILength = dimension_size(VertexIndices,1);
polygon{
VILength+1,
#for(VertexIndex,0,VILength-1,1)
#local Index = VertexIndices[VertexIndex];
#local Vertex = Vertices[Index];
#local Vertex = VProject_Plane(Vertex, Normal);
#local Vertex = vtransform(Vertex, Transform);
#if (VertexIndex = 0)
#local P0 = Vertex;
#end
Vertex
#end
P0
}
#end
}
object{Voro}
#end
#macro Id2Coord(Id, Xn, Yn)
#local Zi = floor(Id/(Xn*Yn));
#local Yi = floor((Id-Zi*Xn*Yn)/Xn);
#local Xi = Id-(Zi*Xn*Yn)-(Yi*Xn);
(<Xi,Yi,Zi>)
#end
#macro Coord2Id(Coord, Xn, Yn)
(Coord.x+Xn*Coord.y+Xn*Yn*Coord.z)
#end
#declare OutFile = "vorobox.txt";
#declare IncFile = "vorobox.inc";
#declare Stream = seed(7);
#declare N = 40;
#declare N3 = pow(N,3);
#declare arrGrid = array[N3];
// first run
#if(!file_exists(OutFile))
#fopen Voro OutFile write
#for(I,0,N3,1)
#declare Coord = Id2Coord(I,N,N);
#declare Probability = pow(10, 3*(Coord.y+1)/N)/1000;
#if(Prob(Probability, Stream))
#declare Point = VRand(Stream) + Coord;
#write(Voro,I," ",Point.x," ",Point.y," ",Point.z,"\n")
#end
#end
#end
// second run
#if(file_exists(IncFile))
#include IncFile
#debug concat(IncFile,"\n")
#declare Length = dimension_size(Labels,1);
union{
#for(I,0,Length-1,1)
#declare Label = Labels[I];
#declare Coord = Id2Coord(Label, N, N);
#if(Coord.y<Coord.x+5)
#declare Point = Points[I];
#declare Normals = FaceNormal[I];
#declare FaceIndex = FaceVerticesIndex[I];
#declare Vertices = RelVertices[I];
object{VoroPolygons(Normals, FaceIndex, Vertices) scale 0.9
translate Point pigment{rgb 1}}
#end
#end
translate <-N/2,-N/2,-N/2>
rotate <0,-45,0>
}
#end
Post a reply to this message
Attachments:
Download 'vororocks.jpg' (151 KB)
Preview of image 'vororocks.jpg'
|
|
povray