replace malloc of a single DeformVertex with stack variable.
[blender.git] / source / blender / editors / sculpt_paint / paint_vertex.c
1 /*
2  * $Id$
3  *
4  * ***** BEGIN GPL LICENSE BLOCK *****
5  *
6  * This program is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU General Public License
8  * as published by the Free Software Foundation; either version 2
9  * of the License, or (at your option) any later version.
10  *
11  * This program is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with this program; if not, write to the Free Software Foundation,
18  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
19  *
20  * The Original Code is Copyright (C) 2001-2002 by NaN Holding BV.
21  * All rights reserved.
22  *
23  * The Original Code is: all of this file.
24  *
25  * Contributor(s): none yet.
26  *
27  * ***** END GPL LICENSE BLOCK *****
28  */
29
30 /** \file blender/editors/sculpt_paint/paint_vertex.c
31  *  \ingroup edsculpt
32  */
33
34
35 #include <math.h>
36 #include <string.h>
37
38 #ifdef WIN32
39 #include <io.h>
40 #else
41 #include <unistd.h>
42 #endif   
43
44 #include "MEM_guardedalloc.h"
45
46 #include "BLI_blenlib.h"
47 #include "BLI_math.h"
48 #include "BLI_utildefines.h"
49 #include "BLI_ghash.h"
50
51 #include "IMB_imbuf.h"
52 #include "IMB_imbuf_types.h"
53
54 #include "DNA_armature_types.h"
55 #include "DNA_mesh_types.h"
56 #include "DNA_particle_types.h"
57 #include "DNA_scene_types.h"
58 #include "DNA_brush_types.h"
59 #include "DNA_object_types.h"
60 #include "DNA_meshdata_types.h"
61
62 #include "RNA_access.h"
63 #include "RNA_define.h"
64 #include "RNA_enum_types.h"
65
66 #include "BKE_DerivedMesh.h"
67 #include "BKE_armature.h"
68 #include "BKE_action.h"
69 #include "BKE_brush.h"
70 #include "BKE_context.h"
71 #include "BKE_depsgraph.h"
72 #include "BKE_deform.h"
73 #include "BKE_mesh.h"
74 #include "BKE_modifier.h"
75 #include "BKE_object.h"
76 #include "BKE_paint.h"
77 #include "BKE_report.h"
78
79 #include "WM_api.h"
80 #include "WM_types.h"
81
82
83 #include "ED_armature.h"
84 #include "ED_mesh.h"
85 #include "ED_screen.h"
86 #include "ED_view3d.h"
87
88 #include "paint_intern.h"
89
90 /* brush->vertexpaint_tool */
91 #define VP_MIX  0
92 #define VP_ADD  1
93 #define VP_SUB  2
94 #define VP_MUL  3
95 #define VP_BLUR 4
96 #define VP_LIGHTEN      5
97 #define VP_DARKEN       6
98
99 /* polling - retrieve whether cursor should be set or operator should be done */
100
101
102 /* Returns true if vertex paint mode is active */
103 int vertex_paint_mode_poll(bContext *C)
104 {
105         Object *ob = CTX_data_active_object(C);
106
107         return ob && ob->mode == OB_MODE_VERTEX_PAINT && ((Mesh *)ob->data)->totface;
108 }
109
110 int vertex_paint_poll(bContext *C)
111 {
112         if(vertex_paint_mode_poll(C) && 
113            paint_brush(&CTX_data_tool_settings(C)->vpaint->paint)) {
114                 ScrArea *sa= CTX_wm_area(C);
115                 if(sa->spacetype==SPACE_VIEW3D) {
116                         ARegion *ar= CTX_wm_region(C);
117                         if(ar->regiontype==RGN_TYPE_WINDOW)
118                                 return 1;
119                         }
120                 }
121         return 0;
122 }
123
124 int weight_paint_mode_poll(bContext *C)
125 {
126         Object *ob = CTX_data_active_object(C);
127
128         return ob && ob->mode == OB_MODE_WEIGHT_PAINT && ((Mesh *)ob->data)->totface;
129 }
130
131 int weight_paint_poll(bContext *C)
132 {
133         Object *ob= CTX_data_active_object(C);
134         ScrArea *sa;
135
136         if(     (ob != NULL) &&
137                 (ob->mode & OB_MODE_WEIGHT_PAINT) &&
138                 (paint_brush(&CTX_data_tool_settings(C)->wpaint->paint) != NULL) &&
139                 (sa= CTX_wm_area(C)) &&
140                 (sa->spacetype == SPACE_VIEW3D)
141         ) {
142                 ARegion *ar= CTX_wm_region(C);
143                 if(ar->regiontype==RGN_TYPE_WINDOW) {
144                         return 1;
145                 }
146         }
147         return 0;
148 }
149
150 static VPaint *new_vpaint(int wpaint)
151 {
152         VPaint *vp= MEM_callocN(sizeof(VPaint), "VPaint");
153         
154         vp->flag= VP_AREA+VP_SPRAY;
155         
156         if(wpaint)
157                 vp->flag= VP_AREA;
158
159         return vp;
160 }
161
162 static int *get_indexarray(Mesh *me)
163 {
164         return MEM_mallocN(sizeof(int)*(me->totface+1), "vertexpaint");
165 }
166
167
168 /* in contradiction to cpack drawing colors, the MCOL colors (vpaint colors) are per byte! 
169    so not endian sensitive. Mcol = ABGR!!! so be cautious with cpack calls */
170
171 static unsigned int rgba_to_mcol(float r, float g, float b, float a)
172 {
173         int ir, ig, ib, ia;
174         unsigned int col;
175         char *cp;
176         
177         ir= floor(255.0f * r);
178         if(ir<0) ir= 0; else if(ir>255) ir= 255;
179         ig= floor(255.0f * g);
180         if(ig<0) ig= 0; else if(ig>255) ig= 255;
181         ib= floor(255.0f * b);
182         if(ib<0) ib= 0; else if(ib>255) ib= 255;
183         ia= floor(255.0f * a);
184         if(ia<0) ia= 0; else if(ia>255) ia= 255;
185         
186         cp= (char *)&col;
187         cp[0]= ia;
188         cp[1]= ib;
189         cp[2]= ig;
190         cp[3]= ir;
191         
192         return col;
193         
194 }
195
196 unsigned int vpaint_get_current_col(VPaint *vp)
197 {
198         Brush *brush = paint_brush(&vp->paint);
199         return rgba_to_mcol(brush->rgb[0], brush->rgb[1], brush->rgb[2], 1.0f);
200 }
201
202 static void do_shared_vertexcol(Mesh *me)
203 {
204         /* if no mcol: do not do */
205         /* if tface: only the involved faces, otherwise all */
206         MFace *mface;
207         MTFace *tface;
208         int a;
209         short *scolmain, *scol;
210         char *mcol;
211         
212         if(me->mcol==NULL || me->totvert==0 || me->totface==0) return;
213         
214         scolmain= MEM_callocN(4*sizeof(short)*me->totvert, "colmain");
215         
216         tface= me->mtface;
217         mface= me->mface;
218         mcol= (char *)me->mcol;
219         for(a=me->totface; a>0; a--, mface++, mcol+=16) {
220                 if((tface && tface->mode & TF_SHAREDCOL) || (me->editflag & ME_EDIT_PAINT_MASK)==0) {
221                         scol= scolmain+4*mface->v1;
222                         scol[0]++; scol[1]+= mcol[1]; scol[2]+= mcol[2]; scol[3]+= mcol[3];
223                         scol= scolmain+4*mface->v2;
224                         scol[0]++; scol[1]+= mcol[5]; scol[2]+= mcol[6]; scol[3]+= mcol[7];
225                         scol= scolmain+4*mface->v3;
226                         scol[0]++; scol[1]+= mcol[9]; scol[2]+= mcol[10]; scol[3]+= mcol[11];
227                         if(mface->v4) {
228                                 scol= scolmain+4*mface->v4;
229                                 scol[0]++; scol[1]+= mcol[13]; scol[2]+= mcol[14]; scol[3]+= mcol[15];
230                         }
231                 }
232                 if(tface) tface++;
233         }
234         
235         a= me->totvert;
236         scol= scolmain;
237         while(a--) {
238                 if(scol[0]>1) {
239                         scol[1]/= scol[0];
240                         scol[2]/= scol[0];
241                         scol[3]/= scol[0];
242                 }
243                 scol+= 4;
244         }
245         
246         tface= me->mtface;
247         mface= me->mface;
248         mcol= (char *)me->mcol;
249         for(a=me->totface; a>0; a--, mface++, mcol+=16) {
250                 if((tface && tface->mode & TF_SHAREDCOL) || (me->editflag & ME_EDIT_PAINT_MASK)==0) {
251                         scol= scolmain+4*mface->v1;
252                         mcol[1]= scol[1]; mcol[2]= scol[2]; mcol[3]= scol[3];
253                         scol= scolmain+4*mface->v2;
254                         mcol[5]= scol[1]; mcol[6]= scol[2]; mcol[7]= scol[3];
255                         scol= scolmain+4*mface->v3;
256                         mcol[9]= scol[1]; mcol[10]= scol[2]; mcol[11]= scol[3];
257                         if(mface->v4) {
258                                 scol= scolmain+4*mface->v4;
259                                 mcol[13]= scol[1]; mcol[14]= scol[2]; mcol[15]= scol[3];
260                         }
261                 }
262                 if(tface) tface++;
263         }
264
265         MEM_freeN(scolmain);
266 }
267
268 static void make_vertexcol(Object *ob)  /* single ob */
269 {
270         Mesh *me;
271         if(!ob || ob->id.lib) return;
272         me= get_mesh(ob);
273         if(me==NULL) return;
274         if(me->edit_mesh) return;
275
276         /* copies from shadedisplist to mcol */
277         if(!me->mcol) {
278                 CustomData_add_layer(&me->fdata, CD_MCOL, CD_CALLOC, NULL, me->totface);
279                 mesh_update_customdata_pointers(me);
280         }
281
282         //if(shade)
283         //      shadeMeshMCol(scene, ob, me);
284         //else
285
286         memset(me->mcol, 255, 4*sizeof(MCol)*me->totface);
287         
288         DAG_id_tag_update(&me->id, 0);
289         
290 }
291
292 /* mirror_vgroup is set to -1 when invalid */
293 static void wpaint_mirror_vgroup_ensure(Object *ob, int *vgroup_mirror)
294 {
295         bDeformGroup *defgroup= BLI_findlink(&ob->defbase, ob->actdef - 1);
296
297         if(defgroup) {
298                 bDeformGroup *curdef;
299                 int mirrdef;
300                 char name[MAXBONENAME];
301
302                 flip_side_name(name, defgroup->name, FALSE);
303
304                 if(strcmp(name, defgroup->name) != 0) {
305                         for (curdef= ob->defbase.first, mirrdef= 0; curdef; curdef=curdef->next, mirrdef++) {
306                                 if (!strcmp(curdef->name, name)) {
307                                         break;
308                                 }
309                         }
310
311                         if(curdef==NULL) {
312                                 int olddef= ob->actdef; /* tsk, ED_vgroup_add sets the active defgroup */
313                                 curdef= ED_vgroup_add_name(ob, name);
314                                 ob->actdef= olddef;
315                         }
316
317                         /* curdef should never be NULL unless this is
318                          * a  lamp and ED_vgroup_add_name fails */
319                         if(curdef) {
320                                 *vgroup_mirror= mirrdef;
321                                 return;
322                         }
323                 }
324         }
325
326         *vgroup_mirror= -1;
327 }
328
329 static void copy_vpaint_prev(VPaint *vp, unsigned int *mcol, int tot)
330 {
331         if(vp->vpaint_prev) {
332                 MEM_freeN(vp->vpaint_prev);
333                 vp->vpaint_prev= NULL;
334         }
335         vp->tot= tot;   
336         
337         if(mcol==NULL || tot==0) return;
338         
339         vp->vpaint_prev= MEM_mallocN(4*sizeof(int)*tot, "vpaint_prev");
340         memcpy(vp->vpaint_prev, mcol, 4*sizeof(int)*tot);
341         
342 }
343
344 static void copy_wpaint_prev (VPaint *wp, MDeformVert *dverts, int dcount)
345 {
346         if (wp->wpaint_prev) {
347                 free_dverts(wp->wpaint_prev, wp->tot);
348                 wp->wpaint_prev= NULL;
349         }
350         
351         if(dverts && dcount) {
352                 
353                 wp->wpaint_prev = MEM_mallocN (sizeof(MDeformVert)*dcount, "wpaint prev");
354                 wp->tot = dcount;
355                 copy_dverts (wp->wpaint_prev, dverts, dcount);
356         }
357 }
358
359
360 void vpaint_fill(Object *ob, unsigned int paintcol)
361 {
362         Mesh *me;
363         MFace *mf;
364         unsigned int *mcol;
365         int i, selected;
366
367         me= get_mesh(ob);
368         if(me==NULL || me->totface==0) return;
369
370         if(!me->mcol)
371                 make_vertexcol(ob);
372
373         selected= (me->editflag & ME_EDIT_PAINT_MASK);
374
375         mf = me->mface;
376         mcol = (unsigned int*)me->mcol;
377         for (i = 0; i < me->totface; i++, mf++, mcol+=4) {
378                 if (!selected || mf->flag & ME_FACE_SEL) {
379                         mcol[0] = paintcol;
380                         mcol[1] = paintcol;
381                         mcol[2] = paintcol;
382                         mcol[3] = paintcol;
383                 }
384         }
385         
386         DAG_id_tag_update(&me->id, 0);
387 }
388
389
390 /* fills in the selected faces with the current weight and vertex group */
391 void wpaint_fill(VPaint *wp, Object *ob, float paintweight)
392 {
393         Mesh *me;
394         MFace *mface;
395         MDeformWeight *dw, *uw;
396         int *indexar;
397         int index, vgroup;
398         unsigned int faceverts[5]={0,0,0,0,0};
399         unsigned char i;
400         int vgroup_mirror= -1;
401         int selected;
402         /* Radish */
403         int use_vert_sel;
404
405         me= ob->data;
406         if(me==NULL || me->totface==0 || me->dvert==NULL || !me->mface) return;
407         
408         selected= (me->editflag & ME_EDIT_PAINT_MASK);
409         /* Radish */
410         use_vert_sel= (me->editflag & ME_EDIT_VERT_SEL) != 0;
411
412         indexar= get_indexarray(me);
413
414         if(selected) {
415                 for(index=0, mface=me->mface; index<me->totface; index++, mface++) {
416                         if((mface->flag & ME_FACE_SEL)==0)
417                                 indexar[index]= 0;
418                         else
419                                 indexar[index]= index+1;
420                 }
421         }
422         else {
423                 for(index=0; index<me->totface; index++)
424                         indexar[index]= index+1;
425         }
426         
427         vgroup= ob->actdef-1;
428
429         /* if mirror painting, find the other group */          
430         if(me->editflag & ME_EDIT_MIRROR_X) {
431                 wpaint_mirror_vgroup_ensure(ob, &vgroup_mirror);
432         }
433         
434         copy_wpaint_prev(wp, me->dvert, me->totvert);
435         
436         for(index=0; index<me->totface; index++) {
437                 if(indexar[index] && indexar[index]<=me->totface) {
438                         mface= me->mface + (indexar[index]-1);
439                         /* just so we can loop through the verts */
440                         faceverts[0]= mface->v1;
441                         faceverts[1]= mface->v2;
442                         faceverts[2]= mface->v3;
443                         faceverts[3]= mface->v4;
444                         for (i=0; i<3 || faceverts[i]; i++) {
445                                 if(!me->dvert[faceverts[i]].flag) {
446                                         /* Radish */
447                                         if(use_vert_sel && ((me->mvert[faceverts[i]].flag & SELECT) == 0)) {
448                                                 continue;
449                                         }
450
451                                         dw= defvert_verify_index(me->dvert+faceverts[i], vgroup);
452                                         if(dw) {
453                                                 uw= defvert_verify_index(wp->wpaint_prev+faceverts[i], vgroup);
454                                                 uw->weight= dw->weight; /* set the undo weight */
455                                                 dw->weight= paintweight;
456                                                 
457                                                 if(me->editflag & ME_EDIT_MIRROR_X) {   /* x mirror painting */
458                                                         int j= mesh_get_x_mirror_vert(ob, faceverts[i]);
459                                                         if(j>=0) {
460                                                                 /* copy, not paint again */
461                                                                 if(vgroup_mirror != -1) {
462                                                                         dw= defvert_verify_index(me->dvert+j, vgroup_mirror);
463                                                                         uw= defvert_verify_index(wp->wpaint_prev+j, vgroup_mirror);
464                                                                 } else {
465                                                                         dw= defvert_verify_index(me->dvert+j, vgroup);
466                                                                         uw= defvert_verify_index(wp->wpaint_prev+j, vgroup);
467                                                                 }
468                                                                 uw->weight= dw->weight; /* set the undo weight */
469                                                                 dw->weight= paintweight;
470                                                         }
471                                                 }
472                                         }
473                                         (me->dvert+faceverts[i])->flag= 1;
474                                 }
475                         }
476                 }
477         }
478         
479         index=0;
480         while (index<me->totvert) {
481                 (me->dvert+index)->flag= 0;
482                 index++;
483         }
484         
485         MEM_freeN(indexar);
486         copy_wpaint_prev(wp, NULL, 0);
487
488         DAG_id_tag_update(&me->id, 0);
489 }
490
491 /* XXX: should be re-implemented as a vertex/weight paint 'color correct' operator
492  
493 void vpaint_dogamma(Scene *scene)
494 {
495         VPaint *vp= scene->toolsettings->vpaint;
496         Mesh *me;
497         Object *ob;
498         float igam, fac;
499         int a, temp;
500         unsigned char *cp, gamtab[256];
501
502         ob= OBACT;
503         me= get_mesh(ob);
504
505         if(!(ob->mode & OB_MODE_VERTEX_PAINT)) return;
506         if(me==0 || me->mcol==0 || me->totface==0) return;
507
508         igam= 1.0/vp->gamma;
509         for(a=0; a<256; a++) {
510                 
511                 fac= ((float)a)/255.0;
512                 fac= vp->mul*pow( fac, igam);
513                 
514                 temp= 255.9*fac;
515                 
516                 if(temp<=0) gamtab[a]= 0;
517                 else if(temp>=255) gamtab[a]= 255;
518                 else gamtab[a]= temp;
519         }
520
521         a= 4*me->totface;
522         cp= (unsigned char *)me->mcol;
523         while(a--) {
524                 
525                 cp[1]= gamtab[ cp[1] ];
526                 cp[2]= gamtab[ cp[2] ];
527                 cp[3]= gamtab[ cp[3] ];
528                 
529                 cp+= 4;
530         }
531 }
532  */
533
534 static unsigned int mcol_blend(unsigned int col1, unsigned int col2, int fac)
535 {
536         char *cp1, *cp2, *cp;
537         int mfac;
538         unsigned int col=0;
539         
540         if(fac==0) return col1;
541         if(fac>=255) return col2;
542
543         mfac= 255-fac;
544         
545         cp1= (char *)&col1;
546         cp2= (char *)&col2;
547         cp=  (char *)&col;
548         
549         cp[0]= 255;
550         cp[1]= (mfac*cp1[1]+fac*cp2[1])/255;
551         cp[2]= (mfac*cp1[2]+fac*cp2[2])/255;
552         cp[3]= (mfac*cp1[3]+fac*cp2[3])/255;
553         
554         return col;
555 }
556
557 static unsigned int mcol_add(unsigned int col1, unsigned int col2, int fac)
558 {
559         char *cp1, *cp2, *cp;
560         int temp;
561         unsigned int col=0;
562         
563         if(fac==0) return col1;
564         
565         cp1= (char *)&col1;
566         cp2= (char *)&col2;
567         cp=  (char *)&col;
568         
569         cp[0]= 255;
570         temp= cp1[1] + ((fac*cp2[1])/255);
571         if(temp>254) cp[1]= 255; else cp[1]= temp;
572         temp= cp1[2] + ((fac*cp2[2])/255);
573         if(temp>254) cp[2]= 255; else cp[2]= temp;
574         temp= cp1[3] + ((fac*cp2[3])/255);
575         if(temp>254) cp[3]= 255; else cp[3]= temp;
576         
577         return col;
578 }
579
580 static unsigned int mcol_sub(unsigned int col1, unsigned int col2, int fac)
581 {
582         char *cp1, *cp2, *cp;
583         int temp;
584         unsigned int col=0;
585         
586         if(fac==0) return col1;
587         
588         cp1= (char *)&col1;
589         cp2= (char *)&col2;
590         cp=  (char *)&col;
591         
592         cp[0]= 255;
593         temp= cp1[1] - ((fac*cp2[1])/255);
594         if(temp<0) cp[1]= 0; else cp[1]= temp;
595         temp= cp1[2] - ((fac*cp2[2])/255);
596         if(temp<0) cp[2]= 0; else cp[2]= temp;
597         temp= cp1[3] - ((fac*cp2[3])/255);
598         if(temp<0) cp[3]= 0; else cp[3]= temp;
599         
600         return col;
601 }
602
603 static unsigned int mcol_mul(unsigned int col1, unsigned int col2, int fac)
604 {
605         char *cp1, *cp2, *cp;
606         int mfac;
607         unsigned int col=0;
608         
609         if(fac==0) return col1;
610
611         mfac= 255-fac;
612         
613         cp1= (char *)&col1;
614         cp2= (char *)&col2;
615         cp=  (char *)&col;
616         
617         /* first mul, then blend the fac */
618         cp[0]= 255;
619         cp[1]= (mfac*cp1[1] + fac*((cp2[1]*cp1[1])/255)  )/255;
620         cp[2]= (mfac*cp1[2] + fac*((cp2[2]*cp1[2])/255)  )/255;
621         cp[3]= (mfac*cp1[3] + fac*((cp2[3]*cp1[3])/255)  )/255;
622
623         
624         return col;
625 }
626
627 static unsigned int mcol_lighten(unsigned int col1, unsigned int col2, int fac)
628 {
629         char *cp1, *cp2, *cp;
630         int mfac;
631         unsigned int col=0;
632         
633         if(fac==0) return col1;
634         if(fac>=255) return col2;
635
636         mfac= 255-fac;
637         
638         cp1= (char *)&col1;
639         cp2= (char *)&col2;
640         cp=  (char *)&col;
641         
642         /* See if are lighter, if so mix, else dont do anything.
643         if the paint col is darker then the original, then ignore */
644         if (cp1[1]+cp1[2]+cp1[3] > cp2[1]+cp2[2]+cp2[3])
645                 return col1;
646         
647         cp[0]= 255;
648         cp[1]= (mfac*cp1[1]+fac*cp2[1])/255;
649         cp[2]= (mfac*cp1[2]+fac*cp2[2])/255;
650         cp[3]= (mfac*cp1[3]+fac*cp2[3])/255;
651         
652         return col;
653 }
654
655 static unsigned int mcol_darken(unsigned int col1, unsigned int col2, int fac)
656 {
657         char *cp1, *cp2, *cp;
658         int mfac;
659         unsigned int col=0;
660         
661         if(fac==0) return col1;
662         if(fac>=255) return col2;
663
664         mfac= 255-fac;
665         
666         cp1= (char *)&col1;
667         cp2= (char *)&col2;
668         cp=  (char *)&col;
669         
670         /* See if were darker, if so mix, else dont do anything.
671         if the paint col is brighter then the original, then ignore */
672         if (cp1[1]+cp1[2]+cp1[3] < cp2[1]+cp2[2]+cp2[3])
673                 return col1;
674         
675         cp[0]= 255;
676         cp[1]= (mfac*cp1[1]+fac*cp2[1])/255;
677         cp[2]= (mfac*cp1[2]+fac*cp2[2])/255;
678         cp[3]= (mfac*cp1[3]+fac*cp2[3])/255;
679         return col;
680 }
681
682 static void vpaint_blend(VPaint *vp, unsigned int *col, unsigned int *colorig, unsigned int paintcol, int alpha)
683 {
684         Brush *brush = paint_brush(&vp->paint);
685
686         if(brush->vertexpaint_tool==VP_MIX || brush->vertexpaint_tool==VP_BLUR) *col= mcol_blend( *col, paintcol, alpha);
687         else if(brush->vertexpaint_tool==VP_ADD) *col= mcol_add( *col, paintcol, alpha);
688         else if(brush->vertexpaint_tool==VP_SUB) *col= mcol_sub( *col, paintcol, alpha);
689         else if(brush->vertexpaint_tool==VP_MUL) *col= mcol_mul( *col, paintcol, alpha);
690         else if(brush->vertexpaint_tool==VP_LIGHTEN) *col= mcol_lighten( *col, paintcol, alpha);
691         else if(brush->vertexpaint_tool==VP_DARKEN) *col= mcol_darken( *col, paintcol, alpha);
692         
693         /* if no spray, clip color adding with colorig & orig alpha */
694         if((vp->flag & VP_SPRAY)==0) {
695                 unsigned int testcol=0, a;
696                 char *cp, *ct, *co;
697                 
698                 alpha= (int)(255.0f*brush_alpha(brush));
699                 
700                 if(brush->vertexpaint_tool==VP_MIX || brush->vertexpaint_tool==VP_BLUR) testcol= mcol_blend( *colorig, paintcol, alpha);
701                 else if(brush->vertexpaint_tool==VP_ADD) testcol= mcol_add( *colorig, paintcol, alpha);
702                 else if(brush->vertexpaint_tool==VP_SUB) testcol= mcol_sub( *colorig, paintcol, alpha);
703                 else if(brush->vertexpaint_tool==VP_MUL) testcol= mcol_mul( *colorig, paintcol, alpha);
704                 else if(brush->vertexpaint_tool==VP_LIGHTEN)  testcol= mcol_lighten( *colorig, paintcol, alpha);
705                 else if(brush->vertexpaint_tool==VP_DARKEN)   testcol= mcol_darken( *colorig, paintcol, alpha);
706                 
707                 cp= (char *)col;
708                 ct= (char *)&testcol;
709                 co= (char *)colorig;
710                 
711                 for(a=0; a<4; a++) {
712                         if( ct[a]<co[a] ) {
713                                 if( cp[a]<ct[a] ) cp[a]= ct[a];
714                                 else if( cp[a]>co[a] ) cp[a]= co[a];
715                         }
716                         else {
717                                 if( cp[a]<co[a] ) cp[a]= co[a];
718                                 else if( cp[a]>ct[a] ) cp[a]= ct[a];
719                         }
720                 }
721         }
722 }
723
724
725 static int sample_backbuf_area(ViewContext *vc, int *indexar, int totface, int x, int y, float size)
726 {
727         struct ImBuf *ibuf;
728         int a, tot=0, index;
729         
730         /* brecht: disabled this because it obviously failes for
731            brushes with size > 64, why is this here? */
732         /*if(size>64.0) size= 64.0;*/
733         
734         ibuf= view3d_read_backbuf(vc, x-size, y-size, x+size, y+size);
735         if(ibuf) {
736                 unsigned int *rt= ibuf->rect;
737
738                 memset(indexar, 0, sizeof(int)*(totface+1));
739                 
740                 size= ibuf->x*ibuf->y;
741                 while(size--) {
742                                 
743                         if(*rt) {
744                                 index= WM_framebuffer_to_index(*rt);
745                                 if(index>0 && index<=totface)
746                                         indexar[index] = 1;
747                         }
748                 
749                         rt++;
750                 }
751                 
752                 for(a=1; a<=totface; a++) {
753                         if(indexar[a]) indexar[tot++]= a;
754                 }
755
756                 IMB_freeImBuf(ibuf);
757         }
758         
759         return tot;
760 }
761
762 static float calc_vp_alpha_dl(VPaint *vp, ViewContext *vc, float vpimat[][3], float *vert_nor, const float mval[2], float pressure)
763 {
764         Brush *brush = paint_brush(&vp->paint);
765         float fac, fac_2, size, dx, dy;
766         float alpha;
767         int vertco[2];
768         const int radius= brush_size(brush);
769
770         project_int_noclip(vc->ar, vert_nor, vertco);
771         dx= mval[0]-vertco[0];
772         dy= mval[1]-vertco[1];
773         
774         if (brush_use_size_pressure(brush))
775                 size = pressure * radius;
776         else
777                 size = radius;
778         
779         fac_2= dx*dx + dy*dy;
780         if(fac_2 > size*size) return 0.f;
781         fac = sqrtf(fac_2);
782         
783         alpha= brush_alpha(brush) * brush_curve_strength_clamp(brush, fac, size);
784         
785         if (brush_use_alpha_pressure(brush))
786                 alpha *= pressure;
787                 
788         if(vp->flag & VP_NORMALS) {
789                 float *no= vert_nor+3;
790                 
791                 /* transpose ! */
792                 fac= vpimat[2][0]*no[0]+vpimat[2][1]*no[1]+vpimat[2][2]*no[2];
793                 if(fac > 0.0f) {
794                         dx= vpimat[0][0]*no[0]+vpimat[0][1]*no[1]+vpimat[0][2]*no[2];
795                         dy= vpimat[1][0]*no[0]+vpimat[1][1]*no[1]+vpimat[1][2]*no[2];
796                         
797                         alpha*= fac/sqrtf(dx*dx + dy*dy + fac*fac);
798                 }
799                 else return 0.f;
800         }
801         
802         return alpha;
803 }
804
805 static void wpaint_blend(VPaint *wp, MDeformWeight *dw, MDeformWeight *uw, float alpha, float paintval, int flip, int multipaint)
806 {
807         Brush *brush = paint_brush(&wp->paint);
808         int tool = brush->vertexpaint_tool;
809         
810         if(dw==NULL || uw==NULL) return;
811         
812         if (flip) {
813                 switch(tool) {
814                         case VP_MIX:
815                                 paintval = 1.f - paintval; break;
816                         case VP_ADD:
817                                 tool= VP_SUB; break;
818                         case VP_SUB:
819                                 tool= VP_ADD; break;
820                         case VP_LIGHTEN:
821                                 tool= VP_DARKEN; break;
822                         case VP_DARKEN:
823                                 tool= VP_LIGHTEN; break;
824                 }
825         }
826         
827         if(tool==VP_MIX || tool==VP_BLUR)
828                 dw->weight = paintval*alpha + dw->weight*(1.0f-alpha);
829         else if(tool==VP_ADD)
830                 dw->weight += paintval*alpha;
831         else if(tool==VP_SUB) 
832                 dw->weight -= paintval*alpha;
833         else if(tool==VP_MUL) 
834                 /* first mul, then blend the fac */
835                 dw->weight = ((1.0f-alpha) + alpha*paintval)*dw->weight;
836         else if(tool==VP_LIGHTEN) {
837                 if (dw->weight < paintval)
838                         dw->weight = paintval*alpha + dw->weight*(1.0f-alpha);
839         } else if(tool==VP_DARKEN) {
840                 if (dw->weight > paintval)
841                         dw->weight = paintval*alpha + dw->weight*(1.0f-alpha);
842         }
843         /* Radish delay clamping until the end so multi-paint can function when the active group is at the limits */
844         //CLAMP(dw->weight, 0.0f, 1.0f);
845         
846         /* if no spray, clip result with orig weight & orig alpha */
847         if((wp->flag & VP_SPRAY)==0) {
848                 float testw=0.0f;
849                 
850                 alpha= brush_alpha(brush);
851                 if(tool==VP_MIX || tool==VP_BLUR)
852                         testw = paintval*alpha + uw->weight*(1.0f-alpha);
853                 else if(tool==VP_ADD)
854                         testw = uw->weight + paintval*alpha;
855                 else if(tool==VP_SUB) 
856                         testw = uw->weight - paintval*alpha;
857                 else if(tool==VP_MUL) 
858                         /* first mul, then blend the fac */
859                         testw = ((1.0f-alpha) + alpha*paintval)*uw->weight;
860                 else if(tool==VP_LIGHTEN) {
861                         if (uw->weight < paintval)
862                                 testw = paintval*alpha + uw->weight*(1.0f-alpha);
863                         else
864                                 testw = uw->weight;
865                 } else if(tool==VP_DARKEN) {
866                         if (uw->weight > paintval)
867                                 testw = paintval*alpha + uw->weight*(1.0f-alpha);
868                         else
869                                 testw = uw->weight;
870                 }
871                 //CLAMP(testw, 0.0f, 1.0f);
872                 if(!multipaint) {
873                         if( testw<uw->weight ) {
874                                 if(dw->weight < testw) dw->weight= testw;
875                                 else if(dw->weight > uw->weight) dw->weight= uw->weight;
876                         }
877                         else {
878                                 if(dw->weight > testw) dw->weight= testw;
879                                 else if(dw->weight < uw->weight) dw->weight= uw->weight;
880                         }
881                 }
882         }
883         
884 }
885
886 /* ----------------------------------------------------- */
887
888
889 /* sets wp->weight to the closest weight value to vertex */
890 /* note: we cant sample frontbuf, weight colors are interpolated too unpredictable */
891 static int weight_sample_invoke(bContext *C, wmOperator *op, wmEvent *event)
892 {
893         ViewContext vc;
894         Mesh *me;
895         short change= FALSE;
896
897         view3d_set_viewcontext(C, &vc);
898         me= get_mesh(vc.obact);
899
900         if (me && me->dvert && vc.v3d && vc.rv3d) {
901                 int index;
902
903                 view3d_operator_needs_opengl(C);
904
905                 index= view3d_sample_backbuf(&vc, event->mval[0], event->mval[1]);
906
907                 if(index && index<=me->totface) {
908                         DerivedMesh *dm= mesh_get_derived_final(vc.scene, vc.obact, CD_MASK_BAREMESH);
909
910                         if(dm->getVertCo==NULL) {
911                                 BKE_report(op->reports, RPT_WARNING, "The modifier used does not support deformed locations");
912                         }
913                         else {
914                                 MFace *mf= ((MFace *)me->mface) + index-1;
915                                 const int vgroup= vc.obact->actdef - 1;
916                                 ToolSettings *ts= vc.scene->toolsettings;
917                                 float mval_f[2];
918                                 int v_idx_best= -1;
919                                 int fidx;
920                                 float len_best= FLT_MAX;
921
922                                 mval_f[0]= (float)event->mval[0];
923                                 mval_f[1]= (float)event->mval[1];
924
925                                 fidx= mf->v4 ? 3:2;
926                                 do {
927                                         float co[3], sco[3], len;
928                                         const int v_idx= (*(&mf->v1 + fidx));
929                                         dm->getVertCo(dm, v_idx, co);
930                                         project_float_noclip(vc.ar, co, sco);
931                                         len= len_squared_v2v2(mval_f, sco);
932                                         if(len < len_best) {
933                                                 len_best= len;
934                                                 v_idx_best= v_idx;
935                                         }
936                                 } while (fidx--);
937
938                                 if(v_idx_best != -1) { /* should always be valid */
939                                         ts->vgroup_weight= defvert_find_weight(&me->dvert[v_idx_best], vgroup);
940                                         change= TRUE;
941                                 }
942                         }
943                         dm->release(dm);
944                 }
945         }
946
947         if(change) {
948                 /* not really correct since the brush didnt change, but redraws the toolbar */
949                 WM_main_add_notifier(NC_BRUSH|NA_EDITED, NULL); /* ts->wpaint->paint.brush */
950
951                 return OPERATOR_FINISHED;
952         }
953         else {
954                 return OPERATOR_CANCELLED;
955         }
956 }
957
958 void PAINT_OT_weight_sample(wmOperatorType *ot)
959 {
960         /* identifiers */
961         ot->name= "Weight Paint Sample Weight";
962         ot->idname= "PAINT_OT_weight_sample";
963
964         /* api callbacks */
965         ot->invoke= weight_sample_invoke;
966         ot->poll= weight_paint_mode_poll;
967
968         /* flags */
969         ot->flag= OPTYPE_UNDO;
970 }
971
972 /* samples cursor location, and gives menu with vertex groups to activate */
973 static EnumPropertyItem *weight_paint_sample_enum_itemf(bContext *C, PointerRNA *UNUSED(ptr), PropertyRNA *UNUSED(prop), int *free)
974 {
975         if (C) {
976                 wmWindow *win= CTX_wm_window(C);
977                 if(win && win->eventstate) {
978                         ViewContext vc;
979                         Mesh *me;
980
981                         view3d_set_viewcontext(C, &vc);
982                         me= get_mesh(vc.obact);
983
984                         if (me && me->dvert && vc.v3d && vc.rv3d) {
985                                 int index;
986
987                                 view3d_operator_needs_opengl(C);
988
989                                 index= view3d_sample_backbuf(&vc, win->eventstate->x - vc.ar->winrct.xmin, win->eventstate->y - vc.ar->winrct.ymin);
990
991                                 if(index && index<=me->totface) {
992                                         const int totgroup= BLI_countlist(&vc.obact->defbase);
993                                         if(totgroup) {
994                                                 MFace *mf= ((MFace *)me->mface) + index-1;
995                                                 int fidx= mf->v4 ? 3:2;
996                                                 int *groups= MEM_callocN(totgroup*sizeof(int), "groups");
997                                                 int found= FALSE;
998
999                                                 do {
1000                                                         MDeformVert *dvert= me->dvert + (*(&mf->v1 + fidx));
1001                                                         int i= dvert->totweight;
1002                                                         MDeformWeight *dw;
1003                                                         for(dw= dvert->dw; i > 0; dw++, i--) {
1004                                                                 groups[dw->def_nr]= TRUE;
1005                                                                 found= TRUE;
1006                                                         }
1007                                                 } while (fidx--);
1008
1009                                                 if(found==FALSE) {
1010                                                         MEM_freeN(groups);
1011                                                 }
1012                                                 else {
1013                                                         EnumPropertyItem *item= NULL, item_tmp= {0};
1014                                                         int totitem= 0;
1015                                                         int i= 0;
1016                                                         bDeformGroup *dg;
1017                                                         for(dg= vc.obact->defbase.first; dg && i<totgroup; i++, dg= dg->next) {
1018                                                                 if(groups[i]) {
1019                                                                         item_tmp.identifier= item_tmp.name= dg->name;
1020                                                                         item_tmp.value= i;
1021                                                                         RNA_enum_item_add(&item, &totitem, &item_tmp);
1022                                                                 }
1023                                                         }
1024
1025                                                         RNA_enum_item_end(&item, &totitem);
1026                                                         *free= 1;
1027
1028                                                         MEM_freeN(groups);
1029                                                         return item;
1030                                                 }
1031                                         }
1032                                 }
1033                         }
1034                 }
1035         }
1036
1037         return DummyRNA_NULL_items;
1038 }
1039
1040 static int weight_sample_group_exec(bContext *C, wmOperator *op)
1041 {
1042         int type= RNA_enum_get(op->ptr, "group");
1043         ViewContext vc;
1044         view3d_set_viewcontext(C, &vc);
1045
1046         vc.obact->actdef= type + 1;
1047
1048         DAG_id_tag_update(&vc.obact->id, OB_RECALC_DATA);
1049         WM_event_add_notifier(C, NC_OBJECT|ND_DRAW, vc.obact);
1050         return OPERATOR_FINISHED;
1051 }
1052
1053 /* TODO, we could make this a menu into OBJECT_OT_vertex_group_set_active rather than its own operator */
1054 void PAINT_OT_weight_sample_group(wmOperatorType *ot)
1055 {
1056         PropertyRNA *prop= NULL;
1057
1058         /* identifiers */
1059         ot->name= "Weight Paint Sample Group";
1060         ot->idname= "PAINT_OT_weight_sample_group";
1061
1062         /* api callbacks */
1063         ot->exec= weight_sample_group_exec;
1064         ot->invoke= WM_menu_invoke;
1065         ot->poll= weight_paint_mode_poll;
1066
1067         /* flags */
1068         ot->flag= OPTYPE_UNDO;
1069
1070         /* keyingset to use (dynamic enum) */
1071         prop= RNA_def_enum(ot->srna, "group", DummyRNA_DEFAULT_items, 0, "Keying Set", "The Keying Set to use");
1072         RNA_def_enum_funcs(prop, weight_paint_sample_enum_itemf);
1073         ot->prop= prop;
1074 }
1075
1076
1077 #if 0 /* UNUSED */
1078 static void do_weight_paint_auto_normalize(MDeformVert *dvert, 
1079                                            int paint_nr, char *map)
1080 {
1081 //      MDeformWeight *dw = dvert->dw;
1082         float sum=0.0f, fac=0.0f, paintw=0.0f;
1083         int i, tot=0;
1084
1085         if (!map)
1086                 return;
1087
1088         for (i=0; i<dvert->totweight; i++) {
1089                 if (dvert->dw[i].def_nr == paint_nr)
1090                         paintw = dvert->dw[i].weight;
1091
1092                 if (map[dvert->dw[i].def_nr]) {
1093                         tot += 1;
1094                         if (dvert->dw[i].def_nr != paint_nr)
1095                                 sum += dvert->dw[i].weight;
1096                 }
1097         }
1098         
1099         if (!tot || sum <= (1.0f - paintw))
1100                 return;
1101
1102         fac = sum / (1.0f - paintw);
1103         fac = fac==0.0f ? 1.0f : 1.0f / fac;
1104
1105         for (i=0; i<dvert->totweight; i++) {
1106                 if (map[dvert->dw[i].def_nr]) {
1107                         if (dvert->dw[i].def_nr != paint_nr)
1108                                 dvert->dw[i].weight *= fac;
1109                 }
1110         }
1111 }
1112 #endif
1113
1114 /* Radish: the active group should be involved in auto normalize */
1115 static void do_weight_paint_auto_normalize_all_groups(MDeformVert *dvert, char *map)
1116 {
1117 //      MDeformWeight *dw = dvert->dw;
1118         float sum=0.0f, fac=0.0f;
1119         int i, tot=0;
1120
1121         if (!map)
1122                 return;
1123
1124         for (i=0; i<dvert->totweight; i++) {
1125                 if (map[dvert->dw[i].def_nr]) {
1126                         tot += 1;
1127                         sum += dvert->dw[i].weight;
1128                 }
1129         }
1130         
1131         if (!tot || sum == 1.0f)
1132                 return;
1133
1134         fac = sum;
1135         fac = fac==0.0f ? 1.0f : 1.0f / fac;
1136
1137         for (i=0; i<dvert->totweight; i++) {
1138                 if (map[dvert->dw[i].def_nr]) {
1139                         dvert->dw[i].weight *= fac;
1140                 }
1141         }
1142 }
1143 /* Radish */
1144 /*
1145 See if the current deform vertex has a locked group
1146 */
1147 static char has_locked_group(MDeformVert *dvert, char *flags)
1148 {
1149         int i;
1150         for(i = 0; i < dvert->totweight; i++) {
1151                 if(flags[(dvert->dw+i)->def_nr] && (dvert->dw+i)->weight > 0.0f) {
1152                         return TRUE;
1153                 }
1154         }
1155         return FALSE;
1156 }
1157 /* Radish
1158  * gen_lck_flags gets the status of "flag" for each bDeformGroup
1159  *in ob->defbase and returns an array containing them
1160  */
1161 static char *gen_lck_flags(Object* ob, int defbase_len)
1162 {
1163         char is_locked = FALSE;
1164         int i;
1165         //int defbase_len = BLI_countlist(&ob->defbase);
1166         char *flags = MEM_mallocN(defbase_len*sizeof(char), "defflags");
1167         bDeformGroup *defgroup;
1168
1169         for(i = 0, defgroup = ob->defbase.first; i < defbase_len && defgroup; defgroup = defgroup->next, i++) {
1170                 flags[i] = ((defgroup->flag & DG_LOCK_WEIGHT) != 0);
1171                 is_locked |= flags[i];
1172         }
1173         if(is_locked){
1174                 return flags;
1175         }
1176         // don't forget to free it if it is unneeded
1177         MEM_freeN(flags);
1178         return NULL;
1179 }
1180 /* Radish */
1181 static int has_locked_group_selected(int defbase_len, char *defbase_sel, char *flags) {
1182         int i;
1183         for(i = 0; i < defbase_len; i++) {
1184                 if(defbase_sel[i] && flags[i]) {
1185                         return TRUE;
1186                 }
1187         }
1188         return FALSE;
1189 }
1190
1191 /* Radish */
1192 #if 0 /* UNUSED */
1193 static int has_unselected_unlocked_bone_group(int defbase_len, char *defbase_sel, int selected, char *flags, char *bone_groups) {
1194         int i;
1195         if(defbase_len == selected) {
1196                 return FALSE;
1197         }
1198         for(i = 0; i < defbase_len; i++) {
1199                 if(bone_groups[i] && !defbase_sel[i] && !flags[i]) {
1200                         return TRUE;
1201                 }
1202         }
1203         return FALSE;
1204 }
1205 #endif
1206
1207 /* Radish */
1208 static void multipaint_selection(MDeformVert *dvert, float change, char *defbase_sel, int defbase_len) {
1209         int i;
1210         MDeformWeight *dw;
1211         float val;
1212         // make sure they are all at most 1 after the change
1213         for(i = 0; i < defbase_len; i++) {
1214                 if(defbase_sel[i]) {
1215                         dw = defvert_find_index(dvert, i);
1216                         if(dw && dw->weight) {
1217                                 val = dw->weight * change;
1218                                 if(val > 1) {
1219                                         /* Radish TODO: when the change is reduced, you need to recheck the earlier values to make sure they are not 0 (precision error) */
1220                                         change = 1.0f/dw->weight;
1221                                 }
1222                                 // the value should never reach zero while multi-painting if it was nonzero beforehand
1223                                 if(val <= 0) {
1224                                         return;
1225                                 }
1226                         }
1227                 }
1228         }
1229         // apply the valid change
1230         for(i = 0; i < defbase_len; i++) {
1231                 if(defbase_sel[i]) {
1232                         dw = defvert_find_index(dvert, i);
1233                         if(dw && dw->weight) {
1234                                 dw->weight = dw->weight * change;
1235                         }
1236                 }
1237         }
1238 }
1239 /* Radish */
1240 // move all change onto valid, unchanged groups.  If there is change left over, then return it.
1241 // assumes there are valid groups to shift weight onto
1242 static float redistribute_change(MDeformVert *ndv, char *change_status, int changeme, int changeto, char *validmap, float totchange, float total_valid) {
1243         float was_change;
1244         float change;
1245         float oldval;
1246         MDeformWeight *ndw;
1247         int i;
1248         do {
1249                 // assume there is no change until you see one
1250                 was_change = FALSE;
1251                 // change each group by the same amount each time
1252                 change = totchange/total_valid;
1253                 for(i = 0; i < ndv->totweight && total_valid && totchange; i++) {
1254                         ndw = (ndv->dw+i);
1255                         // change only the groups with a valid status
1256                         if(change_status[ndw->def_nr] == changeme) {
1257                                 oldval = ndw->weight;
1258                                 // if auto normalize is active, don't worry about upper bounds
1259                                 if(!validmap && ndw->weight + change > 1) {
1260                                         totchange -= 1-ndw->weight;
1261                                         ndw->weight = 1;
1262                                         // stop the changes to this group
1263                                         change_status[ndw->def_nr] = changeto;
1264                                         total_valid--;
1265                                 } else if(ndw->weight + change < 0) { // check the lower bound
1266                                         totchange -= ndw->weight;
1267                                         ndw->weight = 0;
1268                                         change_status[ndw->def_nr] = changeto;
1269                                         total_valid--;
1270                                 } else {// a perfectly valid change occurred to ndw->weight
1271                                         totchange -= change;
1272                                         ndw->weight += change;
1273                                 }
1274                                 // see if there was a change
1275                                 if(oldval != ndw->weight) {
1276                                         was_change = TRUE;
1277                                 }
1278                         }
1279                 }
1280         // don't go again if there was no change, if there is no valid group, or there is no change left
1281         }while(was_change && total_valid && totchange);
1282         // left overs
1283         return totchange;
1284 }
1285 /* Radish */
1286 // observe the changes made to the weights of groups.
1287 // make sure all locked groups on the vertex have the same deformation
1288 // by moving the changes made to groups onto other unlocked groups
1289 static void enforce_locks(MDeformVert *odv, MDeformVert *ndv, int defbase_len, char *flags, char *bone_groups, char *validmap) {
1290         float totchange = 0.0f;
1291         float totchange_allowed = 0.0f;
1292         float left_over;
1293
1294         int total_valid = 0;
1295         int total_changed = 0;
1296         int i;
1297         MDeformWeight *ndw;
1298         MDeformWeight *odw;
1299         MDeformWeight *ndw2;
1300         MDeformWeight *odw2;
1301         int designatedw = -1;
1302         int designatedw_changed = FALSE;
1303         float storedw;
1304         char *change_status;
1305         char new_weight_has_zero = FALSE;
1306
1307         if(!flags || !has_locked_group(ndv, flags)) {
1308                 return;
1309         }
1310         // record if a group was changed, unlocked and not changed, or locked
1311         change_status = MEM_callocN(sizeof(char)*defbase_len, "unlocked_unchanged");
1312
1313         for(i = 0; i < defbase_len; i++) {
1314                 ndw = defvert_find_index(ndv, i);
1315                 odw = defvert_find_index(odv, i);
1316                 // the weights are zero, so we can assume a lot
1317                 if(!ndw || !odw) {
1318                         if (!flags[i] && bone_groups[i]){
1319                                 defvert_verify_index(odv, i);
1320                                 defvert_verify_index(ndv, i);
1321                                 total_valid++;
1322                                 change_status[i] = 1; // can be altered while redistributing
1323                         }
1324                         continue;
1325                 }
1326                 // locked groups should not be changed
1327                 if(flags[i]) {
1328                         ndw->weight = odw->weight;
1329                 } else if(ndw->weight != odw->weight) { // changed groups are handled here
1330                         totchange += ndw->weight-odw->weight;
1331                         change_status[i] = 2; // was altered already
1332                         total_changed++;
1333                         if(ndw->weight == 0) {
1334                                 new_weight_has_zero = TRUE;
1335                         } else if(designatedw == -1){
1336                                 designatedw = i;
1337                         }
1338                 } // unchanged, unlocked bone groups are handled here
1339                 else if (bone_groups[i]){
1340                         totchange_allowed += ndw->weight;
1341                         total_valid++;
1342                         change_status[i] = 1; // can be altered while redistributing
1343                 }
1344         }
1345         // if there was any change, redistribute it
1346         if(total_changed) {
1347                 // auto normalize will allow weights to temporarily go above 1 in redistribution
1348                 if(validmap && total_changed < 0 && total_valid) {
1349                         totchange_allowed = total_valid;
1350                 }
1351                 // there needs to be change allowed, or you should not bother
1352                 if(totchange_allowed) {
1353                         // the way you modify the unlocked+unchanged groups is different depending
1354                         // on whether or not you are painting the weight(s) up or down
1355                         if(totchange < 0) {
1356                                 totchange_allowed = total_valid - totchange_allowed;
1357                         } else {
1358                                 totchange_allowed *= -1;
1359                         }
1360                         left_over = 0;
1361                         if(fabs(totchange_allowed) < fabs(totchange)) {
1362                                 // this amount goes back onto the changed, unlocked weights
1363                                 left_over = fabs(fabs(totchange)-fabs(totchange_allowed));
1364                                 if(totchange > 0) {
1365                                         left_over *= -1;
1366                                 }
1367                         }else {
1368                                 // all of the change will be permitted
1369                                 totchange_allowed = -totchange;
1370                         }
1371                         // move the weight evenly between the allowed groups, move excess back onto the used groups based on the change
1372                         totchange_allowed = redistribute_change(ndv, change_status, 1, -1, validmap, totchange_allowed, total_valid);
1373                         left_over += totchange_allowed;
1374                         if(left_over) {
1375                                 // more than one nonzero weights were changed with the same ratio, so keep them changed that way!
1376                                 if(total_changed > 1 && !new_weight_has_zero && designatedw >= 0) {
1377                                         // this dw is special, it is used as a base to determine how to change the others
1378                                         ndw = defvert_find_index(ndv, designatedw);
1379                                         odw = defvert_find_index(odv, designatedw);
1380                                         storedw = ndw->weight;
1381                                         for(i = 0; i < ndv->totweight; i++) {
1382                                                 if(change_status[ndw->def_nr] == 2) {
1383                                                         odw2 = (odv->dw+i);
1384                                                         ndw2 = (ndv->dw+i);
1385                                                         if(!designatedw_changed) {
1386                                                                 ndw->weight = (totchange_allowed + odw->weight + odw2->weight)/(1 + ndw2->weight/ndw->weight);
1387                                                                 designatedw_changed = TRUE;
1388                                                         }
1389                                                         ndw2->weight = ndw->weight*ndw2->weight/storedw;
1390                                                 }
1391                                         }
1392                                 }
1393                                 // a weight was changed to zero, only one weight was changed, or designatedw is still -1
1394                                 // put weight back as evenly as possible
1395                                 else {
1396                                         redistribute_change(ndv, change_status, 2, -2, validmap, left_over, total_changed);
1397                                 }
1398                         }
1399                 } else {
1400                         // reset the weights
1401                         for(i = 0; i < ndv->totweight; i++) {
1402                                 (ndv->dw+i)->weight = (odv->dw+i)->weight;
1403                         }
1404                 }
1405         }
1406
1407         MEM_freeN(change_status);
1408 }
1409 /* Radish */
1410 // multi-paint's initial, potential change is computed here based on the user's stroke
1411 static float get_mp_change(MDeformVert *odv, char *defbase_sel, float brush_change) {
1412         float selwsum = 0.0f;
1413         int i;
1414         MDeformWeight *dw;
1415         for(i=0; i < odv->totweight; i++) {
1416                 if(defbase_sel[(dw = (odv->dw+i))->def_nr]) {
1417                         selwsum += dw->weight;
1418                 }
1419         }
1420         if(selwsum && selwsum+brush_change > 0) {
1421                 return (selwsum+brush_change)/selwsum;
1422         }
1423         return 0.0f;
1424 }
1425 /* Radish */
1426 // change the weights back to the wv's weights
1427 // it assumes you already have the correct pointer index
1428 static void reset_to_prev(MDeformVert *wv, MDeformVert *dv) {
1429         int i;
1430         MDeformWeight *d;
1431         MDeformWeight *w;
1432         for(i = 0; i < dv->totweight; i++) {
1433                 d = dv->dw+i;
1434                 w = defvert_find_index(wv, d->def_nr);
1435                 // if there was no w when there is a d, then the old weight was 0
1436                 if(w) {
1437                         d->weight = w->weight;
1438                 } else {
1439                         d->weight = 0;
1440                 }
1441         }
1442 }
1443 /* Radish */
1444 static void clamp_weights(MDeformVert *dvert) {
1445         int i;
1446         for (i = 0; i < dvert->totweight; i++) {
1447                 CLAMP((dvert->dw+i)->weight, 0.0f, 1.0f);
1448         }
1449 }
1450 /* Radish */
1451 /* fresh start to make multi-paint and locking modular */
1452 /* returns TRUE if it thinks you need to reset the weights due to normalizing while multi-painting */
1453 static int apply_mp_lcks_normalize(Mesh *me, int index, MDeformWeight *dw, MDeformWeight *tdw, int defbase_len, float change, float oldChange, float oldw, float neww, char *defbase_sel, int selected, char *bone_groups, char *validmap, char *flags, int multipaint) {
1454         MDeformVert *dvert = me->dvert+index;
1455         MDeformVert dv= {NULL};
1456
1457         dv.dw= MEM_dupallocN(dvert->dw);
1458         dv.flag = dvert->flag;
1459         dv.totweight = dvert->totweight;
1460         // do not multi-paint if a locked group is selected or the active group is locked
1461         // !flags[dw->def_nr] helps if nothing is selected, but active group is locked
1462         if((flags == NULL) || (has_locked_group_selected(defbase_len, defbase_sel, flags) == FALSE && flags[dw->def_nr] == FALSE)) {
1463                 if(multipaint && selected > 1) {
1464                         if(change && change!=1) {
1465                                 multipaint_selection(dvert, change, defbase_sel, defbase_len);
1466                         }
1467                 } else {// this lets users paint normally, but don't let them paint locked groups
1468                         dw->weight = neww;
1469                 }
1470         }
1471         clamp_weights(dvert);
1472
1473         enforce_locks(&dv, dvert, defbase_len, flags, bone_groups, validmap);
1474
1475         do_weight_paint_auto_normalize_all_groups(dvert, validmap);
1476
1477         if(oldChange && multipaint && selected > 1) {
1478                 if(tdw->weight != oldw) {
1479                         if( neww > oldw ) {
1480                                 if(tdw->weight <= oldw) {
1481                                         MEM_freeN(dv.dw);
1482                                         return TRUE;
1483                                 }
1484                         } else {
1485                                 if(tdw->weight >= oldw) {
1486                                         MEM_freeN(dv.dw);
1487                                         return TRUE;
1488                                 }
1489                         }
1490                 }
1491         }
1492         MEM_freeN(dv.dw);
1493         return FALSE;
1494 }
1495
1496 // within the current dvert index, get the dw that is selected and has a weight above 0
1497 // this helps multi-paint
1498 static int get_first_selected_nonzero_weight(MDeformVert *dvert, char *defbase_sel) {
1499         int i;
1500         MDeformWeight *dw;
1501         for(i=0; i< dvert->totweight; i++) {
1502                 dw = dvert->dw+i;
1503                 if(defbase_sel[dw->def_nr] && dw->weight > 0) {
1504                         return i;
1505                 }
1506         }
1507         return -1;
1508 }
1509 /* Radish */
1510 static char *wpaint_make_validmap(Object *ob);
1511
1512 static void do_weight_paint_vertex(VPaint *wp, Object *ob, int index, 
1513                                    float alpha, float paintweight, int flip, 
1514                                    int vgroup_mirror, char *validmap, int multipaint)
1515 {
1516         Mesh *me= ob->data;
1517         //                                              Radish: tdw, tuw
1518         MDeformWeight *dw, *uw, *tdw = NULL, *tuw;
1519         int vgroup= ob->actdef-1;
1520         
1521         /* Radish */
1522         char *flags;
1523         char *bone_groups;
1524         char *defbase_sel;
1525         int selected;
1526         float oldw;
1527         float neww;
1528         float testw=0;
1529         int defbase_len;
1530         float change = 0;
1531         float oldChange = 0;
1532         int i;
1533         MDeformVert dv= {NULL};
1534
1535         // Need to know which groups are bone groups
1536         if(validmap) {
1537                 bone_groups = validmap;
1538         }else {
1539                 bone_groups = wpaint_make_validmap(ob);
1540         }
1541
1542         if(wp->flag & VP_ONLYVGROUP) {
1543                 dw= defvert_find_index(me->dvert+index, vgroup);
1544                 uw= defvert_find_index(wp->wpaint_prev+index, vgroup);
1545         }
1546         else {
1547                 dw= defvert_verify_index(me->dvert+index, vgroup);
1548                 uw= defvert_verify_index(wp->wpaint_prev+index, vgroup);
1549         }
1550         if(dw==NULL || uw==NULL)
1551                 return;
1552         /* Radish */
1553         flags = gen_lck_flags(ob, defbase_len = BLI_countlist(&ob->defbase));
1554         defbase_sel = MEM_mallocN(defbase_len * sizeof(char), "dg_selected_flags");
1555         selected = get_selected_defgroups(ob, defbase_sel, defbase_len);
1556         if(!selected && ob->actdef) {
1557                 selected = 1;
1558         }
1559         
1560         oldw = dw->weight;
1561         wpaint_blend(wp, dw, uw, alpha, paintweight, flip, multipaint && selected >1);
1562         neww = dw->weight;
1563         dw->weight = oldw;
1564         
1565         // setup multi-paint
1566         if(selected > 1 && multipaint) {
1567                 dv.dw= MEM_dupallocN((me->dvert+index)->dw);
1568                 dv.flag = me->dvert[index].flag;
1569                 dv.totweight = (me->dvert+index)->totweight;
1570                 tdw = dw;
1571                 tuw = uw;
1572                 change = get_mp_change(wp->wpaint_prev+index, defbase_sel, neww-oldw);
1573                 if(change) {
1574                         if(!tdw->weight) {
1575                                 i = get_first_selected_nonzero_weight(me->dvert+index, defbase_sel);
1576                                 if(i>=0) {
1577                                         tdw = ((me->dvert+index)->dw+i);
1578                                         tuw = defvert_verify_index(wp->wpaint_prev+index, tdw->def_nr);
1579                                 } else {
1580                                         change = 0;
1581                                 }
1582                         }
1583                         if(change && tuw->weight && tuw->weight * change) {
1584                                 if(tdw->weight != tuw->weight) {
1585                                         oldChange = tdw->weight/tuw->weight;
1586                                         testw = tuw->weight*change;
1587                                         if( testw > tuw->weight ) {
1588                                                 if(change > oldChange) {
1589                                                         // reset the weights and use the new change
1590                                                         reset_to_prev(wp->wpaint_prev+index, me->dvert+index);
1591                                                 } else {
1592                                                         // the old change was more significant,
1593                                                         // so set the change to 0 so that it will not do another multi-paint
1594                                                         change = 0;
1595                                                 }
1596                                         } else {
1597                                                 if(change < oldChange) {
1598                                                         reset_to_prev(wp->wpaint_prev+index, me->dvert+index);
1599                                                 } else {
1600                                                         change = 0;
1601                                                 }
1602                                         }
1603                                 }
1604                         } else {
1605                                 change = 0;
1606                         }
1607                 }
1608         }
1609         /* Radish */
1610         if(apply_mp_lcks_normalize(me, index, dw, tdw, defbase_len, change, oldChange, oldw, neww, defbase_sel, selected, bone_groups, validmap, flags, multipaint)) {
1611                 reset_to_prev(&dv, me->dvert+index);
1612                 change = 0;
1613                 oldChange = 0;
1614         }
1615         if(dv.dw) {
1616                 MEM_freeN(dv.dw);
1617         }
1618         // dvert may have been altered greatly
1619         dw = defvert_find_index(me->dvert+index, vgroup);
1620
1621         if(me->editflag & ME_EDIT_MIRROR_X) {   /* x mirror painting */
1622                 int j= mesh_get_x_mirror_vert(ob, index);
1623                 if(j>=0) {
1624                         /* copy, not paint again */
1625                         if(vgroup_mirror != -1)
1626                                 uw= defvert_verify_index(me->dvert+j, vgroup_mirror);
1627                         else
1628                                 uw= defvert_verify_index(me->dvert+j, vgroup);
1629                         /* Radish */
1630                         //uw->weight= dw->weight;
1631                         /* Radish */
1632                         apply_mp_lcks_normalize(me, j, uw, tdw, defbase_len, change, oldChange, oldw, neww, defbase_sel, selected, bone_groups, validmap, flags, multipaint);
1633                 }
1634         }
1635         /* Radish */
1636         if(flags) {
1637                 MEM_freeN(flags);
1638         }
1639         MEM_freeN(defbase_sel);
1640         if(!validmap) {
1641                 MEM_freeN(bone_groups);
1642         }
1643 }
1644
1645
1646 /* *************** set wpaint operator ****************** */
1647
1648 static int set_wpaint(bContext *C, wmOperator *UNUSED(op))              /* toggle */
1649 {               
1650         Object *ob= CTX_data_active_object(C);
1651         Scene *scene= CTX_data_scene(C);
1652         VPaint *wp= scene->toolsettings->wpaint;
1653         Mesh *me;
1654         
1655         me= get_mesh(ob);
1656         if(ob->id.lib || me==NULL) return OPERATOR_PASS_THROUGH;
1657         
1658         if(ob->mode & OB_MODE_WEIGHT_PAINT) ob->mode &= ~OB_MODE_WEIGHT_PAINT;
1659         else ob->mode |= OB_MODE_WEIGHT_PAINT;
1660         
1661         
1662         /* Weightpaint works by overriding colors in mesh,
1663                 * so need to make sure we recalc on enter and
1664                 * exit (exit needs doing regardless because we
1665                                 * should redeform).
1666                 */
1667         DAG_id_tag_update(&me->id, 0);
1668         
1669         if(ob->mode & OB_MODE_WEIGHT_PAINT) {
1670                 Object *par;
1671                 
1672                 if(wp==NULL)
1673                         wp= scene->toolsettings->wpaint= new_vpaint(1);
1674
1675                 paint_init(&wp->paint, PAINT_CURSOR_WEIGHT_PAINT);
1676                 paint_cursor_start(C, weight_paint_poll);
1677                 
1678                 mesh_octree_table(ob, NULL, NULL, 's');
1679                 
1680                 /* verify if active weight group is also active bone */
1681                 par= modifiers_isDeformedByArmature(ob);
1682                 if(par && (par->mode & OB_MODE_POSE)) {
1683                         bArmature *arm= par->data;
1684
1685                         if(arm->act_bone)
1686                                 ED_vgroup_select_by_name(ob, arm->act_bone->name);
1687                 }
1688         }
1689         else {
1690                 mesh_octree_table(NULL, NULL, NULL, 'e');
1691                 mesh_mirrtopo_table(NULL, 'e');
1692         }
1693         
1694         WM_event_add_notifier(C, NC_SCENE|ND_MODE, scene);
1695         
1696         return OPERATOR_FINISHED;
1697 }
1698
1699 /* for switching to/from mode */
1700 static int paint_poll_test(bContext *C)
1701 {
1702         if(CTX_data_edit_object(C))
1703                 return 0;
1704         if(CTX_data_active_object(C)==NULL)
1705                 return 0;
1706         return 1;
1707 }
1708
1709 void PAINT_OT_weight_paint_toggle(wmOperatorType *ot)
1710 {
1711         
1712         /* identifiers */
1713         ot->name= "Weight Paint Mode";
1714         ot->idname= "PAINT_OT_weight_paint_toggle";
1715         
1716         /* api callbacks */
1717         ot->exec= set_wpaint;
1718         ot->poll= paint_poll_test;
1719         
1720         /* flags */
1721         ot->flag= OPTYPE_REGISTER|OPTYPE_UNDO;
1722         
1723 }
1724
1725 /* ************ weight paint operator ********** */
1726
1727 struct WPaintData {
1728         ViewContext vc;
1729         int *indexar;
1730         int vgroup_mirror;
1731         float *vertexcosnos;
1732         float wpimat[3][3];
1733         
1734         /*variables for auto normalize*/
1735         int auto_normalize;
1736         char *vgroup_validmap; /*stores if vgroups tie to deforming bones or not*/
1737 };
1738
1739 static char *wpaint_make_validmap(Object *ob)
1740 {
1741         bDeformGroup *dg;
1742         ModifierData *md;
1743         char *validmap;
1744         bPose *pose;
1745         bPoseChannel *chan;
1746         ArmatureModifierData *amd;
1747         GHash *gh = BLI_ghash_new(BLI_ghashutil_strhash, BLI_ghashutil_strcmp, "wpaint_make_validmap gh");
1748         int i = 0, step1=1;
1749
1750         /*add all names to a hash table*/
1751         for (dg=ob->defbase.first, i=0; dg; dg=dg->next, i++) {
1752                 BLI_ghash_insert(gh, dg->name, NULL);
1753         }
1754
1755         if (!i)
1756                 return NULL;
1757
1758         validmap = MEM_callocN(i, "wpaint valid map");
1759
1760         /*now loop through the armature modifiers and identify deform bones*/
1761         for (md = ob->modifiers.first; md; md= !md->next && step1 ? (step1=0), modifiers_getVirtualModifierList(ob) : md->next) {
1762                 if (!(md->mode & (eModifierMode_Realtime|eModifierMode_Virtual)))
1763                         continue;
1764
1765                 if (md->type == eModifierType_Armature) 
1766                 {
1767                         amd = (ArmatureModifierData*) md;
1768
1769                         if(amd->object && amd->object->pose) {
1770                                 pose = amd->object->pose;
1771                                 
1772                                 for (chan=pose->chanbase.first; chan; chan=chan->next) {
1773                                         if (chan->bone->flag & BONE_NO_DEFORM)
1774                                                 continue;
1775
1776                                         if (BLI_ghash_haskey(gh, chan->name)) {
1777                                                 BLI_ghash_remove(gh, chan->name, NULL, NULL);
1778                                                 BLI_ghash_insert(gh, chan->name, SET_INT_IN_POINTER(1));
1779                                         }
1780                                 }
1781                         }
1782                 }
1783         }
1784         
1785         /*add all names to a hash table*/
1786         for (dg=ob->defbase.first, i=0; dg; dg=dg->next, i++) {
1787                 if (BLI_ghash_lookup(gh, dg->name) != NULL) {
1788                         validmap[i] = 1;
1789                 }
1790         }
1791
1792         BLI_ghash_free(gh, NULL, NULL);
1793
1794         return validmap;
1795 }
1796
1797 static int wpaint_stroke_test_start(bContext *C, wmOperator *op, wmEvent *UNUSED(event))
1798 {
1799         Scene *scene= CTX_data_scene(C);
1800         struct PaintStroke *stroke = op->customdata;
1801         ToolSettings *ts= CTX_data_tool_settings(C);
1802         VPaint *wp= ts->wpaint;
1803         Object *ob= CTX_data_active_object(C);
1804         struct WPaintData *wpd;
1805         Mesh *me;
1806         float mat[4][4], imat[4][4];
1807         
1808         if(scene->obedit) return OPERATOR_CANCELLED;
1809         
1810         me= get_mesh(ob);
1811         if(me==NULL || me->totface==0) return OPERATOR_PASS_THROUGH;
1812         
1813         /* if nothing was added yet, we make dverts and a vertex deform group */
1814         if (!me->dvert) {
1815                 ED_vgroup_data_create(&me->id);
1816                 WM_event_add_notifier(C, NC_GEOM|ND_DATA, me);
1817         }
1818         
1819         /* make mode data storage */
1820         wpd= MEM_callocN(sizeof(struct WPaintData), "WPaintData");
1821         paint_stroke_set_mode_data(stroke, wpd);
1822         view3d_set_viewcontext(C, &wpd->vc);
1823         wpd->vgroup_mirror= -1;
1824         
1825         /*set up auto-normalize, and generate map for detecting which
1826           vgroups affect deform bones*/
1827         wpd->auto_normalize = ts->auto_normalize;
1828         if (wpd->auto_normalize)
1829                 wpd->vgroup_validmap = wpaint_make_validmap(ob);
1830         
1831         //      if(qual & LR_CTRLKEY) {
1832         //              sample_wpaint(scene, ar, v3d, 0);
1833         //              return;
1834         //      }
1835         //      if(qual & LR_SHIFTKEY) {
1836         //              sample_wpaint(scene, ar, v3d, 1);
1837         //              return;
1838         //      }
1839         
1840         /* ALLOCATIONS! no return after this line */
1841         /* painting on subsurfs should give correct points too, this returns me->totvert amount */
1842         wpd->vertexcosnos= mesh_get_mapped_verts_nors(scene, ob);
1843         wpd->indexar= get_indexarray(me);
1844         copy_wpaint_prev(wp, me->dvert, me->totvert);
1845         
1846         /* this happens on a Bone select, when no vgroup existed yet */
1847         if(ob->actdef<=0) {
1848                 Object *modob;
1849                 if((modob = modifiers_isDeformedByArmature(ob))) {
1850                         Bone *actbone= ((bArmature *)modob->data)->act_bone;
1851                         if(actbone) {
1852                                 bPoseChannel *pchan= get_pose_channel(modob->pose, actbone->name);
1853
1854                                 if(pchan) {
1855                                         bDeformGroup *dg= defgroup_find_name(ob, pchan->name);
1856                                         if(dg==NULL)
1857                                                 dg= ED_vgroup_add_name(ob, pchan->name);        /* sets actdef */
1858                                         else
1859                                                 ob->actdef= 1 + defgroup_find_index(ob, dg);
1860                                 }
1861                         }
1862                 }
1863         }
1864         if(ob->defbase.first==NULL) {
1865                 ED_vgroup_add(ob);
1866         }
1867         
1868         //      if(ob->lay & v3d->lay); else error("Active object is not in this layer");
1869         
1870         /* imat for normals */
1871         mul_m4_m4m4(mat, ob->obmat, wpd->vc.rv3d->viewmat);
1872         invert_m4_m4(imat, mat);
1873         copy_m3_m4(wpd->wpimat, imat);
1874         
1875         /* if mirror painting, find the other group */
1876         if(me->editflag & ME_EDIT_MIRROR_X) {
1877                 wpaint_mirror_vgroup_ensure(ob, &wpd->vgroup_mirror);
1878         }
1879         
1880         return 1;
1881 }
1882
1883 static void wpaint_stroke_update_step(bContext *C, struct PaintStroke *stroke, PointerRNA *itemptr)
1884 {
1885         ToolSettings *ts= CTX_data_tool_settings(C);
1886         VPaint *wp= ts->wpaint;
1887         Brush *brush = paint_brush(&wp->paint);
1888         struct WPaintData *wpd= paint_stroke_mode_data(stroke);
1889         ViewContext *vc;
1890         Object *ob;
1891         Mesh *me;
1892         float mat[4][4];
1893         float paintweight;
1894         int *indexar;
1895         int totindex, index, totw, flip;
1896         float alpha;
1897         float mval[2], pressure;
1898         
1899         /* Radish */
1900         int use_vert_sel;
1901
1902         /* cannot paint if there is no stroke data */
1903         if (wpd == NULL) {
1904                 // XXX: force a redraw here, since even though we can't paint, 
1905                 // at least view won't freeze until stroke ends
1906                 ED_region_tag_redraw(CTX_wm_region(C));
1907                 return;
1908         }
1909                 
1910         vc= &wpd->vc;
1911         ob= vc->obact;
1912         me= ob->data;
1913         indexar= wpd->indexar;
1914         
1915         view3d_operator_needs_opengl(C);
1916                 
1917         /* load projection matrix */
1918         mul_m4_m4m4(mat, ob->obmat, vc->rv3d->persmat);
1919
1920         flip = RNA_boolean_get(itemptr, "pen_flip");
1921         pressure = RNA_float_get(itemptr, "pressure");
1922         RNA_float_get_array(itemptr, "mouse", mval);
1923         mval[0]-= vc->ar->winrct.xmin;
1924         mval[1]-= vc->ar->winrct.ymin;
1925                         
1926         swap_m4m4(wpd->vc.rv3d->persmat, mat);
1927
1928         /* Radish */
1929         use_vert_sel= (me->editflag & ME_EDIT_VERT_SEL) != 0;
1930
1931         /* which faces are involved */
1932         if(wp->flag & VP_AREA) {
1933                 // Ugly hack, to avoid drawing vertex index when getting the face index buffer - campbell
1934                 me->editflag &= ~ME_EDIT_VERT_SEL;
1935                 totindex= sample_backbuf_area(vc, indexar, me->totface, mval[0], mval[1], brush_size(brush));
1936                 me->editflag |= use_vert_sel ? ME_EDIT_VERT_SEL : 0;
1937         }
1938         else {
1939                 indexar[0]= view3d_sample_backbuf(vc, mval[0], mval[1]);
1940                 if(indexar[0]) totindex= 1;
1941                 else totindex= 0;
1942         }
1943                         
1944         if(wp->flag & VP_COLINDEX) {
1945                 for(index=0; index<totindex; index++) {
1946                         if(indexar[index] && indexar[index]<=me->totface) {
1947                                 MFace *mface= ((MFace *)me->mface) + (indexar[index]-1);
1948                                                 
1949                                 if(mface->mat_nr!=ob->actcol-1) {
1950                                         indexar[index]= 0;
1951                                 }
1952                         }
1953                 }
1954         }
1955
1956         if((me->editflag & ME_EDIT_PAINT_MASK) && me->mface) {
1957                 for(index=0; index<totindex; index++) {
1958                         if(indexar[index] && indexar[index]<=me->totface) {
1959                                 MFace *mface= ((MFace *)me->mface) + (indexar[index]-1);
1960                                                 
1961                                 if((mface->flag & ME_FACE_SEL)==0) {
1962                                         indexar[index]= 0;
1963                                 }
1964                         }                                       
1965                 }
1966         }
1967                         
1968         /* make sure each vertex gets treated only once */
1969         /* and calculate filter weight */
1970         totw= 0;
1971         if(brush->vertexpaint_tool==VP_BLUR) 
1972                 paintweight= 0.0f;
1973         else
1974                 paintweight= ts->vgroup_weight;
1975                         
1976         for(index=0; index<totindex; index++) {
1977                 if(indexar[index] && indexar[index]<=me->totface) {
1978                         MFace *mface= me->mface + (indexar[index]-1);
1979
1980                         if(use_vert_sel) {
1981                                 me->dvert[mface->v1].flag = (me->mvert[mface->v1].flag & SELECT);
1982                                 me->dvert[mface->v2].flag = (me->mvert[mface->v2].flag & SELECT);
1983                                 me->dvert[mface->v3].flag = (me->mvert[mface->v3].flag & SELECT);
1984                                 if(mface->v4) me->dvert[mface->v4].flag = (me->mvert[mface->v4].flag & SELECT);
1985                         }
1986                         else {
1987                                 me->dvert[mface->v1].flag= 1;
1988                                 me->dvert[mface->v2].flag= 1;
1989                                 me->dvert[mface->v3].flag= 1;
1990                                 if(mface->v4) me->dvert[mface->v4].flag= 1;
1991                         }
1992                                         
1993                         if(brush->vertexpaint_tool==VP_BLUR) {
1994                                 MDeformWeight *dw, *(*dw_func)(MDeformVert *, const int);
1995                                                 
1996                                 if(wp->flag & VP_ONLYVGROUP)
1997                                         dw_func= (MDeformWeight *(*)(MDeformVert *, const int))defvert_find_index;
1998                                 else
1999                                         dw_func= defvert_verify_index;
2000                                                 
2001                                 dw= dw_func(me->dvert+mface->v1, ob->actdef-1);
2002                                 if(dw) {paintweight+= dw->weight; totw++;}
2003                                 dw= dw_func(me->dvert+mface->v2, ob->actdef-1);
2004                                 if(dw) {paintweight+= dw->weight; totw++;}
2005                                 dw= dw_func(me->dvert+mface->v3, ob->actdef-1);
2006                                 if(dw) {paintweight+= dw->weight; totw++;}
2007                                 if(mface->v4) {
2008                                         dw= dw_func(me->dvert+mface->v4, ob->actdef-1);
2009                                         if(dw) {paintweight+= dw->weight; totw++;}
2010                                 }
2011                         }
2012                 }
2013         }
2014                         
2015         if(brush->vertexpaint_tool==VP_BLUR) 
2016                 paintweight/= (float)totw;
2017                         
2018         for(index=0; index<totindex; index++) {
2019                                 
2020                 if(indexar[index] && indexar[index]<=me->totface) {
2021                         MFace *mface= me->mface + (indexar[index]-1);
2022                                         
2023                         if((me->dvert+mface->v1)->flag) {
2024                                 alpha= calc_vp_alpha_dl(wp, vc, wpd->wpimat, wpd->vertexcosnos+6*mface->v1, mval, pressure);
2025                                 if(alpha) {
2026                                         do_weight_paint_vertex(wp, ob, mface->v1, 
2027                                                 alpha, paintweight, flip, wpd->vgroup_mirror, 
2028                                                 wpd->vgroup_validmap, ts->multipaint);
2029                                 }
2030                                 (me->dvert+mface->v1)->flag= 0;
2031                         }
2032                                         
2033                         if((me->dvert+mface->v2)->flag) {
2034                                 alpha= calc_vp_alpha_dl(wp, vc, wpd->wpimat, wpd->vertexcosnos+6*mface->v2, mval, pressure);
2035                                 if(alpha) {
2036                                         do_weight_paint_vertex(wp, ob, mface->v2, 
2037                                                 alpha, paintweight, flip, wpd->vgroup_mirror, 
2038                                                 wpd->vgroup_validmap, ts->multipaint);
2039                                 }
2040                                 (me->dvert+mface->v2)->flag= 0;
2041                         }
2042                                         
2043                         if((me->dvert+mface->v3)->flag) {
2044                                 alpha= calc_vp_alpha_dl(wp, vc, wpd->wpimat, wpd->vertexcosnos+6*mface->v3, mval, pressure);
2045                                 if(alpha) {
2046                                         do_weight_paint_vertex(wp, ob, mface->v3, 
2047                                                 alpha, paintweight, flip, wpd->vgroup_mirror, 
2048                                                 wpd->vgroup_validmap, ts->multipaint);
2049                                 }
2050                                 (me->dvert+mface->v3)->flag= 0;
2051                         }
2052                                         
2053                         if((me->dvert+mface->v4)->flag) {
2054                                 if(mface->v4) {
2055                                         alpha= calc_vp_alpha_dl(wp, vc, wpd->wpimat, wpd->vertexcosnos+6*mface->v4, mval, pressure);
2056                                         if(alpha) {
2057                                                 do_weight_paint_vertex(wp, ob, mface->v4, 
2058                                                         alpha, paintweight, flip, wpd->vgroup_mirror,
2059                                                         wpd->vgroup_validmap, ts->multipaint);
2060                                         }
2061                                         (me->dvert+mface->v4)->flag= 0;
2062                                 }
2063                         }
2064                 }
2065         }
2066                         
2067         swap_m4m4(vc->rv3d->persmat, mat);
2068                         
2069         DAG_id_tag_update(ob->data, 0);
2070         ED_region_tag_redraw(vc->ar);
2071 }
2072
2073 static void wpaint_stroke_done(bContext *C, struct PaintStroke *stroke)
2074 {
2075         ToolSettings *ts= CTX_data_tool_settings(C);
2076         Object *ob= CTX_data_active_object(C);
2077         struct WPaintData *wpd= paint_stroke_mode_data(stroke);
2078         
2079         if(wpd) {
2080                 if(wpd->vertexcosnos)
2081                         MEM_freeN(wpd->vertexcosnos);
2082                 MEM_freeN(wpd->indexar);
2083                 
2084                 if (wpd->vgroup_validmap)
2085                         MEM_freeN(wpd->vgroup_validmap);
2086                 
2087                 MEM_freeN(wpd);
2088         }
2089         
2090         /* frees prev buffer */
2091         copy_wpaint_prev(ts->wpaint, NULL, 0);
2092         
2093         /* and particles too */
2094         if(ob->particlesystem.first) {
2095                 ParticleSystem *psys;
2096                 int i;
2097                 
2098                 for(psys= ob->particlesystem.first; psys; psys= psys->next) {
2099                         for(i=0; i<PSYS_TOT_VG; i++) {
2100                                 if(psys->vgroup[i]==ob->actdef) {
2101                                         psys->recalc |= PSYS_RECALC_RESET;
2102                                         break;
2103                                 }
2104                         }
2105                 }
2106         }
2107         
2108         DAG_id_tag_update(ob->data, 0);
2109 }
2110
2111
2112 static int wpaint_invoke(bContext *C, wmOperator *op, wmEvent *event)
2113 {
2114         
2115         op->customdata = paint_stroke_new(C, NULL, wpaint_stroke_test_start,
2116                                           wpaint_stroke_update_step,
2117                                           wpaint_stroke_done, event->type);
2118         
2119         /* add modal handler */
2120         WM_event_add_modal_handler(C, op);
2121
2122         op->type->modal(C, op, event);
2123         
2124         return OPERATOR_RUNNING_MODAL;
2125 }
2126
2127 static int wpaint_cancel(bContext *C, wmOperator *op)
2128 {
2129         paint_stroke_cancel(C, op);
2130
2131         return OPERATOR_CANCELLED;
2132 }
2133
2134 void PAINT_OT_weight_paint(wmOperatorType *ot)
2135 {
2136         
2137         /* identifiers */
2138         ot->name= "Weight Paint";
2139         ot->idname= "PAINT_OT_weight_paint";
2140         
2141         /* api callbacks */
2142         ot->invoke= wpaint_invoke;
2143         ot->modal= paint_stroke_modal;
2144         /* ot->exec= vpaint_exec; <-- needs stroke property */
2145         ot->poll= weight_paint_poll;
2146         ot->cancel= wpaint_cancel;
2147         
2148         /* flags */
2149         ot->flag= OPTYPE_REGISTER|OPTYPE_UNDO|OPTYPE_BLOCKING;
2150
2151         RNA_def_collection_runtime(ot->srna, "stroke", &RNA_OperatorStrokeElement, "Stroke", "");
2152 }
2153
2154 static int weight_paint_set_exec(bContext *C, wmOperator *UNUSED(op))
2155 {
2156         struct Scene *scene= CTX_data_scene(C);
2157         Object *obact = CTX_data_active_object(C);
2158
2159         wpaint_fill(scene->toolsettings->wpaint, obact, scene->toolsettings->vgroup_weight);
2160         ED_region_tag_redraw(CTX_wm_region(C)); // XXX - should redraw all 3D views
2161         return OPERATOR_FINISHED;
2162 }
2163
2164 void PAINT_OT_weight_set(wmOperatorType *ot)
2165 {
2166         /* identifiers */
2167         ot->name= "Set Weight";
2168         ot->idname= "PAINT_OT_weight_set";
2169
2170         /* api callbacks */
2171         ot->exec= weight_paint_set_exec;
2172         ot->poll= mask_paint_poll; /* Radish, it was facemask_paint_poll */
2173
2174         /* flags */
2175         ot->flag= OPTYPE_REGISTER|OPTYPE_UNDO;
2176 }
2177
2178 /* ************ set / clear vertex paint mode ********** */
2179
2180
2181 static int set_vpaint(bContext *C, wmOperator *op)              /* toggle */
2182 {       
2183         Object *ob= CTX_data_active_object(C);
2184         Scene *scene= CTX_data_scene(C);
2185         VPaint *vp= scene->toolsettings->vpaint;
2186         Mesh *me;
2187         
2188         me= get_mesh(ob);
2189         
2190         if(me==NULL || object_data_is_libdata(ob)) {
2191                 ob->mode &= ~OB_MODE_VERTEX_PAINT;
2192                 return OPERATOR_PASS_THROUGH;
2193         }
2194         
2195         if(me && me->mcol==NULL) make_vertexcol(ob);
2196         
2197         /* toggle: end vpaint */
2198         if(ob->mode & OB_MODE_VERTEX_PAINT) {
2199                 
2200                 ob->mode &= ~OB_MODE_VERTEX_PAINT;
2201         }
2202         else {
2203                 ob->mode |= OB_MODE_VERTEX_PAINT;
2204                 /* Turn off weight painting */
2205                 if (ob->mode & OB_MODE_WEIGHT_PAINT)
2206                         set_wpaint(C, op);
2207                 
2208                 if(vp==NULL)
2209                         vp= scene->toolsettings->vpaint= new_vpaint(0);
2210                 
2211                 paint_cursor_start(C, vertex_paint_poll);
2212
2213                 paint_init(&vp->paint, PAINT_CURSOR_VERTEX_PAINT);
2214         }
2215         
2216         if (me)
2217                 /* update modifier stack for mapping requirements */
2218                 DAG_id_tag_update(&me->id, 0);
2219         
2220         WM_event_add_notifier(C, NC_SCENE|ND_MODE, scene);
2221         
2222         return OPERATOR_FINISHED;
2223 }
2224
2225 void PAINT_OT_vertex_paint_toggle(wmOperatorType *ot)
2226 {
2227         
2228         /* identifiers */
2229         ot->name= "Vertex Paint Mode";
2230         ot->idname= "PAINT_OT_vertex_paint_toggle";
2231         
2232         /* api callbacks */
2233         ot->exec= set_vpaint;
2234         ot->poll= paint_poll_test;
2235         
2236         /* flags */
2237         ot->flag= OPTYPE_REGISTER|OPTYPE_UNDO;
2238 }
2239
2240
2241
2242 /* ********************** vertex paint operator ******************* */
2243
2244 /* Implementation notes:
2245
2246 Operator->invoke()
2247   - validate context (add mcol)
2248   - create customdata storage
2249   - call paint once (mouse click)
2250   - add modal handler 
2251
2252 Operator->modal()
2253   - for every mousemove, apply vertex paint
2254   - exit on mouse release, free customdata
2255         (return OPERATOR_FINISHED also removes handler and operator)
2256
2257 For future:
2258   - implement a stroke event (or mousemove with past positons)
2259   - revise whether op->customdata should be added in object, in set_vpaint
2260
2261 */
2262
2263 typedef struct VPaintData {
2264         ViewContext vc;
2265         unsigned int paintcol;
2266         int *indexar;
2267         float *vertexcosnos;
2268         float vpimat[3][3];
2269 } VPaintData;
2270
2271 static int vpaint_stroke_test_start(bContext *C, struct wmOperator *op, wmEvent *UNUSED(event))
2272 {
2273         ToolSettings *ts= CTX_data_tool_settings(C);
2274         struct PaintStroke *stroke = op->customdata;
2275         VPaint *vp= ts->vpaint;
2276         struct VPaintData *vpd;
2277         Object *ob= CTX_data_active_object(C);
2278         Mesh *me;
2279         float mat[4][4], imat[4][4];
2280
2281         /* context checks could be a poll() */
2282         me= get_mesh(ob);
2283         if(me==NULL || me->totface==0) return OPERATOR_PASS_THROUGH;
2284         
2285         if(me->mcol==NULL) make_vertexcol(ob);
2286         if(me->mcol==NULL) return OPERATOR_CANCELLED;
2287         
2288         /* make mode data storage */
2289         vpd= MEM_callocN(sizeof(struct VPaintData), "VPaintData");
2290         paint_stroke_set_mode_data(stroke, vpd);
2291         view3d_set_viewcontext(C, &vpd->vc);
2292         
2293         vpd->vertexcosnos= mesh_get_mapped_verts_nors(vpd->vc.scene, ob);
2294         vpd->indexar= get_indexarray(me);
2295         vpd->paintcol= vpaint_get_current_col(vp);
2296         
2297         /* for filtering */
2298         copy_vpaint_prev(vp, (unsigned int *)me->mcol, me->totface);
2299         
2300         /* some old cruft to sort out later */
2301         mul_m4_m4m4(mat, ob->obmat, vpd->vc.rv3d->viewmat);
2302         invert_m4_m4(imat, mat);
2303         copy_m3_m4(vpd->vpimat, imat);
2304
2305         return 1;
2306 }
2307
2308 static void vpaint_paint_face(VPaint *vp, VPaintData *vpd, Object *ob, int index, const float mval[2], float pressure, int UNUSED(flip))
2309 {
2310         ViewContext *vc = &vpd->vc;
2311         Brush *brush = paint_brush(&vp->paint);
2312         Mesh *me = get_mesh(ob);
2313         MFace *mface= ((MFace*)me->mface) + index;
2314         unsigned int *mcol= ((unsigned int*)me->mcol) + 4*index;
2315         unsigned int *mcolorig= ((unsigned int*)vp->vpaint_prev) + 4*index;
2316         float alpha;
2317         int i;
2318         
2319         if((vp->flag & VP_COLINDEX && mface->mat_nr!=ob->actcol-1) ||
2320            ((me->editflag & ME_EDIT_PAINT_MASK) && !(mface->flag & ME_FACE_SEL)))
2321                 return;
2322
2323         if(brush->vertexpaint_tool==VP_BLUR) {
2324                 unsigned int fcol1= mcol_blend( mcol[0], mcol[1], 128);
2325                 if(mface->v4) {
2326                         unsigned int fcol2= mcol_blend( mcol[2], mcol[3], 128);
2327                         vpd->paintcol= mcol_blend( fcol1, fcol2, 128);
2328                 }
2329                 else {
2330                         vpd->paintcol= mcol_blend( mcol[2], fcol1, 170);
2331                 }
2332                 
2333         }
2334
2335         for(i = 0; i < (mface->v4 ? 4 : 3); ++i) {
2336                 alpha= calc_vp_alpha_dl(vp, vc, vpd->vpimat, vpd->vertexcosnos+6*(&mface->v1)[i], mval, pressure);
2337                 if(alpha)
2338                         vpaint_blend(vp, mcol+i, mcolorig+i, vpd->paintcol, (int)(alpha*255.0f));
2339         }
2340 }
2341
2342 static void vpaint_stroke_update_step(bContext *C, struct PaintStroke *stroke, PointerRNA *itemptr)
2343 {
2344         ToolSettings *ts= CTX_data_tool_settings(C);
2345         struct VPaintData *vpd = paint_stroke_mode_data(stroke);
2346         VPaint *vp= ts->vpaint;
2347         Brush *brush = paint_brush(&vp->paint);
2348         ViewContext *vc= &vpd->vc;
2349         Object *ob= vc->obact;
2350         Mesh *me= ob->data;
2351         float mat[4][4];
2352         int *indexar= vpd->indexar;
2353         int totindex, index, flip;
2354         float pressure, mval[2];
2355
2356         RNA_float_get_array(itemptr, "mouse", mval);
2357         flip = RNA_boolean_get(itemptr, "pen_flip");
2358         pressure = RNA_float_get(itemptr, "pressure");
2359                         
2360         view3d_operator_needs_opengl(C);
2361                         
2362         /* load projection matrix */
2363         mul_m4_m4m4(mat, ob->obmat, vc->rv3d->persmat);
2364
2365         mval[0]-= vc->ar->winrct.xmin;
2366         mval[1]-= vc->ar->winrct.ymin;
2367
2368                         
2369         /* which faces are involved */
2370         if(vp->flag & VP_AREA) {
2371                 totindex= sample_backbuf_area(vc, indexar, me->totface, mval[0], mval[1], brush_size(brush));
2372         }
2373         else {
2374                 indexar[0]= view3d_sample_backbuf(vc, mval[0], mval[1]);
2375                 if(indexar[0]) totindex= 1;
2376                 else totindex= 0;
2377         }
2378                         
2379         swap_m4m4(vc->rv3d->persmat, mat);
2380                         
2381         for(index=0; index<totindex; index++) {                         
2382                 if(indexar[index] && indexar[index]<=me->totface)
2383                         vpaint_paint_face(vp, vpd, ob, indexar[index]-1, mval, pressure, flip);
2384         }
2385                                                 
2386         swap_m4m4(vc->rv3d->persmat, mat);
2387
2388         /* was disabled because it is slow, but necessary for blur */
2389         if(brush->vertexpaint_tool == VP_BLUR)
2390                 do_shared_vertexcol(me);
2391                         
2392         ED_region_tag_redraw(vc->ar);
2393                         
2394         DAG_id_tag_update(ob->data, 0);
2395 }
2396
2397 static void vpaint_stroke_done(bContext *C, struct PaintStroke *stroke)
2398 {
2399         ToolSettings *ts= CTX_data_tool_settings(C);
2400         struct VPaintData *vpd= paint_stroke_mode_data(stroke);
2401         
2402         if(vpd->vertexcosnos)
2403                 MEM_freeN(vpd->vertexcosnos);
2404         MEM_freeN(vpd->indexar);
2405         
2406         /* frees prev buffer */
2407         copy_vpaint_prev(ts->vpaint, NULL, 0);
2408         
2409         MEM_freeN(vpd);
2410 }
2411
2412 static int vpaint_invoke(bContext *C, wmOperator *op, wmEvent *event)
2413 {
2414         
2415         op->customdata = paint_stroke_new(C, NULL, vpaint_stroke_test_start,
2416                                           vpaint_stroke_update_step,
2417                                           vpaint_stroke_done, event->type);
2418         
2419         /* add modal handler */
2420         WM_event_add_modal_handler(C, op);
2421
2422         op->type->modal(C, op, event);
2423         
2424         return OPERATOR_RUNNING_MODAL;
2425 }
2426
2427 static int vpaint_cancel(bContext *C, wmOperator *op)
2428 {
2429         paint_stroke_cancel(C, op);
2430
2431         return OPERATOR_CANCELLED;
2432 }
2433
2434 void PAINT_OT_vertex_paint(wmOperatorType *ot)
2435 {
2436         /* identifiers */
2437         ot->name= "Vertex Paint";
2438         ot->idname= "PAINT_OT_vertex_paint";
2439         
2440         /* api callbacks */
2441         ot->invoke= vpaint_invoke;
2442         ot->modal= paint_stroke_modal;
2443         /* ot->exec= vpaint_exec; <-- needs stroke property */
2444         ot->poll= vertex_paint_poll;
2445         ot->cancel= vpaint_cancel;
2446         
2447         /* flags */
2448         ot->flag= OPTYPE_REGISTER|OPTYPE_UNDO|OPTYPE_BLOCKING;
2449
2450         RNA_def_collection_runtime(ot->srna, "stroke", &RNA_OperatorStrokeElement, "Stroke", "");
2451 }
2452
2453 /* ********************** weight from bones operator ******************* */
2454
2455 static int weight_from_bones_poll(bContext *C)
2456 {
2457         Object *ob= CTX_data_active_object(C);
2458
2459         return (ob && (ob->mode & OB_MODE_WEIGHT_PAINT) && modifiers_isDeformedByArmature(ob));
2460 }
2461
2462 static int weight_from_bones_exec(bContext *C, wmOperator *op)
2463 {
2464         Scene *scene= CTX_data_scene(C);
2465         Object *ob= CTX_data_active_object(C);
2466         Object *armob= modifiers_isDeformedByArmature(ob);
2467         Mesh *me= ob->data;
2468         int type= RNA_enum_get(op->ptr, "type");
2469
2470         create_vgroups_from_armature(op->reports, scene, ob, armob, type, (me->editflag & ME_EDIT_MIRROR_X));
2471
2472         DAG_id_tag_update(&me->id, 0);
2473         WM_event_add_notifier(C, NC_GEOM|ND_DATA, me);
2474
2475         return OPERATOR_FINISHED;
2476 }
2477
2478 void PAINT_OT_weight_from_bones(wmOperatorType *ot)
2479 {
2480         static EnumPropertyItem type_items[]= {
2481                 {ARM_GROUPS_AUTO, "AUTOMATIC", 0, "Automatic", "Automatic weights froms bones"},
2482                 {ARM_GROUPS_ENVELOPE, "ENVELOPES", 0, "From Envelopes", "Weights from envelopes with user defined radius"},
2483                 {0, NULL, 0, NULL, NULL}};
2484
2485         /* identifiers */
2486         ot->name= "Weight from Bones";
2487         ot->idname= "PAINT_OT_weight_from_bones";
2488         
2489         /* api callbacks */
2490         ot->exec= weight_from_bones_exec;
2491         ot->invoke= WM_menu_invoke;
2492         ot->poll= weight_from_bones_poll;
2493         
2494         /* flags */
2495         ot->flag= OPTYPE_REGISTER|OPTYPE_UNDO;
2496
2497         /* properties */
2498         ot->prop= RNA_def_enum(ot->srna, "type", type_items, 0, "Type", "Method to use for assigning weights.");
2499 }
2500