svn merge ^/trunk/blender -r43976:43995
[blender-staging.git] / source / blender / bmesh / operators / primitiveops.c
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  *
18  * Contributor(s): Joseph Eagar.
19  *
20  * ***** END GPL LICENSE BLOCK *****
21  */
22
23 #include "MEM_guardedalloc.h"
24
25 #include "DNA_meshdata_types.h"
26 #include "DNA_mesh_types.h"
27 #include "DNA_object_types.h"
28 #include "DNA_scene_types.h"
29
30 #include "BLI_utildefines.h"
31
32 #include "BLI_math.h"
33 #include "BLI_ghash.h"
34 #include "BLI_array.h"
35
36 #include "ED_mesh.h"
37
38 #include "bmesh.h"
39 #include "mesh_intern.h"
40 #include "bmesh_private.h"
41
42 #include <stdio.h>
43 #include <stdlib.h>
44 #include <string.h>
45
46
47 /* ************************ primitives ******************* */
48
49 static float icovert[12][3] = {
50         {0.0f,0.0f,-200.0f}, 
51         {144.72f, -105.144f,-89.443f},
52         {-55.277f, -170.128,-89.443f}, 
53         {-178.885f,0.0f,-89.443f},
54         {-55.277f,170.128f,-89.443f}, 
55         {144.72f,105.144f,-89.443f},
56         {55.277f,-170.128f,89.443f},
57         {-144.72f,-105.144f,89.443f},
58         {-144.72f,105.144f,89.443f},
59         {55.277f,170.128f,89.443f},
60         {178.885f,0.0f,89.443f},
61         {0.0f,0.0f,200.0f}
62 };
63
64 static short icoface[20][3] = {
65         {0,1,2},
66         {1,0,5},
67         {0,2,3},
68         {0,3,4},
69         {0,4,5},
70         {1,5,10},
71         {2,1,6},
72         {3,2,7},
73         {4,3,8},
74         {5,4,9},
75         {1,10,6},
76         {2,6,7},
77         {3,7,8},
78         {4,8,9},
79         {5,9,10},
80         {6,10,11},
81         {7,6,11},
82         {8,7,11},
83         {9,8,11},
84         {10,9,11}
85 };
86
87 // HACK: these can also be found in cmoview.tga.c, but are here so that they can be found by linker
88 // this hack is only used so that scons & mingw + split-sources hack works
89         // ------------------------------- start copied code
90 /* these are not the monkeys you are looking for */
91 static int monkeyo = 4;
92 static int monkeynv = 271;
93 static int monkeynf = 250;
94 static signed char monkeyv[271][3] = {
95 {-71,21,98},{-63,12,88},{-57,7,74},{-82,-3,79},{-82,4,92},
96 {-82,17,100},{-92,21,102},{-101,12,95},{-107,7,83},
97 {-117,31,84},{-109,31,95},{-96,31,102},{-92,42,102},
98 {-101,50,95},{-107,56,83},{-82,66,79},{-82,58,92},
99 {-82,46,100},{-71,42,98},{-63,50,88},{-57,56,74},
100 {-47,31,72},{-55,31,86},{-67,31,97},{-66,31,99},
101 {-70,43,100},{-82,48,103},{-93,43,105},{-98,31,105},
102 {-93,20,105},{-82,31,106},{-82,15,103},{-70,20,100},
103 {-127,55,95},{-127,45,105},{-127,-87,94},{-127,-41,100},
104 {-127,-24,102},{-127,-99,92},{-127,52,77},{-127,73,73},
105 {-127,115,-70},{-127,72,-109},{-127,9,-106},{-127,-49,-45},
106 {-101,-24,72},{-87,-56,73},{-82,-89,73},{-80,-114,68},
107 {-85,-121,67},{-104,-124,71},{-127,-126,74},{-71,-18,68},
108 {-46,-5,69},{-21,19,57},{-17,55,76},{-36,62,80},
109 {-64,77,88},{-86,97,94},{-107,92,97},{-119,63,96},
110 {-106,53,99},{-111,39,98},{-101,12,95},{-79,2,90},
111 {-64,8,86},{-47,24,83},{-45,38,83},{-50,48,85},
112 {-72,56,92},{-95,60,97},{-127,-98,94},{-113,-92,94},
113 {-112,-107,91},{-119,-113,89},{-127,-114,88},{-127,-25,96},
114 {-127,-18,95},{-114,-19,95},{-111,-29,96},{-116,-37,95},
115 {-76,-6,86},{-48,7,80},{-34,26,77},{-32,48,84},
116 {-39,53,93},{-71,70,102},{-87,82,107},{-101,79,109},
117 {-114,55,108},{-111,-13,104},{-100,-57,91},{-95,-90,88},
118 {-93,-105,85},{-97,-117,81},{-106,-119,81},{-127,-121,82},
119 {-127,6,93},{-127,27,98},{-85,61,95},{-106,18,96},
120 {-110,27,97},{-112,-88,94},{-117,-57,96},{-127,-57,96},
121 {-127,-42,95},{-115,-35,100},{-110,-29,102},{-113,-17,100},
122 {-122,-16,100},{-127,-26,106},{-121,-19,104},{-115,-20,104},
123 {-113,-29,106},{-117,-32,103},{-127,-37,103},{-94,-40,71},
124 {-106,-31,91},{-104,-40,91},{-97,-32,71},{-127,-112,88},
125 {-121,-111,88},{-115,-105,91},{-115,-95,93},{-127,-100,84},
126 {-115,-96,85},{-115,-104,82},{-121,-109,81},{-127,-110,81},
127 {-105,28,100},{-103,20,99},{-84,55,97},{-92,54,99},
128 {-73,51,99},{-55,45,89},{-52,37,88},{-53,25,87},
129 {-66,13,92},{-79,8,95},{-98,14,100},{-104,38,100},
130 {-100,48,100},{-97,46,97},{-102,38,97},{-96,16,97},
131 {-79,11,93},{-68,15,90},{-57,27,86},{-56,36,86},
132 {-59,43,87},{-74,50,96},{-91,51,98},{-84,52,96},
133 {-101,22,96},{-102,29,96},{-113,59,78},{-102,85,79},
134 {-84,88,76},{-65,71,71},{-40,58,63},{-25,52,59},
135 {-28,21,48},{-50,0,53},{-71,-12,60},{-127,115,37},
136 {-127,126,-10},{-127,-25,-86},{-127,-59,24},{-127,-125,59},
137 {-127,-103,44},{-127,-73,41},{-127,-62,36},{-18,30,7},
138 {-17,41,-6},{-28,34,-56},{-68,56,-90},{-33,-6,9},
139 {-51,-16,-21},{-45,-1,-55},{-84,7,-85},{-97,-45,52},
140 {-104,-53,33},{-90,-91,49},{-95,-64,50},{-85,-117,51},
141 {-109,-97,47},{-111,-69,46},{-106,-121,56},{-99,-36,55},
142 {-100,-29,60},{-101,-22,64},{-100,-50,21},{-89,-40,-34},
143 {-83,-19,-69},{-69,111,-49},{-69,119,-9},{-69,109,30},
144 {-68,67,55},{-34,52,43},{-46,58,36},{-45,90,7},
145 {-25,72,16},{-25,79,-15},{-45,96,-25},{-45,87,-57},
146 {-25,69,-46},{-48,42,-75},{-65,3,-70},{-22,42,-26},
147 {-75,-22,19},{-72,-25,-27},{-13,52,-30},{-28,-18,-16},
148 {6,-13,-42},{37,7,-55},{46,41,-54},{31,65,-54},
149 {4,61,-40},{3,53,-37},{25,56,-50},{35,37,-52},
150 {28,10,-52},{5,-5,-39},{-21,-9,-17},{-9,46,-28},
151 {-6,39,-37},{-14,-3,-27},{6,0,-47},{25,12,-57},
152 {31,32,-57},{23,46,-56},{4,44,-46},{-19,37,-27},
153 {-20,22,-35},{-30,12,-35},{-22,11,-35},{-19,2,-35},
154 {-23,-2,-35},{-34,0,-9},{-35,-3,-22},{-35,5,-24},
155 {-25,26,-27},{-13,31,-34},{-13,30,-41},{-23,-2,-41},
156 {-18,2,-41},{-21,10,-41},{-29,12,-41},{-19,22,-41},
157 {6,42,-53},{25,44,-62},{34,31,-63},{28,11,-62},
158 {7,0,-54},{-14,-2,-34},{-5,37,-44},{-13,14,-42},
159 {-7,8,-43},{1,16,-47},{-4,22,-45},{3,30,-48},
160 {8,24,-49},{15,27,-50},{12,35,-50},{4,56,-62},
161 {33,60,-70},{48,38,-64},{41,7,-68},{6,-11,-63},
162 {-26,-16,-42},{-17,49,-49},
163 };
164
165 static signed char monkeyf[250][4] = {
166 {27,4,5,26}, {25,4,5,24}, {3,6,5,4}, {1,6,5,2}, {5,6,7,4}, 
167 {3,6,7,2}, {5,8,7,6}, {3,8,7,4}, {7,8,9,6}, 
168 {5,8,9,4}, {7,10,9,8}, {5,10,9,6}, {9,10,11,8}, 
169 {7,10,11,6}, {9,12,11,10}, {7,12,11,8}, {11,6,13,12}, 
170 {5,4,13,12}, {3,-2,13,12}, {-3,-4,13,12}, {-5,-10,13,12}, 
171 {-11,-12,14,12}, {-13,-18,14,13}, {-19,4,5,13}, {10,12,4,4}, 
172 {10,11,9,9}, {8,7,9,9}, {7,5,6,6}, {6,3,4,4}, 
173 {5,1,2,2}, {4,-1,0,0}, {3,-3,-2,-2}, {22,67,68,23}, 
174 {20,65,66,21}, {18,63,64,19}, {16,61,62,17}, {14,59,60,15}, 
175 {12,19,48,57}, {18,19,48,47}, {18,19,48,47}, {18,19,48,47}, 
176 {18,19,48,47}, {18,19,48,47}, {18,19,48,47}, {18,19,48,47}, 
177 {18,19,48,47}, {18,-9,-8,47}, {18,27,45,46}, {26,55,43,44}, 
178 {24,41,42,54}, {22,39,40,23}, {20,37,38,21}, {18,35,36,19}, 
179 {16,33,34,17}, {14,31,32,15}, {12,39,30,13}, {11,48,45,38}, 
180 {8,36,-19,9}, {8,-20,44,47}, {42,45,46,43}, {18,19,40,39}, 
181 {16,17,38,37}, {14,15,36,35}, {32,44,43,33}, {12,33,32,42}, 
182 {19,44,43,42}, {40,41,42,-27}, {8,9,39,-28}, {15,43,42,16}, 
183 {13,43,42,14}, {11,43,42,12}, {9,-30,42,10}, {37,12,38,-32}, 
184 {-33,37,45,46}, {-33,40,41,39}, {38,40,41,37}, {36,40,41,35}, 
185 {34,40,41,33}, {36,39,38,37}, {35,40,39,38}, {1,2,14,21}, 
186 {1,2,40,13}, {1,2,40,39}, {1,24,12,39}, {-34,36,38,11}, 
187 {35,38,36,37}, {-37,8,35,37}, {-11,-12,-45,40}, {-11,-12,39,38}, 
188 {-11,-12,37,36}, {-11,-12,35,34}, {33,34,40,41}, {33,34,38,39}, 
189 {33,34,36,37}, {33,-52,34,35}, {33,37,36,34}, {33,35,34,34}, 
190 {8,7,37,36}, {-32,7,35,46}, {-34,-33,45,46}, {4,-33,43,34}, 
191 {-34,-33,41,42}, {-34,-33,39,40}, {-34,-33,37,38}, {-34,-33,35,36}, 
192 {-34,-33,33,34}, {-34,-33,31,32}, {-34,-4,28,30}, {-5,-34,28,27}, 
193 {-35,-44,36,27}, {26,35,36,45}, {24,25,44,45}, {25,23,44,42}, 
194 {25,24,41,40}, {25,24,39,38}, {25,24,37,36}, {25,24,35,34}, 
195 {25,24,33,32}, {25,24,31,30}, {15,24,29,38}, {25,24,27,26}, 
196 {23,12,37,26}, {11,12,35,36}, {-86,-59,36,-80}, {-60,-61,36,35}, 
197 {-62,-63,36,35}, {-64,-65,36,35}, {-66,-67,36,35}, {-68,-69,36,35}, 
198 {-70,-71,36,35}, {-72,-73,36,35}, {-74,-75,36,35}, {42,43,53,58}, 
199 {40,41,57,56}, {38,39,55,57}, {-81,-80,37,56}, {-83,-82,55,52}, 
200 {-85,-84,51,49}, {-87,-86,48,49}, {47,50,51,48}, {46,48,51,49}, 
201 {43,46,49,44}, {-92,-91,45,42}, {-23,49,50,-20}, {-94,40,48,-24}, 
202 {-96,-22,48,49}, {-97,48,21,-90}, {-100,36,50,23}, {22,49,48,-100}, 
203 {-101,47,46,22}, {21,45,35,25}, {33,34,44,41}, {13,14,28,24}, 
204 {-107,26,30,-106}, {14,46,45,15}, {14,44,43,-110}, {-111,42,23,-110}, 
205 {6,7,45,46}, {45,44,47,46}, {45,46,47,48}, {47,46,49,48}, 
206 {17,49,47,48}, {17,36,46,48}, {35,36,44,45}, {35,36,40,43}, 
207 {35,36,38,39}, {-4,-3,37,35}, {-123,34,33,1}, {-9,-8,-7,-6}, 
208 {-10,-7,32,-125}, {-127,-11,-126,-126}, {-7,-6,5,31}, {4,5,33,30}, 
209 {4,39,33,32}, {4,35,32,38}, {20,21,39,38}, {4,37,38,5}, 
210 {-11,-10,36,3}, {-11,15,14,35}, {13,16,34,34}, {-13,14,13,13}, 
211 {-3,1,30,29}, {-3,28,29,1}, {-2,31,28,-1}, {12,13,27,30}, 
212 {-2,26,12,12}, {35,29,42,36}, {34,35,36,33}, {32,35,36,31}, 
213 {30,35,36,29}, {28,35,36,27}, {26,35,36,25}, {34,39,38,35}, 
214 {32,39,38,33}, {30,39,38,31}, {28,39,38,29}, {26,39,38,27}, 
215 {25,31,32,38}, {-18,-17,45,44}, {-18,17,28,44}, {-24,-20,42,-23}, 
216 {11,35,27,14}, {25,28,39,41}, {37,41,40,38}, {34,40,36,35}, 
217 {32,40,39,33}, {30,39,31,40}, {21,29,39,22}, {-31,37,28,4}, 
218 {-32,33,35,36}, {32,33,34,34}, {18,35,36,48}, {34,25,40,35}, 
219 {24,25,38,39}, {24,25,36,37}, {24,25,34,35}, {24,25,32,33}, 
220 {24,13,41,31}, {17,11,41,35}, {15,16,34,35}, {13,14,34,35}, 
221 {11,12,34,35}, {9,10,34,35}, {7,8,34,35}, {26,25,37,36}, 
222 {35,36,37,38}, {37,36,39,38}, {37,38,39,40}, {25,31,36,39}, 
223 {18,34,35,30}, {17,22,30,33}, {19,29,21,20}, {16,26,29,17}, 
224 {24,29,28,25}, {22,31,28,23}, {20,31,30,21}, {18,31,30,19}, 
225 {16,30,17,17}, {-21,-22,35,34}, {-21,-22,33,32}, {-21,-22,31,30}, 
226 {-21,-22,29,28}, {-21,-22,27,26}, {-28,-22,25,31}, {24,28,29,30}, 
227 {23,24,26,27}, {23,24,25,25}, {-69,-35,-32,27}, {-70,26,25,-66}, 
228 {-68,-67,24,-33}, 
229 };
230
231 #define VERT_MARK       1
232
233 #define EDGE_ORIG       1
234 #define EDGE_MARK       2
235
236 #define FACE_MARK       1
237 #define FACE_NEW        2
238
239 void bmesh_create_grid_exec(BMesh *bm, BMOperator *op)
240 {
241         BMOperator bmop, prevop;
242         BMVert *eve, *preveve;
243         BMEdge *e;
244         float vec[3], mat[4][4], phi, phid, dia = BMO_Get_Float(op, "size");
245         int a, tot = BMO_Get_Int(op, "xsegments"), seg = BMO_Get_Int(op, "ysegments");
246
247         if (tot < 2) tot = 2;
248         if (seg < 2) seg = 2;
249
250         BMO_Get_Mat4(op, "mat", mat);
251
252         /* one segment first: the X axis */
253         phi = 1.0f;
254         phid = 2.0f / ((float)tot - 1);
255         for (a = 0; a < tot; a++) {
256                 vec[0] = dia * phi;
257                 vec[1] = -dia;
258                 vec[2] = 0.0f;
259                 mul_m4_v3(mat, vec);
260
261                 eve = BM_Make_Vert(bm, vec, NULL);
262                 BM_Select(bm, eve, TRUE);
263
264                 if (a) {
265                         e = BM_Make_Edge(bm, preveve, eve, NULL, TRUE);
266                         BMO_SetFlag(bm, e, EDGE_ORIG);
267                 }
268
269                 preveve = eve;
270                 phi -= phid;
271         }
272
273         /* extrude and translate */
274         vec[0] = vec[2] = 0.0f;
275         vec[1] = dia * phid;
276         mul_mat3_m4_v3(mat, vec);
277
278         for (a = 0; a < seg - 1; a++) {
279                 if (a) {
280                         BMO_InitOpf(bm, &bmop, "extrude_edge_only edges=%s", &prevop, "geomout");
281                         BMO_Exec_Op(bm, &bmop);
282                         BMO_Finish_Op(bm, &prevop);
283
284                         BMO_Flag_Buffer(bm, &bmop, "geomout", VERT_MARK, BM_VERT);
285                 }
286                 else {
287                         BMO_InitOpf(bm, &bmop, "extrude_edge_only edges=%fe", EDGE_ORIG);
288                         BMO_Exec_Op(bm, &bmop);
289                         BMO_Flag_Buffer(bm, &bmop, "geomout", VERT_MARK, BM_VERT);
290                 }
291
292                 BMO_CallOpf(bm, "translate vec=%v verts=%s", vec, &bmop, "geomout");
293                 prevop = bmop;
294         }
295
296         if (a)
297                 BMO_Finish_Op(bm, &bmop);
298
299         BMO_Flag_To_Slot(bm, op, "vertout", VERT_MARK, BM_VERT);
300 }
301
302 void bmesh_create_uvsphere_exec(BMesh *bm, BMOperator *op)
303 {
304         BMOperator bmop, prevop;
305         BMVert *eve, *preveve;
306         BMEdge *e;
307         BMIter iter;
308         float vec[3], mat[4][4], cmat[3][3], phi, q[4];
309         float phid, dia = BMO_Get_Float(op, "diameter");
310         int a, seg = BMO_Get_Int(op, "segments"), tot = BMO_Get_Int(op, "revolutions");
311
312         BMO_Get_Mat4(op, "mat", mat);
313
314         phid = 2.0f * (float)M_PI / tot;
315         phi = 0.25f * (float)M_PI;
316
317         /* one segment first */
318         phi = 0;
319         phid /= 2;
320         for (a = 0; a <= tot; a++) {
321                 /* Going in this direction, then edge extruding, makes normals face outward */
322                 vec[0] = -dia * sinf(phi);
323                 vec[1] = 0.0;
324                 vec[2] = dia * cosf(phi);
325                 eve = BM_Make_Vert(bm, vec, NULL);
326                 BMO_SetFlag(bm, eve, VERT_MARK);
327
328                 if (a != 0) {
329                         e = BM_Make_Edge(bm, preveve, eve, NULL, FALSE);
330                         BMO_SetFlag(bm, e, EDGE_ORIG);
331                 }
332
333                 phi+= phid;
334                 preveve = eve;
335         }
336
337         /* extrude and rotate; negative phi to make normals face outward */
338         phi = -M_PI / seg;
339         q[0] = cosf(phi);
340         q[3] = sinf(phi);
341         q[1] = q[2] = 0.0f;
342         quat_to_mat3(cmat, q);
343
344         for (a = 0; a < seg; a++) {
345                 if (a) {
346                         BMO_InitOpf(bm, &bmop, "extrude_edge_only edges=%s", &prevop, "geomout");
347                         BMO_Exec_Op(bm, &bmop);
348                         BMO_Finish_Op(bm, &prevop);
349                 }
350                 else {
351                         BMO_InitOpf(bm, &bmop, "extrude_edge_only edges=%fe", EDGE_ORIG);
352                         BMO_Exec_Op(bm, &bmop);
353                 }
354
355                 BMO_Flag_Buffer(bm, &bmop, "geomout", VERT_MARK, BM_VERT);
356                 BMO_CallOpf(bm, "rotate cent=%v mat=%m3 verts=%s", vec, cmat, &bmop, "geomout");
357                 
358                 prevop = bmop;
359         }
360
361         if (a)
362                 BMO_Finish_Op(bm, &bmop);
363
364         {
365                 float len, len2, vec2[3];
366
367                 len= 2*dia*sinf(phid / 2.0f);
368
369                 /* length of one segment in shortest parallen */
370                 vec[0]= dia*sinf(phid);
371                 vec[1]= 0.0;
372                 vec[2]= dia*cosf(phid);
373
374                 mul_v3_m3v3(vec2, cmat, vec);
375                 len2= len_v3v3(vec, vec2);
376
377                 /* use shortest segment length divided by 3 as merge threshold */
378                 BMO_CallOpf(bm, "removedoubles verts=%fv dist=%f", VERT_MARK, MIN2(len, len2) / 3.0f);
379         }
380
381         /* and now do imat */
382         BM_ITER(eve, &iter, bm, BM_VERTS_OF_MESH, NULL) {
383                 if (BMO_TestFlag(bm, eve, VERT_MARK)) {
384                         mul_m4_v3(mat, eve->co);
385                 }
386         }
387
388         BMO_Flag_To_Slot(bm, op, "vertout", VERT_MARK, BM_VERT);
389 }
390
391 void bmesh_create_icosphere_exec(BMesh *bm, BMOperator *op)
392 {
393         BMVert *eva[12];
394         BMIter liter;
395         BMLoop *l;
396         float vec[3], mat[4][4] /* , phi, phid */;
397         float dia = BMO_Get_Float(op, "diameter");
398         int a, subdiv = BMO_Get_Int(op, "subdivisions");
399
400         BMO_Get_Mat4(op, "mat", mat);
401
402         /* phid = 2.0f * (float)M_PI / subdiv; */ /* UNUSED */
403         /* phi = 0.25f * (float)M_PI; */         /* UNUSED */
404
405         dia /= 200.0f;
406         for (a = 0; a < 12; a++) {
407                 vec[0] = dia * icovert[a][0];
408                 vec[1] = dia * icovert[a][1];
409                 vec[2] = dia * icovert[a][2];
410                 eva[a] = BM_Make_Vert(bm, vec, NULL);
411
412                 mul_m4_v3(mat, eva[a]->co);
413                 BM_Select(bm, eva[a], TRUE);
414         }
415
416         for (a = 0; a < 20; a++) {
417                 BMFace *eftemp;
418                 BMVert *v1, *v2, *v3;
419
420                 v1 = eva[icoface[a][0]];
421                 v2 = eva[icoface[a][1]];
422                 v3 = eva[icoface[a][2]];
423
424                 eftemp = BM_Make_Face_QuadTri(bm, v1, v2, v3, NULL, NULL, FALSE);
425                 
426                 BM_ITER(l, &liter, bm, BM_LOOPS_OF_FACE, eftemp) {
427                         BMO_SetFlag(bm, l->e, EDGE_MARK);
428                 }
429
430                 BMO_SetFlag(bm, eftemp, FACE_MARK);
431         }
432
433         dia *= 200.0f;
434
435         for (a = 1; a < subdiv; a++) {
436                 BMOperator bmop;
437
438                 BMO_InitOpf(bm, &bmop,
439                             "esubd edges=%fe smooth=%f numcuts=%i gridfill=%i beauty=%i",
440                             EDGE_MARK, dia, 1, 1, B_SPHERE);
441                 BMO_Exec_Op(bm, &bmop);
442                 BMO_Flag_Buffer(bm, &bmop, "geomout", VERT_MARK, BM_VERT);
443                 BMO_Flag_Buffer(bm, &bmop, "geomout", EDGE_MARK, BM_EDGE);
444                 BMO_Finish_Op(bm, &bmop);
445         }
446
447         BMO_Flag_To_Slot(bm, op, "vertout", VERT_MARK, BM_VERT);
448 }
449
450 void bmesh_create_monkey_exec(BMesh *bm, BMOperator *op)
451 {
452         BMVert *eve;
453         BMVert **tv = MEM_mallocN(sizeof(*tv)*monkeynv * 2, "tv");
454         float mat[4][4];
455         int i;
456
457         BMO_Get_Mat4(op, "mat", mat);
458
459         for (i = 0; i < monkeynv; i++) {
460                 float v[3];
461
462                 v[0] = (monkeyv[i][0] + 127) / 128.0, v[1] = monkeyv[i][1] / 128.0, v[2] = monkeyv[i][2] / 128.0;
463
464                 tv[i] = BM_Make_Vert(bm, v, NULL);
465                 BMO_SetFlag(bm, tv[i], VERT_MARK);
466
467                 tv[monkeynv + i] = (fabsf(v[0] = -v[0]) < 0.001f) ? tv[i]: (eve = BM_Make_Vert(bm, v, NULL), mul_m4_v3(mat, eve->co), eve);
468                 BMO_SetFlag(bm, tv[monkeynv + i], VERT_MARK);
469
470                 mul_m4_v3(mat, tv[i]->co);
471         }
472
473         for (i = 0; i < monkeynf; i++) {
474                 BM_Make_Face_QuadTri(bm,
475                                      tv[monkeyf[i][0] + i - monkeyo],
476                                      tv[monkeyf[i][1] + i - monkeyo],
477                                      tv[monkeyf[i][2] + i - monkeyo],
478                                      (monkeyf[i][3] != monkeyf[i][2]) ? tv[monkeyf[i][3] + i - monkeyo] : NULL,
479                                      NULL, FALSE);
480
481                 BM_Make_Face_QuadTri(bm,
482                                      tv[monkeynv + monkeyf[i][2] + i - monkeyo],
483                                      tv[monkeynv + monkeyf[i][1] + i - monkeyo],
484                                      tv[monkeynv + monkeyf[i][0] + i - monkeyo],
485                                      (monkeyf[i][3] != monkeyf[i][2]) ? tv[monkeynv + monkeyf[i][3] + i - monkeyo]: NULL,
486                                      NULL, FALSE);
487         }
488
489         MEM_freeN(tv);
490
491         BMO_Flag_To_Slot(bm, op, "vertout", VERT_MARK, BM_VERT);
492 }
493
494
495 void bmesh_create_circle_exec(BMesh *bm, BMOperator *op)
496 {
497         BMVert *v1, *lastv1 = NULL, *cent1, *firstv1 = NULL;
498         float vec[3], mat[4][4], phi, phid;
499         float dia = BMO_Get_Float(op, "diameter");
500         int cap_ends = BMO_Get_Int(op, "cap_ends"), segs = BMO_Get_Int(op, "segments");
501         int cap_tris = BMO_Get_Int(op, "cap_tris");
502         int a;
503         
504         if (!segs)
505                 return;
506         
507         BMO_Get_Mat4(op, "mat", mat);
508
509         phid = 2.0f * (float)M_PI / segs;
510         phi = .25f * (float)M_PI;
511
512         if (cap_ends) {
513                 vec[0] = vec[1] = 0.0f;
514                 vec[2] = 0.0;
515                 mul_m4_v3(mat, vec);
516                 
517                 cent1 = BM_Make_Vert(bm, vec, NULL);
518         }
519
520         for (a = 0; a < segs; a++, phi += phid) {
521                 /* Going this way ends up with normal(s) upward */
522                 vec[0] = -dia * sinf(phi);
523                 vec[1] = dia * cosf(phi);
524                 vec[2] = 0.0f;
525                 mul_m4_v3(mat, vec);
526                 v1 = BM_Make_Vert(bm, vec, NULL);
527
528                 BMO_SetFlag(bm, v1, VERT_MARK);
529                 
530                 if (lastv1)
531                         BM_Make_Edge(bm, v1, lastv1, NULL, FALSE);
532                 
533                 if (a && cap_ends) {
534                         BMFace *f;
535                         
536                         f = BM_Make_Face_QuadTri(bm, cent1, lastv1, v1, NULL, NULL, FALSE);
537                         BMO_SetFlag(bm, f, FACE_NEW);
538                 }
539                 
540                 if (!firstv1)
541                         firstv1 = v1;
542
543                 lastv1 = v1;
544         }
545
546         if (!a)
547                 return;
548
549         BM_Make_Edge(bm, lastv1, firstv1, NULL, FALSE);
550
551         if (cap_ends) {
552                 BMFace *f;
553                 
554                 f = BM_Make_Face_QuadTri(bm, cent1, v1, firstv1, NULL, NULL, FALSE);
555                 BMO_SetFlag(bm, f, FACE_NEW);
556         }
557         
558         if (!cap_tris) {
559                 BMO_CallOpf(bm, "dissolvefaces faces=%ff", FACE_NEW);
560         }
561         
562         BMO_Flag_To_Slot(bm, op, "vertout", VERT_MARK, BM_VERT);
563 }
564
565 void bmesh_create_cone_exec(BMesh *bm, BMOperator *op)
566 {
567         BMVert *v1, *v2, *lastv1 = NULL, *lastv2 = NULL, *cent1, *cent2, *firstv1, *firstv2;
568         float vec[3], mat[4][4], phi, phid;
569         float dia1 = BMO_Get_Float(op, "diameter1");
570         float dia2 = BMO_Get_Float(op, "diameter2");
571         float depth = BMO_Get_Float(op, "depth");
572         int cap_ends = BMO_Get_Int(op, "cap_ends"), segs = BMO_Get_Int(op, "segments");
573         int cap_tris = BMO_Get_Int(op, "cap_tris");
574         int a;
575         
576         if (!segs)
577                 return;
578         
579         BMO_Get_Mat4(op, "mat", mat);
580
581         phid = 2.0f * (float)M_PI / segs;
582         phi = 0.25f * (float)M_PI;
583
584         depth *= 0.5f;
585         if (cap_ends) {
586                 vec[0] = vec[1] = 0.0f;
587                 vec[2] = -depth;
588                 mul_m4_v3(mat, vec);
589                 
590                 cent1 = BM_Make_Vert(bm, vec, NULL);
591
592                 vec[0] = vec[1] = 0.0f;
593                 vec[2] = depth;
594                 mul_m4_v3(mat, vec);
595                 
596                 cent2 = BM_Make_Vert(bm, vec, NULL);
597
598                 BMO_SetFlag(bm, cent1, VERT_MARK);
599                 BMO_SetFlag(bm, cent2, VERT_MARK);
600         }
601
602         for (a = 0; a < segs; a++, phi += phid) {
603                 vec[0] = dia1 * sinf(phi);
604                 vec[1] = dia1 * cosf(phi);
605                 vec[2] = -depth;
606                 mul_m4_v3(mat, vec);
607                 v1 = BM_Make_Vert(bm, vec, NULL);
608
609                 vec[0] = dia2 * sinf(phi);
610                 vec[1] = dia2 * cosf(phi);
611                 vec[2] = depth;
612                 mul_m4_v3(mat, vec);
613                 v2 = BM_Make_Vert(bm, vec, NULL);
614
615                 BMO_SetFlag(bm, v1, VERT_MARK);
616                 BMO_SetFlag(bm, v2, VERT_MARK);
617
618                 if (a) {
619                         if (cap_ends) {
620                                 BMFace *f;
621                                 
622                                 f = BM_Make_Face_QuadTri(bm, cent1, lastv1, v1, NULL, NULL, FALSE);
623                                 BMO_SetFlag(bm, f, FACE_NEW);
624                                 f = BM_Make_Face_QuadTri(bm, cent2, v2, lastv2, NULL, NULL, FALSE);
625                                 BMO_SetFlag(bm, f, FACE_NEW);
626                         }
627                         BM_Make_Face_QuadTri(bm, lastv1, lastv2, v2, v1, NULL, FALSE);
628                 }
629                 else {
630                         firstv1 = v1;
631                         firstv2 = v2;
632                 }
633
634                 lastv1 = v1;
635                 lastv2 = v2;
636         }
637
638         if (!a)
639                 return;
640
641         if (cap_ends) {
642                 BMFace *f;
643                 
644                 f = BM_Make_Face_QuadTri(bm, cent1, v1, firstv1, NULL, NULL, FALSE);
645                 BMO_SetFlag(bm, f, FACE_NEW);
646                 f = BM_Make_Face_QuadTri(bm, cent2, firstv2, v2, NULL, NULL, FALSE);
647                 BMO_SetFlag(bm, f, FACE_NEW);
648         }
649         
650         if (!cap_tris) {
651                 BMO_CallOpf(bm, "dissolvefaces faces=%ff", FACE_NEW);
652         }
653         
654         BM_Make_Face_QuadTri(bm, v1, v2, firstv2, firstv1, NULL, FALSE);
655
656         BMO_CallOpf(bm, "removedoubles verts=%fv dist=%f", VERT_MARK, 0.000001);
657         BMO_Flag_To_Slot(bm, op, "vertout", VERT_MARK, BM_VERT);
658 }
659
660 void bmesh_create_cube_exec(BMesh *bm, BMOperator *op)
661 {
662         BMVert *v1, *v2, *v3, *v4, *v5, *v6, *v7, *v8;
663         float vec[3], mat[4][4], off = BMO_Get_Float(op, "size") / 2.0f;
664
665         BMO_Get_Mat4(op, "mat", mat);
666
667         if (!off) off = 0.5f;
668
669         vec[0] = -off;
670         vec[1] = -off;
671         vec[2] = -off;
672         mul_m4_v3(mat, vec);
673         v1 = BM_Make_Vert(bm, vec, NULL);
674         BMO_SetFlag(bm, v1, VERT_MARK);
675
676         vec[0] = -off;
677         vec[1] = off;
678         vec[2] = -off;
679         mul_m4_v3(mat, vec);
680         v2 = BM_Make_Vert(bm, vec, NULL);
681         BMO_SetFlag(bm, v2, VERT_MARK);
682
683         vec[0] = off;
684         vec[1] = off;
685         vec[2] = -off;
686         mul_m4_v3(mat, vec);
687         v3 = BM_Make_Vert(bm, vec, NULL);
688         BMO_SetFlag(bm, v3, VERT_MARK);
689
690         vec[0] = off;
691         vec[1] = -off;
692         vec[2] = -off;
693         mul_m4_v3(mat, vec);
694         v4 = BM_Make_Vert(bm, vec, NULL);
695         BMO_SetFlag(bm, v4, VERT_MARK);
696
697         vec[0] = -off;
698         vec[1] = -off;
699         vec[2] = off;
700         mul_m4_v3(mat, vec);
701         v5 = BM_Make_Vert(bm, vec, NULL);
702         BMO_SetFlag(bm, v5, VERT_MARK);
703
704         vec[0] = -off;
705         vec[1] = off;
706         vec[2] = off;
707         mul_m4_v3(mat, vec);
708         v6 = BM_Make_Vert(bm, vec, NULL);
709         BMO_SetFlag(bm, v6, VERT_MARK);
710
711         vec[0] = off;
712         vec[1] = off;
713         vec[2] = off;
714         mul_m4_v3(mat, vec);
715         v7 = BM_Make_Vert(bm, vec, NULL);
716         BMO_SetFlag(bm, v7, VERT_MARK);
717
718         vec[0] = off;
719         vec[1] = -off;
720         vec[2] = off;
721         mul_m4_v3(mat, vec);
722         v8 = BM_Make_Vert(bm, vec, NULL);
723         BMO_SetFlag(bm, v8, VERT_MARK);
724
725         /* the four sides */
726         BM_Make_Face_QuadTri(bm, v5, v6, v2, v1, NULL, FALSE);
727         BM_Make_Face_QuadTri(bm, v6, v7, v3, v2, NULL, FALSE);
728         BM_Make_Face_QuadTri(bm, v7, v8, v4, v3, NULL, FALSE);
729         BM_Make_Face_QuadTri(bm, v8, v5, v1, v4, NULL, FALSE);
730         
731         /* top/bottom */
732         BM_Make_Face_QuadTri(bm, v1, v2, v3, v4, NULL, FALSE);
733         BM_Make_Face_QuadTri(bm, v8, v7, v6, v5, NULL, FALSE);
734
735         BMO_Flag_To_Slot(bm, op, "vertout", VERT_MARK, BM_VERT);
736 }