Dyntopo:
[blender.git] / source / blender / blenkernel / intern / customdata.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  * The Original Code is Copyright (C) 2006 Blender Foundation.
19  * All rights reserved.
20  *
21  * The Original Code is: all of this file.
22  *
23  * Contributor(s): Ben Batt <benbatt@gmail.com>
24  *
25  * ***** END GPL LICENSE BLOCK *****
26  *
27  * Implementation of CustomData.
28  *
29  * BKE_customdata.h contains the function prototypes for this file.
30  *
31  */
32
33 /** \file blender/blenkernel/intern/customdata.c
34  *  \ingroup bke
35  */
36  
37
38 #include <math.h>
39 #include <string.h>
40 #include <assert.h>
41
42 #include "MEM_guardedalloc.h"
43
44 #include "DNA_meshdata_types.h"
45 #include "DNA_ID.h"
46
47 #include "BLI_utildefines.h"
48 #include "BLI_string.h"
49 #include "BLI_path_util.h"
50 #include "BLI_math.h"
51 #include "BLI_mempool.h"
52 #include "BLI_alloca.h"
53
54 #include "BLF_translation.h"
55
56 #include "BKE_customdata.h"
57 #include "BKE_customdata_file.h"
58 #include "BKE_global.h"
59 #include "BKE_main.h"
60 #include "BKE_multires.h"
61
62 #include "bmesh.h"
63
64 #include <math.h>
65 #include <string.h>
66
67 /* number of layers to add when growing a CustomData object */
68 #define CUSTOMDATA_GROW 5
69
70 /* ensure typemap size is ok */
71 BLI_STATIC_ASSERT(sizeof(((CustomData *)NULL)->typemap) /
72                   sizeof(((CustomData *)NULL)->typemap[0]) == CD_NUMTYPES,
73                   "size mismatch");
74
75
76 /********************* Layer type information **********************/
77 typedef struct LayerTypeInfo {
78         int size;          /* the memory size of one element of this layer's data */
79
80         /** name of the struct used, for file writing */
81         const char *structname;
82         /** number of structs per element, for file writing */
83         int structnum;
84
85         /**
86          * default layer name.
87          * note! when NULL this is a way to ensure there is only ever one item
88          * see: CustomData_layertype_is_singleton() */
89         const char *defaultname;
90
91         /**
92          * a function to copy count elements of this layer's data
93          * (deep copy if appropriate)
94          * if NULL, memcpy is used
95          */
96         void (*copy)(const void *source, void *dest, int count);
97
98         /**
99          * a function to free any dynamically allocated components of this
100          * layer's data (note the data pointer itself should not be freed)
101          * size should be the size of one element of this layer's data (e.g.
102          * LayerTypeInfo.size)
103          */
104         void (*free)(void *data, int count, int size);
105
106         /**
107          * a function to interpolate between count source elements of this
108          * layer's data and store the result in dest
109          * if weights == NULL or sub_weights == NULL, they should default to 1
110          *
111          * weights gives the weight for each element in sources
112          * sub_weights gives the sub-element weights for each element in sources
113          *    (there should be (sub element count)^2 weights per element)
114          * count gives the number of elements in sources
115          *
116          * \note in some cases \a dest pointer is in \a sources
117          *       so all functions have to take this into account and delay
118          *       applying changes while reading from sources.
119          *       See bug [#32395] - Campbell.
120          */
121         void (*interp)(void **sources, const float *weights, const float *sub_weights,
122                        int count, void *dest);
123
124         /** a function to swap the data in corners of the element */
125         void (*swap)(void *data, const int *corner_indices);
126
127         /**
128          * a function to set a layer's data to default values. if NULL, the
129          * default is assumed to be all zeros */
130         void (*set_default)(void *data, int count);
131
132         /** functions necessary for geometry collapse */
133         bool (*equal)(void *data1, void *data2);
134         void (*multiply)(void *data, float fac);
135         void (*initminmax)(void *min, void *max);
136         void (*add)(void *data1, void *data2);
137         void (*dominmax)(void *data1, void *min, void *max);
138         void (*copyvalue)(void *source, void *dest);
139
140         /** a function to read data from a cdf file */
141         int (*read)(CDataFile *cdf, void *data, int count);
142
143         /** a function to write data to a cdf file */
144         int (*write)(CDataFile *cdf, void *data, int count);
145
146         /** a function to determine file size */
147         size_t (*filesize)(CDataFile *cdf, void *data, int count);
148 } LayerTypeInfo;
149
150 static void layerCopy_mdeformvert(const void *source, void *dest,
151                                   int count)
152 {
153         int i, size = sizeof(MDeformVert);
154
155         memcpy(dest, source, count * size);
156
157         for (i = 0; i < count; ++i) {
158                 MDeformVert *dvert = (MDeformVert *)((char *)dest + i * size);
159
160                 if (dvert->totweight) {
161                         MDeformWeight *dw = MEM_mallocN(dvert->totweight * sizeof(*dw),
162                                                         "layerCopy_mdeformvert dw");
163
164                         memcpy(dw, dvert->dw, dvert->totweight * sizeof(*dw));
165                         dvert->dw = dw;
166                 }
167                 else
168                         dvert->dw = NULL;
169         }
170 }
171
172 static void layerFree_mdeformvert(void *data, int count, int size)
173 {
174         int i;
175
176         for (i = 0; i < count; ++i) {
177                 MDeformVert *dvert = (MDeformVert *)((char *)data + i * size);
178
179                 if (dvert->dw) {
180                         MEM_freeN(dvert->dw);
181                         dvert->dw = NULL;
182                         dvert->totweight = 0;
183                 }
184         }
185 }
186
187 /* copy just zeros in this case */
188 static void layerCopy_bmesh_elem_py_ptr(const void *UNUSED(source), void *dest,
189                                         int count)
190 {
191         int i, size = sizeof(void *);
192
193         for (i = 0; i < count; ++i) {
194                 void **ptr = (void **)((char *)dest + i * size);
195                 *ptr = NULL;
196         }
197 }
198
199 #ifndef WITH_PYTHON
200 void bpy_bm_generic_invalidate(struct BPy_BMGeneric *UNUSED(self))
201 {
202         /* dummy */
203 }
204 #endif
205
206 static void layerFree_bmesh_elem_py_ptr(void *data, int count, int size)
207 {
208         int i;
209
210         for (i = 0; i < count; ++i) {
211                 void **ptr = (void *)((char *)data + i * size);
212                 if (*ptr) {
213                         bpy_bm_generic_invalidate(*ptr);
214                 }
215         }
216 }
217
218 static void layerInterp_mdeformvert(void **sources, const float *weights,
219                                     const float *UNUSED(sub_weights), int count, void *dest)
220 {
221         /* a single linked list of MDeformWeight's
222          * use this to avoid double allocs (which LinkNode would do) */
223         struct MDeformWeight_Link {
224                 struct MDeformWeight_Link *next;
225                 MDeformWeight dw;
226         };
227
228         MDeformVert *dvert = dest;
229         struct MDeformWeight_Link *dest_dwlink = NULL;
230         struct MDeformWeight_Link *node;
231         int i, j, totweight;
232
233         if (count <= 0) return;
234
235         /* build a list of unique def_nrs for dest */
236         totweight = 0;
237         for (i = 0; i < count; ++i) {
238                 MDeformVert *source = sources[i];
239                 float interp_weight = weights ? weights[i] : 1.0f;
240
241                 for (j = 0; j < source->totweight; ++j) {
242                         MDeformWeight *dw = &source->dw[j];
243                         float weight = dw->weight * interp_weight;
244
245                         if (weight == 0.0f)
246                                 continue;
247
248                         for (node = dest_dwlink; node; node = node->next) {
249                                 MDeformWeight *tmp_dw = &node->dw;
250
251                                 if (tmp_dw->def_nr == dw->def_nr) {
252                                         tmp_dw->weight += weight;
253                                         break;
254                                 }
255                         }
256
257                         /* if this def_nr is not in the list, add it */
258                         if (!node) {
259                                 struct MDeformWeight_Link *tmp_dwlink = alloca(sizeof(*tmp_dwlink));
260                                 tmp_dwlink->dw.def_nr = dw->def_nr;
261                                 tmp_dwlink->dw.weight = weight;
262
263                                 /* inline linklist */
264                                 tmp_dwlink->next = dest_dwlink;
265                                 dest_dwlink = tmp_dwlink;
266
267                                 totweight++;
268                         }
269                 }
270         }
271
272         /* delay writing to the destination incase dest is in sources */
273
274         /* now we know how many unique deform weights there are, so realloc */
275         if (dvert->dw && (dvert->totweight == totweight)) {
276                 /* pass (fastpath if we don't need to realloc) */
277         }
278         else {
279                 if (dvert->dw) {
280                         MEM_freeN(dvert->dw);
281                 }
282
283                 if (totweight) {
284                         dvert->dw = MEM_mallocN(sizeof(*dvert->dw) * totweight, __func__);
285                 }
286         }
287
288         if (totweight) {
289                 dvert->totweight = totweight;
290                 for (i = 0, node = dest_dwlink; node; node = node->next, i++) {
291                         dvert->dw[i] = node->dw;
292                 }
293         }
294         else {
295                 memset(dvert, 0, sizeof(*dvert));
296         }
297 }
298
299 static void layerCopy_tface(const void *source, void *dest, int count)
300 {
301         const MTFace *source_tf = (const MTFace *)source;
302         MTFace *dest_tf = (MTFace *)dest;
303         int i;
304
305         for (i = 0; i < count; ++i)
306                 dest_tf[i] = source_tf[i];
307 }
308
309 static void layerInterp_tface(void **sources, const float *weights,
310                               const float *sub_weights, int count, void *dest)
311 {
312         MTFace *tf = dest;
313         int i, j, k;
314         float uv[4][2] = {{0.0f}};
315         const float *sub_weight;
316
317         if (count <= 0) return;
318
319         sub_weight = sub_weights;
320         for (i = 0; i < count; ++i) {
321                 float weight = weights ? weights[i] : 1;
322                 MTFace *src = sources[i];
323
324                 for (j = 0; j < 4; ++j) {
325                         if (sub_weights) {
326                                 for (k = 0; k < 4; ++k, ++sub_weight) {
327                                         madd_v2_v2fl(uv[j], src->uv[k], (*sub_weight) * weight);
328                                 }
329                         }
330                         else {
331                                 madd_v2_v2fl(uv[j], src->uv[j], weight);
332                         }
333                 }
334         }
335
336         /* delay writing to the destination incase dest is in sources */
337         *tf = *(MTFace *)(*sources);
338         memcpy(tf->uv, uv, sizeof(tf->uv));
339 }
340
341 static void layerSwap_tface(void *data, const int *corner_indices)
342 {
343         MTFace *tf = data;
344         float uv[4][2];
345         static const short pin_flags[4] = { TF_PIN1, TF_PIN2, TF_PIN3, TF_PIN4 };
346         static const char sel_flags[4] = { TF_SEL1, TF_SEL2, TF_SEL3, TF_SEL4 };
347         short unwrap = tf->unwrap & ~(TF_PIN1 | TF_PIN2 | TF_PIN3 | TF_PIN4);
348         char flag = tf->flag & ~(TF_SEL1 | TF_SEL2 | TF_SEL3 | TF_SEL4);
349         int j;
350
351         for (j = 0; j < 4; ++j) {
352                 const int source_index = corner_indices[j];
353
354                 copy_v2_v2(uv[j], tf->uv[source_index]);
355
356                 /* swap pinning flags around */
357                 if (tf->unwrap & pin_flags[source_index]) {
358                         unwrap |= pin_flags[j];
359                 }
360
361                 /* swap selection flags around */
362                 if (tf->flag & sel_flags[source_index]) {
363                         flag |= sel_flags[j];
364                 }
365         }
366
367         memcpy(tf->uv, uv, sizeof(tf->uv));
368         tf->unwrap = unwrap;
369         tf->flag = flag;
370 }
371
372 static void layerDefault_tface(void *data, int count)
373 {
374         static MTFace default_tf = {{{0, 0}, {1, 0}, {1, 1}, {0, 1}}, NULL,
375                                         0, 0, TF_DYNAMIC | TF_CONVERTED, 0, 0};
376         MTFace *tf = (MTFace *)data;
377         int i;
378
379         for (i = 0; i < count; i++)
380                 tf[i] = default_tf;
381 }
382
383 static void layerCopy_propFloat(const void *source, void *dest,
384                                 int count)
385 {
386         memcpy(dest, source, sizeof(MFloatProperty) * count);
387 }
388
389 static void layerCopy_propInt(const void *source, void *dest,
390                               int count)
391 {
392         memcpy(dest, source, sizeof(MIntProperty) * count);
393 }
394
395 static void layerCopy_propString(const void *source, void *dest,
396                                  int count)
397 {
398         memcpy(dest, source, sizeof(MStringProperty) * count);
399 }
400
401 static void layerCopy_origspace_face(const void *source, void *dest, int count)
402 {
403         const OrigSpaceFace *source_tf = (const OrigSpaceFace *)source;
404         OrigSpaceFace *dest_tf = (OrigSpaceFace *)dest;
405         int i;
406
407         for (i = 0; i < count; ++i)
408                 dest_tf[i] = source_tf[i];
409 }
410
411 static void layerInterp_origspace_face(void **sources, const float *weights,
412                                        const float *sub_weights, int count, void *dest)
413 {
414         OrigSpaceFace *osf = dest;
415         int i, j, k;
416         float uv[4][2] = {{0.0f}};
417         const float *sub_weight;
418
419         if (count <= 0) return;
420
421         sub_weight = sub_weights;
422         for (i = 0; i < count; ++i) {
423                 float weight = weights ? weights[i] : 1;
424                 OrigSpaceFace *src = sources[i];
425
426                 for (j = 0; j < 4; ++j) {
427                         if (sub_weights) {
428                                 for (k = 0; k < 4; ++k, ++sub_weight) {
429                                         madd_v2_v2fl(uv[j], src->uv[k], (*sub_weight) * weight);
430                                 }
431                         }
432                         else {
433                                 madd_v2_v2fl(uv[j], src->uv[j], weight);
434                         }
435                 }
436         }
437
438         /* delay writing to the destination incase dest is in sources */
439
440 #if 0 /* no need, this ONLY contains UV's */
441         *osf = *(OrigSpaceFace *)(*sources);
442 #endif
443         memcpy(osf->uv, uv, sizeof(osf->uv));
444 }
445
446 static void layerSwap_origspace_face(void *data, const int *corner_indices)
447 {
448         OrigSpaceFace *osf = data;
449         float uv[4][2];
450         int j;
451
452         for (j = 0; j < 4; ++j) {
453                 copy_v2_v2(uv[j], osf->uv[corner_indices[j]]);
454         }
455         memcpy(osf->uv, uv, sizeof(osf->uv));
456 }
457
458 static void layerDefault_origspace_face(void *data, int count)
459 {
460         static OrigSpaceFace default_osf = {{{0, 0}, {1, 0}, {1, 1}, {0, 1}}};
461         OrigSpaceFace *osf = (OrigSpaceFace *)data;
462         int i;
463
464         for (i = 0; i < count; i++)
465                 osf[i] = default_osf;
466 }
467
468 static void layerSwap_mdisps(void *data, const int *ci)
469 {
470         MDisps *s = data;
471         float (*d)[3] = NULL;
472         int corners, cornersize, S;
473
474         if (s->disps) {
475                 int nverts = (ci[1] == 3) ? 4 : 3; /* silly way to know vertex count of face */
476                 corners = multires_mdisp_corners(s);
477                 cornersize = s->totdisp / corners;
478
479                 if (corners != nverts) {
480                         /* happens when face changed vertex count in edit mode
481                          * if it happened, just forgot displacement */
482
483                         MEM_freeN(s->disps);
484                         s->totdisp = (s->totdisp / corners) * nverts;
485                         s->disps = MEM_callocN(s->totdisp * sizeof(float) * 3, "mdisp swap");
486                         return;
487                 }
488
489                 d = MEM_callocN(sizeof(float) * 3 * s->totdisp, "mdisps swap");
490
491                 for (S = 0; S < corners; S++)
492                         memcpy(d + cornersize * S, s->disps + cornersize * ci[S], cornersize * 3 * sizeof(float));
493                 
494                 MEM_freeN(s->disps);
495                 s->disps = d;
496         }
497 }
498
499 static void layerCopy_mdisps(const void *source, void *dest, int count)
500 {
501         int i;
502         const MDisps *s = source;
503         MDisps *d = dest;
504
505         for (i = 0; i < count; ++i) {
506                 if (s[i].disps) {
507                         d[i].disps = MEM_dupallocN(s[i].disps);
508                         d[i].hidden = MEM_dupallocN(s[i].hidden);
509                 }
510                 else {
511                         d[i].disps = NULL;
512                         d[i].hidden = NULL;
513                 }
514
515                 /* still copy even if not in memory, displacement can be external */
516                 d[i].totdisp = s[i].totdisp;
517                 d[i].level = s[i].level;
518         }
519 }
520
521 static void layerFree_mdisps(void *data, int count, int UNUSED(size))
522 {
523         int i;
524         MDisps *d = data;
525
526         for (i = 0; i < count; ++i) {
527                 if (d[i].disps)
528                         MEM_freeN(d[i].disps);
529                 if (d[i].hidden)
530                         MEM_freeN(d[i].hidden);
531                 d[i].disps = NULL;
532                 d[i].hidden = NULL;
533                 d[i].totdisp = 0;
534                 d[i].level = 0;
535         }
536 }
537
538 static int layerRead_mdisps(CDataFile *cdf, void *data, int count)
539 {
540         MDisps *d = data;
541         int i;
542
543         for (i = 0; i < count; ++i) {
544                 if (!d[i].disps)
545                         d[i].disps = MEM_callocN(sizeof(float) * 3 * d[i].totdisp, "mdisps read");
546
547                 if (!cdf_read_data(cdf, d[i].totdisp * 3 * sizeof(float), d[i].disps)) {
548                         printf("failed to read multires displacement %d/%d %d\n", i, count, d[i].totdisp);
549                         return 0;
550                 }
551         }
552
553         return 1;
554 }
555
556 static int layerWrite_mdisps(CDataFile *cdf, void *data, int count)
557 {
558         MDisps *d = data;
559         int i;
560
561         for (i = 0; i < count; ++i) {
562                 if (!cdf_write_data(cdf, d[i].totdisp * 3 * sizeof(float), d[i].disps)) {
563                         printf("failed to write multires displacement %d/%d %d\n", i, count, d[i].totdisp);
564                         return 0;
565                 }
566         }
567
568         return 1;
569 }
570
571 static size_t layerFilesize_mdisps(CDataFile *UNUSED(cdf), void *data, int count)
572 {
573         MDisps *d = data;
574         size_t size = 0;
575         int i;
576
577         for (i = 0; i < count; ++i)
578                 size += d[i].totdisp * 3 * sizeof(float);
579
580         return size;
581 }
582
583 static void layerCopy_grid_paint_mask(const void *source, void *dest, int count)
584 {
585         int i;
586         const GridPaintMask *s = source;
587         GridPaintMask *d = dest;
588
589         for (i = 0; i < count; ++i) {
590                 if (s[i].data) {
591                         d[i].data = MEM_dupallocN(s[i].data);
592                         d[i].level = s[i].level;
593                 }
594                 else {
595                         d[i].data = NULL;
596                         d[i].level = 0;
597                 }
598                 
599         }
600 }
601
602 static void layerFree_grid_paint_mask(void *data, int count, int UNUSED(size))
603 {
604         int i;
605         GridPaintMask *gpm = data;
606
607         for (i = 0; i < count; ++i) {
608                 if (gpm[i].data)
609                         MEM_freeN(gpm[i].data);
610                 gpm[i].data = NULL;
611                 gpm[i].level = 0;
612         }
613 }
614
615 /* --------- */
616 static void layerCopyValue_mloopcol(void *source, void *dest)
617 {
618         MLoopCol *m1 = source, *m2 = dest;
619         
620         m2->r = m1->r;
621         m2->g = m1->g;
622         m2->b = m1->b;
623         m2->a = m1->a;
624 }
625
626 static bool layerEqual_mloopcol(void *data1, void *data2)
627 {
628         MLoopCol *m1 = data1, *m2 = data2;
629         float r, g, b, a;
630
631         r = m1->r - m2->r;
632         g = m1->g - m2->g;
633         b = m1->b - m2->b;
634         a = m1->a - m2->a;
635
636         return r * r + g * g + b * b + a * a < 0.001f;
637 }
638
639 static void layerMultiply_mloopcol(void *data, float fac)
640 {
641         MLoopCol *m = data;
642
643         m->r = (float)m->r * fac;
644         m->g = (float)m->g * fac;
645         m->b = (float)m->b * fac;
646         m->a = (float)m->a * fac;
647 }
648
649 static void layerAdd_mloopcol(void *data1, void *data2)
650 {
651         MLoopCol *m = data1, *m2 = data2;
652
653         m->r += m2->r;
654         m->g += m2->g;
655         m->b += m2->b;
656         m->a += m2->a;
657 }
658
659 static void layerDoMinMax_mloopcol(void *data, void *vmin, void *vmax)
660 {
661         MLoopCol *m = data;
662         MLoopCol *min = vmin, *max = vmax;
663
664         if (m->r < min->r) min->r = m->r;
665         if (m->g < min->g) min->g = m->g;
666         if (m->b < min->b) min->b = m->b;
667         if (m->a < min->a) min->a = m->a;
668         
669         if (m->r > max->r) max->r = m->r;
670         if (m->g > max->g) max->g = m->g;
671         if (m->b > max->b) max->b = m->b;
672         if (m->a > max->a) max->a = m->a;
673 }
674
675 static void layerInitMinMax_mloopcol(void *vmin, void *vmax)
676 {
677         MLoopCol *min = vmin, *max = vmax;
678
679         min->r = 255;
680         min->g = 255;
681         min->b = 255;
682         min->a = 255;
683
684         max->r = 0;
685         max->g = 0;
686         max->b = 0;
687         max->a = 0;
688 }
689
690 static void layerDefault_mloopcol(void *data, int count)
691 {
692         MLoopCol default_mloopcol = {255, 255, 255, 255};
693         MLoopCol *mlcol = (MLoopCol *)data;
694         int i;
695         for (i = 0; i < count; i++)
696                 mlcol[i] = default_mloopcol;
697
698 }
699
700 static void layerInterp_mloopcol(void **sources, const float *weights,
701                                  const float *sub_weights, int count, void *dest)
702 {
703         MLoopCol *mc = dest;
704         int i;
705         const float *sub_weight;
706         struct {
707                 float a;
708                 float r;
709                 float g;
710                 float b;
711         } col;
712         col.a = col.r = col.g = col.b = 0;
713
714         sub_weight = sub_weights;
715         for (i = 0; i < count; ++i) {
716                 float weight = weights ? weights[i] : 1;
717                 MLoopCol *src = sources[i];
718                 if (sub_weights) {
719                         col.r += src->r * (*sub_weight) * weight;
720                         col.g += src->g * (*sub_weight) * weight;
721                         col.b += src->b * (*sub_weight) * weight;
722                         col.a += src->a * (*sub_weight) * weight;
723                         sub_weight++;
724                 }
725                 else {
726                         col.r += src->r * weight;
727                         col.g += src->g * weight;
728                         col.b += src->b * weight;
729                         col.a += src->a * weight;
730                 }
731         }
732         
733         /* Subdivide smooth or fractal can cause problems without clamping
734          * although weights should also not cause this situation */
735         CLAMP(col.a, 0.0f, 255.0f);
736         CLAMP(col.r, 0.0f, 255.0f);
737         CLAMP(col.g, 0.0f, 255.0f);
738         CLAMP(col.b, 0.0f, 255.0f);
739
740         /* delay writing to the destination incase dest is in sources */
741         mc->r = (int)col.r;
742         mc->g = (int)col.g;
743         mc->b = (int)col.b;
744         mc->a = (int)col.a;
745 }
746
747 static void layerCopyValue_mloopuv(void *source, void *dest)
748 {
749         MLoopUV *luv1 = source, *luv2 = dest;
750
751         copy_v2_v2(luv2->uv, luv1->uv);
752 }
753
754 static bool layerEqual_mloopuv(void *data1, void *data2)
755 {
756         MLoopUV *luv1 = data1, *luv2 = data2;
757
758         return len_squared_v2v2(luv1->uv, luv2->uv) < 0.00001f;
759 }
760
761 static void layerMultiply_mloopuv(void *data, float fac)
762 {
763         MLoopUV *luv = data;
764
765         mul_v2_fl(luv->uv, fac);
766 }
767
768 static void layerInitMinMax_mloopuv(void *vmin, void *vmax)
769 {
770         MLoopUV *min = vmin, *max = vmax;
771
772         INIT_MINMAX2(min->uv, max->uv);
773 }
774
775 static void layerDoMinMax_mloopuv(void *data, void *vmin, void *vmax)
776 {
777         MLoopUV *min = vmin, *max = vmax, *luv = data;
778
779         minmax_v2v2_v2(min->uv, max->uv, luv->uv);
780 }
781
782 static void layerAdd_mloopuv(void *data1, void *data2)
783 {
784         MLoopUV *l1 = data1, *l2 = data2;
785
786         add_v2_v2(l1->uv, l2->uv);
787 }
788
789 static void layerInterp_mloopuv(void **sources, const float *weights,
790                                 const float *sub_weights, int count, void *dest)
791 {
792         float uv[2];
793         int i;
794
795         zero_v2(uv);
796
797         if (sub_weights) {
798                 const float *sub_weight = sub_weights;
799                 for (i = 0; i < count; i++) {
800                         float weight = weights ? weights[i] : 1.0f;
801                         MLoopUV *src = sources[i];
802                         madd_v2_v2fl(uv, src->uv, (*sub_weight) * weight);
803                         sub_weight++;
804                 }
805         }
806         else {
807                 for (i = 0; i < count; i++) {
808                         float weight = weights ? weights[i] : 1;
809                         MLoopUV *src = sources[i];
810                         madd_v2_v2fl(uv, src->uv, weight);
811                 }
812         }
813
814         /* delay writing to the destination incase dest is in sources */
815         copy_v2_v2(((MLoopUV *)dest)->uv, uv);
816 }
817
818 /* origspace is almost exact copy of mloopuv's, keep in sync */
819 static void layerCopyValue_mloop_origspace(void *source, void *dest)
820 {
821         OrigSpaceLoop *luv1 = source, *luv2 = dest;
822
823         copy_v2_v2(luv2->uv, luv1->uv);
824 }
825
826 static bool layerEqual_mloop_origspace(void *data1, void *data2)
827 {
828         OrigSpaceLoop *luv1 = data1, *luv2 = data2;
829
830         return len_squared_v2v2(luv1->uv, luv2->uv) < 0.00001f;
831 }
832
833 static void layerMultiply_mloop_origspace(void *data, float fac)
834 {
835         OrigSpaceLoop *luv = data;
836
837         mul_v2_fl(luv->uv, fac);
838 }
839
840 static void layerInitMinMax_mloop_origspace(void *vmin, void *vmax)
841 {
842         OrigSpaceLoop *min = vmin, *max = vmax;
843
844         INIT_MINMAX2(min->uv, max->uv);
845 }
846
847 static void layerDoMinMax_mloop_origspace(void *data, void *vmin, void *vmax)
848 {
849         OrigSpaceLoop *min = vmin, *max = vmax, *luv = data;
850
851         minmax_v2v2_v2(min->uv, max->uv, luv->uv);
852 }
853
854 static void layerAdd_mloop_origspace(void *data1, void *data2)
855 {
856         OrigSpaceLoop *l1 = data1, *l2 = data2;
857
858         add_v2_v2(l1->uv, l2->uv);
859 }
860
861 static void layerInterp_mloop_origspace(void **sources, const float *weights,
862                                         const float *sub_weights, int count, void *dest)
863 {
864         float uv[2];
865         int i;
866
867         zero_v2(uv);
868
869         if (sub_weights) {
870                 const float *sub_weight = sub_weights;
871                 for (i = 0; i < count; i++) {
872                         float weight = weights ? weights[i] : 1.0f;
873                         OrigSpaceLoop *src = sources[i];
874                         madd_v2_v2fl(uv, src->uv, (*sub_weight) * weight);
875                         sub_weight++;
876                 }
877         }
878         else {
879                 for (i = 0; i < count; i++) {
880                         float weight = weights ? weights[i] : 1.0f;
881                         OrigSpaceLoop *src = sources[i];
882                         madd_v2_v2fl(uv, src->uv, weight);
883                 }
884         }
885
886         /* delay writing to the destination incase dest is in sources */
887         copy_v2_v2(((OrigSpaceLoop *)dest)->uv, uv);
888 }
889 /* --- end copy */
890
891 static void layerInterp_mcol(void **sources, const float *weights,
892                              const float *sub_weights, int count, void *dest)
893 {
894         MCol *mc = dest;
895         int i, j, k;
896         struct {
897                 float a;
898                 float r;
899                 float g;
900                 float b;
901         } col[4] = {{0.0f}};
902
903         const float *sub_weight;
904
905         if (count <= 0) return;
906         
907         sub_weight = sub_weights;
908         for (i = 0; i < count; ++i) {
909                 float weight = weights ? weights[i] : 1;
910
911                 for (j = 0; j < 4; ++j) {
912                         if (sub_weights) {
913                                 MCol *src = sources[i];
914                                 for (k = 0; k < 4; ++k, ++sub_weight, ++src) {
915                                         const float w = (*sub_weight) * weight;
916                                         col[j].a += src->a * w;
917                                         col[j].r += src->r * w;
918                                         col[j].g += src->g * w;
919                                         col[j].b += src->b * w;
920                                 }
921                         }
922                         else {
923                                 MCol *src = sources[i];
924                                 col[j].a += src[j].a * weight;
925                                 col[j].r += src[j].r * weight;
926                                 col[j].g += src[j].g * weight;
927                                 col[j].b += src[j].b * weight;
928                         }
929                 }
930         }
931
932         /* delay writing to the destination incase dest is in sources */
933         for (j = 0; j < 4; ++j) {
934                 
935                 /* Subdivide smooth or fractal can cause problems without clamping
936                  * although weights should also not cause this situation */
937                 CLAMP(col[j].a, 0.0f, 255.0f);
938                 CLAMP(col[j].r, 0.0f, 255.0f);
939                 CLAMP(col[j].g, 0.0f, 255.0f);
940                 CLAMP(col[j].b, 0.0f, 255.0f);
941                 
942                 mc[j].a = (int)col[j].a;
943                 mc[j].r = (int)col[j].r;
944                 mc[j].g = (int)col[j].g;
945                 mc[j].b = (int)col[j].b;
946         }
947 }
948
949 static void layerSwap_mcol(void *data, const int *corner_indices)
950 {
951         MCol *mcol = data;
952         MCol col[4];
953         int j;
954
955         for (j = 0; j < 4; ++j)
956                 col[j] = mcol[corner_indices[j]];
957
958         memcpy(mcol, col, sizeof(col));
959 }
960
961 static void layerDefault_mcol(void *data, int count)
962 {
963         static MCol default_mcol = {255, 255, 255, 255};
964         MCol *mcol = (MCol *)data;
965         int i;
966
967         for (i = 0; i < 4 * count; i++) {
968                 mcol[i] = default_mcol;
969         }
970 }
971
972 static void layerDefault_origindex(void *data, int count)
973 {
974         fill_vn_i((int *)data, count, ORIGINDEX_NONE);
975 }
976
977 static void layerInterp_bweight(void **sources, const float *weights,
978                                 const float *UNUSED(sub_weights), int count, void *dest)
979 {
980         float f;
981         float **in = (float **)sources;
982         int i;
983         
984         if (count <= 0) return;
985
986         f = 0.0f;
987
988         if (weights) {
989                 for (i = 0; i < count; ++i) {
990                         f += *in[i] * weights[i];
991                 }
992         }
993         else {
994                 for (i = 0; i < count; ++i) {
995                         f += *in[i];
996                 }
997         }
998
999         /* delay writing to the destination incase dest is in sources */
1000         *((float *)dest) = f;
1001 }
1002
1003 static void layerInterp_shapekey(void **sources, const float *weights,
1004                                  const float *UNUSED(sub_weights), int count, void *dest)
1005 {
1006         float co[3];
1007         float **in = (float **)sources;
1008         int i;
1009
1010         if (count <= 0) return;
1011
1012         zero_v3(co);
1013
1014         if (weights) {
1015                 for (i = 0; i < count; ++i) {
1016                         madd_v3_v3fl(co, in[i], weights[i]);
1017                 }
1018         }
1019         else {
1020                 for (i = 0; i < count; ++i) {
1021                         add_v3_v3(co, in[i]);
1022                 }
1023         }
1024
1025         /* delay writing to the destination incase dest is in sources */
1026         copy_v3_v3((float *)dest, co);
1027 }
1028
1029 static void layerDefault_mvert_skin(void *data, int count)
1030 {
1031         MVertSkin *vs = data;
1032         int i;
1033         
1034         for (i = 0; i < count; i++) {
1035                 copy_v3_fl(vs[i].radius, 0.25f);
1036                 vs[i].flag = 0;
1037         }
1038 }
1039
1040 static void layerDefault_dyntopo_node(void *data, int count)
1041 {
1042         int *indices = data;
1043         int i;
1044
1045         for (i = 0; i < count; i++) {
1046                 indices[i] = DYNTOPO_NODE_NONE;
1047         }
1048 }
1049
1050
1051 static void layerInterp_mvert_skin(void **sources, const float *weights,
1052                                    const float *UNUSED(sub_weights),
1053                                    int count, void *dest)
1054 {
1055         float radius[3], w;
1056         MVertSkin *vs;
1057         int i;
1058
1059         zero_v3(radius);
1060         for (i = 0; i < count; i++) {
1061                 w = weights ? weights[i] : 1.0f;
1062                 vs = sources[i];
1063
1064                 madd_v3_v3fl(radius, vs->radius, w);
1065         }
1066
1067         /* delay writing to the destination incase dest is in sources */
1068         vs = dest;
1069         copy_v3_v3(vs->radius, radius);
1070         vs->flag &= ~MVERT_SKIN_ROOT;
1071 }
1072
1073 static const LayerTypeInfo LAYERTYPEINFO[CD_NUMTYPES] = {
1074         /* 0: CD_MVERT */
1075         {sizeof(MVert), "MVert", 1, NULL, NULL, NULL, NULL, NULL, NULL},
1076         /* 1: CD_MSTICKY */  /* DEPRECATED */
1077         {sizeof(float) * 2, "", 1, NULL, NULL, NULL, NULL, NULL,
1078          NULL},
1079         /* 2: CD_MDEFORMVERT */
1080         {sizeof(MDeformVert), "MDeformVert", 1, NULL, layerCopy_mdeformvert,
1081          layerFree_mdeformvert, layerInterp_mdeformvert, NULL, NULL},
1082         /* 3: CD_MEDGE */
1083         {sizeof(MEdge), "MEdge", 1, NULL, NULL, NULL, NULL, NULL, NULL},
1084         /* 4: CD_MFACE */
1085         {sizeof(MFace), "MFace", 1, NULL, NULL, NULL, NULL, NULL, NULL},
1086         /* 5: CD_MTFACE */
1087         {sizeof(MTFace), "MTFace", 1, N_("UVMap"), layerCopy_tface, NULL,
1088          layerInterp_tface, layerSwap_tface, layerDefault_tface},
1089         /* 6: CD_MCOL */
1090         /* 4 MCol structs per face */
1091         {sizeof(MCol) * 4, "MCol", 4, N_("Col"), NULL, NULL, layerInterp_mcol,
1092          layerSwap_mcol, layerDefault_mcol},
1093         /* 7: CD_ORIGINDEX */
1094         {sizeof(int), "", 0, NULL, NULL, NULL, NULL, NULL, layerDefault_origindex},
1095         /* 8: CD_NORMAL */
1096         /* 3 floats per normal vector */
1097         {sizeof(float) * 3, "vec3f", 1, NULL, NULL, NULL, NULL, NULL, NULL},
1098         /* 9: CD_POLYINDEX (deprecated) */
1099         {sizeof(int), "", 0, NULL, NULL, NULL, NULL, NULL, NULL},
1100         /* 10: CD_PROP_FLT */
1101         {sizeof(MFloatProperty), "MFloatProperty", 1, N_("Float"), layerCopy_propFloat, NULL, NULL, NULL},
1102         /* 11: CD_PROP_INT */
1103         {sizeof(MIntProperty), "MIntProperty", 1, N_("Int"), layerCopy_propInt, NULL, NULL, NULL},
1104         /* 12: CD_PROP_STR */
1105         {sizeof(MStringProperty), "MStringProperty", 1, N_("String"), layerCopy_propString, NULL, NULL, NULL},
1106         /* 13: CD_ORIGSPACE */
1107         {sizeof(OrigSpaceFace), "OrigSpaceFace", 1, N_("UVMap"), layerCopy_origspace_face, NULL,
1108          layerInterp_origspace_face, layerSwap_origspace_face, layerDefault_origspace_face},
1109         /* 14: CD_ORCO */
1110         {sizeof(float) * 3, "", 0, NULL, NULL, NULL, NULL, NULL, NULL},
1111         /* 15: CD_MTEXPOLY */
1112         /* note, when we expose the UV Map / TexFace split to the user, change this back to face Texture */
1113         {sizeof(MTexPoly), "MTexPoly", 1, N_("UVMap") /* "Face Texture" */, NULL, NULL, NULL, NULL, NULL},
1114         /* 16: CD_MLOOPUV */
1115         {sizeof(MLoopUV), "MLoopUV", 1, N_("UVMap"), NULL, NULL, layerInterp_mloopuv, NULL, NULL,
1116          layerEqual_mloopuv, layerMultiply_mloopuv, layerInitMinMax_mloopuv, 
1117          layerAdd_mloopuv, layerDoMinMax_mloopuv, layerCopyValue_mloopuv},
1118         /* 17: CD_MLOOPCOL */
1119         {sizeof(MLoopCol), "MLoopCol", 1, N_("Col"), NULL, NULL, layerInterp_mloopcol, NULL,
1120          layerDefault_mloopcol, layerEqual_mloopcol, layerMultiply_mloopcol, layerInitMinMax_mloopcol, 
1121          layerAdd_mloopcol, layerDoMinMax_mloopcol, layerCopyValue_mloopcol},
1122         /* 18: CD_TANGENT */
1123         {sizeof(float) * 4 * 4, "", 0, NULL, NULL, NULL, NULL, NULL, NULL},
1124         /* 19: CD_MDISPS */
1125         {sizeof(MDisps), "MDisps", 1, NULL, layerCopy_mdisps,
1126          layerFree_mdisps, NULL, layerSwap_mdisps, NULL,
1127          NULL, NULL, NULL, NULL, NULL, NULL, 
1128          layerRead_mdisps, layerWrite_mdisps, layerFilesize_mdisps},
1129         /* 20: CD_PREVIEW_MCOL */
1130         {sizeof(MCol) * 4, "MCol", 4, N_("PreviewCol"), NULL, NULL, layerInterp_mcol,
1131          layerSwap_mcol, layerDefault_mcol},
1132         /* 21: CD_ID_MCOL */
1133         {sizeof(MCol) * 4, "MCol", 4, N_("IDCol"), NULL, NULL, layerInterp_mcol,
1134          layerSwap_mcol, layerDefault_mcol},
1135         /* 22: CD_TEXTURE_MCOL */
1136         {sizeof(MCol) * 4, "MCol", 4, N_("TexturedCol"), NULL, NULL, layerInterp_mcol,
1137          layerSwap_mcol, layerDefault_mcol},
1138         /* 23: CD_CLOTH_ORCO */
1139         {sizeof(float) * 3, "", 0, NULL, NULL, NULL, NULL, NULL, NULL},
1140         /* 24: CD_RECAST */
1141         {sizeof(MRecast), "MRecast", 1, N_("Recast"), NULL, NULL, NULL, NULL},
1142
1143 /* BMESH ONLY */
1144         /* 25: CD_MPOLY */
1145         {sizeof(MPoly), "MPoly", 1, N_("NGon Face"), NULL, NULL, NULL, NULL, NULL},
1146         /* 26: CD_MLOOP */
1147         {sizeof(MLoop), "MLoop", 1, N_("NGon Face-Vertex"), NULL, NULL, NULL, NULL, NULL},
1148         /* 27: CD_SHAPE_KEYINDEX */
1149         {sizeof(int), "", 0, NULL, NULL, NULL, NULL, NULL, NULL},
1150         /* 28: CD_SHAPEKEY */
1151         {sizeof(float) * 3, "", 0, N_("ShapeKey"), NULL, NULL, layerInterp_shapekey},
1152         /* 29: CD_BWEIGHT */
1153         {sizeof(float), "", 0, N_("BevelWeight"), NULL, NULL, layerInterp_bweight},
1154         /* 30: CD_CREASE */
1155         {sizeof(float), "", 0, N_("SubSurfCrease"), NULL, NULL, layerInterp_bweight},
1156         /* 31: CD_ORIGSPACE_MLOOP */
1157         {sizeof(OrigSpaceLoop), "OrigSpaceLoop", 1, N_("OS Loop"), NULL, NULL, layerInterp_mloop_origspace, NULL, NULL,
1158          layerEqual_mloop_origspace, layerMultiply_mloop_origspace, layerInitMinMax_mloop_origspace,
1159          layerAdd_mloop_origspace, layerDoMinMax_mloop_origspace, layerCopyValue_mloop_origspace},
1160         /* 32: CD_PREVIEW_MLOOPCOL */
1161         {sizeof(MLoopCol), "MLoopCol", 1, N_("PreviewLoopCol"), NULL, NULL, layerInterp_mloopcol, NULL,
1162          layerDefault_mloopcol, layerEqual_mloopcol, layerMultiply_mloopcol, layerInitMinMax_mloopcol,
1163          layerAdd_mloopcol, layerDoMinMax_mloopcol, layerCopyValue_mloopcol},
1164         /* 33: CD_BM_ELEM_PYPTR */
1165         {sizeof(void *), "", 1, NULL, layerCopy_bmesh_elem_py_ptr,
1166          layerFree_bmesh_elem_py_ptr, NULL, NULL, NULL},
1167
1168 /* END BMESH ONLY */
1169
1170         /* 34: CD_PAINT_MASK */
1171         {sizeof(float), "", 0, NULL, NULL, NULL, NULL, NULL, NULL},
1172         /* 35: CD_GRID_PAINT_MASK */
1173         {sizeof(GridPaintMask), "GridPaintMask", 1, NULL, layerCopy_grid_paint_mask,
1174          layerFree_grid_paint_mask, NULL, NULL, NULL},
1175         /* 36: CD_SKIN_NODE */
1176         {sizeof(MVertSkin), "MVertSkin", 1, NULL, NULL, NULL,
1177          layerInterp_mvert_skin, NULL, layerDefault_mvert_skin},
1178         /* 37: CD_FREESTYLE_EDGE */
1179         {sizeof(FreestyleEdge), "FreestyleEdge", 1, NULL, NULL, NULL, NULL, NULL, NULL},
1180         /* 38: CD_FREESTYLE_FACE */
1181         {sizeof(FreestyleFace), "FreestyleFace", 1, NULL, NULL, NULL, NULL, NULL, NULL},
1182         /* 39: CD_MLOOPTANGENT */
1183         {sizeof(float[4]), "", 0, NULL, NULL, NULL, NULL, NULL, NULL},
1184         /* 40: CD_TESSLOOPNORMAL */
1185         {sizeof(short[4][3]), "", 0, NULL, NULL, NULL, NULL, NULL, NULL},
1186     /* 41: CD_DYNTOPO_NODE */
1187         {sizeof(int), "", 0, NULL, NULL, NULL, NULL, NULL, layerDefault_dyntopo_node},
1188 };
1189
1190 /* note, numbers are from trunk and need updating for bmesh */
1191
1192 static const char *LAYERTYPENAMES[CD_NUMTYPES] = {
1193         /*   0-4 */ "CDMVert", "CDMSticky", "CDMDeformVert", "CDMEdge", "CDMFace",
1194         /*   5-9 */ "CDMTFace", "CDMCol", "CDOrigIndex", "CDNormal", "CDFlags",
1195         /* 10-14 */ "CDMFloatProperty", "CDMIntProperty", "CDMStringProperty", "CDOrigSpace", "CDOrco",
1196         /* 15-19 */ "CDMTexPoly", "CDMLoopUV", "CDMloopCol", "CDTangent", "CDMDisps",
1197         /* 20-24 */ "CDPreviewMCol", "CDIDMCol", "CDTextureMCol", "CDClothOrco", "CDMRecast",
1198
1199 /* BMESH ONLY */
1200         /* 25-29 */ "CDMPoly", "CDMLoop", "CDShapeKeyIndex", "CDShapeKey", "CDBevelWeight",
1201         /* 30-34 */ "CDSubSurfCrease", "CDOrigSpaceLoop", "CDPreviewLoopCol", "CDBMElemPyPtr", "CDPaintMask",
1202         /* 35-36 */ "CDGridPaintMask", "CDMVertSkin",
1203         /* 37-40 */ "CDFreestyleEdge", "CDFreestyleFace", "CDMLoopTangent", "CDTessLoopNormal",
1204     /* 41 */ "CDDyntopoNode"
1205 };
1206
1207
1208 const CustomDataMask CD_MASK_BAREMESH =
1209     CD_MASK_MVERT | CD_MASK_MEDGE | CD_MASK_MFACE | CD_MASK_MLOOP | CD_MASK_MPOLY | CD_MASK_BWEIGHT;
1210 const CustomDataMask CD_MASK_MESH =
1211     CD_MASK_MVERT | CD_MASK_MEDGE | CD_MASK_MFACE |
1212     CD_MASK_MSTICKY | CD_MASK_MDEFORMVERT | CD_MASK_MTFACE | CD_MASK_MCOL |
1213     CD_MASK_PROP_FLT | CD_MASK_PROP_INT | CD_MASK_PROP_STR | CD_MASK_MDISPS |
1214     CD_MASK_MLOOPUV | CD_MASK_MLOOPCOL | CD_MASK_MPOLY | CD_MASK_MLOOP |
1215     CD_MASK_MTEXPOLY | CD_MASK_RECAST | CD_MASK_PAINT_MASK |
1216     CD_MASK_GRID_PAINT_MASK | CD_MASK_MVERT_SKIN | CD_MASK_FREESTYLE_EDGE | CD_MASK_FREESTYLE_FACE;
1217 const CustomDataMask CD_MASK_EDITMESH =
1218     CD_MASK_MSTICKY | CD_MASK_MDEFORMVERT | CD_MASK_MTFACE | CD_MASK_MLOOPUV |
1219     CD_MASK_MLOOPCOL | CD_MASK_MTEXPOLY | CD_MASK_SHAPE_KEYINDEX |
1220     CD_MASK_MCOL | CD_MASK_PROP_FLT | CD_MASK_PROP_INT | CD_MASK_PROP_STR |
1221     CD_MASK_MDISPS | CD_MASK_SHAPEKEY | CD_MASK_RECAST | CD_MASK_PAINT_MASK |
1222     CD_MASK_GRID_PAINT_MASK | CD_MASK_MVERT_SKIN;
1223 const CustomDataMask CD_MASK_DERIVEDMESH =
1224     CD_MASK_MSTICKY | CD_MASK_MDEFORMVERT | CD_MASK_MTFACE |
1225     CD_MASK_MCOL | CD_MASK_PROP_FLT | CD_MASK_PROP_INT | CD_MASK_CLOTH_ORCO |
1226     CD_MASK_MLOOPUV | CD_MASK_MLOOPCOL | CD_MASK_MTEXPOLY | CD_MASK_PREVIEW_MLOOPCOL |
1227     CD_MASK_PROP_STR | CD_MASK_ORIGSPACE | CD_MASK_ORIGSPACE_MLOOP | CD_MASK_ORCO | CD_MASK_TANGENT |
1228     CD_MASK_PREVIEW_MCOL | CD_MASK_SHAPEKEY | CD_MASK_RECAST |
1229     CD_MASK_ORIGINDEX | CD_MASK_MVERT_SKIN | CD_MASK_FREESTYLE_EDGE | CD_MASK_FREESTYLE_FACE;
1230 const CustomDataMask CD_MASK_BMESH =
1231     CD_MASK_MLOOPUV | CD_MASK_MLOOPCOL | CD_MASK_MTEXPOLY |
1232     CD_MASK_MSTICKY | CD_MASK_MDEFORMVERT | CD_MASK_PROP_FLT | CD_MASK_PROP_INT |
1233     CD_MASK_PROP_STR | CD_MASK_SHAPEKEY | CD_MASK_SHAPE_KEYINDEX | CD_MASK_MDISPS |
1234     CD_MASK_CREASE | CD_MASK_BWEIGHT | CD_MASK_RECAST | CD_MASK_PAINT_MASK |
1235     CD_MASK_GRID_PAINT_MASK | CD_MASK_MVERT_SKIN | CD_MASK_FREESTYLE_EDGE | CD_MASK_FREESTYLE_FACE;
1236 const CustomDataMask CD_MASK_FACECORNERS =  /* XXX Not used anywhere! */
1237     CD_MASK_MTFACE | CD_MASK_MCOL | CD_MASK_MTEXPOLY | CD_MASK_MLOOPUV |
1238     CD_MASK_MLOOPCOL | CD_MASK_NORMAL | CD_MASK_MLOOPTANGENT;
1239 const CustomDataMask CD_MASK_EVERYTHING =
1240     CD_MASK_MVERT | CD_MASK_MSTICKY /* DEPRECATED */ | CD_MASK_MDEFORMVERT | CD_MASK_MEDGE | CD_MASK_MFACE |
1241     CD_MASK_MTFACE | CD_MASK_MCOL | CD_MASK_ORIGINDEX | CD_MASK_NORMAL /* | CD_MASK_POLYINDEX */ | CD_MASK_PROP_FLT |
1242     CD_MASK_PROP_INT | CD_MASK_PROP_STR | CD_MASK_ORIGSPACE | CD_MASK_ORCO | CD_MASK_MTEXPOLY | CD_MASK_MLOOPUV |
1243     CD_MASK_MLOOPCOL | CD_MASK_TANGENT | CD_MASK_MDISPS | CD_MASK_PREVIEW_MCOL | CD_MASK_CLOTH_ORCO | CD_MASK_RECAST |
1244     /* BMESH ONLY START */
1245     CD_MASK_MPOLY | CD_MASK_MLOOP | CD_MASK_SHAPE_KEYINDEX | CD_MASK_SHAPEKEY | CD_MASK_BWEIGHT | CD_MASK_CREASE |
1246     CD_MASK_ORIGSPACE_MLOOP | CD_MASK_PREVIEW_MLOOPCOL | CD_MASK_BM_ELEM_PYPTR |
1247     /* BMESH ONLY END */
1248     CD_MASK_PAINT_MASK | CD_MASK_GRID_PAINT_MASK | CD_MASK_MVERT_SKIN |
1249     CD_MASK_FREESTYLE_EDGE | CD_MASK_FREESTYLE_FACE |
1250     CD_MASK_MLOOPTANGENT | CD_MASK_TESSLOOPNORMAL;
1251
1252 static const LayerTypeInfo *layerType_getInfo(int type)
1253 {
1254         if (type < 0 || type >= CD_NUMTYPES) return NULL;
1255
1256         return &LAYERTYPEINFO[type];
1257 }
1258
1259 static const char *layerType_getName(int type)
1260 {
1261         if (type < 0 || type >= CD_NUMTYPES) return NULL;
1262
1263         return LAYERTYPENAMES[type];
1264 }
1265
1266 void customData_mask_layers__print(CustomDataMask mask)
1267 {
1268         int i;
1269
1270         printf("mask=0x%lx:\n", (long unsigned int)mask);
1271         for (i = 0; i < CD_NUMTYPES; i++) {
1272                 if (mask & CD_TYPE_AS_MASK(i)) {
1273                         printf("  %s\n", layerType_getName(i));
1274                 }
1275         }
1276 }
1277
1278 /********************* CustomData functions *********************/
1279 static void customData_update_offsets(CustomData *data);
1280
1281 static CustomDataLayer *customData_add_layer__internal(CustomData *data, int type, int alloctype, void *layerdata,
1282                                                        int totelem, const char *name);
1283
1284 void CustomData_update_typemap(CustomData *data)
1285 {
1286         int i, lasttype = -1;
1287
1288         for (i = 0; i < CD_NUMTYPES; i++) {
1289                 data->typemap[i] = -1;
1290         }
1291
1292         for (i = 0; i < data->totlayer; i++) {
1293                 if (data->layers[i].type != lasttype) {
1294                         data->typemap[data->layers[i].type] = i;
1295                 }
1296                 lasttype = data->layers[i].type;
1297         }
1298 }
1299
1300 /* currently only used in BLI_assert */
1301 #ifndef NDEBUG
1302 static bool customdata_typemap_is_valid(const CustomData *data)
1303 {
1304         CustomData data_copy = *data;
1305         CustomData_update_typemap(&data_copy);
1306         return (memcmp(data->typemap, data_copy.typemap, sizeof(data->typemap)) == 0);
1307 }
1308 #endif
1309
1310 bool CustomData_merge(const struct CustomData *source, struct CustomData *dest,
1311                       CustomDataMask mask, int alloctype, int totelem)
1312 {
1313         /*const LayerTypeInfo *typeInfo;*/
1314         CustomDataLayer *layer, *newlayer;
1315         void *data;
1316         int i, type, number = 0, lasttype = -1, lastactive = 0, lastrender = 0, lastclone = 0, lastmask = 0, lastflag = 0;
1317         bool changed = false;
1318
1319         for (i = 0; i < source->totlayer; ++i) {
1320                 layer = &source->layers[i];
1321                 /*typeInfo = layerType_getInfo(layer->type);*/ /*UNUSED*/
1322
1323                 type = layer->type;
1324
1325                 if (type != lasttype) {
1326                         number = 0;
1327                         lastactive = layer->active;
1328                         lastrender = layer->active_rnd;
1329                         lastclone = layer->active_clone;
1330                         lastmask = layer->active_mask;
1331                         lasttype = type;
1332                         lastflag = layer->flag;
1333                 }
1334                 else
1335                         number++;
1336
1337                 if (lastflag & CD_FLAG_NOCOPY) continue;
1338                 else if (!(mask & CD_TYPE_AS_MASK(type))) continue;
1339                 else if (CustomData_get_layer_named(dest, type, layer->name)) continue;
1340
1341                 switch (alloctype) {
1342                         case CD_ASSIGN:
1343                         case CD_REFERENCE:
1344                         case CD_DUPLICATE:
1345                                 data = layer->data;
1346                                 break;
1347                         default:
1348                                 data = NULL;
1349                                 break;
1350                 }
1351
1352                 if ((alloctype == CD_ASSIGN) && (lastflag & CD_FLAG_NOFREE))
1353                         newlayer = customData_add_layer__internal(dest, type, CD_REFERENCE,
1354                                                                   data, totelem, layer->name);
1355                 else
1356                         newlayer = customData_add_layer__internal(dest, type, alloctype,
1357                                                                   data, totelem, layer->name);
1358                 
1359                 if (newlayer) {
1360                         newlayer->uid = layer->uid;
1361                         
1362                         newlayer->active = lastactive;
1363                         newlayer->active_rnd = lastrender;
1364                         newlayer->active_clone = lastclone;
1365                         newlayer->active_mask = lastmask;
1366                         newlayer->flag |= lastflag & (CD_FLAG_EXTERNAL | CD_FLAG_IN_MEMORY);
1367                         changed = true;
1368                 }
1369         }
1370
1371         CustomData_update_typemap(dest);
1372         return changed;
1373 }
1374
1375 void CustomData_copy(const struct CustomData *source, struct CustomData *dest,
1376                      CustomDataMask mask, int alloctype, int totelem)
1377 {
1378         CustomData_reset(dest);
1379
1380         if (source->external)
1381                 dest->external = MEM_dupallocN(source->external);
1382
1383         CustomData_merge(source, dest, mask, alloctype, totelem);
1384 }
1385
1386 static void customData_free_layer__internal(CustomDataLayer *layer, int totelem)
1387 {
1388         const LayerTypeInfo *typeInfo;
1389
1390         if (!(layer->flag & CD_FLAG_NOFREE) && layer->data) {
1391                 typeInfo = layerType_getInfo(layer->type);
1392
1393                 if (typeInfo->free)
1394                         typeInfo->free(layer->data, totelem, typeInfo->size);
1395
1396                 if (layer->data)
1397                         MEM_freeN(layer->data);
1398         }
1399 }
1400
1401 static void CustomData_external_free(CustomData *data)
1402 {
1403         if (data->external) {
1404                 MEM_freeN(data->external);
1405                 data->external = NULL;
1406         }
1407 }
1408
1409 void CustomData_reset(CustomData *data)
1410 {
1411         memset(data, 0, sizeof(*data));
1412         fill_vn_i(data->typemap, CD_NUMTYPES, -1);
1413 }
1414
1415 void CustomData_free(CustomData *data, int totelem)
1416 {
1417         int i;
1418
1419         for (i = 0; i < data->totlayer; ++i)
1420                 customData_free_layer__internal(&data->layers[i], totelem);
1421
1422         if (data->layers)
1423                 MEM_freeN(data->layers);
1424         
1425         CustomData_external_free(data);
1426         CustomData_reset(data);
1427 }
1428
1429 static void customData_update_offsets(CustomData *data)
1430 {
1431         const LayerTypeInfo *typeInfo;
1432         int i, offset = 0;
1433
1434         for (i = 0; i < data->totlayer; ++i) {
1435                 typeInfo = layerType_getInfo(data->layers[i].type);
1436
1437                 data->layers[i].offset = offset;
1438                 offset += typeInfo->size;
1439         }
1440
1441         data->totsize = offset;
1442         CustomData_update_typemap(data);
1443 }
1444
1445 /* to use when we're in the middle of modifying layers */
1446 static int CustomData_get_layer_index__notypemap(const CustomData *data, int type)
1447 {
1448         int i;
1449
1450         for (i = 0; i < data->totlayer; ++i)
1451                 if (data->layers[i].type == type)
1452                         return i;
1453
1454         return -1;
1455 }
1456
1457 /* -------------------------------------------------------------------- */
1458 /* index values to access the layers (offset from the layer start) */
1459
1460 int CustomData_get_layer_index(const CustomData *data, int type)
1461 {
1462         BLI_assert(customdata_typemap_is_valid(data));
1463         return data->typemap[type];
1464 }
1465
1466 int CustomData_get_layer_index_n(const struct CustomData *data, int type, int n)
1467 {
1468         int i = CustomData_get_layer_index(data, type);
1469
1470         if (i != -1) {
1471                 BLI_assert(i + n < data->totlayer);
1472                 i = (data->layers[i + n].type == type) ? (i + n) : (-1);
1473         }
1474
1475         return i;
1476 }
1477
1478 int CustomData_get_named_layer_index(const CustomData *data, int type, const char *name)
1479 {
1480         int i;
1481
1482         for (i = 0; i < data->totlayer; ++i)
1483                 if (data->layers[i].type == type)
1484                         if (strcmp(data->layers[i].name, name) == 0)
1485                                 return i;
1486
1487         return -1;
1488 }
1489
1490 int CustomData_get_active_layer_index(const CustomData *data, int type)
1491 {
1492         const int layer_index = data->typemap[type];
1493         BLI_assert(customdata_typemap_is_valid(data));
1494         return (layer_index != -1) ? layer_index + data->layers[layer_index].active : -1;
1495 }
1496
1497 int CustomData_get_render_layer_index(const CustomData *data, int type)
1498 {
1499         const int layer_index = data->typemap[type];
1500         BLI_assert(customdata_typemap_is_valid(data));
1501         return (layer_index != -1) ? layer_index + data->layers[layer_index].active_rnd : -1;
1502 }
1503
1504 int CustomData_get_clone_layer_index(const CustomData *data, int type)
1505 {
1506         const int layer_index = data->typemap[type];
1507         BLI_assert(customdata_typemap_is_valid(data));
1508         return (layer_index != -1) ? layer_index + data->layers[layer_index].active_clone : -1;
1509 }
1510
1511 int CustomData_get_stencil_layer_index(const CustomData *data, int type)
1512 {
1513         const int layer_index = data->typemap[type];
1514         BLI_assert(customdata_typemap_is_valid(data));
1515         return (layer_index != -1) ? layer_index + data->layers[layer_index].active_mask : -1;
1516 }
1517
1518
1519 /* -------------------------------------------------------------------- */
1520 /* index values per layer type */
1521
1522 int CustomData_get_named_layer(const struct CustomData *data, int type, const char *name)
1523 {
1524         const int named_index = CustomData_get_named_layer_index(data, type, name);
1525         const int layer_index = data->typemap[type];
1526         BLI_assert(customdata_typemap_is_valid(data));
1527         return (named_index != -1) ? named_index - layer_index : -1;
1528 }
1529
1530 int CustomData_get_active_layer(const CustomData *data, int type)
1531 {
1532         const int layer_index = data->typemap[type];
1533         BLI_assert(customdata_typemap_is_valid(data));
1534         return (layer_index != -1) ? data->layers[layer_index].active : -1;
1535 }
1536
1537 int CustomData_get_render_layer(const CustomData *data, int type)
1538 {
1539         const int layer_index = data->typemap[type];
1540         BLI_assert(customdata_typemap_is_valid(data));
1541         return (layer_index != -1) ? data->layers[layer_index].active_rnd : -1;
1542 }
1543
1544 int CustomData_get_clone_layer(const CustomData *data, int type)
1545 {
1546         const int layer_index = data->typemap[type];
1547         BLI_assert(customdata_typemap_is_valid(data));
1548         return (layer_index != -1) ? data->layers[layer_index].active_clone : -1;
1549 }
1550
1551 int CustomData_get_stencil_layer(const CustomData *data, int type)
1552 {
1553         const int layer_index = data->typemap[type];
1554         BLI_assert(customdata_typemap_is_valid(data));
1555         return (layer_index != -1) ? data->layers[layer_index].active_mask : -1;
1556 }
1557
1558 void CustomData_set_layer_active(CustomData *data, int type, int n)
1559 {
1560         int i;
1561
1562         for (i = 0; i < data->totlayer; ++i)
1563                 if (data->layers[i].type == type)
1564                         data->layers[i].active = n;
1565 }
1566
1567 void CustomData_set_layer_render(CustomData *data, int type, int n)
1568 {
1569         int i;
1570
1571         for (i = 0; i < data->totlayer; ++i)
1572                 if (data->layers[i].type == type)
1573                         data->layers[i].active_rnd = n;
1574 }
1575
1576 void CustomData_set_layer_clone(CustomData *data, int type, int n)
1577 {
1578         int i;
1579
1580         for (i = 0; i < data->totlayer; ++i)
1581                 if (data->layers[i].type == type)
1582                         data->layers[i].active_clone = n;
1583 }
1584
1585 void CustomData_set_layer_stencil(CustomData *data, int type, int n)
1586 {
1587         int i;
1588
1589         for (i = 0; i < data->totlayer; ++i)
1590                 if (data->layers[i].type == type)
1591                         data->layers[i].active_mask = n;
1592 }
1593
1594 /* for using with an index from CustomData_get_active_layer_index and CustomData_get_render_layer_index */
1595 void CustomData_set_layer_active_index(CustomData *data, int type, int n)
1596 {
1597         int i;
1598
1599         for (i = 0; i < data->totlayer; ++i)
1600                 if (data->layers[i].type == type)
1601                         data->layers[i].active = n - i;
1602 }
1603
1604 void CustomData_set_layer_render_index(CustomData *data, int type, int n)
1605 {
1606         int i;
1607
1608         for (i = 0; i < data->totlayer; ++i)
1609                 if (data->layers[i].type == type)
1610                         data->layers[i].active_rnd = n - i;
1611 }
1612
1613 void CustomData_set_layer_clone_index(CustomData *data, int type, int n)
1614 {
1615         int i;
1616
1617         for (i = 0; i < data->totlayer; ++i)
1618                 if (data->layers[i].type == type)
1619                         data->layers[i].active_clone = n - i;
1620 }
1621
1622 void CustomData_set_layer_stencil_index(CustomData *data, int type, int n)
1623 {
1624         int i;
1625
1626         for (i = 0; i < data->totlayer; ++i)
1627                 if (data->layers[i].type == type)
1628                         data->layers[i].active_mask = n - i;
1629 }
1630
1631 void CustomData_set_layer_flag(struct CustomData *data, int type, int flag)
1632 {
1633         int i;
1634
1635         for (i = 0; i < data->totlayer; ++i)
1636                 if (data->layers[i].type == type)
1637                         data->layers[i].flag |= flag;
1638 }
1639
1640 static int customData_resize(CustomData *data, int amount)
1641 {
1642         CustomDataLayer *tmp = MEM_callocN(sizeof(*tmp) * (data->maxlayer + amount),
1643                                            "CustomData->layers");
1644         if (!tmp) return 0;
1645
1646         data->maxlayer += amount;
1647         if (data->layers) {
1648                 memcpy(tmp, data->layers, sizeof(*tmp) * data->totlayer);
1649                 MEM_freeN(data->layers);
1650         }
1651         data->layers = tmp;
1652
1653         return 1;
1654 }
1655
1656 static CustomDataLayer *customData_add_layer__internal(CustomData *data, int type, int alloctype, void *layerdata,
1657                                                        int totelem, const char *name)
1658 {
1659         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
1660         int size = typeInfo->size * totelem, flag = 0, index = data->totlayer;
1661         void *newlayerdata = NULL;
1662
1663         /* Passing a layerdata to copy from with an alloctype that won't copy is
1664          * most likely a bug */
1665         BLI_assert(!layerdata ||
1666                    (alloctype == CD_ASSIGN) ||
1667                    (alloctype == CD_DUPLICATE) ||
1668                    (alloctype == CD_REFERENCE));
1669
1670         BLI_assert(size >= 0);
1671
1672         if (!typeInfo->defaultname && CustomData_has_layer(data, type))
1673                 return &data->layers[CustomData_get_layer_index(data, type)];
1674
1675         if ((alloctype == CD_ASSIGN) || (alloctype == CD_REFERENCE)) {
1676                 newlayerdata = layerdata;
1677         }
1678         else if (size > 0) {
1679                 newlayerdata = MEM_callocN(size, layerType_getName(type));
1680                 if (!newlayerdata)
1681                         return NULL;
1682         }
1683
1684         if (alloctype == CD_DUPLICATE && layerdata) {
1685                 if (typeInfo->copy)
1686                         typeInfo->copy(layerdata, newlayerdata, totelem);
1687                 else
1688                         memcpy(newlayerdata, layerdata, size);
1689         }
1690         else if (alloctype == CD_DEFAULT) {
1691                 if (typeInfo->set_default)
1692                         typeInfo->set_default((char *)newlayerdata, totelem);
1693         }
1694         else if (alloctype == CD_REFERENCE)
1695                 flag |= CD_FLAG_NOFREE;
1696
1697         if (index >= data->maxlayer) {
1698                 if (!customData_resize(data, CUSTOMDATA_GROW)) {
1699                         if (newlayerdata != layerdata)
1700                                 MEM_freeN(newlayerdata);
1701                         return NULL;
1702                 }
1703         }
1704         
1705         data->totlayer++;
1706
1707         /* keep layers ordered by type */
1708         for (; index > 0 && data->layers[index - 1].type > type; --index)
1709                 data->layers[index] = data->layers[index - 1];
1710
1711         data->layers[index].type = type;
1712         data->layers[index].flag = flag;
1713         data->layers[index].data = newlayerdata;
1714
1715         if (name || (name = DATA_(typeInfo->defaultname))) {
1716                 BLI_strncpy(data->layers[index].name, name, sizeof(data->layers[index].name));
1717                 CustomData_set_layer_unique_name(data, index);
1718         }
1719         else
1720                 data->layers[index].name[0] = '\0';
1721
1722         if (index > 0 && data->layers[index - 1].type == type) {
1723                 data->layers[index].active = data->layers[index - 1].active;
1724                 data->layers[index].active_rnd = data->layers[index - 1].active_rnd;
1725                 data->layers[index].active_clone = data->layers[index - 1].active_clone;
1726                 data->layers[index].active_mask = data->layers[index - 1].active_mask;
1727         }
1728         else {
1729                 data->layers[index].active = 0;
1730                 data->layers[index].active_rnd = 0;
1731                 data->layers[index].active_clone = 0;
1732                 data->layers[index].active_mask = 0;
1733         }
1734         
1735         customData_update_offsets(data);
1736
1737         return &data->layers[index];
1738 }
1739
1740 void *CustomData_add_layer(CustomData *data, int type, int alloctype,
1741                            void *layerdata, int totelem)
1742 {
1743         CustomDataLayer *layer;
1744         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
1745         
1746         layer = customData_add_layer__internal(data, type, alloctype, layerdata,
1747                                                totelem, typeInfo->defaultname);
1748         CustomData_update_typemap(data);
1749
1750         if (layer)
1751                 return layer->data;
1752
1753         return NULL;
1754 }
1755
1756 /*same as above but accepts a name*/
1757 void *CustomData_add_layer_named(CustomData *data, int type, int alloctype,
1758                                  void *layerdata, int totelem, const char *name)
1759 {
1760         CustomDataLayer *layer;
1761         
1762         layer = customData_add_layer__internal(data, type, alloctype, layerdata,
1763                                                totelem, name);
1764         CustomData_update_typemap(data);
1765
1766         if (layer)
1767                 return layer->data;
1768
1769         return NULL;
1770 }
1771
1772
1773 bool CustomData_free_layer(CustomData *data, int type, int totelem, int index)
1774 {
1775         const int n = index - CustomData_get_layer_index(data, type);
1776         int i;
1777         
1778         if (index < 0) return 0;
1779
1780         customData_free_layer__internal(&data->layers[index], totelem);
1781
1782         for (i = index + 1; i < data->totlayer; ++i)
1783                 data->layers[i - 1] = data->layers[i];
1784
1785         data->totlayer--;
1786
1787         /* if layer was last of type in array, set new active layer */
1788         i = CustomData_get_layer_index__notypemap(data, type);
1789
1790         if (i != -1) {
1791                 /* don't decrement zero index */
1792                 const int index_nonzero = n ? n : 1;
1793                 CustomDataLayer *layer;
1794
1795                 for (layer = &data->layers[i]; i < data->totlayer && layer->type == type; i++, layer++) {
1796                         if (layer->active       >= index_nonzero) layer->active--;
1797                         if (layer->active_rnd   >= index_nonzero) layer->active_rnd--;
1798                         if (layer->active_clone >= index_nonzero) layer->active_clone--;
1799                         if (layer->active_mask  >= index_nonzero) layer->active_mask--;
1800                 }
1801         }
1802
1803         if (data->totlayer <= data->maxlayer - CUSTOMDATA_GROW)
1804                 customData_resize(data, -CUSTOMDATA_GROW);
1805
1806         customData_update_offsets(data);
1807
1808         return 1;
1809 }
1810
1811 bool CustomData_free_layer_active(CustomData *data, int type, int totelem)
1812 {
1813         int index = 0;
1814         index = CustomData_get_active_layer_index(data, type);
1815         if (index == -1) return 0;
1816         return CustomData_free_layer(data, type, totelem, index);
1817 }
1818
1819
1820 void CustomData_free_layers(CustomData *data, int type, int totelem)
1821 {
1822         while (CustomData_has_layer(data, type))
1823                 CustomData_free_layer_active(data, type, totelem);
1824 }
1825
1826 bool CustomData_has_layer(const CustomData *data, int type)
1827 {
1828         return (CustomData_get_layer_index(data, type) != -1);
1829 }
1830
1831 int CustomData_number_of_layers(const CustomData *data, int type)
1832 {
1833         int i, number = 0;
1834
1835         for (i = 0; i < data->totlayer; i++)
1836                 if (data->layers[i].type == type)
1837                         number++;
1838         
1839         return number;
1840 }
1841
1842 int CustomData_number_of_layers_typemask(const CustomData *data, CustomDataMask mask)
1843 {
1844         int i, number = 0;
1845
1846         for (i = 0; i < data->totlayer; i++)
1847                 if (mask & CD_TYPE_AS_MASK(data->layers[i].type))
1848                         number++;
1849
1850         return number;
1851 }
1852
1853 void *CustomData_duplicate_referenced_layer(struct CustomData *data, const int type, const int totelem)
1854 {
1855         CustomDataLayer *layer;
1856         int layer_index;
1857
1858         /* get the layer index of the first layer of type */
1859         layer_index = CustomData_get_active_layer_index(data, type);
1860         if (layer_index == -1) return NULL;
1861
1862         layer = &data->layers[layer_index];
1863
1864         if (layer->flag & CD_FLAG_NOFREE) {
1865                 /* MEM_dupallocN won't work in case of complex layers, like e.g.
1866                  * CD_MDEFORMVERT, which has pointers to allocated data...
1867                  * So in case a custom copy function is defined, use it!
1868                  */
1869                 const LayerTypeInfo *typeInfo = layerType_getInfo(layer->type);
1870
1871                 if (typeInfo->copy) {
1872                         char *dest_data = MEM_mallocN(typeInfo->size * totelem, "CD duplicate ref layer");
1873                         typeInfo->copy(layer->data, dest_data, totelem);
1874                         layer->data = dest_data;
1875                 }
1876                 else
1877                         layer->data = MEM_dupallocN(layer->data);
1878
1879                 layer->flag &= ~CD_FLAG_NOFREE;
1880         }
1881
1882         return layer->data;
1883 }
1884
1885 void *CustomData_duplicate_referenced_layer_named(struct CustomData *data,
1886                                                   const int type, const char *name, const int totelem)
1887 {
1888         CustomDataLayer *layer;
1889         int layer_index;
1890
1891         /* get the layer index of the desired layer */
1892         layer_index = CustomData_get_named_layer_index(data, type, name);
1893         if (layer_index == -1) return NULL;
1894
1895         layer = &data->layers[layer_index];
1896
1897         if (layer->flag & CD_FLAG_NOFREE) {
1898                 /* MEM_dupallocN won't work in case of complex layers, like e.g.
1899                  * CD_MDEFORMVERT, which has pointers to allocated data...
1900                  * So in case a custom copy function is defined, use it!
1901                  */
1902                 const LayerTypeInfo *typeInfo = layerType_getInfo(layer->type);
1903
1904                 if (typeInfo->copy) {
1905                         char *dest_data = MEM_mallocN(typeInfo->size * totelem, "CD duplicate ref layer");
1906                         typeInfo->copy(layer->data, dest_data, totelem);
1907                         layer->data = dest_data;
1908                 }
1909                 else
1910                         layer->data = MEM_dupallocN(layer->data);
1911
1912                 layer->flag &= ~CD_FLAG_NOFREE;
1913         }
1914
1915         return layer->data;
1916 }
1917
1918 bool CustomData_is_referenced_layer(struct CustomData *data, int type)
1919 {
1920         CustomDataLayer *layer;
1921         int layer_index;
1922
1923         /* get the layer index of the first layer of type */
1924         layer_index = CustomData_get_active_layer_index(data, type);
1925         if (layer_index == -1) return 0;
1926
1927         layer = &data->layers[layer_index];
1928
1929         return (layer->flag & CD_FLAG_NOFREE) != 0;
1930 }
1931
1932 void CustomData_free_temporary(CustomData *data, int totelem)
1933 {
1934         CustomDataLayer *layer;
1935         int i, j;
1936         bool changed = false;
1937
1938         for (i = 0, j = 0; i < data->totlayer; ++i) {
1939                 layer = &data->layers[i];
1940
1941                 if (i != j)
1942                         data->layers[j] = data->layers[i];
1943
1944                 if ((layer->flag & CD_FLAG_TEMPORARY) == CD_FLAG_TEMPORARY) {
1945                         customData_free_layer__internal(layer, totelem);
1946                         changed = true;
1947                 }
1948                 else
1949                         j++;
1950         }
1951
1952         data->totlayer = j;
1953
1954         if (data->totlayer <= data->maxlayer - CUSTOMDATA_GROW) {
1955                 customData_resize(data, -CUSTOMDATA_GROW);
1956                 changed = true;
1957         }
1958
1959         if (changed) {
1960                 customData_update_offsets(data);
1961         }
1962 }
1963
1964 void CustomData_set_only_copy(const struct CustomData *data,
1965                               CustomDataMask mask)
1966 {
1967         int i;
1968
1969         for (i = 0; i < data->totlayer; ++i)
1970                 if (!(mask & CD_TYPE_AS_MASK(data->layers[i].type)))
1971                         data->layers[i].flag |= CD_FLAG_NOCOPY;
1972 }
1973
1974 void CustomData_copy_elements(int type, void *source, void *dest, int count)
1975 {
1976         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
1977
1978         if (typeInfo->copy)
1979                 typeInfo->copy(source, dest, count);
1980         else
1981                 memcpy(dest, source, typeInfo->size * count);
1982 }
1983
1984 static void CustomData_copy_data_layer(const CustomData *source, CustomData *dest,
1985                                        int src_i, int dest_i,
1986                                        int source_index, int dest_index, int count) {
1987         const LayerTypeInfo *typeInfo;
1988         int src_offset;
1989         int dest_offset;
1990
1991         char *src_data = source->layers[src_i].data;
1992         char *dest_data = dest->layers[dest_i].data;
1993
1994         typeInfo = layerType_getInfo(source->layers[src_i].type);
1995
1996         src_offset = source_index * typeInfo->size;
1997         dest_offset = dest_index * typeInfo->size;
1998
1999         if (!src_data || !dest_data) {
2000                 if (!(src_data == NULL && dest_data == NULL)) {
2001                         printf("%s: warning null data for %s type (%p --> %p), skipping\n",
2002                                    __func__, layerType_getName(source->layers[src_i].type),
2003                                    (void *)src_data, (void *)dest_data);
2004                 }
2005                 return;
2006         }
2007
2008         if (typeInfo->copy)
2009                 typeInfo->copy(src_data + src_offset,
2010                                dest_data + dest_offset,
2011                                count);
2012         else
2013                 memcpy(dest_data + dest_offset,
2014                        src_data + src_offset,
2015                        count * typeInfo->size);
2016 }
2017
2018 void CustomData_copy_data_named(const CustomData *source, CustomData *dest,
2019                                 int source_index, int dest_index, int count)
2020 {
2021         int src_i, dest_i;
2022
2023         /* copies a layer at a time */
2024         for (src_i = 0; src_i < source->totlayer; ++src_i) {
2025
2026                 dest_i = CustomData_get_named_layer_index(dest, source->layers[src_i].type, source->layers[src_i].name);
2027
2028                 /* if we found a matching layer, copy the data */
2029                 if (dest_i != -1) {
2030                         CustomData_copy_data_layer(source, dest, src_i, dest_i, source_index, dest_index, count);
2031                 }
2032         }
2033 }
2034
2035 void CustomData_copy_data(const CustomData *source, CustomData *dest,
2036                           int source_index, int dest_index, int count)
2037 {
2038         int src_i, dest_i;
2039
2040         /* copies a layer at a time */
2041         dest_i = 0;
2042         for (src_i = 0; src_i < source->totlayer; ++src_i) {
2043
2044                 /* find the first dest layer with type >= the source type
2045                  * (this should work because layers are ordered by type)
2046                  */
2047                 while (dest_i < dest->totlayer && dest->layers[dest_i].type < source->layers[src_i].type) {
2048                         dest_i++;
2049                 }
2050
2051                 /* if there are no more dest layers, we're done */
2052                 if (dest_i >= dest->totlayer) return;
2053
2054                 /* if we found a matching layer, copy the data */
2055                 if (dest->layers[dest_i].type == source->layers[src_i].type) {
2056                         CustomData_copy_data_layer(source, dest, src_i, dest_i, source_index, dest_index, count);
2057                         
2058                         /* if there are multiple source & dest layers of the same type,
2059                          * we don't want to copy all source layers to the same dest, so
2060                          * increment dest_i
2061                          */
2062                         dest_i++;
2063                 }
2064         }
2065 }
2066
2067 void CustomData_free_elem(CustomData *data, int index, int count)
2068 {
2069         int i;
2070         const LayerTypeInfo *typeInfo;
2071
2072         for (i = 0; i < data->totlayer; ++i) {
2073                 if (!(data->layers[i].flag & CD_FLAG_NOFREE)) {
2074                         typeInfo = layerType_getInfo(data->layers[i].type);
2075
2076                         if (typeInfo->free) {
2077                                 int offset = typeInfo->size * index;
2078
2079                                 typeInfo->free((char *)data->layers[i].data + offset,
2080                                                count, typeInfo->size);
2081                         }
2082                 }
2083         }
2084 }
2085
2086 #define SOURCE_BUF_SIZE 100
2087
2088 void CustomData_interp(const CustomData *source, CustomData *dest,
2089                        int *src_indices, float *weights, float *sub_weights,
2090                        int count, int dest_index)
2091 {
2092         int src_i, dest_i;
2093         int dest_offset;
2094         int j;
2095         void *source_buf[SOURCE_BUF_SIZE];
2096         void **sources = source_buf;
2097
2098         /* slow fallback in case we're interpolating a ridiculous number of
2099          * elements
2100          */
2101         if (count > SOURCE_BUF_SIZE)
2102                 sources = MEM_callocN(sizeof(*sources) * count,
2103                                       "CustomData_interp sources");
2104
2105         /* interpolates a layer at a time */
2106         dest_i = 0;
2107         for (src_i = 0; src_i < source->totlayer; ++src_i) {
2108                 const LayerTypeInfo *typeInfo = layerType_getInfo(source->layers[src_i].type);
2109                 if (!typeInfo->interp) continue;
2110
2111                 /* find the first dest layer with type >= the source type
2112                  * (this should work because layers are ordered by type)
2113                  */
2114                 while (dest_i < dest->totlayer && dest->layers[dest_i].type < source->layers[src_i].type) {
2115                         dest_i++;
2116                 }
2117
2118                 /* if there are no more dest layers, we're done */
2119                 if (dest_i >= dest->totlayer) break;
2120
2121                 /* if we found a matching layer, copy the data */
2122                 if (dest->layers[dest_i].type == source->layers[src_i].type) {
2123                         void *src_data = source->layers[src_i].data;
2124
2125                         for (j = 0; j < count; ++j) {
2126                                 sources[j] = (char *)src_data + typeInfo->size * src_indices[j];
2127                         }
2128
2129                         dest_offset = dest_index * typeInfo->size;
2130
2131                         typeInfo->interp(sources, weights, sub_weights, count,
2132                                          (char *)dest->layers[dest_i].data + dest_offset);
2133
2134                         /* if there are multiple source & dest layers of the same type,
2135                          * we don't want to copy all source layers to the same dest, so
2136                          * increment dest_i
2137                          */
2138                         dest_i++;
2139                 }
2140         }
2141
2142         if (count > SOURCE_BUF_SIZE) MEM_freeN(sources);
2143 }
2144
2145 void CustomData_swap(struct CustomData *data, int index, const int *corner_indices)
2146 {
2147         const LayerTypeInfo *typeInfo;
2148         int i;
2149
2150         for (i = 0; i < data->totlayer; ++i) {
2151                 typeInfo = layerType_getInfo(data->layers[i].type);
2152
2153                 if (typeInfo->swap) {
2154                         int offset = typeInfo->size * index;
2155
2156                         typeInfo->swap((char *)data->layers[i].data + offset, corner_indices);
2157                 }
2158         }
2159 }
2160
2161 void *CustomData_get(const CustomData *data, int index, int type)
2162 {
2163         int offset;
2164         int layer_index;
2165         
2166         BLI_assert(index >= 0);
2167
2168         /* get the layer index of the active layer of type */
2169         layer_index = CustomData_get_active_layer_index(data, type);
2170         if (layer_index == -1) return NULL;
2171
2172         /* get the offset of the desired element */
2173         offset = layerType_getInfo(type)->size * index;
2174
2175         return (char *)data->layers[layer_index].data + offset;
2176 }
2177
2178 void *CustomData_get_n(const CustomData *data, int type, int index, int n)
2179 {
2180         int layer_index;
2181         int offset;
2182
2183         BLI_assert(index >= 0 && n >= 0);
2184
2185         /* get the layer index of the first layer of type */
2186         layer_index = data->typemap[type];
2187         if (layer_index == -1) return NULL;
2188
2189         offset = layerType_getInfo(type)->size * index;
2190         return (char *)data->layers[layer_index + n].data + offset;
2191 }
2192
2193 void *CustomData_get_layer(const CustomData *data, int type)
2194 {
2195         /* get the layer index of the active layer of type */
2196         int layer_index = CustomData_get_active_layer_index(data, type);
2197         if (layer_index == -1) return NULL;
2198
2199         return data->layers[layer_index].data;
2200 }
2201
2202 void *CustomData_get_layer_n(const CustomData *data, int type, int n)
2203 {
2204         /* get the layer index of the active layer of type */
2205         int layer_index = CustomData_get_layer_index_n(data, type, n);
2206         if (layer_index == -1) return NULL;
2207
2208         return data->layers[layer_index].data;
2209 }
2210
2211 void *CustomData_get_layer_named(const struct CustomData *data, int type,
2212                                  const char *name)
2213 {
2214         int layer_index = CustomData_get_named_layer_index(data, type, name);
2215         if (layer_index == -1) return NULL;
2216
2217         return data->layers[layer_index].data;
2218 }
2219
2220 int CustomData_get_offset(const CustomData *data, int type)
2221 {
2222         /* get the layer index of the active layer of type */
2223         int layer_index = CustomData_get_active_layer_index(data, type);
2224         if (layer_index == -1) return -1;
2225
2226         return data->layers[layer_index].offset;
2227 }
2228
2229 int CustomData_get_n_offset(const CustomData *data, int type, int n)
2230 {
2231         /* get the layer index of the active layer of type */
2232         int layer_index = CustomData_get_layer_index_n(data, type, n);
2233         if (layer_index == -1) return -1;
2234
2235         return data->layers[layer_index].offset;
2236 }
2237
2238 bool CustomData_set_layer_name(const CustomData *data, int type, int n, const char *name)
2239 {
2240         /* get the layer index of the first layer of type */
2241         int layer_index = CustomData_get_layer_index_n(data, type, n);
2242
2243         if (layer_index == -1) return false;
2244         if (!name) return false;
2245         
2246         BLI_strncpy(data->layers[layer_index].name, name, sizeof(data->layers[layer_index].name));
2247         
2248         return true;
2249 }
2250
2251 void *CustomData_set_layer(const CustomData *data, int type, void *ptr)
2252 {
2253         /* get the layer index of the first layer of type */
2254         int layer_index = CustomData_get_active_layer_index(data, type);
2255
2256         if (layer_index == -1) return NULL;
2257
2258         data->layers[layer_index].data = ptr;
2259
2260         return ptr;
2261 }
2262
2263 void *CustomData_set_layer_n(const struct CustomData *data, int type, int n, void *ptr)
2264 {
2265         /* get the layer index of the first layer of type */
2266         int layer_index = CustomData_get_layer_index_n(data, type, n);
2267         if (layer_index == -1) return NULL;
2268
2269         data->layers[layer_index].data = ptr;
2270
2271         return ptr;
2272 }
2273
2274 void CustomData_set(const CustomData *data, int index, int type, void *source)
2275 {
2276         void *dest = CustomData_get(data, index, type);
2277         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2278
2279         if (!dest) return;
2280
2281         if (typeInfo->copy)
2282                 typeInfo->copy(source, dest, 1);
2283         else
2284                 memcpy(dest, source, typeInfo->size);
2285 }
2286
2287 /* BMesh functions */
2288 /* needed to convert to/from different face reps */
2289 void CustomData_to_bmeshpoly(CustomData *fdata, CustomData *pdata, CustomData *ldata,
2290                              int totloop, int totpoly)
2291 {
2292         int i;
2293         for (i = 0; i < fdata->totlayer; i++) {
2294                 if (fdata->layers[i].type == CD_MTFACE) {
2295                         CustomData_add_layer_named(pdata, CD_MTEXPOLY, CD_CALLOC, NULL, totpoly, fdata->layers[i].name);
2296                         CustomData_add_layer_named(ldata, CD_MLOOPUV, CD_CALLOC, NULL, totloop, fdata->layers[i].name);
2297                 }
2298                 else if (fdata->layers[i].type == CD_MCOL) {
2299                         CustomData_add_layer_named(ldata, CD_MLOOPCOL, CD_CALLOC, NULL, totloop, fdata->layers[i].name);
2300                 }
2301                 else if (fdata->layers[i].type == CD_MDISPS) {
2302                         CustomData_add_layer_named(ldata, CD_MDISPS, CD_CALLOC, NULL, totloop, fdata->layers[i].name);
2303                 }
2304                 else if (fdata->layers[i].type == CD_TESSLOOPNORMAL) {
2305                         CustomData_add_layer_named(ldata, CD_NORMAL, CD_CALLOC, NULL, totloop, fdata->layers[i].name);
2306                 }
2307         }
2308 }
2309
2310 void CustomData_from_bmeshpoly(CustomData *fdata, CustomData *pdata, CustomData *ldata, int total)
2311 {
2312         int i;
2313         for (i = 0; i < pdata->totlayer; i++) {
2314                 if (pdata->layers[i].type == CD_MTEXPOLY) {
2315                         CustomData_add_layer_named(fdata, CD_MTFACE, CD_CALLOC, NULL, total, pdata->layers[i].name);
2316                 }
2317         }
2318         for (i = 0; i < ldata->totlayer; i++) {
2319                 if (ldata->layers[i].type == CD_MLOOPCOL) {
2320                         CustomData_add_layer_named(fdata, CD_MCOL, CD_CALLOC, NULL, total, ldata->layers[i].name);
2321                 }
2322                 else if (ldata->layers[i].type == CD_PREVIEW_MLOOPCOL) {
2323                         CustomData_add_layer_named(fdata, CD_PREVIEW_MCOL, CD_CALLOC, NULL, total, ldata->layers[i].name);
2324                 }
2325                 else if (ldata->layers[i].type == CD_ORIGSPACE_MLOOP) {
2326                         CustomData_add_layer_named(fdata, CD_ORIGSPACE, CD_CALLOC, NULL, total, ldata->layers[i].name);
2327                 }
2328                 else if (ldata->layers[i].type == CD_NORMAL) {
2329                         CustomData_add_layer_named(fdata, CD_TESSLOOPNORMAL, CD_CALLOC, NULL, total, ldata->layers[i].name);
2330                 }
2331         }
2332
2333         CustomData_bmesh_update_active_layers(fdata, pdata, ldata);
2334 }
2335
2336 void CustomData_bmesh_update_active_layers(CustomData *fdata, CustomData *pdata, CustomData *ldata)
2337 {
2338         int act;
2339
2340         if (CustomData_has_layer(pdata, CD_MTEXPOLY)) {
2341                 act = CustomData_get_active_layer(pdata, CD_MTEXPOLY);
2342                 CustomData_set_layer_active(ldata, CD_MLOOPUV, act);
2343                 CustomData_set_layer_active(fdata, CD_MTFACE, act);
2344
2345                 act = CustomData_get_render_layer(pdata, CD_MTEXPOLY);
2346                 CustomData_set_layer_render(ldata, CD_MLOOPUV, act);
2347                 CustomData_set_layer_render(fdata, CD_MTFACE, act);
2348
2349                 act = CustomData_get_clone_layer(pdata, CD_MTEXPOLY);
2350                 CustomData_set_layer_clone(ldata, CD_MLOOPUV, act);
2351                 CustomData_set_layer_clone(fdata, CD_MTFACE, act);
2352
2353                 act = CustomData_get_stencil_layer(pdata, CD_MTEXPOLY);
2354                 CustomData_set_layer_stencil(ldata, CD_MLOOPUV, act);
2355                 CustomData_set_layer_stencil(fdata, CD_MTFACE, act);
2356         }
2357
2358         if (CustomData_has_layer(ldata, CD_MLOOPCOL)) {
2359                 act = CustomData_get_active_layer(ldata, CD_MLOOPCOL);
2360                 CustomData_set_layer_active(fdata, CD_MCOL, act);
2361
2362                 act = CustomData_get_render_layer(ldata, CD_MLOOPCOL);
2363                 CustomData_set_layer_render(fdata, CD_MCOL, act);
2364
2365                 act = CustomData_get_clone_layer(ldata, CD_MLOOPCOL);
2366                 CustomData_set_layer_clone(fdata, CD_MCOL, act);
2367
2368                 act = CustomData_get_stencil_layer(ldata, CD_MLOOPCOL);
2369                 CustomData_set_layer_stencil(fdata, CD_MCOL, act);
2370         }
2371 }
2372
2373 /* update active indices for active/render/clone/stencil custom data layers
2374  * based on indices from fdata layers
2375  * used by do_versions in readfile.c when creating pdata and ldata for pre-bmesh
2376  * meshes and needed to preserve active/render/clone/stencil flags set in pre-bmesh files
2377  */
2378 void CustomData_bmesh_do_versions_update_active_layers(CustomData *fdata, CustomData *pdata, CustomData *ldata)
2379 {
2380         int act;
2381
2382         if (CustomData_has_layer(fdata, CD_MTFACE)) {
2383                 act = CustomData_get_active_layer(fdata, CD_MTFACE);
2384                 CustomData_set_layer_active(pdata, CD_MTEXPOLY, act);
2385                 CustomData_set_layer_active(ldata, CD_MLOOPUV, act);
2386
2387                 act = CustomData_get_render_layer(fdata, CD_MTFACE);
2388                 CustomData_set_layer_render(pdata, CD_MTEXPOLY, act);
2389                 CustomData_set_layer_render(ldata, CD_MLOOPUV, act);
2390
2391                 act = CustomData_get_clone_layer(fdata, CD_MTFACE);
2392                 CustomData_set_layer_clone(pdata, CD_MTEXPOLY, act);
2393                 CustomData_set_layer_clone(ldata, CD_MLOOPUV, act);
2394
2395                 act = CustomData_get_stencil_layer(fdata, CD_MTFACE);
2396                 CustomData_set_layer_stencil(pdata, CD_MTEXPOLY, act);
2397                 CustomData_set_layer_stencil(ldata, CD_MLOOPUV, act);
2398         }
2399
2400         if (CustomData_has_layer(fdata, CD_MCOL)) {
2401                 act = CustomData_get_active_layer(fdata, CD_MCOL);
2402                 CustomData_set_layer_active(ldata, CD_MLOOPCOL, act);
2403
2404                 act = CustomData_get_render_layer(fdata, CD_MCOL);
2405                 CustomData_set_layer_render(ldata, CD_MLOOPCOL, act);
2406
2407                 act = CustomData_get_clone_layer(fdata, CD_MCOL);
2408                 CustomData_set_layer_clone(ldata, CD_MLOOPCOL, act);
2409
2410                 act = CustomData_get_stencil_layer(fdata, CD_MCOL);
2411                 CustomData_set_layer_stencil(ldata, CD_MLOOPCOL, act);
2412         }
2413 }
2414
2415 void CustomData_bmesh_init_pool(CustomData *data, int totelem, const char htype)
2416 {
2417         int chunksize;
2418
2419         /* Dispose old pools before calling here to avoid leaks */
2420         BLI_assert(data->pool == NULL);
2421
2422         switch (htype) {
2423                 case BM_VERT: chunksize = bm_mesh_chunksize_default.totvert;  break;
2424                 case BM_EDGE: chunksize = bm_mesh_chunksize_default.totedge;  break;
2425                 case BM_LOOP: chunksize = bm_mesh_chunksize_default.totloop;  break;
2426                 case BM_FACE: chunksize = bm_mesh_chunksize_default.totface;  break;
2427                 default:
2428                         BLI_assert(0);
2429                         chunksize = 512;
2430                         break;
2431         }
2432
2433         /* If there are no layers, no pool is needed just yet */
2434         if (data->totlayer) {
2435                 data->pool = BLI_mempool_create(data->totsize, totelem, chunksize, BLI_MEMPOOL_NOP);
2436         }
2437 }
2438
2439 bool CustomData_bmesh_merge(CustomData *source, CustomData *dest,
2440                             CustomDataMask mask, int alloctype, BMesh *bm, const char htype)
2441 {
2442         BMHeader *h;
2443         BMIter iter;
2444         CustomData destold;
2445         void *tmp;
2446         int iter_type;
2447         int totelem;
2448
2449         if (CustomData_number_of_layers_typemask(source, mask) == 0) {
2450                 return false;
2451         }
2452
2453         /* copy old layer description so that old data can be copied into
2454          * the new allocation */
2455         destold = *dest;
2456         if (destold.layers) {
2457                 destold.layers = MEM_dupallocN(destold.layers);
2458         }
2459
2460         if (CustomData_merge(source, dest, mask, alloctype, 0) == false) {
2461                 if (destold.layers)
2462                         MEM_freeN(destold.layers);
2463                 return false;
2464         }
2465
2466         switch (htype) {
2467                 case BM_VERT:
2468                         iter_type = BM_VERTS_OF_MESH;
2469                         totelem = bm->totvert;
2470                         break;
2471                 case BM_EDGE:
2472                         iter_type = BM_EDGES_OF_MESH;
2473                         totelem = bm->totedge;
2474                         break;
2475                 case BM_LOOP:
2476                         iter_type = BM_LOOPS_OF_FACE;
2477                         totelem = bm->totloop;
2478                         break;
2479                 case BM_FACE:
2480                         iter_type = BM_FACES_OF_MESH;
2481                         totelem = bm->totface;
2482                         break;
2483                 default: /* should never happen */
2484                         BLI_assert(!"invalid type given");
2485                         iter_type = BM_VERTS_OF_MESH;
2486                         totelem = bm->totvert;
2487                         break;
2488         }
2489
2490         dest->pool = NULL;
2491         CustomData_bmesh_init_pool(dest, totelem, htype);
2492
2493         if (iter_type != BM_LOOPS_OF_FACE) {
2494                 /*ensure all current elements follow new customdata layout*/
2495                 BM_ITER_MESH (h, &iter, bm, iter_type) {
2496                         tmp = NULL;
2497                         CustomData_bmesh_copy_data(&destold, dest, h->data, &tmp);
2498                         CustomData_bmesh_free_block(&destold, &h->data);
2499                         h->data = tmp;
2500                 }
2501         }
2502         else {
2503                 BMFace *f;
2504                 BMLoop *l;
2505                 BMIter liter;
2506
2507                 /*ensure all current elements follow new customdata layout*/
2508                 BM_ITER_MESH (f, &iter, bm, BM_FACES_OF_MESH) {
2509                         BM_ITER_ELEM (l, &liter, f, BM_LOOPS_OF_FACE) {
2510                                 tmp = NULL;
2511                                 CustomData_bmesh_copy_data(&destold, dest, l->head.data, &tmp);
2512                                 CustomData_bmesh_free_block(&destold, &l->head.data);
2513                                 l->head.data = tmp;
2514                         }
2515                 }
2516         }
2517
2518         if (destold.pool) BLI_mempool_destroy(destold.pool);
2519         if (destold.layers) MEM_freeN(destold.layers);
2520         return true;
2521 }
2522
2523 void CustomData_bmesh_free_block(CustomData *data, void **block)
2524 {
2525         const LayerTypeInfo *typeInfo;
2526         int i;
2527
2528         if (*block == NULL)
2529                 return;
2530
2531         for (i = 0; i < data->totlayer; ++i) {
2532                 if (!(data->layers[i].flag & CD_FLAG_NOFREE)) {
2533                         typeInfo = layerType_getInfo(data->layers[i].type);
2534
2535                         if (typeInfo->free) {
2536                                 int offset = data->layers[i].offset;
2537                                 typeInfo->free((char *)*block + offset, 1, typeInfo->size);
2538                         }
2539                 }
2540         }
2541
2542         if (data->totsize)
2543                 BLI_mempool_free(data->pool, *block);
2544
2545         *block = NULL;
2546 }
2547
2548 /**
2549  * Same as #CustomData_bmesh_free_block but zero the memory rather then freeing.
2550  */
2551 void CustomData_bmesh_free_block_data(CustomData *data, void **block)
2552 {
2553         const LayerTypeInfo *typeInfo;
2554         int i;
2555
2556         if (*block == NULL)
2557                 return;
2558
2559         for (i = 0; i < data->totlayer; ++i) {
2560                 if (!(data->layers[i].flag & CD_FLAG_NOFREE)) {
2561                         typeInfo = layerType_getInfo(data->layers[i].type);
2562
2563                         if (typeInfo->free) {
2564                                 int offset = data->layers[i].offset;
2565                                 typeInfo->free((char *)*block + offset, 1, typeInfo->size);
2566                         }
2567                 }
2568         }
2569
2570         if (data->totsize)
2571                 memset(*block, 0, data->totsize);
2572 }
2573
2574 static void CustomData_bmesh_alloc_block(CustomData *data, void **block)
2575 {
2576
2577         if (*block)
2578                 CustomData_bmesh_free_block(data, block);
2579
2580         if (data->totsize > 0)
2581                 *block = BLI_mempool_alloc(data->pool);
2582         else
2583                 *block = NULL;
2584 }
2585
2586 void CustomData_bmesh_copy_data(const CustomData *source, CustomData *dest,
2587                                 void *src_block, void **dest_block)
2588 {
2589         const LayerTypeInfo *typeInfo;
2590         int dest_i, src_i;
2591
2592         if (*dest_block == NULL) {
2593                 CustomData_bmesh_alloc_block(dest, dest_block);
2594                 if (*dest_block)
2595                         memset(*dest_block, 0, dest->totsize);
2596         }
2597         
2598         /* copies a layer at a time */
2599         dest_i = 0;
2600         for (src_i = 0; src_i < source->totlayer; ++src_i) {
2601
2602                 /* find the first dest layer with type >= the source type
2603                  * (this should work because layers are ordered by type)
2604                  */
2605                 while (dest_i < dest->totlayer && dest->layers[dest_i].type < source->layers[src_i].type) {
2606                         dest_i++;
2607                 }
2608
2609                 /* if there are no more dest layers, we're done */
2610                 if (dest_i >= dest->totlayer) return;
2611
2612                 /* if we found a matching layer, copy the data */
2613                 if (dest->layers[dest_i].type == source->layers[src_i].type &&
2614                     strcmp(dest->layers[dest_i].name, source->layers[src_i].name) == 0)
2615                 {
2616                         char *src_data = (char *)src_block + source->layers[src_i].offset;
2617                         char *dest_data = (char *)*dest_block + dest->layers[dest_i].offset;
2618
2619                         typeInfo = layerType_getInfo(source->layers[src_i].type);
2620
2621                         if (typeInfo->copy)
2622                                 typeInfo->copy(src_data, dest_data, 1);
2623                         else
2624                                 memcpy(dest_data, src_data, typeInfo->size);
2625
2626                         /* if there are multiple source & dest layers of the same type,
2627                          * we don't want to copy all source layers to the same dest, so
2628                          * increment dest_i
2629                          */
2630                         dest_i++;
2631                 }
2632         }
2633 }
2634
2635 /*Bmesh Custom Data Functions. Should replace editmesh ones with these as well, due to more effecient memory alloc*/
2636 void *CustomData_bmesh_get(const CustomData *data, void *block, int type)
2637 {
2638         int layer_index;
2639         
2640         /* get the layer index of the first layer of type */
2641         layer_index = CustomData_get_active_layer_index(data, type);
2642         if (layer_index == -1) return NULL;
2643
2644         return (char *)block + data->layers[layer_index].offset;
2645 }
2646
2647 void *CustomData_bmesh_get_n(const CustomData *data, void *block, int type, int n)
2648 {
2649         int layer_index;
2650         
2651         /* get the layer index of the first layer of type */
2652         layer_index = CustomData_get_layer_index(data, type);
2653         if (layer_index == -1) return NULL;
2654
2655         return (char *)block + data->layers[layer_index + n].offset;
2656 }
2657
2658 /*gets from the layer at physical index n, note: doesn't check type.*/
2659 void *CustomData_bmesh_get_layer_n(const CustomData *data, void *block, int n)
2660 {
2661         if (n < 0 || n >= data->totlayer) return NULL;
2662
2663         return (char *)block + data->layers[n].offset;
2664 }
2665
2666 bool CustomData_layer_has_math(struct CustomData *data, int layer_n)
2667 {
2668         const LayerTypeInfo *typeInfo = layerType_getInfo(data->layers[layer_n].type);
2669         
2670         if (typeInfo->equal && typeInfo->add && typeInfo->multiply && 
2671             typeInfo->initminmax && typeInfo->dominmax)
2672         {
2673                 return true;
2674         }
2675         
2676         return false;
2677 }
2678
2679 bool CustomData_layer_has_interp(struct CustomData *data, int layer_n)
2680 {
2681         const LayerTypeInfo *typeInfo = layerType_getInfo(data->layers[layer_n].type);
2682
2683         if (typeInfo->interp) {
2684                 return true;
2685         }
2686
2687         return false;
2688 }
2689
2690 bool CustomData_has_math(struct CustomData *data)
2691 {
2692         int i;
2693
2694         /* interpolates a layer at a time */
2695         for (i = 0; i < data->totlayer; ++i) {
2696                 if (CustomData_layer_has_math(data, i)) {
2697                         return true;
2698                 }
2699         }
2700
2701         return false;
2702 }
2703
2704 /* a non bmesh version would have to check layer->data */
2705 bool CustomData_bmesh_has_free(struct CustomData *data)
2706 {
2707         const LayerTypeInfo *typeInfo;
2708         int i;
2709
2710         for (i = 0; i < data->totlayer; ++i) {
2711                 if (!(data->layers[i].flag & CD_FLAG_NOFREE)) {
2712                         typeInfo = layerType_getInfo(data->layers[i].type);
2713                         if (typeInfo->free) {
2714                                 return true;
2715                         }
2716                 }
2717         }
2718         return false;
2719 }
2720
2721 bool CustomData_has_interp(struct CustomData *data)
2722 {
2723         int i;
2724
2725         /* interpolates a layer at a time */
2726         for (i = 0; i < data->totlayer; ++i) {
2727                 if (CustomData_layer_has_interp(data, i)) {
2728                         return true;
2729                 }
2730         }
2731
2732         return false;
2733 }
2734
2735 /* copies the "value" (e.g. mloopuv uv or mloopcol colors) from one block to
2736  * another, while not overwriting anything else (e.g. flags)*/
2737 void CustomData_data_copy_value(int type, void *source, void *dest)
2738 {
2739         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2740
2741         if (!dest) return;
2742
2743         if (typeInfo->copyvalue)
2744                 typeInfo->copyvalue(source, dest);
2745         else
2746                 memcpy(dest, source, typeInfo->size);
2747 }
2748
2749 bool CustomData_data_equals(int type, void *data1, void *data2)
2750 {
2751         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2752
2753         if (typeInfo->equal)
2754                 return typeInfo->equal(data1, data2);
2755         else return !memcmp(data1, data2, typeInfo->size);
2756 }
2757
2758 void CustomData_data_initminmax(int type, void *min, void *max)
2759 {
2760         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2761
2762         if (typeInfo->initminmax)
2763                 typeInfo->initminmax(min, max);
2764 }
2765
2766
2767 void CustomData_data_dominmax(int type, void *data, void *min, void *max)
2768 {
2769         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2770
2771         if (typeInfo->dominmax)
2772                 typeInfo->dominmax(data, min, max);
2773 }
2774
2775
2776 void CustomData_data_multiply(int type, void *data, float fac)
2777 {
2778         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2779
2780         if (typeInfo->multiply)
2781                 typeInfo->multiply(data, fac);
2782 }
2783
2784
2785 void CustomData_data_add(int type, void *data1, void *data2)
2786 {
2787         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2788
2789         if (typeInfo->add)
2790                 typeInfo->add(data1, data2);
2791 }
2792
2793 void CustomData_bmesh_set(const CustomData *data, void *block, int type, void *source)
2794 {
2795         void *dest = CustomData_bmesh_get(data, block, type);
2796         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2797
2798         if (!dest) return;
2799
2800         if (typeInfo->copy)
2801                 typeInfo->copy(source, dest, 1);
2802         else
2803                 memcpy(dest, source, typeInfo->size);
2804 }
2805
2806 void CustomData_bmesh_set_n(CustomData *data, void *block, int type, int n, void *source)
2807 {
2808         void *dest = CustomData_bmesh_get_n(data, block, type, n);
2809         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
2810
2811         if (!dest) return;
2812
2813         if (typeInfo->copy)
2814                 typeInfo->copy(source, dest, 1);
2815         else
2816                 memcpy(dest, source, typeInfo->size);
2817 }
2818
2819 void CustomData_bmesh_set_layer_n(CustomData *data, void *block, int n, void *source)
2820 {
2821         void *dest = CustomData_bmesh_get_layer_n(data, block, n);
2822         const LayerTypeInfo *typeInfo = layerType_getInfo(data->layers[n].type);
2823
2824         if (!dest) return;
2825
2826         if (typeInfo->copy)
2827                 typeInfo->copy(source, dest, 1);
2828         else
2829                 memcpy(dest, source, typeInfo->size);
2830 }
2831
2832 /**
2833  * \param src_blocks must be pointers to the data, offset by layer->offset already.
2834  */
2835 void CustomData_bmesh_interp_n(CustomData *data, void **src_blocks, const float *weights,
2836                                const float *sub_weights, int count, void *dest_block, int n)
2837 {
2838         CustomDataLayer *layer = &data->layers[n];
2839         const LayerTypeInfo *typeInfo = layerType_getInfo(layer->type);
2840
2841         typeInfo->interp(src_blocks, weights, sub_weights, count,
2842                          (char *)dest_block + layer->offset);
2843 }
2844
2845 void CustomData_bmesh_interp(CustomData *data, void **src_blocks, const float *weights,
2846                              const float *sub_weights, int count, void *dest_block)
2847 {
2848         int i, j;
2849         void *source_buf[SOURCE_BUF_SIZE];
2850         void **sources = source_buf;
2851
2852         /* slow fallback in case we're interpolating a ridiculous number of
2853          * elements
2854          */
2855         if (count > SOURCE_BUF_SIZE)
2856                 sources = MEM_callocN(sizeof(*sources) * count,
2857                                       "CustomData_interp sources");
2858
2859         /* interpolates a layer at a time */
2860         for (i = 0; i < data->totlayer; ++i) {
2861                 CustomDataLayer *layer = &data->layers[i];
2862                 const LayerTypeInfo *typeInfo = layerType_getInfo(layer->type);
2863                 if (typeInfo->interp) {
2864                         for (j = 0; j < count; ++j) {
2865                                 sources[j] = (char *)src_blocks[j] + layer->offset;
2866                         }
2867                         CustomData_bmesh_interp_n(data, sources, weights, sub_weights, count, dest_block, i);
2868                 }
2869         }
2870
2871         if (count > SOURCE_BUF_SIZE) MEM_freeN(sources);
2872 }
2873
2874 static void CustomData_bmesh_set_default_n(CustomData *data, void **block, int n)
2875 {
2876         const LayerTypeInfo *typeInfo;
2877         int offset = data->layers[n].offset;
2878
2879         typeInfo = layerType_getInfo(data->layers[n].type);
2880
2881         if (typeInfo->set_default) {
2882                 typeInfo->set_default((char *)*block + offset, 1);
2883         }
2884         else {
2885                 memset((char *)*block + offset, 0, typeInfo->size);
2886         }
2887 }
2888
2889 void CustomData_bmesh_set_default(CustomData *data, void **block)
2890 {
2891         int i;
2892
2893         if (*block == NULL)
2894                 CustomData_bmesh_alloc_block(data, block);
2895
2896         for (i = 0; i < data->totlayer; ++i) {
2897                 CustomData_bmesh_set_default_n(data, block, i);
2898         }
2899 }
2900
2901 /**
2902  * \param use_default_init initializes data which can't be copied,
2903  * typically you'll want to use this if the BM_xxx create function
2904  * is called with BM_CREATE_SKIP_CD flag
2905  */
2906 void CustomData_to_bmesh_block(const CustomData *source, CustomData *dest,
2907                                int src_index, void **dest_block, bool use_default_init)
2908 {
2909         const LayerTypeInfo *typeInfo;
2910         int dest_i, src_i, src_offset;
2911
2912         if (*dest_block == NULL)
2913                 CustomData_bmesh_alloc_block(dest, dest_block);
2914         
2915         /* copies a layer at a time */
2916         dest_i = 0;
2917         for (src_i = 0; src_i < source->totlayer; ++src_i) {
2918
2919                 /* find the first dest layer with type >= the source type
2920                  * (this should work because layers are ordered by type)
2921                  */
2922                 while (dest_i < dest->totlayer && dest->layers[dest_i].type < source->layers[src_i].type) {
2923                         if (use_default_init) {
2924                                 CustomData_bmesh_set_default_n(dest, dest_block, dest_i);
2925                         }
2926                         dest_i++;
2927                 }
2928
2929                 /* if there are no more dest layers, we're done */
2930                 if (dest_i >= dest->totlayer) break;
2931
2932                 /* if we found a matching layer, copy the data */
2933                 if (dest->layers[dest_i].type == source->layers[src_i].type) {
2934                         int offset = dest->layers[dest_i].offset;
2935                         char *src_data = source->layers[src_i].data;
2936                         char *dest_data = (char *)*dest_block + offset;
2937
2938                         typeInfo = layerType_getInfo(dest->layers[dest_i].type);
2939                         src_offset = src_index * typeInfo->size;
2940
2941                         if (typeInfo->copy)
2942                                 typeInfo->copy(src_data + src_offset, dest_data, 1);
2943                         else
2944                                 memcpy(dest_data, src_data + src_offset, typeInfo->size);
2945
2946                         /* if there are multiple source & dest layers of the same type,
2947                          * we don't want to copy all source layers to the same dest, so
2948                          * increment dest_i
2949                          */
2950                         dest_i++;
2951                 }
2952         }
2953
2954         if (use_default_init) {
2955                 while (dest_i < dest->totlayer) {
2956                         CustomData_bmesh_set_default_n(dest, dest_block, dest_i);
2957                         dest_i++;
2958                 }
2959         }
2960 }
2961
2962 void CustomData_from_bmesh_block(const CustomData *source, CustomData *dest,
2963                                  void *src_block, int dest_index)
2964 {
2965         const LayerTypeInfo *typeInfo;
2966         int dest_i, src_i, dest_offset;
2967
2968         /* copies a layer at a time */
2969         dest_i = 0;
2970         for (src_i = 0; src_i < source->totlayer; ++src_i) {
2971
2972                 /* find the first dest layer with type >= the source type
2973                  * (this should work because layers are ordered by type)
2974                  */
2975                 while (dest_i < dest->totlayer && dest->layers[dest_i].type < source->layers[src_i].type) {
2976                         dest_i++;
2977                 }
2978
2979                 /* if there are no more dest layers, we're done */
2980                 if (dest_i >= dest->totlayer) return;
2981
2982                 /* if we found a matching layer, copy the data */
2983                 if (dest->layers[dest_i].type == source->layers[src_i].type) {
2984                         int offset = source->layers[src_i].offset;
2985                         char *src_data = (char *)src_block + offset;
2986                         char *dest_data = dest->layers[dest_i].data;
2987
2988                         typeInfo = layerType_getInfo(dest->layers[dest_i].type);
2989                         dest_offset = dest_index * typeInfo->size;
2990
2991                         if (typeInfo->copy)
2992                                 typeInfo->copy(src_data, dest_data + dest_offset, 1);
2993                         else
2994                                 memcpy(dest_data + dest_offset, src_data, typeInfo->size);
2995
2996                         /* if there are multiple source & dest layers of the same type,
2997                          * we don't want to copy all source layers to the same dest, so
2998                          * increment dest_i
2999                          */
3000                         dest_i++;
3001                 }
3002         }
3003
3004 }
3005
3006 void CustomData_file_write_info(int type, const char **structname, int *structnum)
3007 {
3008         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
3009
3010         *structname = typeInfo->structname;
3011         *structnum = typeInfo->structnum;
3012 }
3013
3014 int CustomData_sizeof(int type)
3015 {
3016         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
3017
3018         return typeInfo->size;
3019 }
3020
3021 const char *CustomData_layertype_name(int type)
3022 {
3023         return layerType_getName(type);
3024 }
3025
3026
3027 /**
3028  * Can only ever be one of these.
3029  */
3030 bool CustomData_layertype_is_singleton(int type)
3031 {
3032         const LayerTypeInfo *typeInfo = layerType_getInfo(type);
3033         return typeInfo->defaultname == NULL;
3034 }
3035
3036 static bool CustomData_is_property_layer(int type)
3037 {
3038         if ((type == CD_PROP_FLT) || (type == CD_PROP_INT) || (type == CD_PROP_STR))
3039                 return true;
3040         return false;
3041 }
3042
3043 static bool cd_layer_find_dupe(CustomData *data, const char *name, int type, int index)
3044 {
3045         int i;
3046         /* see if there is a duplicate */
3047         for (i = 0; i < data->totlayer; i++) {
3048                 if (i != index) {
3049                         CustomDataLayer *layer = &data->layers[i];
3050                         
3051                         if (CustomData_is_property_layer(type)) {
3052                                 if (CustomData_is_property_layer(layer->type) && strcmp(layer->name, name) == 0) {
3053                                         return true;
3054                                 }
3055                         }
3056                         else {
3057                                 if (i != index && layer->type == type && strcmp(layer->name, name) == 0) {
3058                                         return true;
3059                                 }
3060                         }
3061                 }
3062         }
3063         
3064         return false;
3065 }
3066
3067 static bool customdata_unique_check(void *arg, const char *name)
3068 {
3069         struct {CustomData *data; int type; int index; } *data_arg = arg;
3070         return cd_layer_find_dupe(data_arg->data, name, data_arg->type, data_arg->index);
3071 }
3072
3073 void CustomData_set_layer_unique_name(CustomData *data, int index)
3074 {       
3075         CustomDataLayer *nlayer = &data->layers[index];
3076         const LayerTypeInfo *typeInfo = layerType_getInfo(nlayer->type);
3077
3078         struct {CustomData *data; int type; int index; } data_arg;
3079         data_arg.data = data;
3080         data_arg.type = nlayer->type;
3081         data_arg.index = index;
3082
3083         if (!typeInfo->defaultname)
3084                 return;
3085
3086         BLI_uniquename_cb(customdata_unique_check, &data_arg, DATA_(typeInfo->defaultname), '.', nlayer->name,
3087                           sizeof(nlayer->name));
3088 }
3089
3090 void CustomData_validate_layer_name(const CustomData *data, int type, const char *name, char *outname)
3091 {
3092         int index = -1;
3093
3094         /* if a layer name was given, try to find that layer */
3095         if (name[0])
3096                 index = CustomData_get_named_layer_index(data, type, name);
3097
3098         if (index == -1) {
3099                 /* either no layer was specified, or the layer we want has been
3100                  * deleted, so assign the active layer to name
3101                  */
3102                 index = CustomData_get_active_layer_index(data, type);
3103                 BLI_strncpy(outname, data->layers[index].name, MAX_CUSTOMDATA_LAYER_NAME);
3104         }
3105         else {
3106                 BLI_strncpy(outname, name, MAX_CUSTOMDATA_LAYER_NAME);
3107         }
3108 }
3109
3110 bool CustomData_verify_versions(struct CustomData *data, int index)
3111 {
3112         const LayerTypeInfo *typeInfo;
3113         CustomDataLayer *layer = &data->layers[index];
3114         bool keeplayer = true;
3115         int i;
3116
3117         if (layer->type >= CD_NUMTYPES) {
3118                 keeplayer = false; /* unknown layer type from future version */
3119         }
3120         else {
3121                 typeInfo = layerType_getInfo(layer->type);
3122
3123                 if (!typeInfo->defaultname && (index > 0) &&
3124                     data->layers[index - 1].type == layer->type)
3125                 {
3126                         keeplayer = false; /* multiple layers of which we only support one */
3127                 }
3128         }
3129
3130         if (!keeplayer) {
3131                 for (i = index + 1; i < data->totlayer; ++i)
3132                         data->layers[i - 1] = data->layers[i];
3133                 data->totlayer--;
3134         }
3135
3136         return keeplayer;
3137 }
3138
3139 /****************************** External Files *******************************/
3140
3141 static void customdata_external_filename(char filename[FILE_MAX], ID *id, CustomDataExternal *external)
3142 {
3143         BLI_strncpy(filename, external->filename, FILE_MAX);
3144         BLI_path_abs(filename, ID_BLEND_PATH(G.main, id));
3145 }
3146
3147 void CustomData_external_reload(CustomData *data, ID *UNUSED(id), CustomDataMask mask, int totelem)
3148 {
3149         CustomDataLayer *layer;
3150         const LayerTypeInfo *typeInfo;
3151         int i;
3152
3153         for (i = 0; i < data->totlayer; i++) {
3154                 layer = &data->layers[i];
3155                 typeInfo = layerType_getInfo(layer->type);
3156
3157                 if (!(mask & CD_TYPE_AS_MASK(layer->type))) {
3158                         /* pass */
3159                 }
3160                 else if ((layer->flag & CD_FLAG_EXTERNAL) && (layer->flag & CD_FLAG_IN_MEMORY)) {
3161                         if (typeInfo->free)
3162                                 typeInfo->free(layer->data, totelem, typeInfo->size);
3163                         layer->flag &= ~CD_FLAG_IN_MEMORY;
3164                 }
3165         }
3166 }
3167
3168 void CustomData_external_read(CustomData *data, ID *id, CustomDataMask mask, int totelem)
3169 {
3170         CustomDataExternal *external = data->external;
3171         CustomDataLayer *layer;
3172         CDataFile *cdf;
3173         CDataFileLayer *blay;
3174         char filename[FILE_MAX];
3175         const LayerTypeInfo *typeInfo;
3176         int i, update = 0;
3177
3178         if (!external)
3179                 return;
3180         
3181         for (i = 0; i < data->totlayer; i++) {
3182                 layer = &data->layers[i];
3183                 typeInfo = layerType_getInfo(layer->type);
3184
3185                 if (!(mask & CD_TYPE_AS_MASK(layer->type))) {
3186                         /* pass */
3187                 }
3188                 else if (layer->flag & CD_FLAG_IN_MEMORY) {
3189                         /* pass */
3190                 }
3191                 else if ((layer->flag & CD_FLAG_EXTERNAL) && typeInfo->read) {
3192                         update = 1;
3193                 }
3194         }
3195
3196         if (!update)
3197                 return;
3198
3199         customdata_external_filename(filename, id, external);
3200
3201         cdf = cdf_create(CDF_TYPE_MESH);
3202         if (!cdf_read_open(cdf, filename)) {
3203                 fprintf(stderr, "Failed to read %s layer from %s.\n", layerType_getName(layer->type), filename);
3204                 return;
3205         }
3206
3207         for (i = 0; i < data->totlayer; i++) {
3208                 layer = &data->layers[i];
3209                 typeInfo = layerType_getInfo(layer->type);
3210
3211                 if (!(mask & CD_TYPE_AS_MASK(layer->type))) {
3212                         /* pass */
3213                 }
3214                 else if (layer->flag & CD_FLAG_IN_MEMORY) {
3215                         /* pass */
3216                 }
3217                 else if ((layer->flag & CD_FLAG_EXTERNAL) && typeInfo->read) {
3218                         blay = cdf_layer_find(cdf, layer->type, layer->name);
3219
3220                         if (blay) {
3221                                 if (cdf_read_layer(cdf, blay)) {
3222                                         if (typeInfo->read(cdf, layer->data, totelem)) {
3223                                                 /* pass */
3224                                         }
3225                                         else {
3226                                                 break;
3227                                         }
3228                                         layer->flag |= CD_FLAG_IN_MEMORY;
3229                                 }
3230                                 else
3231                                         break;
3232                         }
3233                 }
3234         }
3235
3236         cdf_read_close(cdf);
3237         cdf_free(cdf);
3238 }
3239
3240 void CustomData_external_write(CustomData *data, ID *id, CustomDataMask mask, int totelem, int free)
3241 {
3242         CustomDataExternal *external = data->external;
3243         CustomDataLayer *layer;
3244         CDataFile *cdf;
3245         CDataFileLayer *blay;
3246         const LayerTypeInfo *typeInfo;
3247         int i, update = 0;
3248         char filename[FILE_MAX];
3249
3250         if (!external)
3251                 return;
3252
3253         /* test if there is anything to write */
3254         for (i = 0; i < data->totlayer; i++) {
3255                 layer = &data->layers[i];
3256                 typeInfo = layerType_getInfo(layer->type);
3257
3258                 if (!(mask & CD_TYPE_AS_MASK(layer->type))) {
3259                         /* pass */
3260                 }
3261                 else if ((layer->flag & CD_FLAG_EXTERNAL) && typeInfo->write) {
3262                         update = 1;
3263                 }
3264         }
3265
3266         if (!update)
3267                 return;
3268
3269         /* make sure data is read before we try to write */
3270         CustomData_external_read(data, id, mask, totelem);
3271         customdata_external_filename(filename, id, external);
3272
3273         cdf = cdf_create(CDF_TYPE_MESH);
3274
3275         for (i = 0; i < data->totlayer; i++) {
3276                 layer = &data->layers[i];
3277                 typeInfo = layerType_getInfo(layer->type);
3278
3279                 if ((layer->flag & CD_FLAG_EXTERNAL) && typeInfo->filesize) {
3280                         if (layer->flag & CD_FLAG_IN_MEMORY) {
3281                                 cdf_layer_add(cdf, layer->type, layer->name,
3282                                               typeInfo->filesize(cdf, layer->data, totelem));
3283                         }
3284                         else {
3285                                 cdf_free(cdf);
3286                                 return; /* read failed for a layer! */
3287                         }
3288                 }
3289         }
3290
3291         if (!cdf_write_open(cdf, filename)) {
3292                 fprintf(stderr, "Failed to open %s for writing.\n", filename);
3293                 return;
3294         }
3295
3296         for (i = 0; i < data->totlayer; i++) {
3297                 layer = &data->layers[i];
3298                 typeInfo = layerType_getInfo(layer->type);
3299
3300                 if ((layer->flag & CD_FLAG_EXTERNAL) && typeInfo->write) {
3301                         blay = cdf_layer_find(cdf, layer->type, layer->name);
3302
3303                         if (cdf_write_layer(cdf, blay)) {
3304                                 if (typeInfo->write(cdf, layer->data, totelem)) {
3305                                         /* pass */
3306                                 }
3307                                 else {
3308                                         break;
3309                                 }
3310                         }
3311                         else {
3312                                 break;
3313                         }
3314                 }
3315         }
3316
3317         if (i != data->totlayer) {
3318                 fprintf(stderr, "Failed to write data to %s.\n", filename);
3319                 cdf_free(cdf);
3320                 return;
3321         }
3322
3323         for (i = 0; i < data->totlayer; i++) {