fix for type mismatch with SWAP() macro.
[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 callocstructN(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 = ob_get_keyblock(obedit);
1237
1238                 if (actkey) {
1239                         // XXX strcpy(G.editModeTitleExtra, "(Key) ");
1240                         undo_editmode_clear();
1241                         key_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         else return OPERATOR_CANCELLED;
1584
1585         cu->lastsel = NULL;
1586
1587         nu = editnurb->first;
1588         while (nu) {
1589                 next = nu->next;
1590
1591                 /* is entire nurb selected */
1592                 bp = nu->bp;
1593                 a = nu->pntsu * nu->pntsv;
1594                 while (a) {
1595                         a--;
1596                         if (bp->f1 & flag) ;
1597                         else break;
1598                         bp++;
1599                 }
1600                 if (a == 0) {
1601                         BLI_remlink(editnurb, nu);
1602                         keyIndex_delNurb(cu->editnurb, nu);
1603                         BKE_nurb_free(nu); nu = NULL;
1604                 }
1605                 else {
1606                         /* is nurb in U direction selected */
1607                         newv = nu->pntsv;
1608                         bp = nu->bp;
1609                         for (b = 0; b < nu->pntsv; b++) {
1610                                 sel = 0;
1611                                 for (a = 0; a < nu->pntsu; a++, bp++) {
1612                                         if (bp->f1 & flag) sel++;
1613                                 }
1614                                 if (sel == nu->pntsu) {
1615                                         newv--;
1616                                 }
1617                                 else if (sel >= 1) {
1618                                         /* don't delete */
1619                                         break;
1620                                 }
1621                         }
1622                         if (newv != nu->pntsv && b == nu->pntsv) {
1623                                 /* delete */
1624                                 bp = nu->bp;
1625                                 bpn = newbp = (BPoint *)MEM_mallocN(newv * nu->pntsu * sizeof(BPoint), "deleteNurb");
1626                                 for (b = 0; b < nu->pntsv; b++) {
1627                                         if ((bp->f1 & flag) == 0) {
1628                                                 memcpy(bpn, bp, nu->pntsu * sizeof(BPoint));
1629                                                 keyIndex_updateBP(cu->editnurb, bp, bpn, nu->pntsu);
1630                                                 bpn += nu->pntsu;
1631                                         }
1632                                         else {
1633                                                 keyIndex_delBP(cu->editnurb, bp);
1634                                         }
1635                                         bp += nu->pntsu;
1636                                 }
1637                                 nu->pntsv = newv;
1638                                 MEM_freeN(nu->bp);
1639                                 nu->bp = newbp;
1640                                 BKE_nurb_order_clamp_v(nu);
1641
1642                                 BKE_nurb_knot_calc_v(nu);
1643                         }
1644                         else {
1645                                 /* is the nurb in V direction selected */
1646                                 newu = nu->pntsu;
1647                                 for (a = 0; a < nu->pntsu; a++) {
1648                                         bp = &nu->bp[a];
1649                                         sel = 0;
1650                                         for (b = 0; b < nu->pntsv; b++, bp += nu->pntsu) {
1651                                                 if (bp->f1 & flag) sel++;
1652                                         }
1653                                         if (sel == nu->pntsv) {
1654                                                 newu--;
1655                                         }
1656                                         else if (sel >= 1) {
1657                                                 /* don't delete */
1658                                                 break;
1659                                         }
1660                                 }
1661                                 if (newu != nu->pntsu && a == nu->pntsu) {
1662                                         /* delete */
1663                                         bp = nu->bp;
1664                                         bpn = newbp = (BPoint *)MEM_mallocN(newu * nu->pntsv * sizeof(BPoint), "deleteNurb");
1665                                         for (b = 0; b < nu->pntsv; b++) {
1666                                                 for (a = 0; a < nu->pntsu; a++, bp++) {
1667                                                         if ((bp->f1 & flag) == 0) {
1668                                                                 *bpn = *bp;
1669                                                                 keyIndex_updateBP(cu->editnurb, bp, bpn, 1);
1670                                                                 bpn++;
1671                                                         }
1672                                                         else {
1673                                                                 keyIndex_delBP(cu->editnurb, bp);
1674                                                         }
1675                                                 }
1676                                         }
1677                                         MEM_freeN(nu->bp);
1678                                         nu->bp = newbp;
1679                                         if (newu == 1 && nu->pntsv > 1) {    /* make a U spline */
1680                                                 nu->pntsu = nu->pntsv;
1681                                                 nu->pntsv = 1;
1682                                                 SWAP(short, nu->orderu, nu->orderv);
1683                                                 BKE_nurb_order_clamp_u(nu);
1684                                                 if (nu->knotsv) MEM_freeN(nu->knotsv);
1685                                                 nu->knotsv = NULL;
1686                                         }
1687                                         else {
1688                                                 nu->pntsu = newu;
1689                                                 BKE_nurb_order_clamp_u(nu);
1690                                         }
1691                                         BKE_nurb_knot_calc_u(nu);
1692                                 }
1693                         }
1694                 }
1695                 nu = next;
1696         }
1697
1698         if (ED_curve_updateAnimPaths(obedit->data))
1699                 WM_event_add_notifier(C, NC_OBJECT | ND_KEYS, obedit);
1700
1701         return OPERATOR_FINISHED;
1702 }
1703
1704 /* only for OB_SURF */
1705 static short extrudeflagNurb(EditNurb *editnurb, int flag)
1706 {
1707         Nurb *nu;
1708         BPoint *bp, *bpn, *newbp;
1709         int ok = 0, a, u, v, len;
1710
1711         nu = editnurb->nurbs.first;
1712         while (nu) {
1713
1714                 if (nu->pntsv == 1) {
1715                         bp = nu->bp;
1716                         a = nu->pntsu;
1717                         while (a) {
1718                                 if (bp->f1 & flag) ;
1719                                 else break;
1720                                 bp++;
1721                                 a--;
1722                         }
1723                         if (a == 0) {
1724                                 ok = 1;
1725                                 newbp = (BPoint *)MEM_mallocN(2 * nu->pntsu * sizeof(BPoint), "extrudeNurb1");
1726                                 ED_curve_bpcpy(editnurb, newbp, nu->bp, nu->pntsu);
1727                                 bp = newbp + nu->pntsu;
1728                                 ED_curve_bpcpy(editnurb, bp, nu->bp, nu->pntsu);
1729                                 MEM_freeN(nu->bp);
1730                                 nu->bp = newbp;
1731                                 a = nu->pntsu;
1732                                 while (a--) {
1733                                         select_bpoint(bp, SELECT, flag, HIDDEN);
1734                                         select_bpoint(newbp, DESELECT, flag, HIDDEN);
1735                                         bp++; 
1736                                         newbp++;
1737                                 }
1738
1739                                 nu->pntsv = 2;
1740                                 nu->orderv = 2;
1741                                 BKE_nurb_knot_calc_v(nu);
1742                         }
1743                 }
1744                 else {
1745                         /* which row or column is selected */
1746
1747                         if (isNurbselUV(nu, &u, &v, flag) ) {
1748
1749                                 /* deselect all */
1750                                 bp = nu->bp;
1751                                 a = nu->pntsu * nu->pntsv;
1752                                 while (a--) {
1753                                         select_bpoint(bp, DESELECT, flag, HIDDEN);
1754                                         bp++;
1755                                 }
1756
1757                                 if (u == 0 || u == nu->pntsv - 1) {      /* row in u-direction selected */
1758                                         ok = 1;
1759                                         newbp = (BPoint *)MEM_mallocN(nu->pntsu * (nu->pntsv + 1) *
1760                                                                       sizeof(BPoint), "extrudeNurb1");
1761                                         if (u == 0) {
1762                                                 len = nu->pntsv * nu->pntsu;
1763                                                 ED_curve_bpcpy(editnurb, newbp + nu->pntsu, nu->bp, len);
1764                                                 ED_curve_bpcpy(editnurb, newbp, nu->bp, nu->pntsu);
1765                                                 bp = newbp;
1766                                         }
1767                                         else {
1768                                                 len = nu->pntsv * nu->pntsu;
1769                                                 ED_curve_bpcpy(editnurb, newbp, nu->bp, len);
1770                                                 ED_curve_bpcpy(editnurb, newbp + len, &nu->bp[len - nu->pntsu], nu->pntsu);
1771                                                 bp = newbp + len;
1772                                         }
1773
1774                                         a = nu->pntsu;
1775                                         while (a--) {
1776                                                 select_bpoint(bp, SELECT, flag, HIDDEN);
1777                                                 bp++;
1778                                         }
1779
1780                                         MEM_freeN(nu->bp);
1781                                         nu->bp = newbp;
1782                                         nu->pntsv++;
1783                                         BKE_nurb_knot_calc_v(nu);
1784                                 }
1785                                 else if (v == 0 || v == nu->pntsu - 1) {     /* column in v-direction selected */
1786                                         ok = 1;
1787                                         bpn = newbp = (BPoint *)MEM_mallocN((nu->pntsu + 1) * nu->pntsv * sizeof(BPoint), "extrudeNurb1");
1788                                         bp = nu->bp;
1789
1790                                         for (a = 0; a < nu->pntsv; a++) {
1791                                                 if (v == 0) {
1792                                                         *bpn = *bp;
1793                                                         bpn->f1 |= flag;
1794                                                         bpn++;
1795                                                 }
1796                                                 ED_curve_bpcpy(editnurb, bpn, bp, nu->pntsu);
1797                                                 bp += nu->pntsu;
1798                                                 bpn += nu->pntsu;
1799                                                 if (v == nu->pntsu - 1) {
1800                                                         *bpn = *(bp - 1);
1801                                                         bpn->f1 |= flag;
1802                                                         bpn++;
1803                                                 }
1804                                         }
1805
1806                                         MEM_freeN(nu->bp);
1807                                         nu->bp = newbp;
1808                                         nu->pntsu++;
1809                                         BKE_nurb_knot_calc_u(nu);
1810                                 }
1811                         }
1812                 }
1813                 nu = nu->next;
1814         }
1815
1816         return ok;
1817 }
1818
1819 static void adduplicateflagNurb(Object *obedit, short flag)
1820 {
1821         ListBase *editnurb = object_editcurve_get(obedit);
1822         Nurb *nu, *newnu;
1823         BezTriple *bezt, *bezt1;
1824         BPoint *bp, *bp1;
1825         Curve *cu = (Curve *)obedit->data;
1826         int a, b, starta, enda, newu, newv;
1827         char *usel;
1828
1829         cu->lastsel = NULL;
1830
1831         nu = editnurb->last;
1832         while (nu) {
1833                 if (nu->type == CU_BEZIER) {
1834                         bezt = nu->bezt;
1835                         for (a = 0; a < nu->pntsu; a++) {
1836                                 enda = -1;
1837                                 starta = a;
1838                                 while ( (bezt->f1 & flag) || (bezt->f2 & flag) || (bezt->f3 & flag) ) {
1839                                         select_beztriple(bezt, DESELECT, flag, HIDDEN);
1840                                         enda = a;
1841                                         if (a >= nu->pntsu - 1) break;
1842                                         a++;
1843                                         bezt++;
1844                                 }
1845                                 if (enda >= starta) {
1846                                         newnu = (Nurb *)MEM_mallocN(sizeof(Nurb), "adduplicateN");
1847                                         memcpy(newnu, nu, sizeof(Nurb));
1848                                         BLI_addtail(editnurb, newnu);
1849                                         set_actNurb(obedit, newnu);
1850                                         newnu->pntsu = enda - starta + 1;
1851                                         newnu->bezt = (BezTriple *)MEM_mallocN((enda - starta + 1) * sizeof(BezTriple), "adduplicateN");
1852                                         memcpy(newnu->bezt, &nu->bezt[starta], newnu->pntsu * sizeof(BezTriple));
1853
1854                                         b = newnu->pntsu;
1855                                         bezt1 = newnu->bezt;
1856                                         while (b--) {
1857                                                 select_beztriple(bezt1, SELECT, flag, HIDDEN);
1858                                                 bezt1++;
1859                                         }
1860
1861                                         if (nu->flagu & CU_NURB_CYCLIC) {
1862                                                 if (starta != 0 || enda != nu->pntsu - 1) {
1863                                                         newnu->flagu &= ~CU_NURB_CYCLIC;
1864                                                 }
1865                                         }
1866                                 }
1867                                 bezt++;
1868                         }
1869                 }
1870                 else if (nu->pntsv == 1) {    /* because UV Nurb has a different method for dupli */
1871                         bp = nu->bp;
1872                         for (a = 0; a < nu->pntsu; a++) {
1873                                 enda = -1;
1874                                 starta = a;
1875                                 while (bp->f1 & flag) {
1876                                         select_bpoint(bp, DESELECT, flag, HIDDEN);
1877                                         enda = a;
1878                                         if (a >= nu->pntsu - 1) break;
1879                                         a++;
1880                                         bp++;
1881                                 }
1882                                 if (enda >= starta) {
1883                                         newnu = (Nurb *)MEM_mallocN(sizeof(Nurb), "adduplicateN3");
1884                                         memcpy(newnu, nu, sizeof(Nurb));
1885                                         set_actNurb(obedit, newnu);
1886                                         BLI_addtail(editnurb, newnu);
1887                                         newnu->pntsu = enda - starta + 1;
1888                                         newnu->bp = (BPoint *)MEM_mallocN((enda - starta + 1) * sizeof(BPoint), "adduplicateN4");
1889                                         memcpy(newnu->bp, &nu->bp[starta], newnu->pntsu * sizeof(BPoint));
1890
1891                                         b = newnu->pntsu;
1892                                         bp1 = newnu->bp;
1893                                         while (b--) {
1894                                                 select_bpoint(bp1, SELECT, flag, HIDDEN);
1895                                                 bp1++;
1896                                         }
1897
1898                                         if (nu->flagu & CU_NURB_CYCLIC) {
1899                                                 if (starta != 0 || enda != nu->pntsu - 1) {
1900                                                         newnu->flagu &= ~CU_NURB_CYCLIC;
1901                                                 }
1902                                         }
1903
1904                                         /* knots */
1905                                         newnu->knotsu = NULL;
1906                                         BKE_nurb_knot_calc_u(newnu);
1907                                 }
1908                                 bp++;
1909                         }
1910                 }
1911                 else {
1912                         /* a rectangular area in nurb has to be selected */
1913                         if (isNurbsel(nu)) {
1914                                 usel = MEM_callocN(nu->pntsu, "adduplicateN4");
1915                                 bp = nu->bp;
1916                                 for (a = 0; a < nu->pntsv; a++) {
1917                                         for (b = 0; b < nu->pntsu; b++, bp++) {
1918                                                 if (bp->f1 & flag) usel[b]++;
1919                                         }
1920                                 }
1921                                 newu = 0;
1922                                 newv = 0;
1923                                 for (a = 0; a < nu->pntsu; a++) {
1924                                         if (usel[a]) {
1925                                                 if (newv == 0 || usel[a] == newv) {
1926                                                         newv = usel[a];
1927                                                         newu++;
1928                                                 }
1929                                                 else {
1930                                                         newv = 0;
1931                                                         break;
1932                                                 }
1933                                         }
1934                                 }
1935                                 if (newu == 0 || newv == 0) {
1936                                         if (G.debug & G_DEBUG)
1937                                                 printf("Can't duplicate Nurb\n");
1938                                 }
1939                                 else {
1940
1941                                         if (newu == 1) SWAP(int, newu, newv);
1942
1943                                         newnu = (Nurb *)MEM_mallocN(sizeof(Nurb), "adduplicateN5");
1944                                         memcpy(newnu, nu, sizeof(Nurb));
1945                                         BLI_addtail(editnurb, newnu);
1946                                         set_actNurb(obedit, newnu);
1947                                         newnu->pntsu = newu;
1948                                         newnu->pntsv = newv;
1949                                         newnu->bp = (BPoint *)MEM_mallocN(newu * newv * sizeof(BPoint), "adduplicateN6");
1950                                         BKE_nurb_order_clamp_u(newnu);
1951                                         BKE_nurb_order_clamp_v(newnu);
1952                                         
1953                                         newnu->knotsu = newnu->knotsv = NULL;
1954                                         
1955                                         bp = newnu->bp;
1956                                         bp1 = nu->bp;
1957                                         for (a = 0; a < nu->pntsv; a++) {
1958                                                 for (b = 0; b < nu->pntsu; b++, bp1++) {
1959                                                         if (bp1->f1 & flag) {
1960                                                                 memcpy(bp, bp1, sizeof(BPoint));
1961                                                                 select_bpoint(bp1, DESELECT, flag, HIDDEN);
1962                                                                 bp++;
1963                                                         }
1964                                                 }
1965                                         }
1966                                         if (BKE_nurb_check_valid_u(newnu)) {
1967                                                 if (nu->pntsu == newnu->pntsu && nu->knotsu) {
1968                                                         newnu->knotsu = MEM_dupallocN(nu->knotsu);
1969                                                 }
1970                                                 else {
1971                                                         BKE_nurb_knot_calc_u(newnu);
1972                                                 }
1973                                         }
1974                                         if (BKE_nurb_check_valid_v(newnu)) {
1975                                                 if (nu->pntsv == newnu->pntsv && nu->knotsv) {
1976                                                         newnu->knotsv = MEM_dupallocN(nu->knotsv);
1977                                                 }
1978                                                 else {
1979                                                         BKE_nurb_knot_calc_v(newnu);
1980                                                 }
1981                                         }
1982                                 }
1983                                 MEM_freeN(usel);
1984                         }
1985                 }
1986
1987                 nu = nu->prev;
1988         }
1989         
1990         /* actnu changed */
1991 }
1992
1993 /**************** switch direction operator ***************/
1994
1995 static int switch_direction_exec(bContext *C, wmOperator *UNUSED(op))
1996 {
1997         Object *obedit = CTX_data_edit_object(C);
1998         Curve *cu = (Curve *)obedit->data;
1999         EditNurb *editnurb = cu->editnurb;
2000         Nurb *nu;
2001
2002         for (nu = editnurb->nurbs.first; nu; nu = nu->next)
2003                 if (isNurbsel(nu)) {
2004                         BKE_nurb_direction_switch(nu);
2005                         keyData_switchDirectionNurb(cu, nu);
2006                 }
2007
2008         if (ED_curve_updateAnimPaths(obedit->data))
2009                 WM_event_add_notifier(C, NC_OBJECT | ND_KEYS, obedit);
2010
2011         DAG_id_tag_update(obedit->data, 0);
2012         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2013
2014         return OPERATOR_FINISHED;
2015 }
2016
2017 void CURVE_OT_switch_direction(wmOperatorType *ot)
2018 {
2019         /* identifiers */
2020         ot->name = "Switch Direction";
2021         ot->description = "Switch direction of selected splines";
2022         ot->idname = "CURVE_OT_switch_direction";
2023         
2024         /* api callbacks */
2025         ot->exec = switch_direction_exec;
2026         ot->poll = ED_operator_editsurfcurve;
2027
2028         /* flags */
2029         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2030 }
2031
2032 /****************** set weight operator *******************/
2033
2034 static int set_goal_weight_exec(bContext *C, wmOperator *op)
2035 {
2036         Object *obedit = CTX_data_edit_object(C);
2037         ListBase *editnurb = object_editcurve_get(obedit);
2038         Nurb *nu;
2039         BezTriple *bezt;
2040         BPoint *bp;
2041         float weight = RNA_float_get(op->ptr, "weight");
2042         int a;
2043                                 
2044         for (nu = editnurb->first; nu; nu = nu->next) {
2045                 if (nu->bezt) {
2046                         for (bezt = nu->bezt, a = 0; a < nu->pntsu; a++, bezt++) {
2047                                 if (bezt->f2 & SELECT)
2048                                         bezt->weight = weight;
2049                         }
2050                 }
2051                 else if (nu->bp) {
2052                         for (bp = nu->bp, a = 0; a < nu->pntsu * nu->pntsv; a++, bp++) {
2053                                 if (bp->f1 & SELECT)
2054                                         bp->weight = weight;
2055                         }
2056                 }
2057         }       
2058
2059         DAG_id_tag_update(obedit->data, 0);
2060         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2061
2062         return OPERATOR_FINISHED;
2063 }
2064
2065 void CURVE_OT_spline_weight_set(wmOperatorType *ot)
2066 {
2067         /* identifiers */
2068         ot->name = "Set Goal Weight";
2069         ot->description = "Set softbody goal weight for selected points";
2070         ot->idname = "CURVE_OT_spline_weight_set";
2071         
2072         /* api callbacks */
2073         ot->exec = set_goal_weight_exec;
2074         ot->invoke = WM_operator_props_popup;
2075         ot->poll = ED_operator_editsurfcurve;
2076
2077         /* flags */
2078         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2079
2080         /* properties */
2081         RNA_def_float_factor(ot->srna, "weight", 1.0f, 0.0f, 1.0f, "Weight", "", 0.0f, 1.0f);
2082 }
2083
2084 /******************* set radius operator ******************/
2085
2086 static int set_radius_exec(bContext *C, wmOperator *op)
2087 {
2088         Object *obedit = CTX_data_edit_object(C);
2089         ListBase *editnurb = object_editcurve_get(obedit);
2090         Nurb *nu;
2091         BezTriple *bezt;
2092         BPoint *bp;
2093         float radius = RNA_float_get(op->ptr, "radius");
2094         int a;
2095         
2096         for (nu = editnurb->first; nu; nu = nu->next) {
2097                 if (nu->bezt) {
2098                         for (bezt = nu->bezt, a = 0; a < nu->pntsu; a++, bezt++) {
2099                                 if (bezt->f2 & SELECT)
2100                                         bezt->radius = radius;
2101                         }
2102                 }
2103                 else if (nu->bp) {
2104                         for (bp = nu->bp, a = 0; a < nu->pntsu * nu->pntsv; a++, bp++) {
2105                                 if (bp->f1 & SELECT)
2106                                         bp->radius = radius;
2107                         }
2108                 }
2109         }       
2110
2111         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2112         DAG_id_tag_update(obedit->data, 0);
2113
2114         return OPERATOR_FINISHED;
2115 }
2116
2117 void CURVE_OT_radius_set(wmOperatorType *ot)
2118 {
2119         /* identifiers */
2120         ot->name = "Set Curve Radius";
2121         ot->description = "Set per-point radius which is used for bevel tapering";
2122         ot->idname = "CURVE_OT_radius_set";
2123         
2124         /* api callbacks */
2125         ot->exec = set_radius_exec;
2126         ot->invoke = WM_operator_props_popup;
2127         ot->poll = ED_operator_editsurfcurve;
2128
2129         /* flags */
2130         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2131
2132         /* properties */
2133         RNA_def_float(ot->srna, "radius", 1.0f, 0.0f, FLT_MAX, "Radius", "", 0.0001f, 10.0f);
2134 }
2135
2136 /********************* smooth operator ********************/
2137
2138 static int smooth_exec(bContext *C, wmOperator *UNUSED(op))
2139 {
2140         Object *obedit = CTX_data_edit_object(C);
2141         ListBase *editnurb = object_editcurve_get(obedit);
2142         Nurb *nu;
2143         BezTriple *bezt, *beztOrig;
2144         BPoint *bp, *bpOrig;
2145         float val, newval, offset;
2146         int a, i, change = 0;
2147         
2148         for (nu = editnurb->first; nu; nu = nu->next) {
2149                 if (nu->bezt) {
2150                         change = 0;
2151                         beztOrig = MEM_dupallocN(nu->bezt);
2152                         for (bezt = &nu->bezt[1], a = 1; a < nu->pntsu - 1; a++, bezt++) {
2153                                 if (bezt->f2 & SELECT) {
2154                                         for (i = 0; i < 3; i++) {
2155                                                 val = bezt->vec[1][i];
2156                                                 newval = ((beztOrig + (a - 1))->vec[1][i] * 0.5f) + ((beztOrig + (a + 1))->vec[1][i] * 0.5f);
2157                                                 offset = (val * ((1.0f / 6.0f) * 5.0f)) + (newval * (1.0f / 6.0f)) - val;
2158                                                 /* offset handles */
2159                                                 bezt->vec[1][i] += offset;
2160                                                 bezt->vec[0][i] += offset;
2161                                                 bezt->vec[2][i] += offset;
2162                                         }
2163                                         change = 1;
2164                                 }
2165                         }
2166                         MEM_freeN(beztOrig);
2167                         if (change)
2168                                 BKE_nurb_handles_calc(nu);
2169                 }
2170                 else if (nu->bp) {
2171                         bpOrig = MEM_dupallocN(nu->bp);
2172                         /* Same as above, keep these the same! */
2173                         for (bp = &nu->bp[1], a = 1; a < nu->pntsu - 1; a++, bp++) {
2174                                 if (bp->f1 & SELECT) {
2175                                         for (i = 0; i < 3; i++) {
2176                                                 val = bp->vec[i];
2177                                                 newval = ((bpOrig + (a - 1))->vec[i] * 0.5f) + ((bpOrig + (a + 1))->vec[i] * 0.5f);
2178                                                 offset = (val * ((1.0f / 6.0f) * 5.0f)) + (newval * (1.0f / 6.0f)) - val;
2179                                         
2180                                                 bp->vec[i] += offset;
2181                                         }
2182                                 }
2183                         }
2184                         MEM_freeN(bpOrig);
2185                 }
2186         }
2187
2188         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2189         DAG_id_tag_update(obedit->data, 0);
2190
2191         return OPERATOR_FINISHED;
2192 }
2193
2194 void CURVE_OT_smooth(wmOperatorType *ot)
2195 {
2196         /* identifiers */
2197         ot->name = "Smooth";
2198         ot->description = "Flatten angles of selected points";
2199         ot->idname = "CURVE_OT_smooth";
2200         
2201         /* api callbacks */
2202         ot->exec = smooth_exec;
2203         ot->poll = ED_operator_editsurfcurve;
2204
2205         /* flags */
2206         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2207 }
2208
2209 /**************** smooth curve radius operator *************/
2210
2211 /* TODO, make smoothing distance based */
2212 static int smooth_radius_exec(bContext *C, wmOperator *UNUSED(op))
2213 {
2214         Object *obedit = CTX_data_edit_object(C);
2215         ListBase *editnurb = object_editcurve_get(obedit);
2216         Nurb *nu;
2217         BezTriple *bezt;
2218         BPoint *bp;
2219         int a;
2220         
2221         /* use for smoothing */
2222         int last_sel;
2223         int start_sel, end_sel; /* selection indices, inclusive */
2224         float start_rad, end_rad, fac, range;
2225         
2226         for (nu = editnurb->first; nu; nu = nu->next) {
2227                 if (nu->bezt) {
2228                         
2229                         for (last_sel = 0; last_sel < nu->pntsu; last_sel++) {
2230                                 /* loop over selection segments of a curve, smooth each */
2231                                 
2232                                 /* Start BezTriple code, this is duplicated below for points, make sure these functions stay in sync */
2233                                 start_sel = -1;
2234                                 for (bezt = &nu->bezt[last_sel], a = last_sel; a < nu->pntsu; a++, bezt++) {
2235                                         if (bezt->f2 & SELECT) {
2236                                                 start_sel = a;
2237                                                 break;
2238                                         }
2239                                 }
2240                                 /* in case there are no other selected verts */
2241                                 end_sel = start_sel;
2242                                 for (bezt = &nu->bezt[start_sel + 1], a = start_sel + 1; a < nu->pntsu; a++, bezt++) {
2243                                         if ((bezt->f2 & SELECT) == 0) {
2244                                                 break;
2245                                         }
2246                                         end_sel = a;
2247                                 }
2248                                 
2249                                 if (start_sel == -1) {
2250                                         last_sel = nu->pntsu; /* next... */
2251                                 }
2252                                 else {
2253                                         last_sel = end_sel; /* before we modify it */
2254                                         
2255                                         /* now blend between start and end sel */
2256                                         start_rad = end_rad = -1.0;
2257                                         
2258                                         if (start_sel == end_sel) {
2259                                                 /* simple, only 1 point selected */
2260                                                 if (start_sel > 0)                         start_rad = nu->bezt[start_sel - 1].radius;
2261                                                 if (end_sel != -1 && end_sel < nu->pntsu)  end_rad   = nu->bezt[start_sel + 1].radius;
2262                                                 
2263                                                 if      (start_rad >= 0.0f && end_rad >= 0.0f)  nu->bezt[start_sel].radius = (start_rad + end_rad) / 2.0f;
2264                                                 else if (start_rad >= 0.0f)                     nu->bezt[start_sel].radius = start_rad;
2265                                                 else if (end_rad >= 0.0f)                       nu->bezt[start_sel].radius = end_rad;
2266                                         }
2267                                         else {
2268                                                 /* if endpoints selected, then use them */
2269                                                 if (start_sel == 0) {
2270                                                         start_rad = nu->bezt[start_sel].radius;
2271                                                         start_sel++; /* we don't want to edit the selected endpoint */
2272                                                 }
2273                                                 else {
2274                                                         start_rad = nu->bezt[start_sel - 1].radius;
2275                                                 }
2276                                                 if (end_sel == nu->pntsu - 1) {
2277                                                         end_rad = nu->bezt[end_sel].radius;
2278                                                         end_sel--; /* we don't want to edit the selected endpoint */
2279                                                 }
2280                                                 else {
2281                                                         end_rad = nu->bezt[end_sel + 1].radius;
2282                                                 }
2283                                                 
2284                                                 /* Now Blend between the points */
2285                                                 range = (float)(end_sel - start_sel) + 2.0f;
2286                                                 for (bezt = &nu->bezt[start_sel], a = start_sel; a <= end_sel; a++, bezt++) {
2287                                                         fac = (float)(1 + a - start_sel) / range;
2288                                                         bezt->radius = start_rad * (1.0f - fac) + end_rad * fac;
2289                                                 }
2290                                         }
2291                                 }
2292                         }
2293                 }
2294                 else if (nu->bp) {
2295                         /* Same as above, keep these the same! */
2296                         for (last_sel = 0; last_sel < nu->pntsu; last_sel++) {
2297                                 /* loop over selection segments of a curve, smooth each */
2298                                 
2299                                 /* Start BezTriple code, this is duplicated below for points, make sure these functions stay in sync */
2300                                 start_sel = -1;
2301                                 for (bp = &nu->bp[last_sel], a = last_sel; a < nu->pntsu; a++, bp++) {
2302                                         if (bp->f1 & SELECT) {
2303                                                 start_sel = a;
2304                                                 break;
2305                                         }
2306                                 }
2307                                 /* in case there are no other selected verts */
2308                                 end_sel = start_sel;
2309                                 for (bp = &nu->bp[start_sel + 1], a = start_sel + 1; a < nu->pntsu; a++, bp++) {
2310                                         if ((bp->f1 & SELECT) == 0) {
2311                                                 break;
2312                                         }
2313                                         end_sel = a;
2314                                 }
2315                                 
2316                                 if (start_sel == -1) {
2317                                         last_sel = nu->pntsu; /* next... */
2318                                 }
2319                                 else {
2320                                         last_sel = end_sel; /* before we modify it */
2321                                         
2322                                         /* now blend between start and end sel */
2323                                         start_rad = end_rad = -1.0;
2324                                         
2325                                         if (start_sel == end_sel) {
2326                                                 /* simple, only 1 point selected */
2327                                                 if (start_sel > 0) start_rad = nu->bp[start_sel - 1].radius;
2328                                                 if (end_sel != -1 && end_sel < nu->pntsu) end_rad = nu->bp[start_sel + 1].radius;
2329                                                 
2330                                                 if      (start_rad >= 0.0f && end_rad >= 0.0f)  nu->bp[start_sel].radius = (start_rad + end_rad) / 2;
2331                                                 else if (start_rad >= 0.0f)                     nu->bp[start_sel].radius = start_rad;
2332                                                 else if (end_rad >= 0.0f)                       nu->bp[start_sel].radius = end_rad;
2333                                         }
2334                                         else {
2335                                                 /* if endpoints selected, then use them */
2336                                                 if (start_sel == 0) {
2337                                                         start_rad = nu->bp[start_sel].radius;
2338                                                         start_sel++; /* we don't want to edit the selected endpoint */
2339                                                 }
2340                                                 else {
2341                                                         start_rad = nu->bp[start_sel - 1].radius;
2342                                                 }
2343                                                 if (end_sel == nu->pntsu - 1) {
2344                                                         end_rad = nu->bp[end_sel].radius;
2345                                                         end_sel--; /* we don't want to edit the selected endpoint */
2346                                                 }
2347                                                 else {
2348                                                         end_rad = nu->bp[end_sel + 1].radius;
2349                                                 }
2350                                                 
2351                                                 /* Now Blend between the points */
2352                                                 range = (float)(end_sel - start_sel) + 2.0f;
2353                                                 for (bp = &nu->bp[start_sel], a = start_sel; a <= end_sel; a++, bp++) {
2354                                                         fac = (float)(1 + a - start_sel) / range;
2355                                                         bp->radius = start_rad * (1.0f - fac) + end_rad * fac;
2356                                                 }
2357                                         }
2358                                 }
2359                         }
2360                 }
2361         }
2362
2363         WM_event_add_notifier(C, NC_GEOM | ND_DATA, obedit->data);
2364         DAG_id_tag_update(obedit->data, 0);
2365
2366         return OPERATOR_FINISHED;
2367 }
2368
2369 void CURVE_OT_smooth_radius(wmOperatorType *ot)
2370 {
2371         /* identifiers */
2372         ot->name = "Smooth Curve Radius";
2373         ot->description = "Flatten radii of selected points";
2374         ot->idname = "CURVE_OT_smooth_radius";
2375         
2376         /* api clastbacks */
2377         ot->exec = smooth_radius_exec;
2378         ot->poll = ED_operator_editsurfcurve;
2379         
2380         /* flags */
2381         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2382 }
2383
2384 /***************** selection utility *************************/
2385
2386 /* next == 1 -> select next         */
2387 /* next == -1 -> select previous    */
2388 /* cont == 1 -> select continuously */
2389 /* selstatus, inverts behavior          */
2390 static void select_adjacent_cp(ListBase *editnurb, short next, short cont, short selstatus)
2391 {
2392         Nurb *nu;
2393         BezTriple *bezt;
2394         BPoint *bp;
2395         int a;
2396         short lastsel = 0;
2397         
2398         if (next == 0) return;
2399         
2400         for (nu = editnurb->first; nu; nu = nu->next) {
2401                 lastsel = 0;
2402                 if (nu->type == CU_BEZIER) {                    
2403                         a = nu->pntsu;
2404                         bezt = nu->bezt;
2405                         if (next < 0) bezt = &nu->bezt[a - 1];
2406                         while (a--) {
2407                                 if (a - abs(next) < 0) break;
2408                                 if ((lastsel == 0) && (bezt->hide == 0) && ((bezt->f2 & SELECT) || (selstatus == 0))) {
2409                                         bezt += next;
2410                                         if (!(bezt->f2 & SELECT) || (selstatus == 0)) {
2411                                                 short sel = select_beztriple(bezt, selstatus, 1, VISIBLE);
2412                                                 if ((sel == 1) && (cont == 0)) lastsel = 1;
2413                                         }
2414                                 }
2415                                 else {
2416                                         bezt += next;
2417                                         lastsel = 0;
2418                                 }
2419                                 /* move around in zigzag way so that we go through each */                              
2420                                 bezt -= (next - next / abs(next));
2421                         }
2422                 }
2423                 else {
2424                         a = nu->pntsu * nu->pntsv;
2425                         bp = nu->bp;
2426                         if (next < 0) bp = &nu->bp[a - 1];
2427                         while (a--) {
2428                                 if (a - abs(next) < 0) break;
2429                                 if ((lastsel == 0) && (bp->hide == 0) && ((bp->f1 & SELECT) || (selstatus == 0))) {
2430                                         bp += next;
2431                                         if (!(bp->f1 & SELECT) || (selstatus == 0)) {
2432                                                 short sel = select_bpoint(bp, selstatus, 1, VISIBLE);
2433                                                 if ((sel == 1) && (cont == 0)) lastsel = 1;
2434                                         }                       
2435                                 }
2436                                 else {
2437                                         bp += next;
2438                                         lastsel = 0;
2439                                 }
2440                                 /* move around in zigzag way so that we go through each */
2441                                 bp -= (next - next / abs(next));
2442                         }
2443                 }
2444         }
2445 }
2446
2447 /**************** select start/end operators **************/
2448
2449 /* (de)selects first or last of visible part of each Nurb depending on selFirst     */
2450 /* selFirst: defines the end of which to select                                     */
2451 /* doswap: defines if selection state of each first/last control point is swapped   */
2452 /* selstatus: selection status in case doswap is false                              */
2453 void selectend_nurb(Object *obedit, short selfirst, short doswap, short selstatus)
2454 {
2455         ListBase *editnurb = object_editcurve_get(obedit);
2456         Nurb *nu;
2457         BPoint *bp;
2458         BezTriple *bezt;
2459         Curve *cu;
2460         int a;
2461
2462         if (obedit == NULL) return;
2463
2464         cu = (Curve *)obedit->data;
2465         cu->lastsel = NULL;
2466
2467         for (nu = editnurb->first; nu; nu = nu->next) {
2468                 if (nu->type == CU_BEZIER) {
2469                         a = nu->pntsu;
2470                         
2471                         /* which point? */
2472                         if (selfirst == 0) { /* select last */
2473                                 bezt = &nu->bezt[a - 1];
2474                         }
2475                         else { /* select first */
2476                                 bezt = nu->bezt;
2477                         }
2478                         
2479                         while (a--) {
2480                                 short sel;
2481                                 if (doswap) sel = swap_selection_beztriple(bezt);
2482                                 else sel = select_beztriple(bezt, selstatus, 1, VISIBLE);
2483                                 
2484                                 if (sel == 1) break;
2485                         }
2486                 }
2487                 else {
2488                         a = nu->pntsu * nu->pntsv;
2489                         
2490                         /* which point? */
2491                         if (selfirst == 0) { /* select last */
2492                                 bp = &nu->bp[a - 1];
2493                         }
2494                         else { /* select first */
2495                                 bp = nu->bp;
2496                         }
2497
2498                         while (a--) {
2499                                 if (bp->hide == 0) {
2500                                         short sel;
2501                                         if (doswap) sel = swap_selection_bpoint(bp);
2502                                         else sel = select_bpoint(bp, selstatus, 1, VISIBLE);
2503                                         
2504                                         if (sel == 1) break;
2505                                 }
2506                         }
2507                 }
2508         }
2509 }
2510
2511 static int de_select_first_exec(bContext *C, wmOperator *UNUSED(op))
2512 {
2513         Object *obedit = CTX_data_edit_object(C);
2514
2515         selectend_nurb(obedit, FIRST, 1, DESELECT);
2516         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2517
2518         return OPERATOR_FINISHED;
2519 }
2520
2521 void CURVE_OT_de_select_first(wmOperatorType *ot)
2522 {
2523         /* identifiers */
2524         ot->name = "(De)select First";
2525         ot->idname = "CURVE_OT_de_select_first";
2526         ot->description = "(De)select first of visible part of each Nurb";
2527         
2528         /* api cfirstbacks */
2529         ot->exec = de_select_first_exec;
2530         ot->poll = ED_operator_editcurve;
2531         
2532         /* flags */
2533         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2534 }
2535
2536 static int de_select_last_exec(bContext *C, wmOperator *UNUSED(op))
2537 {
2538         Object *obedit = CTX_data_edit_object(C);
2539
2540         selectend_nurb(obedit, LAST, 1, DESELECT);
2541         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2542
2543         return OPERATOR_FINISHED;
2544 }
2545
2546 void CURVE_OT_de_select_last(wmOperatorType *ot)
2547 {
2548         /* identifiers */
2549         ot->name = "(De)select Last";
2550         ot->idname = "CURVE_OT_de_select_last";
2551         ot->description = "(De)select last of visible part of each Nurb";
2552         
2553         /* api clastbacks */
2554         ot->exec = de_select_last_exec;
2555         ot->poll = ED_operator_editcurve;
2556         
2557         /* flags */
2558         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2559 }
2560
2561 /******************* de select all operator ***************/
2562
2563 static short nurb_has_selected_cps(ListBase *editnurb)
2564 {
2565         Nurb *nu;
2566         BezTriple *bezt;
2567         BPoint *bp;
2568         int a;
2569
2570         for (nu = editnurb->first; nu; nu = nu->next) {
2571                 if (nu->type == CU_BEZIER) {
2572                         a = nu->pntsu;
2573                         bezt = nu->bezt;
2574                         while (a--) {
2575                                 if (bezt->hide == 0) {
2576                                         if ((bezt->f1 & SELECT) ||
2577                                             (bezt->f2 & SELECT) ||
2578                                             (bezt->f3 & SELECT))
2579                                         {
2580                                                 return 1;
2581                                         }
2582                                 }
2583                                 bezt++;
2584                         }
2585                 }
2586                 else {
2587                         a = nu->pntsu * nu->pntsv;
2588                         bp = nu->bp;
2589                         while (a--) {
2590                                 if ((bp->hide == 0) && (bp->f1 & SELECT)) return 1;
2591                                 bp++;
2592                         }
2593                 }
2594         }
2595         
2596         return 0;
2597 }
2598
2599 static int de_select_all_exec(bContext *C, wmOperator *op)
2600 {
2601         Object *obedit = CTX_data_edit_object(C);
2602         ListBase *editnurb = object_editcurve_get(obedit);
2603         int action = RNA_enum_get(op->ptr, "action");
2604
2605         if (action == SEL_TOGGLE) {
2606                 action = SEL_SELECT;
2607                 if (nurb_has_selected_cps(editnurb))
2608                         action = SEL_DESELECT;
2609         }
2610
2611         switch (action) {
2612                 case SEL_SELECT:
2613                         CU_select_all(obedit);
2614                         break;
2615                 case SEL_DESELECT:
2616                         CU_deselect_all(obedit);
2617                         break;
2618                 case SEL_INVERT:
2619                         CU_select_swap(obedit);
2620                         break;
2621         }
2622         
2623         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2624
2625         return OPERATOR_FINISHED;
2626 }
2627
2628 void CURVE_OT_select_all(wmOperatorType *ot)
2629 {
2630         /* identifiers */
2631         ot->name = "(De)select All";
2632         ot->idname = "CURVE_OT_select_all";
2633         ot->description = "(De)select all control points";
2634         
2635         /* api callbacks */
2636         ot->exec = de_select_all_exec;
2637         ot->poll = ED_operator_editsurfcurve;
2638         
2639         /* flags */
2640         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2641
2642         /* properties */
2643         WM_operator_properties_select_all(ot);
2644 }
2645
2646 /********************** hide operator *********************/
2647
2648 static int hide_exec(bContext *C, wmOperator *op)
2649 {
2650         Object *obedit = CTX_data_edit_object(C);
2651         Curve *cu = obedit->data;
2652         ListBase *editnurb = object_editcurve_get(obedit);
2653         Nurb *nu;
2654         BPoint *bp;
2655         BezTriple *bezt;
2656         int a, sel, invert = RNA_boolean_get(op->ptr, "unselected");
2657
2658         for (nu = editnurb->first; nu; nu = nu->next) {
2659                 if (nu->type == CU_BEZIER) {
2660                         bezt = nu->bezt;
2661                         a = nu->pntsu;
2662                         sel = 0;
2663                         while (a--) {
2664                                 if (invert == 0 && BEZSELECTED_HIDDENHANDLES(cu, bezt)) {
2665                                         select_beztriple(bezt, DESELECT, 1, HIDDEN);
2666                                         bezt->hide = 1;
2667                                 }
2668                                 else if (invert && !BEZSELECTED_HIDDENHANDLES(cu, bezt)) {
2669                                         select_beztriple(bezt, DESELECT, 1, HIDDEN);
2670                                         bezt->hide = 1;
2671                                 }
2672                                 if (bezt->hide) sel++;
2673                                 bezt++;
2674                         }
2675                         if (sel == nu->pntsu) nu->hide = 1;
2676                 }
2677                 else {
2678                         bp = nu->bp;
2679                         a = nu->pntsu * nu->pntsv;
2680                         sel = 0;
2681                         while (a--) {
2682                                 if (invert == 0 && (bp->f1 & SELECT)) {
2683                                         select_bpoint(bp, DESELECT, 1, HIDDEN);
2684                                         bp->hide = 1;
2685                                 }
2686                                 else if (invert && (bp->f1 & SELECT) == 0) {
2687                                         select_bpoint(bp, DESELECT, 1, HIDDEN);
2688                                         bp->hide = 1;
2689                                 }
2690                                 if (bp->hide) sel++;
2691                                 bp++;
2692                         }
2693                         if (sel == nu->pntsu * nu->pntsv) nu->hide = 1;
2694                 }
2695         }
2696
2697         DAG_id_tag_update(obedit->data, 0);
2698         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2699
2700         return OPERATOR_FINISHED;       
2701 }
2702
2703 void CURVE_OT_hide(wmOperatorType *ot)
2704 {
2705         /* identifiers */
2706         ot->name = "Hide Selected";
2707         ot->idname = "CURVE_OT_hide";
2708         ot->description = "Hide (un)selected control points";
2709         
2710         /* api callbacks */
2711         ot->exec = hide_exec;
2712         ot->poll = ED_operator_editsurfcurve;
2713         
2714         /* flags */
2715         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2716         
2717         /* props */
2718         RNA_def_boolean(ot->srna, "unselected", 0, "Unselected", "Hide unselected rather than selected");
2719 }
2720
2721 /********************** reveal operator *********************/
2722
2723 static int reveal_exec(bContext *C, wmOperator *UNUSED(op))
2724 {
2725         Object *obedit = CTX_data_edit_object(C);
2726         ListBase *editnurb = object_editcurve_get(obedit);
2727         Nurb *nu;
2728         BPoint *bp;
2729         BezTriple *bezt;
2730         int a;
2731
2732         for (nu = editnurb->first; nu; nu = nu->next) {
2733                 nu->hide = 0;
2734                 if (nu->type == CU_BEZIER) {
2735                         bezt = nu->bezt;
2736                         a = nu->pntsu;
2737                         while (a--) {
2738                                 if (bezt->hide) {
2739                                         select_beztriple(bezt, SELECT, 1, HIDDEN);
2740                                         bezt->hide = 0;
2741                                 }
2742                                 bezt++;
2743                         }
2744                 }
2745                 else {
2746                         bp = nu->bp;
2747                         a = nu->pntsu * nu->pntsv;
2748                         while (a--) {
2749                                 if (bp->hide) {
2750                                         select_bpoint(bp, SELECT, 1, HIDDEN);
2751                                         bp->hide = 0;
2752                                 }
2753                                 bp++;
2754                         }
2755                 }
2756         }
2757
2758         DAG_id_tag_update(obedit->data, 0);
2759         WM_event_add_notifier(C, NC_GEOM | ND_SELECT, obedit->data);
2760
2761         return OPERATOR_FINISHED;       
2762 }
2763
2764 void CURVE_OT_reveal(wmOperatorType *ot)
2765 {
2766         /* identifiers */
2767         ot->name = "Reveal Hidden";
2768         ot->idname = "CURVE_OT_reveal";
2769         ot->description = "Show again hidden control points";
2770         
2771         /* api callbacks */
2772         ot->exec = reveal_exec;
2773         ot->poll = ED_operator_editsurfcurve;
2774         
2775         /* flags */
2776         ot->flag = OPTYPE_REGISTER | OPTYPE_UNDO;
2777 }
2778
2779 /********************** subdivide operator *********************/
2780
2781 /** Divide the line segments associated with the currently selected
2782  * curve nodes (Bezier or NURB). If there are no valid segment
2783  * selections within the current selection, nothing happens.
2784  */
2785 static void subdividenurb(Object *obedit, int number_cuts)
2786 {
2787         Curve *cu = obedit->data;
2788         EditNurb *editnurb = cu->editnurb;
2789         Nurb *nu;
2790         BezTriple *prevbezt, *bezt, *beztnew, *beztn;
2791         BPoint *bp, *prevbp, *bpnew, *bpn;
2792         float vec[15];
2793         int a, b, sel, amount, *usel, *vsel, i;
2794         float factor;
2795
2796         // printf("*** subdivideNurb: entering subdivide\n");
2797
2798         for (nu = editnurb->nurbs.first; nu; nu = nu->next) {
2799                 amount = 0;
2800                 if (nu->type == CU_BEZIER) {
2801                         /*
2802                          * Insert a point into a 2D Bezier curve.
2803                          * Endpoints are preserved. Otherwise, all selected and inserted points are
2804                          * newly created. Old points are discarded.
2805                          */
2806                         /* count */
2807                         if (nu->flagu & CU_NURB_CYCLIC) {
2808                                 a = nu->pntsu;
2809                                 bezt = nu->bezt;
2810                                 prevbezt = bezt + (a - 1);
2811                         }
2812                         else {
2813                                 a = nu->pntsu - 1;
2814                                 prevbezt = nu->bezt;
2815                                 bezt = prevbezt + 1;
2816                         }
2817                         while (a--) {
2818                                 if (BEZSELECTED_HIDDENHANDLES(cu, prevbezt) && BEZSELECTED_HIDDENHANDLES(cu, bezt) ) amount += number_cuts;
2819                                 prevbezt = bezt;
2820                                 bezt++;
2821                         }
2822
2823                         if (amount) {
2824                                 /* insert */
2825                                 beztnew = (BezTriple *)MEM_mallocN((amount + nu->pntsu) * sizeof(BezTriple), "subdivNurb");
2826                                 beztn = beztnew;
2827                                 if (nu->flagu & CU_NURB_CYCLIC) {
2828                                         a = nu->pntsu;
2829                                         bezt = nu->bezt;
2830                                         prevbezt = bezt + (a - 1);
2831                                 }
2832                                 else {
2833                                         a = nu->pntsu - 1;
2834                                         prevbezt = nu->bezt;
2835                                         bezt = prevbezt + 1;
2836                                 }
2837                                 while (a--) {
2838                                         memcpy(beztn, prevbezt, sizeof(BezTriple));
2839                                         keyIndex_updateBezt(editnurb, prevbezt, beztn, 1);
2840                                         beztn++;
2841
2842                                         if (BEZSELECTED_HIDDENHANDLES(cu, prevbezt) && BEZSELECTED_HIDDENHANDLES(cu, bezt) ) {
2843                                                 float prevvec[3][3];
2844
2845                                                 memcpy(prevvec, prevbezt->vec, sizeof(float) * 9);
2846
2847                                                 for (i = 0; i < number_cuts; i++) {
2848                                                         factor = 1.0f / (number_cuts + 1 - i);
2849
2850          &nbs