Merge branch 'master' into blender2.8
[blender.git] / source / blender / modifiers / intern / MOD_weightvgproximity.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) 2011 by Bastien Montagne.
19  * All rights reserved.
20  *
21  * Contributor(s): None yet.
22  *
23  * ***** END GPL LICENSE BLOCK *****
24  *
25  */
26
27 /** \file blender/modifiers/intern/MOD_weightvgproximity.c
28  *  \ingroup modifiers
29  */
30
31 #include "BLI_utildefines.h"
32 #include "BLI_ghash.h"
33 #include "BLI_listbase.h"
34 #include "BLI_math.h"
35 #include "BLI_rand.h"
36 #include "BLI_task.h"
37
38 #include "DNA_mesh_types.h"
39 #include "DNA_meshdata_types.h"
40 #include "DNA_modifier_types.h"
41 #include "DNA_object_types.h"
42
43 #include "BKE_cdderivedmesh.h"
44 #include "BKE_curve.h"
45 #include "BKE_customdata.h"
46 #include "BKE_deform.h"
47 #include "BKE_library.h"
48 #include "BKE_library_query.h"
49 #include "BKE_mesh.h"
50 #include "BKE_modifier.h"
51 #include "BKE_texture.h"          /* Texture masking. */
52
53 #include "DEG_depsgraph_build.h"
54
55 #include "MEM_guardedalloc.h"
56
57 #include "MOD_weightvg_util.h"
58 #include "MOD_modifiertypes.h"
59 #include "MOD_util.h"
60
61 //#define USE_TIMEIT
62
63 #ifdef USE_TIMEIT
64 #  include "PIL_time.h"
65 #  include "PIL_time_utildefines.h"
66 #endif
67
68 /**************************************
69  * Util functions.                    *
70  **************************************/
71
72 /* Util macro. */
73 #define OUT_OF_MEMORY() ((void)printf("WeightVGProximity: Out of memory.\n"))
74
75 typedef struct Vert2GeomData {
76         /* Read-only data */
77         float (*v_cos)[3];
78
79         const SpaceTransform *loc2trgt;
80
81         BVHTreeFromMesh *treeData[3];
82
83         /* Write data, but not needing locking (two different threads will never write same index). */
84         float *dist[3];
85 } Vert2GeomData;
86
87 /* Data which is localized to each computed chunk (i.e. thread-safe, and with continous subset of index range). */
88 typedef struct Vert2GeomDataChunk {
89         /* Read-only data */
90         float last_hit_co[3][3];
91         bool is_init[3];
92 } Vert2GeomDataChunk;
93
94 /**
95  * Callback used by BLI_task 'for loop' helper.
96  */
97 static void vert2geom_task_cb_ex(
98         void *__restrict userdata,
99         const int iter,
100         const ParallelRangeTLS *__restrict tls)
101 {
102         Vert2GeomData *data = userdata;
103         Vert2GeomDataChunk *data_chunk = tls->userdata_chunk;
104
105         float tmp_co[3];
106         int i;
107
108         /* Convert the vertex to tree coordinates. */
109         copy_v3_v3(tmp_co, data->v_cos[iter]);
110         BLI_space_transform_apply(data->loc2trgt, tmp_co);
111
112         for (i = 0; i < ARRAY_SIZE(data->dist); i++) {
113                 if (data->dist[i]) {
114                         BVHTreeNearest nearest = {0};
115
116                         /* Note that we use local proximity heuristics (to reduce the nearest search).
117                          *
118                          * If we already had an hit before in same chunk of tasks (i.e. previous vertex by index),
119                          * we assume this vertex is going to have a close hit to that other vertex, so we can initiate
120                          * the "nearest.dist" with the expected value to that last hit.
121                          * This will lead in pruning of the search tree.
122                          */
123                         nearest.dist_sq = data_chunk->is_init[i] ? len_squared_v3v3(tmp_co, data_chunk->last_hit_co[i]) : FLT_MAX;
124                         nearest.index = -1;
125
126                         /* Compute and store result. If invalid (-1 idx), keep FLT_MAX dist. */
127                         BLI_bvhtree_find_nearest(data->treeData[i]->tree, tmp_co, &nearest,
128                                                  data->treeData[i]->nearest_callback, data->treeData[i]);
129                         data->dist[i][iter] = sqrtf(nearest.dist_sq);
130
131                         if (nearest.index != -1) {
132                                 copy_v3_v3(data_chunk->last_hit_co[i], nearest.co);
133                                 data_chunk->is_init[i] = true;
134                         }
135                 }
136         }
137 }
138
139 /**
140  * Find nearest vertex and/or edge and/or face, for each vertex (adapted from shrinkwrap.c).
141  */
142 static void get_vert2geom_distance(int numVerts, float (*v_cos)[3],
143                                    float *dist_v, float *dist_e, float *dist_f,
144                                    Mesh *target, const SpaceTransform *loc2trgt)
145 {
146         Vert2GeomData data = {0};
147         Vert2GeomDataChunk data_chunk = {{{0}}};
148
149         BVHTreeFromMesh treeData_v = {NULL};
150         BVHTreeFromMesh treeData_e = {NULL};
151         BVHTreeFromMesh treeData_f = {NULL};
152
153         /* XXX TODO horrible, but simpler for now, bvhtree needs some love first! */
154         DerivedMesh *target_dm = CDDM_from_mesh(target);
155
156         if (dist_v) {
157                 /* Create a bvh-tree of the given target's verts. */
158                 bvhtree_from_mesh_get(&treeData_v, target_dm, BVHTREE_FROM_VERTS, 2);
159                 if (treeData_v.tree == NULL) {
160                         OUT_OF_MEMORY();
161                         return;
162                 }
163         }
164         if (dist_e) {
165                 /* Create a bvh-tree of the given target's edges. */
166                 bvhtree_from_mesh_get(&treeData_e, target_dm, BVHTREE_FROM_EDGES, 2);
167                 if (treeData_e.tree == NULL) {
168                         OUT_OF_MEMORY();
169                         return;
170                 }
171         }
172         if (dist_f) {
173                 /* Create a bvh-tree of the given target's faces. */
174                 bvhtree_from_mesh_get(&treeData_f, target_dm, BVHTREE_FROM_LOOPTRI, 2);
175                 if (treeData_f.tree == NULL) {
176                         OUT_OF_MEMORY();
177                         return;
178                 }
179         }
180
181         data.v_cos = v_cos;
182         data.loc2trgt = loc2trgt;
183         data.treeData[0] = &treeData_v;
184         data.treeData[1] = &treeData_e;
185         data.treeData[2] = &treeData_f;
186         data.dist[0] = dist_v;
187         data.dist[1] = dist_e;
188         data.dist[2] = dist_f;
189
190         ParallelRangeSettings settings;
191         BLI_parallel_range_settings_defaults(&settings);
192         settings.use_threading = (numVerts > 10000);
193         settings.userdata_chunk = &data_chunk;
194         settings.userdata_chunk_size = sizeof(data_chunk);
195         BLI_task_parallel_range(
196                 0, numVerts,
197                 &data,
198                 vert2geom_task_cb_ex,
199                 &settings);
200
201         if (dist_v)
202                 free_bvhtree_from_mesh(&treeData_v);
203         if (dist_e)
204                 free_bvhtree_from_mesh(&treeData_e);
205         if (dist_f)
206                 free_bvhtree_from_mesh(&treeData_f);
207
208         target_dm->release(target_dm);
209 }
210
211 /**
212  * Returns the real distance between a vertex and another reference object.
213  * Note that it works in final world space (i.e. with constraints etc. applied).
214  */
215 static void get_vert2ob_distance(int numVerts, float (*v_cos)[3], float *dist,
216                                  Object *ob, Object *obr)
217 {
218         /* Vertex and ref object coordinates. */
219         float v_wco[3];
220         unsigned int i = numVerts;
221
222         while (i-- > 0) {
223                 /* Get world-coordinates of the vertex (constraints and anim included). */
224                 mul_v3_m4v3(v_wco, ob->obmat, v_cos[i]);
225                 /* Return distance between both coordinates. */
226                 dist[i] = len_v3v3(v_wco, obr->obmat[3]);
227         }
228 }
229
230 /**
231  * Returns the real distance between an object and another reference object.
232  * Note that it works in final world space (i.e. with constraints etc. applied).
233  */
234 static float get_ob2ob_distance(const Object *ob, const Object *obr)
235 {
236         return len_v3v3(ob->obmat[3], obr->obmat[3]); 
237 }
238
239 /**
240  * Maps distances to weights, with an optional "smoothing" mapping.
241  */
242 static void do_map(Object *ob, float *weights, const int nidx, const float min_d, const float max_d, short mode)
243 {
244         const float range_inv = 1.0f / (max_d - min_d); /* invert since multiplication is faster */
245         unsigned int i = nidx;
246         if (max_d == min_d) {
247                 while (i-- > 0) {
248                         weights[i] = (weights[i] >= max_d) ? 1.0f : 0.0f; /* "Step" behavior... */
249                 }
250         }
251         else if (max_d > min_d) {
252                 while (i-- > 0) {
253                         if     (weights[i] >= max_d) weights[i] = 1.0f;  /* most likely case first */
254                         else if (weights[i] <= min_d) weights[i] = 0.0f;
255                         else weights[i] = (weights[i] - min_d) * range_inv;
256                 }
257         }
258         else {
259                 while (i-- > 0) {
260                         if     (weights[i] <= max_d) weights[i] = 1.0f;  /* most likely case first */
261                         else if (weights[i] >= min_d) weights[i] = 0.0f;
262                         else weights[i] = (weights[i] - min_d) * range_inv;
263                 }
264         }
265
266         if (!ELEM(mode, MOD_WVG_MAPPING_NONE, MOD_WVG_MAPPING_CURVE)) {
267                 RNG *rng = NULL;
268
269                 if (mode == MOD_WVG_MAPPING_RANDOM) {
270                         rng = BLI_rng_new_srandom(BLI_ghashutil_strhash(ob->id.name + 2));
271                 }
272
273                 weightvg_do_map(nidx, weights, mode, NULL, rng);
274
275                 if (rng) {
276                         BLI_rng_free(rng);
277                 }
278         }
279 }
280
281 /**************************************
282  * Modifiers functions.               *
283  **************************************/
284 static void initData(ModifierData *md)
285 {
286         WeightVGProximityModifierData *wmd = (WeightVGProximityModifierData *) md;
287
288         wmd->proximity_mode       = MOD_WVG_PROXIMITY_OBJECT;
289         wmd->proximity_flags      = MOD_WVG_PROXIMITY_GEOM_VERTS;
290
291         wmd->falloff_type         = MOD_WVG_MAPPING_NONE;
292
293         wmd->mask_constant        = 1.0f;
294         wmd->mask_tex_use_channel = MOD_WVG_MASK_TEX_USE_INT; /* Use intensity by default. */
295         wmd->mask_tex_mapping     = MOD_DISP_MAP_LOCAL;
296         wmd->max_dist             = 1.0f; /* vert arbitrary distance, but don't use 0 */
297 }
298
299 static CustomDataMask requiredDataMask(Object *UNUSED(ob), ModifierData *md)
300 {
301         WeightVGProximityModifierData *wmd = (WeightVGProximityModifierData *) md;
302         CustomDataMask dataMask = 0;
303
304         /* We need vertex groups! */
305         dataMask |= CD_MASK_MDEFORMVERT;
306
307         /* Ask for UV coordinates if we need them. */
308         if (wmd->mask_tex_mapping == MOD_DISP_MAP_UV)
309                 dataMask |= CD_MASK_MTFACE;
310
311         /* No need to ask for CD_PREVIEW_MLOOPCOL... */
312
313         return dataMask;
314 }
315
316 static bool dependsOnTime(ModifierData *md)
317 {
318         WeightVGProximityModifierData *wmd = (WeightVGProximityModifierData *) md;
319
320         if (wmd->mask_texture)
321                 return BKE_texture_dependsOnTime(wmd->mask_texture);
322         return 0;
323 }
324
325 static void foreachObjectLink(ModifierData *md, Object *ob, ObjectWalkFunc walk, void *userData)
326 {
327         WeightVGProximityModifierData *wmd = (WeightVGProximityModifierData *) md;
328         walk(userData, ob, &wmd->proximity_ob_target, IDWALK_CB_NOP);
329         walk(userData, ob, &wmd->mask_tex_map_obj, IDWALK_CB_NOP);
330 }
331
332 static void foreachIDLink(ModifierData *md, Object *ob, IDWalkFunc walk, void *userData)
333 {
334         WeightVGProximityModifierData *wmd = (WeightVGProximityModifierData *) md;
335
336         walk(userData, ob, (ID **)&wmd->mask_texture, IDWALK_CB_USER);
337
338         foreachObjectLink(md, ob, (ObjectWalkFunc)walk, userData);
339 }
340
341 static void foreachTexLink(ModifierData *md, Object *ob, TexWalkFunc walk, void *userData)
342 {
343         walk(userData, ob, md, "mask_texture");
344 }
345
346 static void updateDepsgraph(ModifierData *md, const ModifierUpdateDepsgraphContext *ctx)
347 {
348         WeightVGProximityModifierData *wmd = (WeightVGProximityModifierData *)md;
349         if (wmd->proximity_ob_target != NULL) {
350                 DEG_add_object_relation(ctx->node, wmd->proximity_ob_target, DEG_OB_COMP_TRANSFORM, "WeightVGProximity Modifier");
351                 DEG_add_object_relation(ctx->node, wmd->proximity_ob_target, DEG_OB_COMP_GEOMETRY, "WeightVGProximity Modifier");
352         }
353         if (wmd->mask_tex_map_obj != NULL && wmd->mask_tex_mapping == MOD_DISP_MAP_OBJECT) {
354                 DEG_add_object_relation(ctx->node, wmd->mask_tex_map_obj, DEG_OB_COMP_TRANSFORM, "WeightVGProximity Modifier");
355                 DEG_add_object_relation(ctx->node, wmd->mask_tex_map_obj, DEG_OB_COMP_GEOMETRY, "WeightVGProximity Modifier");
356         }
357         if (wmd->mask_tex_mapping == MOD_DISP_MAP_GLOBAL) {
358                 DEG_add_object_relation(ctx->node, ctx->object, DEG_OB_COMP_TRANSFORM, "WeightVGProximity Modifier");
359                 DEG_add_object_relation(ctx->node, ctx->object, DEG_OB_COMP_GEOMETRY, "WeightVGProximity Modifier");
360         }
361 }
362
363 static bool isDisabled(ModifierData *md, int UNUSED(useRenderParams))
364 {
365         WeightVGProximityModifierData *wmd = (WeightVGProximityModifierData *) md;
366         /* If no vertex group, bypass. */
367         if (wmd->defgrp_name[0] == '\0') return 1;
368         /* If no target object, bypass. */
369         return (wmd->proximity_ob_target == NULL);
370 }
371
372 static Mesh *applyModifier(ModifierData *md, const ModifierEvalContext *ctx, Mesh *mesh)
373 {
374         BLI_assert(mesh != NULL);
375
376         WeightVGProximityModifierData *wmd = (WeightVGProximityModifierData *) md;
377         MDeformVert *dvert = NULL;
378         MDeformWeight **dw, **tdw;
379         float (*v_cos)[3] = NULL; /* The vertices coordinates. */
380         Object *ob = ctx->object;
381         Object *obr = NULL; /* Our target object. */
382         int defgrp_index;
383         float *tw = NULL;
384         float *org_w = NULL;
385         float *new_w = NULL;
386         int *tidx, *indices = NULL;
387         int numIdx = 0;
388         int i;
389         /* Flags. */
390 #if 0
391         const bool do_prev = (wmd->modifier.mode & eModifierMode_DoWeightPreview) != 0;
392 #endif
393
394 #ifdef USE_TIMEIT
395         TIMEIT_START(perf);
396 #endif
397
398         /* Get number of verts. */
399         const int numVerts = mesh->totvert;
400
401         /* Check if we can just return the original mesh.
402          * Must have verts and therefore verts assigned to vgroups to do anything useful!
403          */
404         if ((numVerts == 0) || BLI_listbase_is_empty(&ctx->object->defbase)) {
405                 return mesh;
406         }
407
408         /* Get our target object. */
409         obr = wmd->proximity_ob_target;
410         if (obr == NULL) {
411                 return mesh;
412         }
413
414         /* Get vgroup idx from its name. */
415         defgrp_index = defgroup_name_index(ob, wmd->defgrp_name);
416         if (defgrp_index == -1) {
417                 return mesh;
418         }
419
420         const bool has_mdef = CustomData_has_layer(&mesh->vdata, CD_MDEFORMVERT);
421         /* If no vertices were ever added to an object's vgroup, dvert might be NULL. */
422         /* As this modifier never add vertices to vgroup, just return. */
423         if (!has_mdef) {
424                 return mesh;
425         }
426
427         Mesh *result;
428         BKE_id_copy_ex(
429                     NULL, &mesh->id, (ID **)&result,
430                     LIB_ID_CREATE_NO_MAIN |
431                     LIB_ID_CREATE_NO_USER_REFCOUNT |
432                     LIB_ID_CREATE_NO_DEG_TAG|
433                     LIB_ID_COPY_NO_PREVIEW,
434                     false);
435
436         if (has_mdef) {
437                 dvert = CustomData_get_layer(&result->vdata, CD_MDEFORMVERT);
438         }
439         else {
440                 /* Add a valid data layer! */
441                 dvert = CustomData_add_layer(&result->vdata, CD_MDEFORMVERT, CD_CALLOC, NULL, numVerts);
442         }
443         /* Ultimate security check. */
444         if (!dvert) {
445                 BKE_id_free(NULL, result);
446                 return mesh;
447         }
448
449         /* Find out which vertices to work on (all vertices in vgroup), and get their relevant weight. */
450         tidx = MEM_malloc_arrayN(numVerts, sizeof(int), "WeightVGProximity Modifier, tidx");
451         tw = MEM_malloc_arrayN(numVerts, sizeof(float), "WeightVGProximity Modifier, tw");
452         tdw = MEM_malloc_arrayN(numVerts, sizeof(MDeformWeight *), "WeightVGProximity Modifier, tdw");
453         for (i = 0; i < numVerts; i++) {
454                 MDeformWeight *_dw = defvert_find_index(&dvert[i], defgrp_index);
455                 if (_dw) {
456                         tidx[numIdx] = i;
457                         tw[numIdx] = _dw->weight;
458                         tdw[numIdx++] = _dw;
459                 }
460         }
461         /* If no vertices found, return org data! */
462         if (numIdx == 0) {
463                 MEM_freeN(tidx);
464                 MEM_freeN(tw);
465                 MEM_freeN(tdw);
466                 BKE_id_free(NULL, result);
467                 return mesh;
468         }
469         if (numIdx != numVerts) {
470                 indices = MEM_malloc_arrayN(numIdx, sizeof(int), "WeightVGProximity Modifier, indices");
471                 memcpy(indices, tidx, sizeof(int) * numIdx);
472                 org_w = MEM_malloc_arrayN(numIdx, sizeof(float), "WeightVGProximity Modifier, org_w");
473                 memcpy(org_w, tw, sizeof(float) * numIdx);
474                 dw = MEM_malloc_arrayN(numIdx, sizeof(MDeformWeight *), "WeightVGProximity Modifier, dw");
475                 memcpy(dw, tdw, sizeof(MDeformWeight *) * numIdx);
476                 MEM_freeN(tw);
477                 MEM_freeN(tdw);
478         }
479         else {
480                 org_w = tw;
481                 dw = tdw;
482         }
483         new_w = MEM_malloc_arrayN(numIdx, sizeof(float), "WeightVGProximity Modifier, new_w");
484         MEM_freeN(tidx);
485
486         /* Get our vertex coordinates. */
487         if (numIdx != numVerts) {
488                 float (*tv_cos)[3] = BKE_mesh_vertexCos_get(result, NULL);
489                 v_cos = MEM_malloc_arrayN(numIdx, sizeof(float[3]), "WeightVGProximity Modifier, v_cos");
490                 for (i = 0; i < numIdx; i++) {
491                         copy_v3_v3(v_cos[i], tv_cos[indices[i]]);
492                 }
493                 MEM_freeN(tv_cos);
494         }
495         else {
496                 v_cos = BKE_mesh_vertexCos_get(result, NULL);
497         }
498
499         /* Compute wanted distances. */
500         if (wmd->proximity_mode == MOD_WVG_PROXIMITY_OBJECT) {
501                 const float dist = get_ob2ob_distance(ob, obr);
502                 for (i = 0; i < numIdx; i++) {
503                         new_w[i] = dist;
504                 }
505         }
506         else if (wmd->proximity_mode == MOD_WVG_PROXIMITY_GEOMETRY) {
507                 const bool use_trgt_verts = (wmd->proximity_flags & MOD_WVG_PROXIMITY_GEOM_VERTS) != 0;
508                 const bool use_trgt_edges = (wmd->proximity_flags & MOD_WVG_PROXIMITY_GEOM_EDGES) != 0;
509                 const bool use_trgt_faces = (wmd->proximity_flags & MOD_WVG_PROXIMITY_GEOM_FACES) != 0;
510
511                 if (use_trgt_verts || use_trgt_edges || use_trgt_faces) {
512                         Mesh *target_mesh = get_mesh_eval_for_modifier(obr, ctx->flag);
513
514                         /* We must check that we do have a valid target_mesh! */
515                         if (target_mesh != NULL) {
516                                 SpaceTransform loc2trgt;
517                                 float *dists_v = use_trgt_verts ? MEM_malloc_arrayN(numIdx, sizeof(float), "dists_v") : NULL;
518                                 float *dists_e = use_trgt_edges ? MEM_malloc_arrayN(numIdx, sizeof(float), "dists_e") : NULL;
519                                 float *dists_f = use_trgt_faces ? MEM_malloc_arrayN(numIdx, sizeof(float), "dists_f") : NULL;
520
521                                 BLI_SPACE_TRANSFORM_SETUP(&loc2trgt, ob, obr);
522                                 get_vert2geom_distance(numIdx, v_cos, dists_v, dists_e, dists_f,
523                                                        target_mesh, &loc2trgt);
524                                 for (i = 0; i < numIdx; i++) {
525                                         new_w[i] = dists_v ? dists_v[i] : FLT_MAX;
526                                         if (dists_e)
527                                                 new_w[i] = min_ff(dists_e[i], new_w[i]);
528                                         if (dists_f)
529                                                 new_w[i] = min_ff(dists_f[i], new_w[i]);
530                                 }
531
532                                 MEM_SAFE_FREE(dists_v);
533                                 MEM_SAFE_FREE(dists_e);
534                                 MEM_SAFE_FREE(dists_f);
535                         }
536                         /* Else, fall back to default obj2vert behavior. */
537                         else {
538                                 get_vert2ob_distance(numIdx, v_cos, new_w, ob, obr);
539                         }
540                 }
541                 else {
542                         get_vert2ob_distance(numIdx, v_cos, new_w, ob, obr);
543                 }
544         }
545
546         /* Map distances to weights. */
547         do_map(ob, new_w, numIdx, wmd->min_dist, wmd->max_dist, wmd->falloff_type);
548
549         /* Do masking. */
550         weightvg_do_mask(numIdx, indices, org_w, new_w, ob, result, wmd->mask_constant,
551                          wmd->mask_defgrp_name, wmd->modifier.scene, wmd->mask_texture,
552                          wmd->mask_tex_use_channel, wmd->mask_tex_mapping,
553                          wmd->mask_tex_map_obj, wmd->mask_tex_uvlayer_name);
554
555         /* Update vgroup. Note we never add nor remove vertices from vgroup here. */
556         weightvg_update_vg(dvert, defgrp_index, dw, numIdx, indices, org_w, false, 0.0f, false, 0.0f);
557
558         /* If weight preview enabled... */
559 #if 0 /* XXX Currently done in mod stack :/ */
560         if (do_prev)
561                 DM_update_weight_mcol(ob, dm, 0, org_w, numIdx, indices);
562 #endif
563
564         /* Freeing stuff. */
565         MEM_freeN(org_w);
566         MEM_freeN(new_w);
567         MEM_freeN(dw);
568         MEM_freeN(v_cos);
569         MEM_SAFE_FREE(indices);
570
571 #ifdef USE_TIMEIT
572         TIMEIT_END(perf);
573 #endif
574
575         /* Return the vgroup-modified mesh. */
576         return result;
577 }
578
579
580 ModifierTypeInfo modifierType_WeightVGProximity = {
581         /* name */              "VertexWeightProximity",
582         /* structName */        "WeightVGProximityModifierData",
583         /* structSize */        sizeof(WeightVGProximityModifierData),
584         /* type */              eModifierTypeType_NonGeometrical,
585         /* flags */             eModifierTypeFlag_AcceptsMesh |
586                                 eModifierTypeFlag_SupportsMapping |
587                                 eModifierTypeFlag_SupportsEditmode |
588                                 eModifierTypeFlag_UsesPreview,
589
590         /* copyData */          modifier_copyData_generic,
591
592         /* deformVerts_DM */    NULL,
593         /* deformMatrices_DM */ NULL,
594         /* deformVertsEM_DM */  NULL,
595         /* deformMatricesEM_DM*/NULL,
596         /* applyModifier_DM */  NULL,
597         /* applyModifierEM_DM */NULL,
598
599         /* deformVerts */       NULL,
600         /* deformMatrices */    NULL,
601         /* deformVertsEM */     NULL,
602         /* deformMatricesEM */  NULL,
603         /* applyModifier */     applyModifier,
604         /* applyModifierEM */   NULL,
605
606         /* initData */          initData,
607         /* requiredDataMask */  requiredDataMask,
608         /* freeData */          NULL,
609         /* isDisabled */        isDisabled,
610         /* updateDepsgraph */   updateDepsgraph,
611         /* dependsOnTime */     dependsOnTime,
612         /* dependsOnNormals */  NULL,
613         /* foreachObjectLink */ foreachObjectLink,
614         /* foreachIDLink */     foreachIDLink,
615         /* foreachTexLink */    foreachTexLink,
616 };