336ff5e1023f6b458a7721fc7296ae450d53463a
[blender.git] / source / blender / editors / curve / editcurve.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/curve/editcurve.c
29  *  \ingroup edcurve
30  */
31
32 #include <math.h>
33 #include <string.h>
34
35 #ifndef WIN32
36 #include <unistd.h>
37 #else
38 #include <io.h>
39 #endif
40 #include <stdlib.h>
41
42 #include "DNA_key_types.h"
43 #include "DNA_object_types.h"
44 #include "DNA_scene_types.h"
45 #include "DNA_anim_types.h"
46
47 #include "MEM_guardedalloc.h"
48
49 #include "BLI_blenlib.h"
50 #include "BLI_math.h"
51 #include "BLI_dynstr.h"
52 #include "BLI_rand.h"
53 #include "BLI_utildefines.h"
54 #include "BLI_ghash.h"
55
56 #include "BKE_context.h"
57 #include "BKE_curve.h"
58 #include "BKE_depsgraph.h"
59 #include "BKE_fcurve.h"
60 #include "BKE_global.h"
61 #include "BKE_key.h"
62 #include "BKE_library.h"
63 #include "BKE_main.h"
64 #include "BKE_report.h"
65 #include "BKE_animsys.h"
66 #include "BKE_action.h"
67
68 #include "WM_api.h"
69 #include "WM_types.h"
70
71 #include "ED_keyframes_edit.h"
72 #include "ED_object.h"
73 #include "ED_screen.h"
74 #include "ED_transform.h"
75 #include "ED_types.h"
76 #include "ED_util.h"
77 #include "ED_view3d.h"
78 #include "ED_curve.h"
79
80 #include "curve_intern.h"
81
82 #include "UI_interface.h"
83 #include "UI_resources.h"
84
85 #include "RNA_access.h"
86 #include "RNA_define.h"
87 #include "RNA_enum_types.h"
88
89 /* Undo stuff */
90 typedef struct {
91         ListBase nubase;
92         void *lastsel;
93         GHash *undoIndex;
94         ListBase fcurves, drivers;
95         int actnu;
96 } UndoCurve;
97
98 /* Definitions needed for shape keys */
99 typedef struct {
100         void *orig_cv;
101         int key_index, nu_index, pt_index;
102         int switched;
103         Nurb *orig_nu;
104 } CVKeyIndex;
105
106 void selectend_nurb(Object *obedit, short selfirst, short doswap, short selstatus);
107 static void select_adjacent_cp(ListBase *editnurb, short next, short cont, short selstatus);
108
109 /* still need to eradicate a few :( */
110 #define CALLOC_STRUCT_N(x, y, name) (x *)MEM_callocN((y) * sizeof(x), name)
111
112 static float nurbcircle[8][2] = {
113         {0.0, -1.0}, {-1.0, -1.0}, {-1.0, 0.0}, {-1.0,  1.0},
114         {0.0,  1.0}, { 1.0,  1.0}, { 1.0, 0.0}, { 1.0, -1.0}
115 };
116
117 ListBase *object_editcurve_get(Object *ob)
118 {
119         if (ob && ELEM(ob->type, OB_CURVE, OB_SURF)) {
120                 Curve *cu = ob->data;
121                 return &cu->editnurb->nurbs;
122         }
123         return NULL;
124 }
125
126 /* this replaces the active flag used in uv/face mode */
127 static void set_actNurb(Object *obedit, Nurb *nu)
128 {
129         Curve *cu = obedit->data;
130         
131         if (nu == NULL)
132                 cu->actnu = -1;
133         else {
134                 ListBase *nurbs = BKE_curve_editNurbs_get(cu);
135                 cu->actnu = BLI_findindex(nurbs, nu);
136         }
137 }
138
139 static Nurb *get_actNurb(Object *obedit)
140 {
141         Curve *cu = obedit->data;
142         ListBase *nurbs = BKE_curve_editNurbs_get(cu);
143
144         return BLI_findlink(nurbs, cu->actnu);
145 }
146
147 /* ******************* SELECTION FUNCTIONS ********************* */
148
149 #define HIDDEN          1
150 #define VISIBLE         0
151
152 #define FIRST           1
153 #define LAST            0
154
155
156 /* returns 1 in case (de)selection was successful */
157 static short select_beztriple(BezTriple *bezt, short selstatus, short flag, short hidden)
158 {       
159         if (bezt) {
160                 if ((bezt->hide == 0) || (hidden == 1)) {
161                         if (selstatus == 1) { /* selects */
162                                 bezt->f1 |= flag;
163                                 bezt->f2 |= flag;
164                                 bezt->f3 |= flag;
165                                 return 1;
166                         }
167                         else { /* deselects */
168                                 bezt->f1 &= ~flag; 
169                                 bezt->f2 &= ~flag; 
170                                 bezt->f3 &= ~flag; 
171                                 return 1;
172                         }
173                 }
174         }
175         
176         return 0;
177 }
178
179 /* returns 1 in case (de)selection was successful */
180 static short select_bpoint(BPoint *bp, short selstatus, short flag, short hidden) 
181 {       
182         if (bp) {
183                 if ((bp->hide == 0) || (hidden == 1)) {
184                         if (selstatus == 1) {
185                                 bp->f1 |= flag;
186                                 return 1;
187                         }
188                         else {
189                                 bp->f1 &= ~flag;
190                                 return 1;
191                         }
192                 }
193         }
194
195         return 0;
196 }
197
198 static short swap_selection_beztriple(BezTriple *bezt)
199 {
200         if (bezt->f2 & SELECT)
201                 return select_beztriple(bezt, DESELECT, 1, VISIBLE);
202         else
203                 return select_beztriple(bezt, SELECT, 1, VISIBLE);
204 }
205
206 static short swap_selection_bpoint(BPoint *bp)
207 {
208         if (bp->f1 & SELECT)
209                 return select_bpoint(bp, DESELECT, 1, VISIBLE);
210         else
211                 return select_bpoint(bp, SELECT, 1, VISIBLE);
212 }
213
214 int isNurbsel(Nurb *nu)
215 {
216         BezTriple *bezt;
217         BPoint *bp;
218         int a;
219
220         if (nu->type == CU_BEZIER) {
221                 bezt = nu->bezt;
222                 a = nu->pntsu;
223                 while (a--) {
224                         if ( (bezt->f1 & SELECT) || (bezt->f2 & SELECT) || (bezt->f3 & SELECT) ) return 1;
225                         bezt++;
226                 }
227         }
228         else {
229                 bp = nu->bp;
230                 a = nu->pntsu * nu->pntsv;
231                 while (a--) {
232                         if (bp->f1 & SELECT) return 1;
233                         bp++;
234                 }
235         }
236         return 0;
237 }
238
239 static int isNurbsel_count(Curve *cu, Nurb *nu)
240 {
241         BezTriple *bezt;
242         BPoint *bp;
243         int a, sel = 0;
244
245         if (nu->type == CU_BEZIER) {
246                 bezt = nu->bezt;
247                 a = nu->pntsu;
248                 while (a--) {
249                         if (BEZSELECTED_HIDDENHANDLES(cu, bezt)) sel++;
250                         bezt++;
251                 }
252         }
253         else {
254                 bp = nu->bp;
255                 a = nu->pntsu * nu->pntsv;
256                 while (a--) {
257                         if ( (bp->f1 & SELECT) ) sel++;
258                         bp++;
259                 }
260         }
261         return sel;
262 }
263
264 /* ******************* PRINTS ********************* */
265
266 void printknots(Object *obedit)
267 {
268         ListBase *editnurb = object_editcurve_get(obedit);
269         Nurb *nu;
270         int a, num;
271
272         for (nu = editnurb->first; nu; nu = nu->next) {
273                 if (isNurbsel(nu) && nu->type == CU_NURBS) {
274                         if (nu->knotsu) {
275                                 num = KNOTSU(nu);
276                                 for (a = 0; a < num; a++) printf("knotu %d: %f\n", a, nu->knotsu[a]);
277                         }
278                         if (nu->knotsv) {
279                                 num = KNOTSV(nu);
280                                 for (a = 0; a < num; a++) printf("knotv %d: %f\n", a, nu->knotsv[a]);
281                         }
282                 }
283         }
284 }
285
286 /* ********************* Shape keys *************** */
287
288 static CVKeyIndex *init_cvKeyIndex(void *cv, int key_index, int nu_index, int pt_index, Nurb *orig_nu)
289 {
290         CVKeyIndex *cvIndex = MEM_callocN(sizeof(CVKeyIndex), "init_cvKeyIndex");
291
292         cvIndex->orig_cv = cv;
293         cvIndex->key_index = key_index;
294         cvIndex->nu_index = nu_index;
295         cvIndex->pt_index = pt_index;
296         cvIndex->switched = 0;
297         cvIndex->orig_nu = orig_nu;
298
299         return cvIndex;
300 }
301
302 static void init_editNurb_keyIndex(EditNurb *editnurb, ListBase *origBase)
303 {
304         Nurb *nu = editnurb->nurbs.first;
305         Nurb *orignu = origBase->first;
306         GHash *gh;
307         BezTriple *bezt, *origbezt;
308         BPoint *bp, *origbp;
309         CVKeyIndex *keyIndex;
310         int a, key_index = 0, nu_index = 0, pt_index = 0;
311
312         if (editnurb->keyindex) return;
313
314         gh = BLI_ghash_ptr_new("editNurb keyIndex");
315
316         while (orignu) {
317                 if (orignu->bezt) {
318                         a = orignu->pntsu;
319                         bezt = nu->bezt;
320                         origbezt = orignu->bezt;
321                         pt_index = 0;
322                         while (a--) {
323                                 keyIndex = init_cvKeyIndex(origbezt, key_index, nu_index, pt_index, orignu);
324                                 BLI_ghash_insert(gh, bezt, keyIndex);
325                                 key_index += 12;
326                                 bezt++;
327                                 origbezt++;
328                                 pt_index++;
329                         }
330                 }
331                 else {
332                         a = orignu->pntsu * orignu->pntsv;
333                         bp = nu->bp;
334                         origbp = orignu->bp;
335                         pt_index = 0;
336                         while (a--) {
337                                 keyIndex = init_cvKeyIndex(origbp, key_index, nu_index, pt_index, orignu);
338                                 BLI_ghash_insert(gh, bp, keyIndex);
339                                 key_index += 4;
340                                 bp++;
341                                 origbp++;
342                                 pt_index++;
343                         }
344                 }
345
346                 nu = nu->next;
347                 orignu = orignu->next;
348                 nu_index++;
349         }
350
351         editnurb->keyindex = gh;
352 }
353
354 static CVKeyIndex *getCVKeyIndex(EditNurb *editnurb, void *cv)
355 {
356         return BLI_ghash_lookup(editnurb->keyindex, cv);
357 }
358
359 static CVKeyIndex *popCVKeyIndex(EditNurb *editnurb, void *cv)
360 {
361         return BLI_ghash_pop(editnurb->keyindex, cv, NULL);
362 }
363
364 static BezTriple *getKeyIndexOrig_bezt(EditNurb *editnurb, BezTriple *bezt)
365 {
366         CVKeyIndex *index = getCVKeyIndex(editnurb, bezt);
367
368         if (!index) {
369                 return NULL;
370         }
371
372         return (BezTriple *)index->orig_cv;
373 }
374
375 static BPoint *getKeyIndexOrig_bp(EditNurb *editnurb, BPoint *bp)
376 {
377         CVKeyIndex *index = getCVKeyIndex(editnurb, bp);
378
379         if (!index) {
380                 return NULL;
381         }
382
383         return (BPoint *)index->orig_cv;
384 }
385
386 static int getKeyIndexOrig_keyIndex(EditNurb *editnurb, void *cv)
387 {
388         CVKeyIndex *index = getCVKeyIndex(editnurb, cv);
389
390         if (!index) {
391                 return -1;
392         }
393
394         return index->key_index;
395 }
396
397 static void keyIndex_delCV(EditNurb *editnurb, void *cv)
398 {
399         if (!editnurb->keyindex) {
400                 return;
401         }
402
403         BLI_ghash_remove(editnurb->keyindex, cv, NULL, (GHashValFreeFP)MEM_freeN);
404 }
405
406 static void keyIndex_delBezt(EditNurb *editnurb, BezTriple *bezt)
407 {
408         keyIndex_delCV(editnurb, bezt);
409 }
410
411 static void keyIndex_delBP(EditNurb *editnurb, BPoint *bp)
412 {
413         keyIndex_delCV(editnurb, bp);
414 }
415
416 static void keyIndex_delNurb(EditNurb *editnurb, Nurb *nu)
417 {
418         int a;
419
420         if (!editnurb->keyindex) {
421                 return;
422         }
423
424         if (nu->bezt) {
425                 BezTriple *bezt = nu->bezt;
426                 a = nu->pntsu;
427
428                 while (a--) {
429                         BLI_ghash_remove(editnurb->keyindex, bezt, NULL, (GHashValFreeFP)MEM_freeN);
430                         bezt++;
431                 }
432         }
433         else {
434                 BPoint *bp = nu->bp;
435                 a = nu->pntsu * nu->pntsv;
436
437                 while (a--) {
438                         BLI_ghash_remove(editnurb->keyindex, bp, NULL, (GHashValFreeFP)MEM_freeN);
439                         bp++;
440                 }
441         }
442 }
443
444 static void keyIndex_delNurbList(EditNurb *editnurb, ListBase *nubase)
445 {
446         Nurb *nu = nubase->first;
447
448         while (nu) {
449                 keyIndex_delNurb(editnurb, nu);
450
451                 nu = nu->next;
452         }
453 }
454
455 static void keyIndex_updateCV(EditNurb *editnurb, char *cv,
456                               char *newcv, int count, int size)
457 {
458         int i;
459         CVKeyIndex *index;
460
461         if (editnurb->keyindex == NULL) {
462                 /* No shape keys - updating not needed */
463                 return;
464         }
465
466         for (i = 0; i < count; i++) {
467                 index = popCVKeyIndex(editnurb, cv);
468
469                 if (index) {
470                         BLI_ghash_insert(editnurb->keyindex, newcv, index);
471                 }
472
473                 newcv += size;
474                 cv += size;
475         }
476 }
477
478 static void keyIndex_updateBezt(EditNurb *editnurb, BezTriple *bezt,
479                                 BezTriple *newbezt, int count)
480 {
481         keyIndex_updateCV(editnurb, (char *)bezt, (char *)newbezt, count, sizeof(BezTriple));
482 }
483
484 static void keyIndex_updateBP(EditNurb *editnurb, BPoint *bp,
485                               BPoint *newbp, int count)
486 {
487         keyIndex_updateCV(editnurb, (char *)bp, (char *)newbp, count, sizeof(BPoint));
488 }
489
490 static void keyIndex_updateNurb(EditNurb *editnurb, Nurb *nu, Nurb *newnu)
491 {
492         if (nu->bezt) {
493                 keyIndex_updateBezt(editnurb, nu->bezt, newnu->bezt, newnu->pntsu);
494         }
495         else {
496                 keyIndex_updateBP(editnurb, nu->bp, newnu->bp, newnu->pntsu * newnu->pntsv);
497         }
498 }
499
500 static void keyIndex_swap(EditNurb *editnurb, void *a, void *b)
501 {
502         CVKeyIndex *index1 = popCVKeyIndex(editnurb, a);
503         CVKeyIndex *index2 = popCVKeyIndex(editnurb, b);
504
505         if (index2) BLI_ghash_insert(editnurb->keyindex, a, index2);
506         if (index1) BLI_ghash_insert(editnurb->keyindex, b, index1);
507 }
508
509 static void keyIndex_switchDirection(EditNurb *editnurb, Nurb *nu)
510 {
511         int a;
512         CVKeyIndex *index1, *index2;
513
514         if (nu->bezt) {
515                 BezTriple *bezt1, *bezt2;
516
517                 a = nu->pntsu;
518
519                 bezt1 = nu->bezt;
520                 bezt2 = bezt1 + (a - 1);
521
522                 if (a & 1) ++a;
523
524                 a /= 2;
525
526                 while (a--) {
527                         index1 = getCVKeyIndex(editnurb, bezt1);
528                         index2 = getCVKeyIndex(editnurb, bezt2);
529
530                         if (index1) index1->switched = !index1->switched;
531
532                         if (bezt1 != bezt2) {
533                                 keyIndex_swap(editnurb, bezt1, bezt2);
534
535                                 if (index2) index2->switched = !index2->switched;
536                         }
537
538                         bezt1++;
539                         bezt2--;
540                 }
541         }
542         else {
543                 BPoint *bp1, *bp2;
544
545                 if (nu->pntsv == 1) {
546                         a = nu->pntsu;
547                         bp1 = nu->bp;
548                         bp2 = bp1 + (a - 1);
549                         a /= 2;
550                         while (bp1 != bp2 && a > 0) {
551                                 index1 = getCVKeyIndex(editnurb, bp1);
552                                 index2 = getCVKeyIndex(editnurb, bp2);
553
554                                 if (index1) index1->switched = !index1->switched;
555
556                                 if (bp1 != bp2) {
557                                         if (index2) index2->switched = !index2->switched;
558
559                                         keyIndex_swap(editnurb, bp1, bp2);
560                                 }
561
562                                 a--;
563                                 bp1++;
564                                 bp2--;
565                         }
566                 }
567                 else {
568                         int b;
569
570                         for (b = 0; b < nu->pntsv; b++) {
571
572                                 bp1 = &nu->bp[b * nu->pntsu];
573                                 a = nu->pntsu;
574                                 bp2 = bp1 + (a - 1);
575                                 a /= 2;
576
577                                 while (bp1 != bp2 && a > 0) {
578                                         index1 = getCVKeyIndex(editnurb, bp1);
579                                         index2 = getCVKeyIndex(editnurb, bp2);
580
581                                         if (index1) index1->switched = !index1->switched;
582
583                                         if (bp1 != bp2) {
584                                                 if (index2) index2->switched = !index2->switched;
585
586                                                 keyIndex_swap(editnurb, bp1, bp2);
587                                         }
588
589                                         a--;
590                                         bp1++;
591                                         bp2--;
592                                 }
593                         }
594
595                 }
596         }
597 }
598
599 static void switch_keys_direction(Curve *cu, Nurb *actnu)
600 {
601         KeyBlock *currkey;
602         EditNurb *editnurb = cu->editnurb;
603         ListBase *nubase = &editnurb->nurbs;
604         Nurb *nu;
605         float *fp;
606         int a;
607
608         currkey = cu->key->block.first;
609         while (currkey) {
610                 fp = currkey->data;
611
612                 nu = nubase->first;
613                 while (nu) {
614                         if (nu->bezt) {
615                                 BezTriple *bezt = nu->bezt;
616                                 a = nu->pntsu;
617                                 if (nu == actnu) {
618                                         while (a--) {
619                                                 if (getKeyIndexOrig_bezt(editnurb, bezt)) {
620                                                         swap_v3_v3(fp, fp + 6);
621                                                         *(fp + 9) = -*(fp + 9);
622                                                         fp += 12;
623                                                 }
624                                                 bezt++;
625                                         }
626                                 }
627                                 else fp += a * 12;
628                         }
629                         else {
630                                 BPoint *bp = nu->bp;
631                                 a = nu->pntsu * nu->pntsv;
632                                 if (nu == actnu) {
633                                         while (a--) {
634                                                 if (getKeyIndexOrig_bp(editnurb, bp)) {
635                                                         *(fp + 3) = -*(fp + 3);
636                                                         fp += 4;
637                                                 }
638                                                 bp++;
639                                         }
640                                 }
641                                 else fp += a * 4;
642                         }
643
644                         nu = nu->next;
645                 }
646
647                 currkey = currkey->next;
648         }
649 }
650
651 static void keyData_switchDirectionNurb(Curve *cu, Nurb *nu)
652 {
653         EditNurb *editnurb = cu->editnurb;
654
655         if (!editnurb->keyindex) {
656                 /* no shape keys - nothing to do */
657                 return;
658         }
659
660         keyIndex_switchDirection(editnurb, nu);
661         if (cu->key)
662                 switch_keys_direction(cu, nu);
663 }
664
665 static GHash *dupli_keyIndexHash(GHash *keyindex)
666 {
667         GHash *gh;
668         GHashIterator *hashIter;
669
670         gh = BLI_ghash_ptr_new("dupli_keyIndex gh");
671
672         for (hashIter = BLI_ghashIterator_new(keyindex);
673              !BLI_ghashIterator_isDone(hashIter);
674              BLI_ghashIterator_step(hashIter))
675         {
676                 void *cv = BLI_ghashIterator_getKey(hashIter);
677                 CVKeyIndex *index = BLI_ghashIterator_getValue(hashIter);
678                 CVKeyIndex *newIndex = MEM_callocN(sizeof(CVKeyIndex), "dupli_keyIndexHash index");
679
680                 memcpy(newIndex, index, sizeof(CVKeyIndex));
681
682                 BLI_ghash_insert(gh, cv, newIndex);
683         }
684
685         BLI_ghashIterator_free(hashIter);
686
687         return gh;
688 }
689
690 static void key_to_bezt(float *key, BezTriple *basebezt, BezTriple *bezt)
691 {
692         memcpy(bezt, basebezt, sizeof(BezTriple));
693         memcpy(bezt->vec, key, sizeof(float) * 9);
694         bezt->alfa = key[9];
695 }
696
697 static void bezt_to_key(BezTriple *bezt, float *key)
698 {
699         memcpy(key, bezt->vec, sizeof(float) * 9);
700         key[9] = bezt->alfa;
701 }
702
703 static void calc_keyHandles(ListBase *nurb, float *key)
704 {
705         Nurb *nu;
706         int a;
707         float *fp = key;
708         BezTriple *bezt;
709
710         nu = nurb->first;
711         while (nu) {
712                 if (nu->bezt) {
713                         BezTriple *prevp, *nextp;
714                         BezTriple cur, prev, next;
715                         float *startfp, *prevfp, *nextfp;
716
717                         bezt = nu->bezt;
718                         a = nu->pntsu;
719                         startfp = fp;
720
721                         if (nu->flagu & CU_NURB_CYCLIC) {
722                                 prevp = bezt + (a - 1);
723                                 prevfp = fp + (12 * (a - 1));
724                         }
725                         else {
726                                 prevp = NULL;
727                                 prevfp = NULL;
728                         }
729
730                         nextp = bezt + 1;
731                         nextfp = fp + 12;
732
733                         while (a--) {
734                                 key_to_bezt(fp, bezt, &cur);
735
736                                 if (nextp) key_to_bezt(nextfp, nextp, &next);
737                                 if (prevp) key_to_bezt(prevfp, prevp, &prev);
738
739                                 BKE_nurb_handle_calc(&cur, prevp ? &prev : NULL, nextp ? &next : NULL, 0);
740                                 bezt_to_key(&cur, fp);
741
742                                 prevp = bezt;
743                                 prevfp = fp;
744                                 if (a == 1) {
745                                         if (nu->flagu & CU_NURB_CYCLIC) {
746                                                 nextp = nu->bezt;
747                                                 nextfp = startfp;
748                                         }
749                                         else {
750                                                 nextp = NULL;
751                                                 nextfp = NULL;
752                                         }
753                                 }
754                                 else {
755                                         nextp++;
756                                         nextfp += 12;
757                                 }
758
759                                 bezt++;
760                                 fp += 12;
761                         }
762                 }
763                 else {
764                         a = nu->pntsu * nu->pntsv;
765                         fp += a * 4;
766                 }
767
768                 nu = nu->next;
769         }
770 }
771
772 static void calc_shapeKeys(Object *obedit)
773 {
774         Curve *cu = (Curve *)obedit->data;
775
776         /* are there keys? */
777         if (cu->key) {
778                 int a, i;
779                 EditNurb *editnurb = cu->editnurb;
780                 KeyBlock *currkey;
781                 KeyBlock *actkey = BLI_findlink(&cu->key->block, editnurb->shapenr - 1);
782                 BezTriple *bezt, *oldbezt;
783                 BPoint *bp, *oldbp;
784                 Nurb *nu;
785                 int totvert = BKE_nurbList_verts_count(&editnurb->nurbs);
786
787                 float (*ofs)[3] = NULL;
788                 float *oldkey, *newkey, *ofp;
789
790                 /* editing the base key should update others */
791                 if (cu->key->type == KEY_RELATIVE) {
792                         int act_is_basis = 0;
793                         /* find if this key is a basis for any others */
794                         for (currkey = cu->key->block.first; currkey; currkey = currkey->next) {
795                                 if (editnurb->shapenr - 1 == currkey->relative) {
796                                         act_is_basis = 1;
797                                         break;
798                                 }
799                         }
800
801                         if (act_is_basis) { /* active key is a base */
802                                 int totvec = 0;
803
804                                 /* Calculate needed memory to store offset */
805                                 nu = editnurb->nurbs.first;
806                                 while (nu) {
807                                         if (nu->bezt) {
808                                                 /* Three vects to store handles and one for alfa */
809                                                 totvec += nu->pntsu * 4;
810                                         }
811                                         else {
812                                                 totvec += 2 * nu->pntsu * nu->pntsv;
813                                         }
814
815                                         nu = nu->next;
816                                 }
817
818                                 ofs = MEM_callocN(sizeof(float) * 3 * totvec,  "currkey->data");
819                                 nu = editnurb->nurbs.first;
820                                 i = 0;
821                                 while (nu) {
822                                         if (nu->bezt) {
823                                                 bezt = nu->bezt;
824                                                 a = nu->pntsu;
825                                                 while (a--) {
826                                                         oldbezt = getKeyIndexOrig_bezt(editnurb, bezt);
827
828                                                         if (oldbezt) {
829                                                                 int j;
830                                                                 for (j = 0; j < 3; ++j) {
831                                                                         sub_v3_v3v3(ofs[i], bezt->vec[j], oldbezt->vec[j]);
832                                                                         i++;
833                                                                 }
834                                                                 ofs[i++][0] = bezt->alfa - oldbezt->alfa;
835                                                         }
836                                                         else {
837                                                                 i += 4;
838                                                         }
839                                                         bezt++;
840                                                 }
841                                         }
842                                         else {
843                                                 bp = nu->bp;
844                                                 a = nu->pntsu * nu->pntsv;
845                                                 while (a--) {
846                                                         oldbp = getKeyIndexOrig_bp(editnurb, bp);
847                                                         if (oldbp) {
848                                                                 sub_v3_v3v3(ofs[i], bp->vec, oldbp->vec);
849                                                                 ofs[i + 1][0] = bp->alfa - oldbp->alfa;
850                                                         }
851                                                         i += 2;
852                                                         bp++;
853                                                 }
854                                         }
855
856                                         nu = nu->next;
857                                 }
858                         }
859                 }
860
861                 currkey = cu->key->block.first;
862                 while (currkey) {
863                         int apply_offset = (ofs && (currkey != actkey) && (editnurb->shapenr - 1 == currkey->relative));
864
865                         float *fp = newkey = MEM_callocN(cu->key->elemsize * totvert,  "currkey->data");
866                         ofp = oldkey = currkey->data;
867
868                         nu = editnurb->nurbs.first;
869                         i = 0;
870                         while (nu) {
871                                 if (currkey == actkey) {
872                                         int restore = actkey != cu->key->refkey;
873
874                                         if (nu->bezt) {
875                                                 bezt = nu->bezt;
876                                                 a = nu->pntsu;
877                                                 while (a--) {
878                                                         int j;
879                                                         oldbezt = getKeyIndexOrig_bezt(editnurb, bezt);
880
881                                                         for (j = 0; j < 3; ++j, ++i) {
882                                                                 copy_v3_v3(fp, bezt->vec[j]);
883
884                                                                 if (restore && oldbezt) {
885                                                                         copy_v3_v3(bezt->vec[j], oldbezt->vec[j]);
886                                                                 }
887
888                                                                 fp += 3;
889                                                         }
890                                                         fp[0] = bezt->alfa;
891
892                                                         if (restore && oldbezt) {
893                                                                 bezt->alfa = oldbezt->alfa;
894                                                         }
895
896                                                         fp += 3; ++i; /* alphas */
897                                                         bezt++;
898                                                 }
899                                         }
900                                         else {
901                                                 bp = nu->bp;
902                                                 a = nu->pntsu * nu->pntsv;
903                                                 while (a--) {
904                                                         oldbp = getKeyIndexOrig_bp(editnurb, bp);
905
906                                                         copy_v3_v3(fp, bp->vec);
907
908                                                         fp[3] = bp->alfa;
909
910                                                         if (restore && oldbp) {
911                                                                 copy_v3_v3(bp->vec, oldbp->vec);
912                                                                 bp->alfa = oldbp->alfa;
913                                                         }
914
915                                                         fp += 4;
916                                                         bp++;
917                                                         i += 2;
918                                                 }
919                                         }
920                                 }
921                                 else {
922                                         int index;
923                                         float *curofp;
924
925                                         if (oldkey) {
926                                                 if (nu->bezt) {
927                                                         bezt = nu->bezt;
928                                                         a = nu->pntsu;
929
930                                                         while (a--) {
931                                                                 index = getKeyIndexOrig_keyIndex(editnurb, bezt);
932                                                                 if (index >= 0) {
933                                                                         int j;
934                                                                         curofp = ofp + index;
935
936                                                                         for (j = 0; j < 3; ++j, ++i) {
937                                                                                 copy_v3_v3(fp, curofp);
938
939                                                                                 if (apply_offset) {
940                                                                                         add_v3_v3(fp, ofs[i]);
941                                                                                 }
942
943                                                                                 fp += 3; curofp += 3;
944                                                                         }
945                                                                         fp[0] = curofp[0];
946
947                                                                         if (apply_offset) {
948                                                                                 /* apply alfa offsets */
949                                                                                 add_v3_v3(fp, ofs[i]);
950                                                                                 i++;
951                                                                         }
952
953                                                                         fp += 3; /* alphas */
954                                                                 }
955                                                                 else {
956                                                                         int j;
957                                                                         for (j = 0; j < 3; ++j, ++i) {
958                                                                                 copy_v3_v3(fp, bezt->vec[j]);
959                                                                                 fp += 3;
960                                                                         }
961                                                                         fp[0] = bezt->alfa;
962
963                                                                         fp += 3; /* alphas */
964                                                                 }
965                                                                 bezt++;
966                                                         }
967                                                 }
968                                                 else {
969                                                         bp = nu->bp;
970                                                         a = nu->pntsu * nu->pntsv;
971                                                         while (a--) {
972                                                                 index = getKeyIndexOrig_keyIndex(editnurb, bp);
973
974                                                                 if (index >= 0) {
975                                                                         curofp = ofp + index;
976                                                                         copy_v3_v3(fp, curofp);
977                                                                         fp[3] = curofp[3];
978
979                                                                         if (apply_offset) {
980                                                                                 add_v3_v3(fp, ofs[i]);
981                                                                                 fp[3] += ofs[i + 1][0];
982                                                                         }
983                                                                 }
984                                                                 else {
985                                                                         copy_v3_v3(fp, bp->vec);
986                                                                         fp[3] = bp->alfa;
987                                                                 }
988
989                                                                 fp += 4;
990                                                                 bp++;
991                                                                 i += 2;
992                                                         }
993                                                 }
994                                         }
995                                 }
996
997                                 nu = nu->next;
998                         }
999
1000                         if (apply_offset) {
1001                                 /* handles could become malicious after offsets applying */
1002                                 calc_keyHandles(&editnurb->nurbs, newkey);
1003                         }
1004
1005                         currkey->totelem = totvert;
1006                         if (currkey->data) MEM_freeN(currkey->data);
1007                         currkey->data = newkey;
1008
1009                         currkey = currkey->next;
1010                 }
1011
1012                 if (ofs) MEM_freeN(ofs);
1013         }
1014 }
1015
1016 /* ********************* Amimation data *************** */
1017
1018 static int curve_is_animated(Curve *cu)
1019 {
1020         AnimData *ad = BKE_animdata_from_id(&cu->id);
1021
1022         return ad && (ad->action || ad->drivers.first);
1023 }
1024
1025 static void fcurve_path_rename(AnimData *ad, char *orig_rna_path, char *rna_path, ListBase *orig_curves, ListBase *curves)
1026 {
1027         FCurve *fcu, *nfcu, *nextfcu;
1028         int len = strlen(orig_rna_path);
1029
1030         fcu = orig_curves->first;
1031         while (fcu) {
1032                 nextfcu = fcu->next;
1033                 if (!strncmp(fcu->rna_path, orig_rna_path, len)) {
1034                         char *spath, *suffix = fcu->rna_path + len;
1035                         nfcu = copy_fcurve(fcu);
1036                         spath = nfcu->rna_path;
1037                         nfcu->rna_path = BLI_sprintfN("%s%s", rna_path, suffix);
1038                         BLI_addtail(curves, nfcu);
1039
1040                         if (fcu->grp) {
1041                                 action_groups_remove_channel(ad->action, fcu);
1042                                 action_groups_add_channel(ad->action, fcu->grp, nfcu);
1043                         }
1044                         else if (ad->action && &ad->action->curves == orig_curves)
1045                                 BLI_remlink(&ad->action->curves, fcu);
1046                         else
1047                                 BLI_remlink(&ad->drivers, fcu);
1048
1049                         free_fcurve(fcu);
1050
1051                         MEM_freeN(spath);
1052                 }
1053                 fcu = nextfcu;
1054         }
1055 }
1056
1057 static void fcurve_remove(AnimData *ad, ListBase *orig_curves, FCurve *fcu)
1058 {
1059         if (orig_curves == &ad->drivers) BLI_remlink(&ad->drivers, fcu);
1060         else action_groups_remove_channel(ad->action, fcu);
1061
1062         free_fcurve(fcu);
1063 }
1064
1065 static void curve_rename_fcurves(Curve *cu, ListBase *orig_curves)
1066 {
1067         int nu_index = 0, a, pt_index;
1068         EditNurb *editnurb = cu->editnurb;
1069         Nurb *nu = editnurb->nurbs.first;
1070         CVKeyIndex *keyIndex;
1071         char rna_path[64], orig_rna_path[64];
1072         AnimData *ad = BKE_animdata_from_id(&cu->id);
1073         ListBase curves = {NULL, NULL};
1074         FCurve *fcu, *next;
1075
1076         while (nu) {
1077                 if (nu->bezt) {
1078                         BezTriple *bezt = nu->bezt;
1079                         a = nu->pntsu;
1080                         pt_index = 0;
1081
1082                         while (a--) {
1083                                 keyIndex = getCVKeyIndex(editnurb, bezt);
1084                                 if (keyIndex) {
1085                                         BLI_snprintf(rna_path, sizeof(rna_path), "splines[%d].bezier_points[%d]", nu_index, pt_index);
1086                                         BLI_snprintf(orig_rna_path, sizeof(orig_rna_path), "splines[%d].bezier_points[%d]", keyIndex->nu_index, keyIndex->pt_index);
1087
1088                                         if (keyIndex->switched) {
1089                                                 char handle_path[64], orig_handle_path[64];
1090                                                 BLI_snprintf(orig_handle_path, sizeof(orig_rna_path), "%s.handle_left", orig_rna_path);
1091                                                 BLI_snprintf(handle_path, sizeof(rna_path), "%s.handle_right", rna_path);
1092                                                 fcurve_path_rename(ad, orig_handle_path, handle_path, orig_curves, &curves);
1093
1094                                                 BLI_snprintf(orig_handle_path, sizeof(orig_rna_path), "%s.handle_right", orig_rna_path);
1095                                                 BLI_snprintf(handle_path, sizeof(rna_path), "%s.handle_left", rna_path);
1096                                                 fcurve_path_rename(ad, orig_handle_path, handle_path, orig_curves, &curves);
1097                                         }
1098
1099                                         fcurve_path_rename(ad, orig_rna_path, rna_path, orig_curves, &curves);
1100
1101                                         keyIndex->nu_index = nu_index;
1102                                         keyIndex->pt_index = pt_index;
1103                                 }
1104
1105                                 bezt++;
1106                                 pt_index++;
1107                         }
1108                 }
1109                 else {
1110                         BPoint *bp = nu->bp;
1111                         a = nu->pntsu * nu->pntsv;
1112                         pt_index = 0;
1113
1114                         while (a--) {
1115                                 keyIndex = getCVKeyIndex(editnurb, bp);
1116                                 if (keyIndex) {
1117                                         BLI_snprintf(rna_path, sizeof(rna_path), "splines[%d].points[%d]", nu_index, pt_index);
1118                                         BLI_snprintf(orig_rna_path, sizeof(orig_rna_path), "splines[%d].points[%d]", keyIndex->nu_index, keyIndex->pt_index);
1119                                         fcurve_path_rename(ad, orig_rna_path, rna_path, orig_curves, &curves);
1120
1121                                         keyIndex->nu_index = nu_index;
1122                                         keyIndex->pt_index = pt_index;
1123                                 }
1124
1125                                 bp++;
1126                                 pt_index++;
1127                         }
1128                 }
1129                 nu = nu->next;
1130                 nu_index++;
1131         }
1132
1133         /* remove paths for removed control points
1134          * need this to make further step with copying non-cv related curves copying
1135          * not touching cv's f-curves */
1136         for (fcu = orig_curves->first; fcu; fcu = next) {
1137                 next = fcu->next;
1138
1139                 if (!strncmp(fcu->rna_path, "splines", 7)) {
1140                         char *ch = strchr(fcu->rna_path, '.');
1141
1142                         if (ch && (!strncmp(ch, ".bezier_points", 14) || !strncmp(ch, ".points", 7)))
1143                                 fcurve_remove(ad, orig_curves, fcu);
1144                 }
1145         }
1146
1147         nu_index = 0;
1148         nu = editnurb->nurbs.first;
1149         while (nu) {
1150                 keyIndex = NULL;
1151                 if (nu->pntsu) {
1152                         if (nu->bezt) keyIndex = getCVKeyIndex(editnurb, &nu->bezt[0]);
1153                         else keyIndex = getCVKeyIndex(editnurb, &nu->bp[0]);
1154                 }
1155
1156                 if (keyIndex) {
1157                         BLI_snprintf(rna_path, sizeof(rna_path), "splines[%d]", nu_index);
1158                         BLI_snprintf(orig_rna_path, sizeof(orig_rna_path), "splines[%d]", keyIndex->nu_index);
1159                         fcurve_path_rename(ad, orig_rna_path, rna_path, orig_curves, &curves);
1160                 }
1161
1162                 nu_index++;
1163                 nu = nu->next;
1164         }
1165
1166         /* the remainders in orig_curves can be copied back (like follow path) */
1167         /* (if it's not path to spline) */
1168         for (fcu = orig_curves->first; fcu; fcu = next) {
1169                 next = fcu->next;
1170
1171                 if (!strncmp(fcu->rna_path, "splines", 7)) fcurve_remove(ad, orig_curves, fcu);
1172                 else BLI_addtail(&curves, fcu);
1173         }
1174
1175         *orig_curves = curves;
1176 }
1177
1178 /* return 0 if animation data wasn't changed, 1 otherwise */
1179 int ED_curve_updateAnimPaths(Curve *cu)
1180 {
1181         AnimData *ad = BKE_animdata_from_id(&cu->id);
1182
1183         if (!curve_is_animated(cu)) return 0;
1184
1185         if (ad->action)
1186                 curve_rename_fcurves(cu, &ad->action->curves);
1187
1188         curve_rename_fcurves(cu, &ad->drivers);
1189
1190         return 1;
1191 }
1192
1193 /* ********************* LOAD and MAKE *************** */
1194
1195 /* load editNurb in object */
1196 void load_editNurb(Object *obedit)
1197 {
1198         ListBase *editnurb = object_editcurve_get(obedit);
1199
1200         if (obedit == NULL) return;
1201
1202         if (ELEM(obedit->type, OB_CURVE, OB_SURF)) {
1203                 Curve *cu = obedit->data;
1204                 Nurb *nu, *newnu;
1205                 ListBase newnurb = {NULL, NULL}, oldnurb = cu->nurb;
1206
1207                 for (nu = editnurb->first; nu; nu = nu->next) {
1208                         newnu = BKE_nurb_duplicate(nu);
1209                         BLI_addtail(&newnurb, newnu);
1210
1211                         if (nu->type == CU_NURBS) {
1212                                 BKE_nurb_order_clamp_u(nu);
1213                         }
1214                 }
1215
1216                 cu->nurb = newnurb;
1217
1218                 calc_shapeKeys(obedit);
1219                 ED_curve_updateAnimPaths(obedit->data);
1220
1221                 BKE_nurbList_free(&oldnurb);
1222         }
1223 }
1224
1225 /* make copy in cu->editnurb */
1226 void make_editNurb(Object *obedit)
1227 {
1228         Curve *cu = (Curve *)obedit->data;
1229         EditNurb *editnurb = cu->editnurb;
1230         Nurb *nu, *newnu, *nu_act = NULL;
1231         KeyBlock *actkey;
1232
1233         set_actNurb(obedit, NULL);
1234
1235         if (ELEM(obedit->type, OB_CURVE, OB_SURF)) {
1236                 actkey = BKE_keyblock_from_object(obedit);
1237
1238                 if (actkey) {
1239                         // XXX strcpy(G.editModeTitleExtra, "(Key) ");
1240                         undo_editmode_clear();
1241                         BKE_key_convert_to_curve(actkey, cu, &cu->nurb);
1242                 }
1243
1244                 if (editnurb) {
1245                         BKE_nurbList_free(&editnurb->nurbs);
1246                         BKE_curve_editNurb_keyIndex_free(editnurb);
1247                         editnurb->keyindex = NULL;
1248                 }
1249                 else {
1250                         editnurb = MEM_callocN(sizeof(EditNurb), "editnurb");
1251                         cu->editnurb = editnurb;
1252                 }
1253
1254                 nu = cu->nurb.first;
1255                 cu->lastsel = NULL;   /* for select row */
1256
1257                 while (nu) {
1258                         newnu = BKE_nurb_duplicate(nu);
1259                         BKE_nurb_test2D(newnu); // after join, or any other creation of curve
1260                         BLI_addtail(&editnurb->nurbs, newnu);
1261
1262                         if (nu_act == NULL && isNurbsel(nu)) {
1263                                 nu_act = newnu;
1264                                 set_actNurb(obedit, newnu);
1265                         }
1266
1267                         nu = nu->next;
1268                 }
1269
1270                 if (actkey)
1271                         editnurb->shapenr = obedit->shapenr;
1272
1273                 /* animation could be added in editmode even if there was no animdata i
1274                  * object mode hence we always need CVs index be created */
1275                 init_editNurb_keyIndex(editnurb, &cu->nurb);
1276         }
1277 }
1278
1279 void free_editNurb(Object *obedit)
1280 {
1281         Curve *cu = obedit->data;
1282
1283         BKE_curve_editNurb_free(cu);
1284 }
1285
1286 void CU_deselect_all(Object *obedit)
1287 {
1288         ListBase *editnurb = object_editcurve_get(obedit);
1289
1290         if (editnurb) {
1291                 selectend_nurb(obedit, FIRST, 0, DESELECT); /* set first control points as unselected */
1292                 select_adjacent_cp(editnurb, 1, 1, DESELECT); /* cascade selection */
1293         }
1294 }
1295
1296 void CU_select_all(Object *obedit)
1297 {
1298         ListBase *editnurb = object_editcurve_get(obedit);
1299
1300         if (editnurb) {
1301                 selectend_nurb(obedit, FIRST, 0, SELECT); /* set first control points as unselected */
1302                 select_adjacent_cp(editnurb, 1, 1, SELECT); /* cascade selection */
1303         }
1304 }
1305
1306 void CU_select_swap(Object *obedit)
1307 {
1308         ListBase *editnurb = object_editcurve_get(obedit);
1309
1310         if (editnurb) {
1311                 Curve *cu = obedit->data;
1312                 Nurb *nu;
1313                 BPoint *bp;
1314                 BezTriple *bezt;
1315                 int a;
1316
1317                 cu->lastsel = NULL;
1318
1319                 for (nu = editnurb->first; nu; nu = nu->next) {
1320                         if (nu->type == CU_BEZIER) {
1321                                 bezt = nu->bezt;
1322                                 a = nu->pntsu;
1323                                 while (a--) {
1324                                         if (bezt->hide == 0) {
1325                                                 bezt->f2 ^= SELECT; /* always do the center point */
1326                                                 if ((cu->drawflag & CU_HIDE_HANDLES) == 0) {
1327                                                         bezt->f1 ^= SELECT;
1328                                                         bezt->f3 ^= SELECT;
1329                                                 }
1330                                         }
1331                                         bezt++;
1332                                 }
1333                         }
1334                         else {
1335                                 bp = nu->bp;
1336                                 a = nu->pntsu * nu->pntsv;
1337                                 while (a--) {
1338                                         swap_selection_bpoint(bp);
1339                                         bp++;
1340                                 }
1341                         }
1342                 }
1343         }
1344 }
1345
1346 /******************** separate operator ***********************/
1347
1348 static int separate_exec(bContext *C, wmOperator *op)
1349 {
1350         Main *bmain = CTX_data_main(C);
1351         Scene *scene = CTX_data_scene(C);
1352         Nurb *nu, *nu1;
1353         Object *oldob, *newob;
1354         Base *oldbase, *newbase;
1355         Curve *oldcu, *newcu;
1356         EditNurb *oldedit, *newedit;
1357
1358         oldbase = CTX_data_active_base(C);
1359         oldob = oldbase->object;
1360         oldcu = oldob->data;
1361         oldedit = oldcu->editnurb;
1362
1363         if (oldcu->key) {
1364                 BKE_report(op->reports, RPT_ERROR, "Can't separate a curve with vertex keys");
1365                 return OPERATOR_CANCELLED;
1366         }
1367
1368         WM_cursor_wait(1);
1369         
1370         /* 1. duplicate the object and data */
1371         newbase = ED_object_add_duplicate(bmain, scene, oldbase, 0); /* 0 = fully linked */
1372         DAG_scene_sort(bmain, scene);
1373
1374         ED_base_object_select(newbase, BA_DESELECT);
1375         newob = newbase->object;
1376
1377         newcu = newob->data = BKE_curve_copy(oldcu);
1378         newcu->editnurb = NULL;
1379         oldcu->id.us--; /* because new curve is a copy: reduce user count */
1380
1381         /* 2. put new object in editmode and clear it */
1382         make_editNurb(newob);
1383         newedit = newcu->editnurb;
1384         BKE_nurbList_free(&newedit->nurbs);
1385         BKE_curve_editNurb_keyIndex_free(newedit);
1386
1387         /* 3. move over parts from old object */
1388         for (nu = oldedit->nurbs.first; nu; nu = nu1) {
1389                 nu1 = nu->next;
1390
1391                 if (isNurbsel(nu)) {
1392                         BLI_remlink(&oldedit->nurbs, nu);
1393                         BLI_addtail(&newedit->nurbs, nu);
1394                 }
1395         }
1396
1397         /* 4. put old object out of editmode */
1398         load_editNurb(newob);
1399         free_editNurb(newob);
1400
1401         DAG_id_tag_update(&oldob->id, OB_RECALC_DATA);  /* this is the original one */
1402         DAG_id_tag_update(&newob->id, OB_RECALC_DATA);  /* this is the separated one */
1403
1404         WM_event_add_notifier(C, NC_GEOM | ND_DATA, oldob->data);
1405
1406         WM_cursor_wait(0);
1407
1408         return OPERATOR_FINISHED;
1409 }
1410
1411 void CURVE_OT_separate(wmOperatorType *ot)
1412 {
1413         /* identifiers */
1414         ot->name = "Separate";
1415         ot->idname = "CURVE_OT_separate";
1416         ot->description = "Separate (partly) selected curves or surfaces into a new object";
1417         
1418         /* api callbacks */
1419         ot->exec = separate_exec;
1420         ot->poll = ED_operator_editsurfcurve;
1421         
1422         /* flags */
1423         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
1424 }
1425
1426 /* ******************* FLAGS ********************* */
1427
1428 static short isNurbselUV(Nurb *nu, int *u, int *v, int flag)
1429 {
1430         /* return (u != -1): 1 row in u-direction selected. U has value between 0-pntsv
1431          * return (v != -1): 1 column in v-direction selected. V has value between 0-pntsu
1432          */
1433         BPoint *bp;
1434         int a, b, sel;
1435
1436         *u = *v = -1;
1437
1438         bp = nu->bp;
1439         for (b = 0; b < nu->pntsv; b++) {
1440                 sel = 0;
1441                 for (a = 0; a < nu->pntsu; a++, bp++) {
1442                         if (bp->f1 & flag) sel++;
1443                 }
1444                 if (sel == nu->pntsu) {
1445                         if (*u == -1) *u = b;
1446                         else return 0;
1447                 }
1448                 else if (sel > 1) return 0;  /* because sel == 1 is still ok */
1449         }
1450
1451         for (a = 0; a < nu->pntsu; a++) {
1452                 sel = 0;
1453                 bp = &nu->bp[a];
1454                 for (b = 0; b < nu->pntsv; b++, bp += nu->pntsu) {
1455                         if (bp->f1 & flag) sel++;
1456                 }
1457                 if (sel == nu->pntsv) {
1458                         if (*v == -1) *v = a;
1459                         else return 0;
1460                 }
1461                 else if (sel > 1) return 0;
1462         }
1463
1464         if (*u == -1 && *v > -1) return 1;
1465         if (*v == -1 && *u > -1) return 1;
1466         return 0;
1467 }
1468
1469 static void setflagsNurb(ListBase *editnurb, short flag)
1470 {
1471         Nurb *nu;
1472         BezTriple *bezt;
1473         BPoint *bp;
1474         int a;
1475
1476         for (nu = editnurb->first; nu; nu = nu->next) {
1477                 if (nu->type == CU_BEZIER) {
1478                         a = nu->pntsu;
1479                         bezt = nu->bezt;
1480                         while (a--) {
1481                                 bezt->f1 = bezt->f2 = bezt->f3 = flag;
1482                                 bezt++;
1483                         }
1484                 }
1485                 else {
1486                         a = nu->pntsu * nu->pntsv;
1487                         bp = nu->bp;
1488                         while (a--) {
1489                                 bp->f1 = flag;
1490                                 bp++;
1491                         }
1492                 }
1493         }
1494 }
1495
1496 static void rotateflagNurb(ListBase *editnurb, short flag, float *cent, float rotmat[][3])
1497 {
1498         /* all verts with (flag & 'flag') rotate */
1499         Nurb *nu;
1500         BPoint *bp;
1501         int a;
1502
1503         for (nu = editnurb->first; nu; nu = nu->next) {
1504                 if (nu->type == CU_NURBS) {
1505                         bp = nu->bp;
1506                         a = nu->pntsu * nu->pntsv;
1507
1508                         while (a--) {
1509                                 if (bp->f1 & flag) {
1510                                         sub_v3_v3(bp->vec, cent);
1511                                         mul_m3_v3(rotmat, bp->vec);
1512                                         add_v3_v3(bp->vec, cent);
1513                                 }
1514                                 bp++;
1515                         }
1516                 }
1517         }
1518 }
1519
1520 static void translateflagNurb(ListBase *editnurb, short flag, const float vec[3])
1521 {
1522         /* all verts with ('flag' & flag) translate */
1523         Nurb *nu;
1524         BezTriple *bezt;
1525         BPoint *bp;
1526         int a;
1527
1528         for (nu = editnurb->first; nu; nu = nu->next) {
1529                 if (nu->type == CU_BEZIER) {
1530                         a = nu->pntsu;
1531                         bezt = nu->bezt;
1532                         while (a--) {
1533                                 if (bezt->f1 & flag) add_v3_v3(bezt->vec[0], vec);
1534                                 if (bezt->f2 & flag) add_v3_v3(bezt->vec[1], vec);
1535                                 if (bezt->f3 & flag) add_v3_v3(bezt->vec[2], vec);
1536                                 bezt++;
1537                         }
1538                 }
1539                 else {
1540                         a = nu->pntsu * nu->pntsv;
1541                         bp = nu->bp;
1542                         while (a--) {
1543                                 if (bp->f1 & flag) add_v3_v3(bp->vec, vec);
1544                                 bp++;
1545                         }
1546                 }
1547
1548                 BKE_nurb_test2D(nu);
1549         }
1550 }
1551
1552 static void weightflagNurb(ListBase *editnurb, short flag, float w)
1553 {
1554         Nurb *nu;
1555         BPoint *bp;
1556         int a;
1557
1558         for (nu = editnurb->first; nu; nu = nu->next) {
1559                 if (nu->type == CU_NURBS) {
1560                         a = nu->pntsu * nu->pntsv;
1561                         bp = nu->bp;
1562                         while (a--) {
1563                                 if (bp->f1 & flag) {
1564                                         /* a mode used to exist for replace/multiple but is was unused */
1565                                         bp->vec[3] *= w;
1566                                 }
1567                                 bp++;
1568                         }
1569                 }
1570         }
1571 }
1572
1573 static int deleteflagNurb(bContext *C, wmOperator *UNUSED(op), int flag)
1574 {
1575         Object *obedit = CTX_data_edit_object(C);
1576         Curve *cu = obedit->data;
1577         ListBase *editnurb = object_editcurve_get(obedit);
1578         Nurb *nu, *next;
1579         BPoint *bp, *bpn, *newbp;
1580         int a, b, newu, newv, sel;
1581
1582         if (obedit->type != OB_SURF) {
1583                 return OPERATOR_CANCELLED;
1584         }
1585
1586         cu->lastsel = NULL;
1587
1588         nu = editnurb->first;
1589         while (nu) {
1590                 next = nu->next;
1591
1592                 /* is entire nurb selected */
1593                 bp = nu->bp;
1594                 a = nu->pntsu * nu->pntsv;
1595                 while (a) {
1596                         a--;
1597                         if (bp->f1 & flag) {
1598                                 /* pass */
1599                         }
1600                         else {
1601                                 break;
1602                         }
1603                         bp++;
1604                 }
1605                 if (a == 0) {
1606                         BLI_remlink(editnurb, nu);
1607                         keyIndex_delNurb(cu->editnurb, nu);
1608                         BKE_nurb_free(nu); nu = NULL;
1609                 }
1610                 else {
1611                         /* is nurb in U direction selected */
1612                         newv = nu->pntsv;
1613                         bp = nu->bp;
1614                         for (b = 0; b < nu->pntsv; b++) {
1615                                 sel = 0;
1616                                 for (a = 0; a < nu->pntsu; a++, bp++) {
1617                                         if (bp->f1 & flag) sel++;
1618                                 }
1619                                 if (sel == nu->pntsu) {
1620                                         newv--;
1621                                 }
1622                                 else if (sel >= 1) {
1623                                         /* don't delete */
1624                                         break;
1625                                 }
1626                         }
1627                         if (newv != nu->pntsv && b == nu->pntsv) {
1628                                 /* delete */
1629                                 bp = nu->bp;
1630                                 bpn = newbp = (BPoint *)MEM_mallocN(newv * nu->pntsu * sizeof(BPoint), "deleteNurb");
1631                                 for (b = 0; b < nu->pntsv; b++) {
1632                                         if ((bp->f1 & flag) == 0) {
1633                                                 memcpy(bpn, bp, nu->pntsu * sizeof(BPoint));
1634                                                 keyIndex_updateBP(cu->editnurb, bp, bpn, nu->pntsu);
1635                                                 bpn += nu->pntsu;
1636                                         }
1637                                         else {
1638                                                 keyIndex_delBP(cu->editnurb, bp);
1639                                         }
1640                                         bp += nu->pntsu;
1641                                 }
1642                                 nu->pntsv = newv;
1643                                 MEM_freeN(nu->bp);
1644                                 nu->bp = newbp;
1645                                 BKE_nurb_order_clamp_v(nu);
1646
1647                                 BKE_nurb_knot_calc_v(nu);
1648                         }
1649                         else {
1650                                 /* is the nurb in V direction selected */
1651                                 newu = nu->pntsu;
1652                                 for (a = 0; a < nu->pntsu; a++) {
1653                                         bp = &nu->bp[a];
1654                                         sel = 0;
1655                                         for (b = 0; b < nu->pntsv; b++, bp += nu->pntsu) {
1656                                                 if (bp->f1 & flag) sel++;
1657                                         }
1658                                         if (sel == nu->pntsv) {
1659                                                 newu--;
1660                                         }
1661                                         else if (sel >= 1) {
1662                                                 /* don't delete */
1663                                                 break;
1664                                         }
1665                                 }
1666                                 if (newu != nu->pntsu && a == nu->pntsu) {
1667                                         /* delete */
1668                                         bp = nu->bp;
1669                                         bpn = newbp = (BPoint *)MEM_mallocN(newu * nu->pntsv * sizeof(BPoint), "deleteNurb");
1670                                         for (b = 0; b < nu->pntsv; b++) {
1671                                                 for (a = 0; a < nu->pntsu; a++, bp++) {
1672                                                         if ((bp->f1 & flag) == 0) {
1673                                                                 *bpn = *bp;
1674                                                                 keyIndex_updateBP(cu->editnurb, bp, bpn, 1);
1675                                                                 bpn++;
1676                                                         }
1677                                                         else {
1678                                                                 keyIndex_delBP(cu->editnurb, bp);
1679                                                         }
1680                                                 }
1681                                         }
1682                                         MEM_freeN(nu->bp);
1683                                         nu->bp = newbp;
1684                                         if (newu == 1 && nu->pntsv > 1) {    /* make a U spline */
1685                                                 nu->pntsu = nu->pntsv;
1686                                                 nu->pntsv = 1;
1687                                                 SWAP(short, nu->orderu, nu->orderv);
1688                                                 BKE_nurb_order_clamp_u(nu);
1689                                                 if (nu->knotsv) MEM_freeN(nu->knotsv);
1690                                                 nu->knotsv = NULL;
1691                                         }
1692                                         else {
1693                                                 nu->pntsu = newu;
1694                                                 BKE_nurb_order_clamp_u(nu);
1695                                         }
1696                                         BKE_nurb_knot_calc_u(nu);
1697                                 }
1698                         }
1699                 }
1700                 nu = next;
1701         }
1702
1703         if (ED_curve_updateAnimPaths(obedit->data))
1704                 WM_event_add_notifier(C, NC_OBJECT | ND_KEYS, obedit);
1705
1706         return OPERATOR_FINISHED;
1707 }
1708
1709 /* only for OB_SURF */
1710 static short extrudeflagNurb(EditNurb *editnurb, int flag)
1711 {
1712         Nurb *nu;
1713         BPoint *bp, *bpn, *newbp;
1714         int ok = 0, a, u, v, len;
1715
1716         nu = editnurb->nurbs.first;
1717         while (nu) {
1718
1719                 if (nu->pntsv == 1) {
1720                         bp = nu->bp;
1721                         a = nu->pntsu;
1722                         while (a) {
1723                                 if (bp->f1 & flag) {
1724                                         /* pass */
1725                                 }
1726                                 else {
1727                                         break;
1728                                 }
1729                                 bp++;
1730                                 a--;
1731                         }
1732                         if (a == 0) {
1733                                 ok = 1;
1734                                 newbp = (BPoint *)MEM_mallocN(2 * nu->pntsu * sizeof(BPoint), "extrudeNurb1");
1735                                 ED_curve_bpcpy(editnurb, newbp, nu->bp, nu->pntsu);
1736                                 bp = newbp + nu->pntsu;
1737                                 ED_curve_bpcpy(editnurb, bp, nu->bp, nu->pntsu);
1738                                 MEM_freeN(nu->bp);
1739                                 nu->bp = newbp;
1740                                 a = nu->pntsu;
1741                                 while (a--) {
1742                                         select_bpoint(bp, SELECT, flag, HIDDEN);
1743                                         select_bpoint(newbp, DESELECT, flag, HIDDEN);
1744                                         bp++; 
1745                                         newbp++;
1746                                 }
1747
1748                                 nu->pntsv = 2;
1749                                 nu->orderv = 2;
1750                                 BKE_nurb_knot_calc_v(nu);
1751                         }
1752                 }
1753                 else {
1754                         /* which row or column is selected */
1755
1756                         if (isNurbselUV(nu, &u, &v, flag) ) {
1757
1758                                 /* deselect all */
1759                                 bp = nu->bp;
1760                                 a = nu->pntsu * nu->pntsv;
1761                                 while (a--) {
1762                                         select_bpoint(bp, DESELECT, flag, HIDDEN);
1763                                         bp++;
1764                                 }
1765
1766                                 if (u == 0 || u == nu->pntsv - 1) {      /* row in u-direction selected */
1767                                         ok = 1;
1768                                         newbp = (BPoint *)MEM_mallocN(nu->pntsu * (nu->pntsv + 1) *
1769                                                                       sizeof(BPoint), "extrudeNurb1");
1770                                         if (u == 0) {
1771                                                 len = nu->pntsv * nu->pntsu;
1772                                                 ED_curve_bpcpy(editnurb, newbp + nu->pntsu, nu->bp, len);
1773                                                 ED_curve_bpcpy(editnurb, newbp, nu->bp, nu->pntsu);
1774                                                 bp = newbp;
1775                                         }
1776                                         else {
1777                                                 len = nu->pntsv * nu->pntsu;
1778                                                 ED_curve_bpcpy(editnurb, newbp, nu->bp, len);
1779                                                 ED_curve_bpcpy(editnurb, newbp + len, &nu->bp[len - nu->pntsu], nu->pntsu);
1780                                                 bp = newbp + len;
1781                                         }
1782
1783                                         a = nu->pntsu;
1784                                         while (a--) {
1785                                                 select_bpoint(bp, SELECT, flag, HIDDEN);
1786                                                 bp++;
1787                                         }
1788
1789                                         MEM_freeN(nu->bp);
1790                                         nu->bp = newbp;
1791                                         nu->pntsv++;
1792                                         BKE_nurb_knot_calc_v(nu);
1793                                 }
1794                                 else if (v == 0 || v == nu->pntsu - 1) {     /* column in v-direction selected */
1795                                         ok = 1;
1796                                         bpn = newbp = (BPoint *)MEM_mallocN((nu->pntsu + 1) * nu->pntsv * sizeof(BPoint), "extrudeNurb1");
1797                                         bp = nu->bp;
1798
1799                                         for (a = 0; a < nu->pntsv; a++) {
1800                                                 if (v == 0) {
1801                                                         *bpn = *bp;
1802                                                         bpn->f1 |= flag;
1803                                                         bpn++;
1804                                                 }
1805                                                 ED_curve_bpcpy(editnurb, bpn, bp, nu->pntsu);
1806                                                 bp += nu->pntsu;
1807                                                 bpn += nu->pntsu;
1808                                                 if (v == nu->pntsu - 1) {
1809                                                         *bpn = *(bp - 1);
1810                                                         bpn->f1 |= flag;
1811                                                         bpn++;
1812                                                 }
1813                                         }
1814
1815                                         MEM_freeN(nu->bp);
1816                                         nu->bp = newbp;
1817                                         nu->pntsu++;
1818                                         BKE_nurb_knot_calc_u(nu);
1819                                 }
1820                         }
1821                 }
1822                 nu = nu->next;
1823         }
1824
1825         return ok;
1826 }
1827
1828 static void adduplicateflagNurb(Object *obedit, short flag)
1829 {
1830         ListBase *editnurb = object_editcurve_get(obedit);
1831         Nurb *nu, *newnu;
1832         BezTriple *bezt, *bezt1;
1833         BPoint *bp, *bp1;
1834         Curve *cu = (Curve *)obedit->data;
1835         int a, b, starta, enda, newu, newv;
1836         char *usel;
1837
1838         cu->lastsel = NULL;
1839
1840         nu = editnurb->last;
1841         while (nu) {
1842                 if (nu->type == CU_BEZIER) {
1843                         bezt = nu->bezt;
1844                         for (a = 0; a < nu->pntsu; a++) {
1845                                 enda = -1;
1846                                 starta = a;
1847                                 while ( (bezt->f1 & flag) || (bezt->f2 & flag) || (bezt->f3 & flag) ) {
1848                                         select_beztriple(bezt, DESELECT, flag, HIDDEN);
1849                                         enda = a;
1850                                         if (a >= nu->pntsu - 1) break;
1851                                         a++;
1852                                         bezt++;
1853                                 }
1854                                 if (enda >= starta) {
1855                                         newnu = (Nurb *)MEM_mallocN(sizeof(Nurb), "adduplicateN");
1856                                         memcpy(newnu, nu, sizeof(Nurb));
1857                                         BLI_addtail(editnurb, newnu);
1858                                         set_actNurb(obedit, newnu);
1859                                         newnu->pntsu = enda - starta + 1;
1860                                         newnu->bezt = (BezTriple *)MEM_mallocN((enda - starta + 1) * sizeof(BezTriple), "adduplicateN");
1861                                         memcpy(newnu->bezt, &nu->bezt[starta], newnu->pntsu * sizeof(BezTriple));
1862
1863                                         b = newnu->pntsu;
1864                                         bezt1 = newnu->bezt;
1865                                         while (b--) {
1866                                                 select_beztriple(bezt1, SELECT, flag, HIDDEN);
1867                                                 bezt1++;
1868                                         }
1869
1870                                         if (nu->flagu & CU_NURB_CYCLIC) {
1871                                                 if (starta != 0 || enda != nu->pntsu - 1) {
1872                                                         newnu->flagu &= ~CU_NURB_CYCLIC;
1873                                                 }
1874                                         }
1875                                 }
1876                                 bezt++;
1877                         }
1878                 }
1879                 else if (nu->pntsv == 1) {    /* because UV Nurb has a different method for dupli */
1880                         bp = nu->bp;
1881                         for (a = 0; a < nu->pntsu; a++) {
1882                                 enda = -1;
1883                                 starta = a;
1884                                 while (bp->f1 & flag) {
1885                                         select_bpoint(bp, DESELECT, flag, HIDDEN);
1886                                         enda = a;
1887                                         if (a >= nu->pntsu - 1) break;
1888                                         a++;
1889                                         bp++;
1890                                 }
1891                                 if (enda >= starta) {
1892                                         newnu = (Nurb *)MEM_mallocN(sizeof(Nurb), "adduplicateN3");
1893                                         memcpy(newnu, nu, sizeof(Nurb));
1894                                         set_actNurb(obedit, newnu);
1895                                         BLI_addtail(editnurb, newnu);
1896                                         newnu->pntsu = enda - starta + 1;
1897                                         newnu->bp = (BPoint *)MEM_mallocN((enda - starta + 1) * sizeof(BPoint), "adduplicateN4");
1898                                         memcpy(newnu->bp, &nu->bp[starta], newnu->pntsu * sizeof(BPoint));
1899
1900                                         b = newnu->pntsu;
1901                                         bp1 = newnu->bp;
1902                                         while (b--) {
1903                                                 select_bpoint(bp1, SELECT, flag, HIDDEN);
1904                                                 bp1++;
1905                                         }
1906
1907                                         if (nu->flagu & CU_NURB_CYCLIC) {
1908                                                 if (starta != 0 || enda != nu->pntsu - 1) {
1909                                                         newnu->flagu &= ~CU_NURB_CYCLIC;
1910                                                 }
1911                                         }
1912
1913                                         /* knots */
1914                                         newnu->knotsu = NULL;
1915                                         BKE_nurb_knot_calc_u(newnu);
1916                                 }
1917                                 bp++;
1918                         }
1919                 }
1920                 else {
1921                         /* a rectangular area in nurb has to be selected */
1922                         if (isNurbsel(nu)) {
1923                                 usel = MEM_callocN(nu->pntsu, "adduplicateN4");
1924                                 bp = nu->bp;
1925                                 for (a = 0; a < nu->pntsv; a++) {
1926                                         for (b = 0; b < nu->pntsu; b++, bp++) {
1927                                                 if (bp->f1 & flag) usel[b]++;
1928                                         }
1929                                 }
1930                                 newu = 0;
1931                                 newv = 0;
1932                                 for (a = 0; a < nu->pntsu; a++) {
1933                                         if (usel[a]) {
1934                                                 if (newv == 0 || usel[a] == newv) {
1935                                                         newv = usel[a];
1936                                                         newu++;
1937                                                 }
1938                                                 else {
1939                                                         newv = 0;
1940                                                         break;
1941                                                 }
1942                                         }
1943                                 }
1944                                 if (newu == 0 || newv == 0) {
1945                                         if (G.debug & G_DEBUG)
1946                                                 printf("Can't duplicate Nurb\n");
1947                                 }
1948                                 else {
1949
1950                                         if (newu == 1) SWAP(int, newu, newv);
1951
1952                                         newnu = (Nurb *)MEM_mallocN(sizeof(Nurb), "adduplicateN5");
1953                                         memcpy(newnu, nu, sizeof(Nurb));
1954                                         BLI_addtail(editnurb, newnu);
1955                                         set_actNurb(obedit, newnu);
1956                                         newnu->pntsu = newu;
1957                                         newnu->pntsv = newv;
1958                                         newnu->bp = (BPoint *)MEM_mallocN(newu * newv * sizeof(BPoint), "adduplicateN6");
1959                                         BKE_nurb_order_clamp_u(newnu);
1960                                         BKE_nurb_order_clamp_v(newnu);
1961                                         
1962                                         newnu->knotsu = newnu->knotsv = NULL;
1963                                         
1964                                         bp = newnu->bp;
1965                                         bp1 = nu->bp;
1966                                         for (a = 0; a < nu->pntsv; a++) {
1967                                                 for (b = 0; b < nu->pntsu; b++, bp1++) {
1968                                                         if (bp1->f1 & flag) {
1969                                                                 memcpy(bp, bp1, sizeof(BPoint));
1970                                                                 select_bpoint(bp1, DESELECT, flag, HIDDEN);
1971                                                                 bp++;
1972                                                         }
1973                                                 }
1974                                         }
1975                                         if (BKE_nurb_check_valid_u(newnu)) {
1976                                                 if (nu->pntsu == newnu->pntsu && nu->knotsu) {
1977                                                         newnu->knotsu = MEM_dupallocN(nu->knotsu);
1978                                                 }
1979                                                 else {
1980                                                         BKE_nurb_knot_calc_u(newnu);
1981                                                 }
1982                                         }
1983                                         if (BKE_nurb_check_valid_v(newnu)) {
1984                                                 if (nu->pntsv == newnu->pntsv && nu->knotsv) {
1985                                                         newnu->knotsv = MEM_dupallocN(nu->knotsv);
1986                                                 }
1987                                                 else {
1988                                                         BKE_nurb_knot_calc_v(newnu);
1989                                                 }
1990                                         }
1991                                 }
1992                                 MEM_freeN(usel);
1993                         }
1994                 }
1995
1996                 nu = nu->prev;
1997         }
1998         
1999         /* actnu changed */
2000 }
2001
2002 /**************** switch direction operator ***************/
2003
2004 static int switch_direction_exec(bContext *C, wmOperator *UNUSED(op))
2005 {
2006         Object *obedit = CTX_data_edit_object(C);
2007         Curve *cu = (Curve *)obedit->data;
2008         EditNurb *editnurb = cu->editnurb;
2009         Nurb *nu;
2010
2011         for (nu = editnurb->nurbs.first; nu; nu = nu->next)
2012                 if (isNurbsel(nu)) {
2013                         BKE_nurb_direction_switch(nu);
2014                         keyData_switchDirectionNurb(cu, nu);
2015                 }
2016
2017         if (ED_curve_updateAnimPaths(obedit->data))
2018                 WM_event_add_notifier(C, NC_OBJECT | ND_KEYS, obedit);
2019
2020         DAG_id_tag_update(obedit->data, 0);
2021         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2022
2023         return OPERATOR_FINISHED;
2024 }
2025
2026 void CURVE_OT_switch_direction(wmOperatorType *ot)
2027 {
2028         /* identifiers */
2029         ot->name = "Switch Direction";
2030         ot->description = "Switch direction of selected splines";
2031         ot->idname = "CURVE_OT_switch_direction";
2032         
2033         /* api callbacks */
2034         ot->exec = switch_direction_exec;
2035         ot->poll = ED_operator_editsurfcurve;
2036
2037         /* flags */
2038         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2039 }
2040
2041 /****************** set weight operator *******************/
2042
2043 static int set_goal_weight_exec(bContext *C, wmOperator *op)
2044 {
2045         Object *obedit = CTX_data_edit_object(C);
2046         ListBase *editnurb = object_editcurve_get(obedit);
2047         Nurb *nu;
2048         BezTriple *bezt;
2049         BPoint *bp;
2050         float weight = RNA_float_get(op->ptr, "weight");
2051         int a;
2052                                 
2053         for (nu = editnurb->first; nu; nu = nu->next) {
2054                 if (nu->bezt) {
2055                         for (bezt = nu->bezt, a = 0; a < nu->pntsu; a++, bezt++) {
2056                                 if (bezt->f2 & SELECT)
2057                                         bezt->weight = weight;
2058                         }
2059                 }
2060                 else if (nu->bp) {
2061                         for (bp = nu->bp, a = 0; a < nu->pntsu * nu->pntsv; a++, bp++) {
2062                                 if (bp->f1 & SELECT)
2063                                         bp->weight = weight;
2064                         }
2065                 }
2066         }
2067
2068         DAG_id_tag_update(obedit->data, 0);
2069         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2070
2071         return OPERATOR_FINISHED;
2072 }
2073
2074 void CURVE_OT_spline_weight_set(wmOperatorType *ot)
2075 {
2076         /* identifiers */
2077         ot->name = "Set Goal Weight";
2078         ot->description = "Set softbody goal weight for selected points";
2079         ot->idname = "CURVE_OT_spline_weight_set";
2080         
2081         /* api callbacks */
2082         ot->exec = set_goal_weight_exec;
2083         ot->invoke = WM_operator_props_popup;
2084         ot->poll = ED_operator_editsurfcurve;
2085
2086         /* flags */
2087         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2088
2089         /* properties */
2090         RNA_def_float_factor(ot->srna, "weight", 1.0f, 0.0f, 1.0f, "Weight", "", 0.0f, 1.0f);
2091 }
2092
2093 /******************* set radius operator ******************/
2094
2095 static int set_radius_exec(bContext *C, wmOperator *op)
2096 {
2097         Object *obedit = CTX_data_edit_object(C);
2098         ListBase *editnurb = object_editcurve_get(obedit);
2099         Nurb *nu;
2100         BezTriple *bezt;
2101         BPoint *bp;
2102         float radius = RNA_float_get(op->ptr, "radius");
2103         int a;
2104         
2105         for (nu = editnurb->first; nu; nu = nu->next) {
2106                 if (nu->bezt) {
2107                         for (bezt = nu->bezt, a = 0; a < nu->pntsu; a++, bezt++) {
2108                                 if (bezt->f2 & SELECT)
2109                                         bezt->radius = radius;
2110                         }
2111                 }
2112                 else if (nu->bp) {
2113                         for (bp = nu->bp, a = 0; a < nu->pntsu * nu->pntsv; a++, bp++) {
2114                                 if (bp->f1 & SELECT)
2115                                         bp->radius = radius;
2116                         }
2117                 }
2118         }
2119
2120         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2121         DAG_id_tag_update(obedit->data, 0);
2122
2123         return OPERATOR_FINISHED;
2124 }
2125
2126 void CURVE_OT_radius_set(wmOperatorType *ot)
2127 {
2128         /* identifiers */
2129         ot->name = "Set Curve Radius";
2130         ot->description = "Set per-point radius which is used for bevel tapering";
2131         ot->idname = "CURVE_OT_radius_set";
2132         
2133         /* api callbacks */
2134         ot->exec = set_radius_exec;
2135         ot->invoke = WM_operator_props_popup;
2136         ot->poll = ED_operator_editsurfcurve;
2137
2138         /* flags */
2139         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2140
2141         /* properties */
2142         RNA_def_float(ot->srna, "radius", 1.0f, 0.0f, FLT_MAX, "Radius", "", 0.0001f, 10.0f);
2143 }
2144
2145 /********************* smooth operator ********************/
2146
2147 static int smooth_exec(bContext *C, wmOperator *UNUSED(op))
2148 {
2149         Object *obedit = CTX_data_edit_object(C);
2150         ListBase *editnurb = object_editcurve_get(obedit);
2151         Nurb *nu;
2152         BezTriple *bezt, *beztOrig;
2153         BPoint *bp, *bpOrig;
2154         float val, newval, offset;
2155         int a, i, change = 0;
2156         
2157         for (nu = editnurb->first; nu; nu = nu->next) {
2158                 if (nu->bezt) {
2159                         change = 0;
2160                         beztOrig = MEM_dupallocN(nu->bezt);
2161                         for (bezt = &nu->bezt[1], a = 1; a < nu->pntsu - 1; a++, bezt++) {
2162                                 if (bezt->f2 & SELECT) {
2163                                         for (i = 0; i < 3; i++) {
2164                                                 val = bezt->vec[1][i];
2165                                                 newval = ((beztOrig + (a - 1))->vec[1][i] * 0.5f) + ((beztOrig + (a + 1))->vec[1][i] * 0.5f);
2166                                                 offset = (val * ((1.0f / 6.0f) * 5.0f)) + (newval * (1.0f / 6.0f)) - val;
2167                                                 /* offset handles */
2168                                                 bezt->vec[1][i] += offset;
2169                                                 bezt->vec[0][i] += offset;
2170                                                 bezt->vec[2][i] += offset;
2171                                         }
2172                                         change = 1;
2173                                 }
2174                         }
2175                         MEM_freeN(beztOrig);
2176                         if (change)
2177                                 BKE_nurb_handles_calc(nu);
2178                 }
2179                 else if (nu->bp) {
2180                         bpOrig = MEM_dupallocN(nu->bp);
2181                         /* Same as above, keep these the same! */
2182                         for (bp = &nu->bp[1], a = 1; a < nu->pntsu - 1; a++, bp++) {
2183                                 if (bp->f1 & SELECT) {
2184                                         for (i = 0; i < 3; i++) {
2185                                                 val = bp->vec[i];
2186                                                 newval = ((bpOrig + (a - 1))->vec[i] * 0.5f) + ((bpOrig + (a + 1))->vec[i] * 0.5f);
2187                                                 offset = (val * ((1.0f / 6.0f) * 5.0f)) + (newval * (1.0f / 6.0f)) - val;
2188                                         
2189                                                 bp->vec[i] += offset;
2190                                         }
2191                                 }
2192                         }
2193                         MEM_freeN(bpOrig);
2194                 }
2195         }
2196
2197         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2198         DAG_id_tag_update(obedit->data, 0);
2199
2200         return OPERATOR_FINISHED;
2201 }
2202
2203 void CURVE_OT_smooth(wmOperatorType *ot)
2204 {
2205         /* identifiers */
2206         ot->name = "Smooth";
2207         ot->description = "Flatten angles of selected points";
2208         ot->idname = "CURVE_OT_smooth";
2209         
2210         /* api callbacks */
2211         ot->exec = smooth_exec;
2212         ot->poll = ED_operator_editsurfcurve;
2213
2214         /* flags */
2215         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2216 }
2217
2218 /**************** smooth curve radius operator *************/
2219
2220 /* TODO, make smoothing distance based */
2221 static int smooth_radius_exec(bContext *C, wmOperator *UNUSED(op))
2222 {
2223         Object *obedit = CTX_data_edit_object(C);
2224         ListBase *editnurb = object_editcurve_get(obedit);
2225         Nurb *nu;
2226         BezTriple *bezt;
2227         BPoint *bp;
2228         int a;
2229         
2230         /* use for smoothing */
2231         int last_sel;
2232         int start_sel, end_sel; /* selection indices, inclusive */
2233         float start_rad, end_rad, fac, range;
2234         
2235         for (nu = editnurb->first; nu; nu = nu->next) {
2236                 if (nu->bezt) {
2237                         
2238                         for (last_sel = 0; last_sel < nu->pntsu; last_sel++) {
2239                                 /* loop over selection segments of a curve, smooth each */
2240                                 
2241                                 /* Start BezTriple code, this is duplicated below for points, make sure these functions stay in sync */
2242                                 start_sel = -1;
2243                                 for (bezt = &nu->bezt[last_sel], a = last_sel; a < nu->pntsu; a++, bezt++) {
2244                                         if (bezt->f2 & SELECT) {
2245                                                 start_sel = a;
2246                                                 break;
2247                                         }
2248                                 }
2249                                 /* in case there are no other selected verts */
2250                                 end_sel = start_sel;
2251                                 for (bezt = &nu->bezt[start_sel + 1], a = start_sel + 1; a < nu->pntsu; a++, bezt++) {
2252                                         if ((bezt->f2 & SELECT) == 0) {
2253                                                 break;
2254                                         }
2255                                         end_sel = a;
2256                                 }
2257                                 
2258                                 if (start_sel == -1) {
2259                                         last_sel = nu->pntsu; /* next... */
2260                                 }
2261                                 else {
2262                                         last_sel = end_sel; /* before we modify it */
2263                                         
2264                                         /* now blend between start and end sel */
2265                                         start_rad = end_rad = -1.0;
2266                                         
2267                                         if (start_sel == end_sel) {
2268                                                 /* simple, only 1 point selected */
2269                                                 if (start_sel > 0)                         start_rad = nu->bezt[start_sel - 1].radius;
2270                                                 if (end_sel != -1 && end_sel < nu->pntsu)  end_rad   = nu->bezt[start_sel + 1].radius;
2271                                                 
2272                                                 if      (start_rad >= 0.0f && end_rad >= 0.0f)  nu->bezt[start_sel].radius = (start_rad + end_rad) / 2.0f;
2273                                                 else if (start_rad >= 0.0f)                     nu->bezt[start_sel].radius = start_rad;
2274                                                 else if (end_rad >= 0.0f)                       nu->bezt[start_sel].radius = end_rad;
2275                                         }
2276                                         else {
2277                                                 /* if endpoints selected, then use them */
2278                                                 if (start_sel == 0) {
2279                                                         start_rad = nu->bezt[start_sel].radius;
2280                                                         start_sel++; /* we don't want to edit the selected endpoint */
2281                                                 }
2282                                                 else {
2283                                                         start_rad = nu->bezt[start_sel - 1].radius;
2284                                                 }
2285                                                 if (end_sel == nu->pntsu - 1) {
2286                                                         end_rad = nu->bezt[end_sel].radius;
2287                                                         end_sel--; /* we don't want to edit the selected endpoint */
2288                                                 }
2289                                                 else {
2290                                                         end_rad = nu->bezt[end_sel + 1].radius;
2291                                                 }
2292                                                 
2293                                                 /* Now Blend between the points */
2294                                                 range = (float)(end_sel - start_sel) + 2.0f;
2295                                                 for (bezt = &nu->bezt[start_sel], a = start_sel; a <= end_sel; a++, bezt++) {
2296                                                         fac = (float)(1 + a - start_sel) / range;
2297                                                         bezt->radius = start_rad * (1.0f - fac) + end_rad * fac;
2298                                                 }
2299                                         }
2300                                 }
2301                         }
2302                 }
2303                 else if (nu->bp) {
2304                         /* Same as above, keep these the same! */
2305                         for (last_sel = 0; last_sel < nu->pntsu; last_sel++) {
2306                                 /* loop over selection segments of a curve, smooth each */
2307                                 
2308                                 /* Start BezTriple code, this is duplicated below for points, make sure these functions stay in sync */
2309                                 start_sel = -1;
2310                                 for (bp = &nu->bp[last_sel], a = last_sel; a < nu->pntsu; a++, bp++) {
2311                                         if (bp->f1 & SELECT) {
2312                                                 start_sel = a;
2313                                                 break;
2314                                         }
2315                                 }
2316                                 /* in case there are no other selected verts */
2317                                 end_sel = start_sel;
2318                                 for (bp = &nu->bp[start_sel + 1], a = start_sel + 1; a < nu->pntsu; a++, bp++) {
2319                                         if ((bp->f1 & SELECT) == 0) {
2320                                                 break;
2321                                         }
2322                                         end_sel = a;
2323                                 }
2324                                 
2325                                 if (start_sel == -1) {
2326                                         last_sel = nu->pntsu; /* next... */
2327                                 }
2328                                 else {
2329                                         last_sel = end_sel; /* before we modify it */
2330                                         
2331                                         /* now blend between start and end sel */
2332                                         start_rad = end_rad = -1.0;
2333                                         
2334                                         if (start_sel == end_sel) {
2335                                                 /* simple, only 1 point selected */
2336                                                 if (start_sel > 0) start_rad = nu->bp[start_sel - 1].radius;
2337                                                 if (end_sel != -1 && end_sel < nu->pntsu) end_rad = nu->bp[start_sel + 1].radius;
2338                                                 
2339                                                 if      (start_rad >= 0.0f && end_rad >= 0.0f)  nu->bp[start_sel].radius = (start_rad + end_rad) / 2;
2340                                                 else if (start_rad >= 0.0f)                     nu->bp[start_sel].radius = start_rad;
2341                                                 else if (end_rad >= 0.0f)                       nu->bp[start_sel].radius = end_rad;
2342                                         }
2343                                         else {
2344                                                 /* if endpoints selected, then use them */
2345                                                 if (start_sel == 0) {
2346                                                         start_rad = nu->bp[start_sel].radius;
2347                                                         start_sel++; /* we don't want to edit the selected endpoint */
2348                                                 }
2349                                                 else {
2350                                                         start_rad = nu->bp[start_sel - 1].radius;
2351                                                 }
2352                                                 if (end_sel == nu->pntsu - 1) {
2353                                                         end_rad = nu->bp[end_sel].radius;
2354                                                         end_sel--; /* we don't want to edit the selected endpoint */
2355                                                 }
2356                                                 else {
2357                                                         end_rad = nu->bp[end_sel + 1].radius;
2358                                                 }
2359                                                 
2360                                                 /* Now Blend between the points */
2361                                                 range = (float)(end_sel - start_sel) + 2.0f;
2362                                                 for (bp = &nu->bp[start_sel], a = start_sel; a <= end_sel; a++, bp++) {
2363                                                         fac = (float)(1 + a - start_sel) / range;
2364                                                         bp->radius = start_rad * (1.0f - fac) + end_rad * fac;
2365                                                 }
2366                                         }
2367                                 }
2368                         }
2369                 }
2370         }
2371
2372         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2373         DAG_id_tag_update(obedit->data, 0);
2374
2375         return OPERATOR_FINISHED;
2376 }
2377
2378 void CURVE_OT_smooth_radius(wmOperatorType *ot)
2379 {
2380         /* identifiers */
2381         ot->name = "Smooth Curve Radius";
2382         ot->description = "Flatten radii of selected points";
2383         ot->idname = "CURVE_OT_smooth_radius";
2384         
2385         /* api clastbacks */
2386         ot->exec = smooth_radius_exec;
2387         ot->poll = ED_operator_editsurfcurve;
2388         
2389         /* flags */
2390         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2391 }
2392
2393 /***************** selection utility *************************/
2394
2395 /* next == 1 -> select next         */
2396 /* next == -1 -> select previous    */
2397 /* cont == 1 -> select continuously */
2398 /* selstatus, inverts behavior          */
2399 static void select_adjacent_cp(ListBase *editnurb, short next, short cont, short selstatus)
2400 {
2401         Nurb *nu;
2402         BezTriple *bezt;
2403         BPoint *bp;
2404         int a;
2405         short lastsel = 0;
2406         
2407         if (next == 0) return;
2408         
2409         for (nu = editnurb->first; nu; nu = nu->next) {
2410                 lastsel = 0;
2411                 if (nu->type == CU_BEZIER) {
2412                         a = nu->pntsu;
2413                         bezt = nu->bezt;
2414                         if (next < 0) bezt = &nu->bezt[a - 1];
2415                         while (a--) {
2416                                 if (a - abs(next) < 0) break;
2417                                 if ((lastsel == 0) && (bezt->hide == 0) && ((bezt->f2 & SELECT) || (selstatus == 0))) {
2418                                         bezt += next;
2419                                         if (!(bezt->f2 & SELECT) || (selstatus == 0)) {
2420                                                 short sel = select_beztriple(bezt, selstatus, 1, VISIBLE);
2421                                                 if ((sel == 1) && (cont == 0)) lastsel = 1;
2422                                         }
2423                                 }
2424                                 else {
2425                                         bezt += next;
2426                                         lastsel = 0;
2427                                 }
2428                                 /* move around in zigzag way so that we go through each */
2429                                 bezt -= (next - next / abs(next));
2430                         }
2431                 }
2432                 else {
2433                         a = nu->pntsu * nu->pntsv;
2434                         bp = nu->bp;
2435                         if (next < 0) bp = &nu->bp[a - 1];
2436                         while (a--) {
2437                                 if (a - abs(next) < 0) break;
2438                                 if ((lastsel == 0) && (bp->hide == 0) && ((bp->f1 & SELECT) || (selstatus == 0))) {
2439                                         bp += next;
2440                                         if (!(bp->f1 & SELECT) || (selstatus == 0)) {
2441                                                 short sel = select_bpoint(bp, selstatus, 1, VISIBLE);
2442                                                 if ((sel == 1) && (cont == 0)) lastsel = 1;
2443                                         }
2444                                 }
2445                                 else {
2446                                         bp += next;
2447                                         lastsel = 0;
2448                                 }
2449                                 /* move around in zigzag way so that we go through each */
2450                                 bp -= (next - next / abs(next));
2451                         }
2452                 }
2453         }
2454 }
2455
2456 /**************** select start/end operators **************/
2457
2458 /* (de)selects first or last of visible part of each Nurb depending on selFirst     */
2459 /* selFirst: defines the end of which to select                                     */
2460 /* doswap: defines if selection state of each first/last control point is swapped   */
2461 /* selstatus: selection status in case doswap is false                              */
2462 void selectend_nurb(Object *obedit, short selfirst, short doswap, short selstatus)
2463 {
2464         ListBase *editnurb = object_editcurve_get(obedit);
2465         Nurb *nu;
2466         BPoint *bp;
2467         BezTriple *bezt;
2468         Curve *cu;
2469         int a;
2470
2471         if (obedit == NULL) return;
2472
2473         cu = (Curve *)obedit->data;
2474         cu->lastsel = NULL;
2475
2476         for (nu = editnurb->first; nu; nu = nu->next) {
2477                 if (nu->type == CU_BEZIER) {
2478                         a = nu->pntsu;
2479                         
2480                         /* which point? */
2481                         if (selfirst == 0) { /* select last */
2482                                 bezt = &nu->bezt[a - 1];
2483                         }
2484                         else { /* select first */
2485                                 bezt = nu->bezt;
2486                         }
2487                         
2488                         while (a--) {
2489                                 short sel;
2490                                 if (doswap) sel = swap_selection_beztriple(bezt);
2491                                 else sel = select_beztriple(bezt, selstatus, 1, VISIBLE);
2492                                 
2493                                 if (sel == 1) break;
2494                         }
2495                 }
2496                 else {
2497                         a = nu->pntsu * nu->pntsv;
2498                         
2499                         /* which point? */
2500                         if (selfirst == 0) { /* select last */
2501                                 bp = &nu->bp[a - 1];
2502                         }
2503                         else { /* select first */
2504                                 bp = nu->bp;
2505                         }
2506
2507                         while (a--) {
2508                                 if (bp->hide == 0) {
2509                                         short sel;
2510                                         if (doswap) sel = swap_selection_bpoint(bp);
2511                                         else sel = select_bpoint(bp, selstatus, 1, VISIBLE);
2512                                         
2513                                         if (sel == 1) break;
2514                                 }
2515                         }
2516                 }
2517         }
2518 }
2519
2520 static int de_select_first_exec(bContext *C, wmOperator *UNUSED(op))
2521 {
2522         Object *obedit = CTX_data_edit_object(C);
2523
2524         selectend_nurb(obedit, FIRST, 1, DESELECT);
2525         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2526
2527         return OPERATOR_FINISHED;
2528 }
2529
2530 void CURVE_OT_de_select_first(wmOperatorType *ot)
2531 {
2532         /* identifiers */
2533         ot->name = "(De)select First";
2534         ot->idname = "CURVE_OT_de_select_first";
2535         ot->description = "(De)select first of visible part of each Nurb";
2536         
2537         /* api cfirstbacks */
2538         ot->exec = de_select_first_exec;
2539         ot->poll = ED_operator_editcurve;
2540         
2541         /* flags */
2542         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2543 }
2544
2545 static int de_select_last_exec(bContext *C, wmOperator *UNUSED(op))
2546 {
2547         Object *obedit = CTX_data_edit_object(C);
2548
2549         selectend_nurb(obedit, LAST, 1, DESELECT);
2550         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2551
2552         return OPERATOR_FINISHED;
2553 }
2554
2555 void CURVE_OT_de_select_last(wmOperatorType *ot)
2556 {
2557         /* identifiers */
2558         ot->name = "(De)select Last";
2559         ot->idname = "CURVE_OT_de_select_last";
2560         ot->description = "(De)select last of visible part of each Nurb";
2561         
2562         /* api clastbacks */
2563         ot->exec = de_select_last_exec;
2564         ot->poll = ED_operator_editcurve;
2565         
2566         /* flags */
2567         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2568 }
2569
2570 /******************* de select all operator ***************/
2571
2572 static short nurb_has_selected_cps(ListBase *editnurb)
2573 {
2574         Nurb *nu;
2575         BezTriple *bezt;
2576         BPoint *bp;
2577         int a;
2578
2579         for (nu = editnurb->first; nu; nu = nu->next) {
2580                 if (nu->type == CU_BEZIER) {
2581                         a = nu->pntsu;
2582                         bezt = nu->bezt;
2583                         while (a--) {
2584                                 if (bezt->hide == 0) {
2585                                         if ((bezt->f1 & SELECT) ||
2586                                             (bezt->f2 & SELECT) ||
2587                                             (bezt->f3 & SELECT))
2588                                         {
2589                                                 return 1;
2590                                         }
2591                                 }
2592                                 bezt++;
2593                         }
2594                 }
2595                 else {
2596                         a = nu->pntsu * nu->pntsv;
2597                         bp = nu->bp;
2598                         while (a--) {
2599                                 if ((bp->hide == 0) && (bp->f1 & SELECT)) return 1;
2600                                 bp++;
2601                         }
2602                 }
2603         }
2604         
2605         return 0;
2606 }
2607
2608 static int de_select_all_exec(bContext *C, wmOperator *op)
2609 {
2610         Object *obedit = CTX_data_edit_object(C);
2611         ListBase *editnurb = object_editcurve_get(obedit);
2612         int action = RNA_enum_get(op->ptr, "action");
2613
2614         if (action == SEL_TOGGLE) {
2615                 action = SEL_SELECT;
2616                 if (nurb_has_selected_cps(editnurb))
2617                         action = SEL_DESELECT;
2618         }
2619
2620         switch (action) {
2621                 case SEL_SELECT:
2622                         CU_select_all(obedit);
2623                         break;
2624                 case SEL_DESELECT:
2625                         CU_deselect_all(obedit);
2626                         break;
2627                 case SEL_INVERT:
2628                         CU_select_swap(obedit);
2629                         break;
2630         }
2631         
2632         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2633
2634         return OPERATOR_FINISHED;
2635 }
2636
2637 void CURVE_OT_select_all(wmOperatorType *ot)
2638 {
2639         /* identifiers */
2640         ot->name = "(De)select All";
2641         ot->idname = "CURVE_OT_select_all";
2642         ot->description = "(De)select all control points";
2643         
2644         /* api callbacks */
2645         ot->exec = de_select_all_exec;
2646         ot->poll = ED_operator_editsurfcurve;
2647         
2648         /* flags */
2649         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2650
2651         /* properties */
2652         WM_operator_properties_select_all(ot);
2653 }
2654
2655 /********************** hide operator *********************/
2656
2657 static int hide_exec(bContext *C, wmOperator *op)
2658 {
2659         Object *obedit = CTX_data_edit_object(C);
2660         Curve *cu = obedit->data;
2661         ListBase *editnurb = object_editcurve_get(obedit);
2662         Nurb *nu;
2663         BPoint *bp;
2664         BezTriple *bezt;
2665         int a, sel, invert = RNA_boolean_get(op->ptr, "unselected");
2666
2667         for (nu = editnurb->first; nu; nu = nu->next) {
2668                 if (nu->type == CU_BEZIER) {
2669                         bezt = nu->bezt;
2670                         a = nu->pntsu;
2671                         sel = 0;
2672                         while (a--) {
2673                                 if (invert == 0 && BEZSELECTED_HIDDENHANDLES(cu, bezt)) {
2674                                         select_beztriple(bezt, DESELECT, 1, HIDDEN);
2675                                         bezt->hide = 1;
2676                                 }
2677                                 else if (invert && !BEZSELECTED_HIDDENHANDLES(cu, bezt)) {
2678                                         select_beztriple(bezt, DESELECT, 1, HIDDEN);
2679                                         bezt->hide = 1;
2680                                 }
2681                                 if (bezt->hide) sel++;
2682                                 bezt++;
2683                         }
2684                         if (sel == nu->pntsu) nu->hide = 1;
2685                 }
2686                 else {
2687                         bp = nu->bp;
2688                         a = nu->pntsu * nu->pntsv;
2689                         sel = 0;
2690                         while (a--) {
2691                                 if (invert == 0 && (bp->f1 & SELECT)) {
2692                                         select_bpoint(bp, DESELECT, 1, HIDDEN);
2693                                         bp->hide = 1;
2694                                 }
2695                                 else if (invert && (bp->f1 & SELECT) == 0) {
2696                                         select_bpoint(bp, DESELECT, 1, HIDDEN);
2697                                         bp->hide = 1;
2698                                 }
2699                                 if (bp->hide) sel++;
2700                                 bp++;
2701                         }
2702                         if (sel == nu->pntsu * nu->pntsv) nu->hide = 1;
2703                 }
2704         }
2705
2706         DAG_id_tag_update(obedit->data, 0);
2707         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2708
2709         return OPERATOR_FINISHED;
2710 }
2711
2712 void CURVE_OT_hide(wmOperatorType *ot)
2713 {
2714         /* identifiers */
2715         ot->name = "Hide Selected";
2716         ot->idname = "CURVE_OT_hide";
2717         ot->description = "Hide (un)selected control points";
2718         
2719         /* api callbacks */
2720         ot->exec = hide_exec;
2721         ot->poll = ED_operator_editsurfcurve;
2722         
2723         /* flags */
2724         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2725         
2726         /* props */
2727         RNA_def_boolean(ot->srna, "unselected", 0, "Unselected", "Hide unselected rather than selected");
2728 }
2729
2730 /********************** reveal operator *********************/
2731
2732 static int reveal_exec(bContext *C, wmOperator *UNUSED(op))
2733 {
2734         Object *obedit = CTX_data_edit_object(C);
2735         ListBase *editnurb = object_editcurve_get(obedit);
2736         Nurb *nu;
2737         BPoint *bp;
2738         BezTriple *bezt;
2739         int a;
2740
2741         for (nu = editnurb->first; nu; nu = nu->next) {
2742                 nu->hide = 0;
2743                 if (nu->type == CU_BEZIER) {
2744                         bezt = nu->bezt;
2745                         a = nu->pntsu;
2746                         while (a--) {
2747                                 if (bezt->hide) {
2748                                         select_beztriple(bezt, SELECT, 1, HIDDEN);
2749                                         bezt->hide = 0;
2750                                 }
2751                                 bezt++;
2752                         }
2753                 }
2754                 else {
2755                         bp = nu->bp;
2756                         a = nu->pntsu * nu->pntsv;
2757                         while (a--) {
2758                                 if (bp->hide) {
2759                                         select_bpoint(bp, SELECT, 1, HIDDEN);
2760                                         bp->hide = 0;
2761                                 }
2762                                 bp++;
2763                         }
2764                 }
2765         }
2766
2767         DAG_id_tag_update(obedit->data, 0);
2768         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2769
2770         return OPERATOR_FINISHED;
2771 }
2772
2773 void CURVE_OT_reveal(wmOperatorType *ot)
2774 {
2775         /* identifiers */
2776         ot->name = "Reveal Hidden";
2777         ot->idname = "CURVE_OT_reveal";
2778         ot->description = "Show again hidden control points";
2779         
2780         /* api callbacks */
2781         ot->exec = reveal_exec;
2782         ot->poll = ED_operator_editsurfcurve;
2783         
2784         /* flags */
2785         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2786 }
2787
2788 /********************** subdivide operator *********************/
2789
2790 /** Divide the line segments associated with the currently selected
2791  * curve nodes (Bezier or NURB). If there are no valid segment
2792  * selections within the current selection, nothing happens.
2793  */
2794 static void subdividenurb(Object *obedit, int number_cuts)
2795 {
2796         Curve *cu = obedit->data;
2797         EditNurb *editnurb = cu->editnurb;
2798         Nurb *nu;
2799         BezTriple *prevbezt, *bezt, *beztnew, *beztn;
2800         BPoint *bp, *prevbp, *bpnew, *bpn;
2801         float vec[15];
2802         int a, b, sel, amount, *usel, *vsel, i;
2803         float factor;
2804
2805         // printf("*** subdivideNurb: entering subdivide\n");
2806
2807         for (nu = editnurb->nurbs.first; nu; nu = nu->next) {
2808                 amount = 0;
2809                 if (nu->type == CU_BEZIER) {
2810                         /*
2811                          * Insert a point into a 2D Bezier curve.
2812                          * Endpoints are preserved. Otherwise, all selected and inserted points are
2813                          * newly created. Old points are discarded.
2814                          */
2815                         /* count */
2816                         if (nu->flagu & CU_NURB_CYCLIC) {
2817                                 a = nu->pntsu;
2818                                 bezt = nu->bezt;
2819                                 prevbezt = bezt + (a - 1);
2820                         }
2821                         else {
2822                                 a = nu->pntsu - 1;
2823                                 prevbezt = nu->bezt;
2824                                 bezt = prevbezt + 1;
2825                         }
2826                         while (a--) {
2827                                 if (BEZSELECTED_HIDDENHANDLES(cu, prevbezt) && BEZSELECTED_HIDDENHANDLES(cu, bezt) ) amount += number_cuts;
2828                                 prevbezt = bezt;
2829                                 bezt++;
2830                         }
2831
2832                         if (amount) {
2833                                 /* insert */
2834                                 beztnew = (BezTriple *)MEM_mallocN((amount + nu->pntsu) * sizeof(BezTriple), "subdivNurb");
2835                                 beztn = beztnew;
2836                                 if (nu->flagu & CU_NURB_CYCLIC) {
2837                                         a = nu->pntsu;
2838                                         bezt = nu->bezt;
2839                                         prevbezt = bezt + (a - 1);
2840                                 }
2841                                 else {
2842                                         a = nu->pntsu - 1;
2843                                         prevbezt = nu->bezt;
2844                                         bezt = prevbezt + 1;
2845                                 }
2846                                 while (a--) {
2847                                         memcpy(beztn, prevbezt, sizeof(BezTriple));
2848                                         keyIndex_updateBezt(editnurb, prevbezt, beztn, 1);
2849