fix T37411: Transform mouse constraint could fail in some situations.
[blender-staging.git] / source / blender / editors / transform / transform_constraints.c
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  *
18  * The Original Code is Copyright (C) 2001-2002 by NaN Holding BV.
19  * All rights reserved.
20  *
21  * The Original Code is: all of this file.
22  *
23  * Contributor(s): none yet.
24  *
25  * ***** END GPL LICENSE BLOCK *****
26  */
27
28 /** \file blender/editors/transform/transform_constraints.c
29  *  \ingroup edtransform
30  */
31
32 #include <stdlib.h>
33 #include <stdio.h>
34 #include <string.h>
35 #include <math.h>
36
37 #ifndef WIN32
38 #  include <unistd.h>
39 #else
40 #  include <io.h>
41 #endif
42
43 #include "DNA_object_types.h"
44 #include "DNA_scene_types.h"
45 #include "DNA_screen_types.h"
46 #include "DNA_space_types.h"
47 #include "DNA_view3d_types.h"
48
49 #include "BIF_gl.h"
50 #include "BIF_glutil.h"
51
52 #include "BLI_math.h"
53 #include "BLI_utildefines.h"
54 #include "BLI_string.h"
55
56 #include "BKE_context.h"
57
58 #include "ED_image.h"
59 #include "ED_view3d.h"
60
61 #include "BLF_translation.h"
62
63 #include "UI_resources.h"
64
65 #include "transform.h"
66
67 static void drawObjectConstraint(TransInfo *t);
68
69 /* ************************** CONSTRAINTS ************************* */
70 static void constraintAutoValues(TransInfo *t, float vec[3])
71 {
72         int mode = t->con.mode;
73         if (mode & CON_APPLY) {
74                 float nval = (t->flag & T_NULL_ONE) ? 1.0f : 0.0f;
75
76                 if ((mode & CON_AXIS0) == 0) {
77                         vec[0] = nval;
78                 }
79                 if ((mode & CON_AXIS1) == 0) {
80                         vec[1] = nval;
81                 }
82                 if ((mode & CON_AXIS2) == 0) {
83                         vec[2] = nval;
84                 }
85         }
86 }
87
88 void constraintNumInput(TransInfo *t, float vec[3])
89 {
90         int mode = t->con.mode;
91         if (mode & CON_APPLY) {
92                 float nval = (t->flag & T_NULL_ONE) ? 1.0f : 0.0f;
93
94                 if (getConstraintSpaceDimension(t) == 2) {
95                         int axis = mode & (CON_AXIS0 | CON_AXIS1 | CON_AXIS2);
96                         if (axis == (CON_AXIS0 | CON_AXIS1)) {
97                                 /* vec[0] = vec[0]; */ /* same */
98                                 /* vec[1] = vec[1]; */ /* same */
99                                 vec[2] = nval;
100                         }
101                         else if (axis == (CON_AXIS1 | CON_AXIS2)) {
102                                 vec[2] = vec[1];
103                                 vec[1] = vec[0];
104                                 vec[0] = nval;
105                         }
106                         else if (axis == (CON_AXIS0 | CON_AXIS2)) {
107                                 /* vec[0] = vec[0]; */  /* same */
108                                 vec[2] = vec[1];
109                                 vec[1] = nval;
110                         }
111                 }
112                 else if (getConstraintSpaceDimension(t) == 1) {
113                         if (mode & CON_AXIS0) {
114                                 /* vec[0] = vec[0]; */ /* same */
115                                 vec[1] = nval;
116                                 vec[2] = nval;
117                         }
118                         else if (mode & CON_AXIS1) {
119                                 vec[1] = vec[0];
120                                 vec[0] = nval;
121                                 vec[2] = nval;
122                         }
123                         else if (mode & CON_AXIS2) {
124                                 vec[2] = vec[0];
125                                 vec[0] = nval;
126                                 vec[1] = nval;
127                         }
128                 }
129         }
130 }
131
132 static void postConstraintChecks(TransInfo *t, float vec[3], float pvec[3])
133 {
134         int i = 0;
135
136         mul_m3_v3(t->con.imtx, vec);
137
138         snapGridIncrement(t, vec);
139
140         if (t->num.flag & T_NULL_ONE) {
141                 if (!(t->con.mode & CON_AXIS0))
142                         vec[0] = 1.0f;
143
144                 if (!(t->con.mode & CON_AXIS1))
145                         vec[1] = 1.0f;
146
147                 if (!(t->con.mode & CON_AXIS2))
148                         vec[2] = 1.0f;
149         }
150
151         if (hasNumInput(&t->num)) {
152                 applyNumInput(&t->num, vec);
153                 constraintNumInput(t, vec);
154                 removeAspectRatio(t, vec);
155         }
156
157         /* autovalues is operator param, use that directly but not if snapping is forced */
158         if (t->flag & T_AUTOVALUES && (t->tsnap.status & SNAP_FORCED) == 0) {
159                 mul_v3_m3v3(vec, t->con.imtx, t->auto_values);
160                 constraintAutoValues(t, vec);
161                 /* inverse transformation at the end */
162         }
163
164         if (t->con.mode & CON_AXIS0) {
165                 pvec[i++] = vec[0];
166         }
167         if (t->con.mode & CON_AXIS1) {
168                 pvec[i++] = vec[1];
169         }
170         if (t->con.mode & CON_AXIS2) {
171                 pvec[i++] = vec[2];
172         }
173
174         mul_m3_v3(t->con.mtx, vec);
175 }
176
177 static void viewAxisCorrectCenter(TransInfo *t, float t_con_center[3])
178 {
179         if (t->spacetype == SPACE_VIEW3D) {
180                 // View3D *v3d = t->sa->spacedata.first;
181                 const float min_dist = 1.0f;  /* v3d->near; */
182                 float dir[3];
183                 float l;
184
185                 sub_v3_v3v3(dir, t_con_center, t->viewinv[3]);
186                 if (dot_v3v3(dir, t->viewinv[2]) < 0.0f) {
187                         negate_v3(dir);
188                 }
189                 project_v3_v3v3(dir, dir, t->viewinv[2]);
190
191                 l = len_v3(dir);
192
193                 if (l < min_dist) {
194                         float diff[3];
195                         normalize_v3_v3(diff, t->viewinv[2]);
196                         mul_v3_fl(diff, min_dist - l);
197
198                         sub_v3_v3(t_con_center, diff);
199                 }
200         }
201 }
202
203 static void axisProjection(TransInfo *t, const float axis[3], const float in[3], float out[3])
204 {
205         float norm[3], vec[3], factor, angle;
206         float t_con_center[3];
207
208         if (is_zero_v3(in)) {
209                 return;
210         }
211
212         copy_v3_v3(t_con_center, t->con.center);
213
214         /* checks for center being too close to the view center */
215         viewAxisCorrectCenter(t, t_con_center);
216         
217         angle = fabsf(angle_v3v3(axis, t->viewinv[2]));
218         if (angle > (float)M_PI / 2.0f) {
219                 angle = (float)M_PI - angle;
220         }
221         angle = RAD2DEGF(angle);
222
223         /* For when view is parallel to constraint... will cause NaNs otherwise
224          * So we take vertical motion in 3D space and apply it to the
225          * constraint axis. Nice for camera grab + MMB */
226         if (angle < 5.0f) {
227                 project_v3_v3v3(vec, in, t->viewinv[1]);
228                 factor = dot_v3v3(t->viewinv[1], vec) * 2.0f;
229                 /* since camera distance is quite relative, use quadratic relationship. holding shift can compensate */
230                 if (factor < 0.0f) factor *= -factor;
231                 else factor *= factor;
232
233                 copy_v3_v3(out, axis);
234                 normalize_v3(out);
235                 mul_v3_fl(out, -factor);  /* -factor makes move down going backwards */
236         }
237         else {
238                 float v[3], i1[3], i2[3];
239                 float v2[3], v4[3];
240                 float norm_center[3];
241                 float plane[3];
242
243                 getViewVector(t, t_con_center, norm_center);
244                 cross_v3_v3v3(plane, norm_center, axis);
245
246                 project_v3_v3v3(vec, in, plane);
247                 sub_v3_v3v3(vec, in, vec);
248                 
249                 add_v3_v3v3(v, vec, t_con_center);
250                 getViewVector(t, v, norm);
251
252                 /* give arbitrary large value if projection is impossible */
253                 factor = dot_v3v3(axis, norm);
254                 if (1.0f - fabsf(factor) < 0.0002f) {
255                         copy_v3_v3(out, axis);
256                         if (factor > 0) {
257                                 mul_v3_fl(out, 1000000000.0f);
258                         }
259                         else {
260                                 mul_v3_fl(out, -1000000000.0f);
261                         }
262                 }
263                 else {
264                         add_v3_v3v3(v2, t_con_center, axis);
265                         add_v3_v3v3(v4, v, norm);
266                         
267                         isect_line_line_v3(t_con_center, v2, v, v4, i1, i2);
268                         
269                         sub_v3_v3v3(v, i2, v);
270         
271                         sub_v3_v3v3(out, i1, t_con_center);
272
273                         /* possible some values become nan when
274                          * viewpoint and object are both zero */
275                         if (!finite(out[0])) out[0] = 0.0f;
276                         if (!finite(out[1])) out[1] = 0.0f;
277                         if (!finite(out[2])) out[2] = 0.0f;
278                 }
279         }
280 }
281
282 static void planeProjection(TransInfo *t, const float in[3], float out[3])
283 {
284         float vec[3], factor, norm[3];
285
286         add_v3_v3v3(vec, in, t->con.center);
287         getViewVector(t, vec, norm);
288
289         sub_v3_v3v3(vec, out, in);
290
291         factor = dot_v3v3(vec, norm);
292         if (fabsf(factor) <= 0.001f) {
293                 return; /* prevent divide by zero */
294         }
295         factor = dot_v3v3(vec, vec) / factor;
296
297         copy_v3_v3(vec, norm);
298         mul_v3_fl(vec, factor);
299
300         add_v3_v3v3(out, in, vec);
301 }
302
303 /*
304  * Generic callback for constant spatial constraints applied to linear motion
305  *
306  * The IN vector in projected into the constrained space and then further
307  * projected along the view vector.
308  * (in perspective mode, the view vector is relative to the position on screen)
309  *
310  */
311
312 static void applyAxisConstraintVec(TransInfo *t, TransData *td, const float in[3], float out[3], float pvec[3])
313 {
314         copy_v3_v3(out, in);
315         if (!td && t->con.mode & CON_APPLY) {
316                 mul_m3_v3(t->con.pmtx, out);
317
318                 // With snap, a projection is alright, no need to correct for view alignment
319                 if (!(t->tsnap.mode != SCE_SNAP_MODE_INCREMENT && activeSnap(t))) {
320                         if (getConstraintSpaceDimension(t) == 2) {
321                                 if (out[0] != 0.0f || out[1] != 0.0f || out[2] != 0.0f) {
322                                         planeProjection(t, in, out);
323                                 }
324                         }
325                         else if (getConstraintSpaceDimension(t) == 1) {
326                                 float c[3];
327
328                                 if (t->con.mode & CON_AXIS0) {
329                                         copy_v3_v3(c, t->con.mtx[0]);
330                                 }
331                                 else if (t->con.mode & CON_AXIS1) {
332                                         copy_v3_v3(c, t->con.mtx[1]);
333                                 }
334                                 else if (t->con.mode & CON_AXIS2) {
335                                         copy_v3_v3(c, t->con.mtx[2]);
336                                 }
337                                 axisProjection(t, c, in, out);
338                         }
339                 }
340                 postConstraintChecks(t, out, pvec);
341         }
342 }
343
344 /*
345  * Generic callback for object based spatial constraints applied to linear motion
346  *
347  * At first, the following is applied to the first data in the array
348  * The IN vector in projected into the constrained space and then further
349  * projected along the view vector.
350  * (in perspective mode, the view vector is relative to the position on screen)
351  *
352  * Further down, that vector is mapped to each data's space.
353  */
354
355 static void applyObjectConstraintVec(TransInfo *t, TransData *td, const float in[3], float out[3], float pvec[3])
356 {
357         copy_v3_v3(out, in);
358         if (t->con.mode & CON_APPLY) {
359                 if (!td) {
360                         mul_m3_v3(t->con.pmtx, out);
361                         if (getConstraintSpaceDimension(t) == 2) {
362                                 if (out[0] != 0.0f || out[1] != 0.0f || out[2] != 0.0f) {
363                                         planeProjection(t, in, out);
364                                 }
365                         }
366                         else if (getConstraintSpaceDimension(t) == 1) {
367                                 float c[3];
368
369                                 if (t->con.mode & CON_AXIS0) {
370                                         copy_v3_v3(c, t->con.mtx[0]);
371                                 }
372                                 else if (t->con.mode & CON_AXIS1) {
373                                         copy_v3_v3(c, t->con.mtx[1]);
374                                 }
375                                 else if (t->con.mode & CON_AXIS2) {
376                                         copy_v3_v3(c, t->con.mtx[2]);
377                                 }
378                                 axisProjection(t, c, in, out);
379                         }
380                         postConstraintChecks(t, out, pvec);
381                         copy_v3_v3(out, pvec);
382                 }
383                 else {
384                         int i = 0;
385
386                         out[0] = out[1] = out[2] = 0.0f;
387                         if (t->con.mode & CON_AXIS0) {
388                                 out[0] = in[i++];
389                         }
390                         if (t->con.mode & CON_AXIS1) {
391                                 out[1] = in[i++];
392                         }
393                         if (t->con.mode & CON_AXIS2) {
394                                 out[2] = in[i++];
395                         }
396
397                         mul_m3_v3(td->axismtx, out);
398                         if (t->flag & T_EDIT) {
399                                 mul_m3_v3(t->obedit_mat, out);
400                         }
401                 }
402         }
403 }
404
405 /*
406  * Generic callback for constant spatial constraints applied to resize motion
407  */
408
409 static void applyAxisConstraintSize(TransInfo *t, TransData *td, float smat[3][3])
410 {
411         if (!td && t->con.mode & CON_APPLY) {
412                 float tmat[3][3];
413
414                 if (!(t->con.mode & CON_AXIS0)) {
415                         smat[0][0] = 1.0f;
416                 }
417                 if (!(t->con.mode & CON_AXIS1)) {
418                         smat[1][1] = 1.0f;
419                 }
420                 if (!(t->con.mode & CON_AXIS2)) {
421                         smat[2][2] = 1.0f;
422                 }
423
424                 mul_m3_m3m3(tmat, smat, t->con.imtx);
425                 mul_m3_m3m3(smat, t->con.mtx, tmat);
426         }
427 }
428
429 /*
430  * Callback for object based spatial constraints applied to resize motion
431  */
432
433 static void applyObjectConstraintSize(TransInfo *t, TransData *td, float smat[3][3])
434 {
435         if (td && t->con.mode & CON_APPLY) {
436                 float tmat[3][3];
437                 float imat[3][3];
438
439                 invert_m3_m3(imat, td->axismtx);
440
441                 if (!(t->con.mode & CON_AXIS0)) {
442                         smat[0][0] = 1.0f;
443                 }
444                 if (!(t->con.mode & CON_AXIS1)) {
445                         smat[1][1] = 1.0f;
446                 }
447                 if (!(t->con.mode & CON_AXIS2)) {
448                         smat[2][2] = 1.0f;
449                 }
450
451                 mul_m3_m3m3(tmat, smat, imat);
452                 if (t->flag & T_EDIT) {
453                         mul_m3_m3m3(smat, t->obedit_mat, smat);
454                 }
455                 mul_m3_m3m3(smat, td->axismtx, tmat);
456         }
457 }
458
459 /*
460  * Generic callback for constant spatial constraints applied to rotations
461  *
462  * The rotation axis is copied into VEC.
463  *
464  * In the case of single axis constraints, the rotation axis is directly the one constrained to.
465  * For planar constraints (2 axis), the rotation axis is the normal of the plane.
466  *
467  * The following only applies when CON_NOFLIP is not set.
468  * The vector is then modified to always point away from the screen (in global space)
469  * This insures that the rotation is always logically following the mouse.
470  * (ie: not doing counterclockwise rotations when the mouse moves clockwise).
471  */
472
473 static void applyAxisConstraintRot(TransInfo *t, TransData *td, float vec[3], float *angle)
474 {
475         if (!td && t->con.mode & CON_APPLY) {
476                 int mode = t->con.mode & (CON_AXIS0 | CON_AXIS1 | CON_AXIS2);
477
478                 switch (mode) {
479                         case CON_AXIS0:
480                         case (CON_AXIS1 | CON_AXIS2):
481                                 copy_v3_v3(vec, t->con.mtx[0]);
482                                 break;
483                         case CON_AXIS1:
484                         case (CON_AXIS0 | CON_AXIS2):
485                                 copy_v3_v3(vec, t->con.mtx[1]);
486                                 break;
487                         case CON_AXIS2:
488                         case (CON_AXIS0 | CON_AXIS1):
489                                 copy_v3_v3(vec, t->con.mtx[2]);
490                                 break;
491                 }
492                 /* don't flip axis if asked to or if num input */
493                 if (angle && (mode & CON_NOFLIP) == 0 && hasNumInput(&t->num) == 0) {
494                         if (dot_v3v3(vec, t->viewinv[2]) > 0.0f) {
495                                 *angle = -(*angle);
496                         }
497                 }
498         }
499 }
500
501 /*
502  * Callback for object based spatial constraints applied to rotations
503  *
504  * The rotation axis is copied into VEC.
505  *
506  * In the case of single axis constraints, the rotation axis is directly the one constrained to.
507  * For planar constraints (2 axis), the rotation axis is the normal of the plane.
508  *
509  * The following only applies when CON_NOFLIP is not set.
510  * The vector is then modified to always point away from the screen (in global space)
511  * This insures that the rotation is always logically following the mouse.
512  * (ie: not doing counterclockwise rotations when the mouse moves clockwise).
513  */
514
515 static void applyObjectConstraintRot(TransInfo *t, TransData *td, float vec[3], float *angle)
516 {
517         if (t->con.mode & CON_APPLY) {
518                 int mode = t->con.mode & (CON_AXIS0 | CON_AXIS1 | CON_AXIS2);
519                 float tmp_axismtx[3][3];
520                 float (*axismtx)[3];
521
522                 /* on setup call, use first object */
523                 if (td == NULL) {
524                         td = t->data;
525                 }
526
527                 if (t->flag & T_EDIT) {
528                         mul_m3_m3m3(tmp_axismtx, t->obedit_mat, td->axismtx);
529                         axismtx = tmp_axismtx;
530                 }
531                 else {
532                         axismtx = td->axismtx;
533                 }
534
535                 switch (mode) {
536                         case CON_AXIS0:
537                         case (CON_AXIS1 | CON_AXIS2):
538                                 copy_v3_v3(vec, axismtx[0]);
539                                 break;
540                         case CON_AXIS1:
541                         case (CON_AXIS0 | CON_AXIS2):
542                                 copy_v3_v3(vec, axismtx[1]);
543                                 break;
544                         case CON_AXIS2:
545                         case (CON_AXIS0 | CON_AXIS1):
546                                 copy_v3_v3(vec, axismtx[2]);
547                                 break;
548                 }
549                 if (angle && (mode & CON_NOFLIP) == 0 && hasNumInput(&t->num) == 0) {
550                         if (dot_v3v3(vec, t->viewinv[2]) > 0.0f) {
551                                 *angle = -(*angle);
552                         }
553                 }
554         }
555 }
556
557 /*--------------------- INTERNAL SETUP CALLS ------------------*/
558
559 void setConstraint(TransInfo *t, float space[3][3], int mode, const char text[])
560 {
561         BLI_strncpy(t->con.text + 1, text, sizeof(t->con.text) - 1);
562         copy_m3_m3(t->con.mtx, space);
563         t->con.mode = mode;
564         getConstraintMatrix(t);
565
566         startConstraint(t);
567
568         t->con.drawExtra = NULL;
569         t->con.applyVec = applyAxisConstraintVec;
570         t->con.applySize = applyAxisConstraintSize;
571         t->con.applyRot = applyAxisConstraintRot;
572         t->redraw = TREDRAW_HARD;
573 }
574
575 /* applies individual td->axismtx constraints */
576 void setAxisMatrixConstraint(TransInfo *t, int mode, const char text[])
577 {
578         if (t->total == 1) {
579                 setConstraint(t, t->data->axismtx, mode, text);
580         }
581         else {
582                 BLI_strncpy(t->con.text + 1, text, sizeof(t->con.text) - 1);
583                 copy_m3_m3(t->con.mtx, t->data->axismtx);
584                 t->con.mode = mode;
585                 getConstraintMatrix(t);
586
587                 startConstraint(t);
588
589                 t->con.drawExtra = drawObjectConstraint;
590                 t->con.applyVec = applyObjectConstraintVec;
591                 t->con.applySize = applyObjectConstraintSize;
592                 t->con.applyRot = applyObjectConstraintRot;
593                 t->redraw = TREDRAW_HARD;
594         }
595 }
596
597 void setLocalConstraint(TransInfo *t, int mode, const char text[])
598 {
599         /* edit-mode now allows local transforms too */
600         if (t->flag & T_EDIT) {
601                 setConstraint(t, t->obedit_mat, mode, text);
602         }
603         else {
604                 setAxisMatrixConstraint(t, mode, text);
605         }
606 }
607
608 /*
609  * Set the constraint according to the user defined orientation
610  *
611  * ftext is a format string passed to BLI_snprintf. It will add the name of
612  * the orientation where %s is (logically).
613  */
614 void setUserConstraint(TransInfo *t, short orientation, int mode, const char ftext[])
615 {
616         char text[40];
617
618         switch (orientation) {
619                 case V3D_MANIP_GLOBAL:
620                 {
621                         float mtx[3][3] = MAT3_UNITY;
622                         BLI_snprintf(text, sizeof(text), ftext, IFACE_("global"));
623                         setConstraint(t, mtx, mode, text);
624                         break;
625                 }
626                 case V3D_MANIP_LOCAL:
627                         BLI_snprintf(text, sizeof(text), ftext, IFACE_("local"));
628                         setLocalConstraint(t, mode, text);
629                         break;
630                 case V3D_MANIP_NORMAL:
631                         BLI_snprintf(text, sizeof(text), ftext, IFACE_("normal"));
632                         if (checkUseAxisMatrix(t)) {
633                                 setAxisMatrixConstraint(t, mode, text);
634                         }
635                         else {
636                                 setConstraint(t, t->spacemtx, mode, text);
637                         }
638                         break;
639                 case V3D_MANIP_VIEW:
640                         BLI_snprintf(text, sizeof(text), ftext, IFACE_("view"));
641                         setConstraint(t, t->spacemtx, mode, text);
642                         break;
643                 case V3D_MANIP_GIMBAL:
644                         BLI_snprintf(text, sizeof(text), ftext, IFACE_("gimbal"));
645                         setConstraint(t, t->spacemtx, mode, text);
646                         break;
647                 default: /* V3D_MANIP_CUSTOM */
648                         BLI_snprintf(text, sizeof(text), ftext, t->spacename);
649                         setConstraint(t, t->spacemtx, mode, text);
650                         break;
651         }
652
653         t->con.orientation = orientation;
654
655         t->con.mode |= CON_USER;
656 }
657
658 /*----------------- DRAWING CONSTRAINTS -------------------*/
659
660 void drawConstraint(TransInfo *t)
661 {
662         TransCon *tc = &(t->con);
663
664         if (!ELEM3(t->spacetype, SPACE_VIEW3D, SPACE_IMAGE, SPACE_NODE))
665                 return;
666         if (!(tc->mode & CON_APPLY))
667                 return;
668         if (t->flag & T_USES_MANIPULATOR)
669                 return;
670         if (t->flag & T_NO_CONSTRAINT)
671                 return;
672
673         /* nasty exception for Z constraint in camera view */
674         // TRANSFORM_FIX_ME
675 //      if ((t->flag & T_OBJECT) && G.vd->camera==OBACT && G.vd->persp==V3D_CAMOB)
676 //              return;
677
678         if (tc->drawExtra) {
679                 tc->drawExtra(t);
680         }
681         else {
682                 if (tc->mode & CON_SELECT) {
683                         float vec[3];
684                         char col2[3] = {255, 255, 255};
685                         int depth_test_enabled;
686
687                         convertViewVec(t, vec, (t->mval[0] - t->con.imval[0]), (t->mval[1] - t->con.imval[1]));
688                         add_v3_v3(vec, tc->center);
689
690                         drawLine(t, tc->center, tc->mtx[0], 'X', 0);
691                         drawLine(t, tc->center, tc->mtx[1], 'Y', 0);
692                         drawLine(t, tc->center, tc->mtx[2], 'Z', 0);
693
694                         glColor3ubv((GLubyte *)col2);
695
696                         depth_test_enabled = glIsEnabled(GL_DEPTH_TEST);
697                         if (depth_test_enabled)
698                                 glDisable(GL_DEPTH_TEST);
699
700                         setlinestyle(1);
701                         glBegin(GL_LINE_STRIP);
702                         glVertex3fv(tc->center);
703                         glVertex3fv(vec);
704                         glEnd();
705                         setlinestyle(0);
706
707                         if (depth_test_enabled)
708                                 glEnable(GL_DEPTH_TEST);
709                 }
710
711                 if (tc->mode & CON_AXIS0) {
712                         drawLine(t, tc->center, tc->mtx[0], 'X', DRAWLIGHT);
713                 }
714                 if (tc->mode & CON_AXIS1) {
715                         drawLine(t, tc->center, tc->mtx[1], 'Y', DRAWLIGHT);
716                 }
717                 if (tc->mode & CON_AXIS2) {
718                         drawLine(t, tc->center, tc->mtx[2], 'Z', DRAWLIGHT);
719                 }
720         }
721 }
722
723 /* called from drawview.c, as an extra per-window draw option */
724 void drawPropCircle(const struct bContext *C, TransInfo *t)
725 {
726         if (t->flag & T_PROP_EDIT) {
727                 RegionView3D *rv3d = CTX_wm_region_view3d(C);
728                 float tmat[4][4], imat[4][4];
729                 float center[3];
730                 int depth_test_enabled;
731
732                 UI_ThemeColor(TH_GRID);
733
734                 if (t->spacetype == SPACE_VIEW3D && rv3d != NULL) {
735                         copy_m4_m4(tmat, rv3d->viewmat);
736                         invert_m4_m4(imat, tmat);
737                 }
738                 else {
739                         unit_m4(tmat);
740                         unit_m4(imat);
741                 }
742
743                 glPushMatrix();
744
745                 copy_v3_v3(center, t->center);
746
747                 if ((t->spacetype == SPACE_VIEW3D) && t->obedit) {
748                         mul_m4_v3(t->obedit->obmat, center); /* because t->center is in local space */
749                 }
750                 else if (t->spacetype == SPACE_IMAGE) {
751                         float aspx, aspy;
752
753                         if (t->options & CTX_MASK) {
754                                 /* untested - mask aspect is TODO */
755                                 ED_space_image_get_aspect(t->sa->spacedata.first, &aspx, &aspy);
756                         }
757                         else {
758                                 ED_space_image_get_uv_aspect(t->sa->spacedata.first, &aspx, &aspy);
759                         }
760                         glScalef(1.0f / aspx, 1.0f / aspy, 1.0);
761                 }
762
763                 depth_test_enabled = glIsEnabled(GL_DEPTH_TEST);
764                 if (depth_test_enabled)
765                         glDisable(GL_DEPTH_TEST);
766
767                 set_inverted_drawing(1);
768                 drawcircball(GL_LINE_LOOP, center, t->prop_size, imat);
769                 set_inverted_drawing(0);
770
771                 if (depth_test_enabled)
772                         glEnable(GL_DEPTH_TEST);
773
774                 glPopMatrix();
775         }
776 }
777
778 static void drawObjectConstraint(TransInfo *t)
779 {
780         /* Draw the first one lighter because that's the one who controls the others.
781          * Meaning the transformation is projected on that one and just copied on the others
782          * constraint space.
783          * In a nutshell, the object with light axis is controlled by the user and the others follow.
784          * Without drawing the first light, users have little clue what they are doing.
785          */
786         short options = DRAWLIGHT;
787         TransData *td = t->data;
788         int i;
789         float tmp_axismtx[3][3];
790
791         for (i = 0; i < t->total; i++, td++) {
792                 float co[3];
793                 float (*axismtx)[3];
794
795                 if (t->flag & T_PROP_EDIT) {
796                         /* we're sorted, so skip the rest */
797                         if (td->factor == 0.0f) {
798                                 break;
799                         }
800                 }
801
802                 if (t->flag & T_OBJECT) {
803                         copy_v3_v3(co, td->ob->obmat[3]);
804                         axismtx = td->axismtx;
805                 }
806                 else if (t->flag & T_EDIT) {
807                         mul_v3_m4v3(co, t->obedit->obmat, td->center);
808
809                         mul_m3_m3m3(tmp_axismtx, t->obedit_mat, td->axismtx);
810                         axismtx = tmp_axismtx;
811                 }
812                 else if (t->flag & T_POSE) {
813                         mul_v3_m4v3(co, t->poseobj->obmat, td->center);
814                         axismtx = td->axismtx;
815                 }
816                 else {
817                         copy_v3_v3(co, td->center);
818                         axismtx = td->axismtx;
819                 }
820
821                 if (t->con.mode & CON_AXIS0) {
822                         drawLine(t, co, axismtx[0], 'X', options);
823                 }
824                 if (t->con.mode & CON_AXIS1) {
825                         drawLine(t, co, axismtx[1], 'Y', options);
826                 }
827                 if (t->con.mode & CON_AXIS2) {
828                         drawLine(t, co, axismtx[2], 'Z', options);
829                 }
830                 options &= ~DRAWLIGHT;
831         }
832 }
833
834 /*--------------------- START / STOP CONSTRAINTS ---------------------- */
835
836 void startConstraint(TransInfo *t)
837 {
838         t->con.mode |= CON_APPLY;
839         *t->con.text = ' ';
840         t->num.idx_max = min_ii(getConstraintSpaceDimension(t) - 1, t->idx_max);
841 }
842
843 void stopConstraint(TransInfo *t)
844 {
845         t->con.mode &= ~(CON_APPLY | CON_SELECT);
846         *t->con.text = '\0';
847         t->num.idx_max = t->idx_max;
848 }
849
850 void getConstraintMatrix(TransInfo *t)
851 {
852         float mat[3][3];
853         invert_m3_m3(t->con.imtx, t->con.mtx);
854         unit_m3(t->con.pmtx);
855
856         if (!(t->con.mode & CON_AXIS0)) {
857                 t->con.pmtx[0][0]       =
858                     t->con.pmtx[0][1]   =
859                     t->con.pmtx[0][2]   = 0.0f;
860         }
861
862         if (!(t->con.mode & CON_AXIS1)) {
863                 t->con.pmtx[1][0]       =
864                     t->con.pmtx[1][1]   =
865                     t->con.pmtx[1][2]   = 0.0f;
866         }
867
868         if (!(t->con.mode & CON_AXIS2)) {
869                 t->con.pmtx[2][0]       =
870                     t->con.pmtx[2][1]   =
871                     t->con.pmtx[2][2]   = 0.0f;
872         }
873
874         mul_m3_m3m3(mat, t->con.pmtx, t->con.imtx);
875         mul_m3_m3m3(t->con.pmtx, t->con.mtx, mat);
876 }
877
878 /*------------------------- MMB Select -------------------------------*/
879
880 void initSelectConstraint(TransInfo *t, float mtx[3][3])
881 {
882         copy_m3_m3(t->con.mtx, mtx);
883         t->con.mode |= CON_APPLY;
884         t->con.mode |= CON_SELECT;
885
886         setNearestAxis(t);
887         t->con.drawExtra = NULL;
888         t->con.applyVec = applyAxisConstraintVec;
889         t->con.applySize = applyAxisConstraintSize;
890         t->con.applyRot = applyAxisConstraintRot;
891 }
892
893 void selectConstraint(TransInfo *t)
894 {
895         if (t->con.mode & CON_SELECT) {
896                 setNearestAxis(t);
897                 startConstraint(t);
898         }
899 }
900
901 void postSelectConstraint(TransInfo *t)
902 {
903         if (!(t->con.mode & CON_SELECT))
904                 return;
905
906         t->con.mode &= ~CON_AXIS0;
907         t->con.mode &= ~CON_AXIS1;
908         t->con.mode &= ~CON_AXIS2;
909         t->con.mode &= ~CON_SELECT;
910
911         setNearestAxis(t);
912
913         startConstraint(t);
914         t->redraw = TREDRAW_HARD;
915 }
916
917 static void setNearestAxis2d(TransInfo *t)
918 {
919         /* no correction needed... just use whichever one is lower */
920         if (abs(t->mval[0] - t->con.imval[0]) < abs(t->mval[1] - t->con.imval[1]) ) {
921                 t->con.mode |= CON_AXIS1;
922                 BLI_strncpy(t->con.text, IFACE_(" along Y axis"), sizeof(t->con.text));
923         }
924         else {
925                 t->con.mode |= CON_AXIS0;
926                 BLI_strncpy(t->con.text, IFACE_(" along X axis"), sizeof(t->con.text));
927         }
928 }
929
930 static void setNearestAxis3d(TransInfo *t)
931 {
932         float zfac;
933         float mvec[3], proj[3];
934         float len[3];
935         int i;
936
937         /* calculate mouse movement */
938         mvec[0] = (float)(t->mval[0] - t->con.imval[0]);
939         mvec[1] = (float)(t->mval[1] - t->con.imval[1]);
940         mvec[2] = 0.0f;
941
942         /* we need to correct axis length for the current zoomlevel of view,
943          * this to prevent projected values to be clipped behind the camera
944          * and to overflow the short integers.
945          * The formula used is a bit stupid, just a simplification of the subtraction
946          * of two 2D points 30 pixels apart (that's the last factor in the formula) after
947          * projecting them with ED_view3d_win_to_delta and then get the length of that vector.
948          */
949         zfac = mul_project_m4_v3_zfac(t->persmat, t->center);
950         zfac = len_v3(t->persinv[0]) * 2.0f / t->ar->winx * zfac * 30.0f;
951
952         for (i = 0; i < 3; i++) {
953                 float axis[3], axis_2d[2];
954
955                 copy_v3_v3(axis, t->con.mtx[i]);
956
957                 mul_v3_fl(axis, zfac);
958                 /* now we can project to get window coordinate */
959                 add_v3_v3(axis, t->con.center);
960                 projectFloatView(t, axis, axis_2d);
961
962                 sub_v2_v2v2(axis, axis_2d, t->center2d);
963                 axis[2] = 0.0f;
964
965                 if (normalize_v3(axis) != 0.0f) {
966                         project_v3_v3v3(proj, mvec, axis);
967                         sub_v3_v3v3(axis, mvec, proj);
968                         len[i] = normalize_v3(axis);
969                 }
970                 else {
971                         len[i] = 10000000000.0f;
972                 }
973         }
974
975         if (len[0] <= len[1] && len[0] <= len[2]) {
976                 if (t->modifiers & MOD_CONSTRAINT_PLANE) {
977                         t->con.mode |= (CON_AXIS1 | CON_AXIS2);
978                         BLI_snprintf(t->con.text, sizeof(t->con.text), IFACE_(" locking %s X axis"), t->spacename);
979                 }
980                 else {
981                         t->con.mode |= CON_AXIS0;
982                         BLI_snprintf(t->con.text, sizeof(t->con.text), IFACE_(" along %s X axis"), t->spacename);
983                 }
984         }
985         else if (len[1] <= len[0] && len[1] <= len[2]) {
986                 if (t->modifiers & MOD_CONSTRAINT_PLANE) {
987                         t->con.mode |= (CON_AXIS0 | CON_AXIS2);
988                         BLI_snprintf(t->con.text, sizeof(t->con.text), IFACE_(" locking %s Y axis"), t->spacename);
989                 }
990                 else {
991                         t->con.mode |= CON_AXIS1;
992                         BLI_snprintf(t->con.text, sizeof(t->con.text), IFACE_(" along %s Y axis"), t->spacename);
993                 }
994         }
995         else if (len[2] <= len[1] && len[2] <= len[0]) {
996                 if (t->modifiers & MOD_CONSTRAINT_PLANE) {
997                         t->con.mode |= (CON_AXIS0 | CON_AXIS1);
998                         BLI_snprintf(t->con.text, sizeof(t->con.text), IFACE_(" locking %s Z axis"), t->spacename);
999                 }
1000                 else {
1001                         t->con.mode |= CON_AXIS2;
1002                         BLI_snprintf(t->con.text, sizeof(t->con.text), IFACE_(" along %s Z axis"), t->spacename);
1003                 }
1004         }
1005 }
1006
1007 void setNearestAxis(TransInfo *t)
1008 {
1009         /* clear any prior constraint flags */
1010         t->con.mode &= ~CON_AXIS0;
1011         t->con.mode &= ~CON_AXIS1;
1012         t->con.mode &= ~CON_AXIS2;
1013
1014         /* constraint setting - depends on spacetype */
1015         if (t->spacetype == SPACE_VIEW3D) {
1016                 /* 3d-view */
1017                 setNearestAxis3d(t);
1018         }
1019         else {
1020                 /* assume that this means a 2D-Editor */
1021                 setNearestAxis2d(t);
1022         }
1023
1024         getConstraintMatrix(t);
1025 }
1026
1027 /*-------------- HELPER FUNCTIONS ----------------*/
1028
1029 char constraintModeToChar(TransInfo *t)
1030 {
1031         if ((t->con.mode & CON_APPLY) == 0) {
1032                 return '\0';
1033         }
1034         switch (t->con.mode & (CON_AXIS0 | CON_AXIS1 | CON_AXIS2)) {
1035                 case (CON_AXIS0):
1036                 case (CON_AXIS1 | CON_AXIS2):
1037                         return 'X';
1038                 case (CON_AXIS1):
1039                 case (CON_AXIS0 | CON_AXIS2):
1040                         return 'Y';
1041                 case (CON_AXIS2):
1042                 case (CON_AXIS0 | CON_AXIS1):
1043                         return 'Z';
1044                 default:
1045                         return '\0';
1046         }
1047 }
1048
1049
1050 bool isLockConstraint(TransInfo *t)
1051 {
1052         int mode = t->con.mode;
1053
1054         if ((mode & (CON_AXIS0 | CON_AXIS1)) == (CON_AXIS0 | CON_AXIS1))
1055                 return true;
1056
1057         if ((mode & (CON_AXIS1 | CON_AXIS2)) == (CON_AXIS1 | CON_AXIS2))
1058                 return true;
1059
1060         if ((mode & (CON_AXIS0 | CON_AXIS2)) == (CON_AXIS0 | CON_AXIS2))
1061                 return true;
1062
1063         return false;
1064 }
1065
1066 /*
1067  * Returns the dimension of the constraint space.
1068  *
1069  * For that reason, the flags always needs to be set to properly evaluate here,
1070  * even if they aren't actually used in the callback function. (Which could happen
1071  * for weird constraints not yet designed. Along a path for example.)
1072  */
1073
1074 int getConstraintSpaceDimension(TransInfo *t)
1075 {
1076         int n = 0;
1077
1078         if (t->con.mode & CON_AXIS0)
1079                 n++;
1080
1081         if (t->con.mode & CON_AXIS1)
1082                 n++;
1083
1084         if (t->con.mode & CON_AXIS2)
1085                 n++;
1086
1087         return n;
1088 /*
1089  * Someone willing to do it cryptically could do the following instead:
1090  *
1091  * return t->con & (CON_AXIS0|CON_AXIS1|CON_AXIS2);
1092  *
1093  * Based on the assumptions that the axis flags are one after the other and start at 1
1094  */
1095 }