Code

- try to use more forward declarations for less dependencies on display/curve.h
[inkscape.git] / src / sp-item-transform.cpp
index 2ad7758a46262720130862ae94502e9f73cb1ae0..08c4d167b2c1672be331110a0df3e0a736618ace 100644 (file)
@@ -80,8 +80,8 @@ void sp_item_move_rel(SPItem *item, NR::translate const &tr)
 }
 
 /*
-** Returns the matrix you need to apply to an object with given bbox and strokewidth to
-scale/move it to the new box x0/y0/x1/y1. Takes into account the "scale stroke"
+** Returns the matrix you need to apply to an object with given visual bbox and strokewidth to
+scale/move it to the new visual bbox x0/y0/x1/y1. Takes into account the "scale stroke"
 preference value passed to it. Has to solve a quadratic equation to make sure
 the goal is met exactly and the stroke scaling is obeyed.
 */
@@ -97,9 +97,9 @@ get_scale_transform_with_stroke (NR::Rect &bbox_param, gdouble strokewidth, bool
     NR::Matrix scale = NR::Matrix (NR::scale (1, 1)); // scale component
     NR::Matrix unbudge = NR::Matrix (NR::translate (0, 0)); // move component to compensate for the drift caused by stroke width change
 
-    gdouble w0 = bbox.extent(NR::X);
+    gdouble w0 = bbox.extent(NR::X); // will return a value >= 0, as required further down the road
     gdouble h0 = bbox.extent(NR::Y);
-    gdouble w1 = x1 - x0;
+    gdouble w1 = x1 - x0; // can have any sign
     gdouble h1 = y1 - y0;
     gdouble r0 = strokewidth;
 
@@ -114,43 +114,42 @@ get_scale_transform_with_stroke (NR::Rect &bbox_param, gdouble strokewidth, bool
         return (p2o * direct * o2n); // can't solve the equation: one of the dimensions is equal to stroke width, so return the straightforward scaler
     }
 
-    // Flip when the width or height changes sign
-    int flip_x = ((w1 < 0) == (w0 < 0)) ? 1 : -1;
-    int flip_y = ((h1 < 0) == (h0 < 0)) ? 1 : -1;
+    int flip_x = (w1 > 0) ? 1 : -1;
+    int flip_y = (h1 > 0) ? 1 : -1;
     
-    // w1 and h1 can be negative, but if so then e.g. w1-r0 won't make sense
-    // therefore we should use fabs() all over the place
-    gdouble ratio_x = (fabs(w1) - fabs(r0)) / (fabs(w0) - fabs(r0));
-    gdouble ratio_y = (fabs(h1) - fabs(r0)) / (fabs(h0) - fabs(r0));    
+    // w1 and h1 will be negative when mirroring, but if so then e.g. w1-r0 won't make sense
+    // Therefore we will use the absolute values from this point on
+    w1 = fabs(w1);
+    h1 = fabs(h1);
+    r0 = fabs(r0);
+    // w0 and h0 will always be positive due to the definition extent()
+
+    gdouble ratio_x = (w1 - r0) / (w0 - r0);
+    gdouble ratio_y = (h1 - r0) / (h0 - r0);
     
-    NR::Matrix direct_constant_r = NR::Matrix (NR::scale(flip_x * ratio_x, flip_y*ratio_y));
+    NR::Matrix direct_constant_r = NR::Matrix (NR::scale(flip_x * ratio_x, flip_y * ratio_y));
 
     if (transform_stroke && r0 != 0 && r0 != NR_HUGE) { // there's stroke, and we need to scale it
         // These coefficients are obtained from the assumption that scaling applies to the
         // non-stroked "shape proper" and that stroke scale is scaled by the expansion of that
-        // matrix
-        // In fact, we're trying to solve this equation:
-        // r1 = r0 * sqrt (((w1-r0)/(w0-r0))*((h1-r1)/(h0-r0)))
-        // To make sense of this, the operant of the sqrt() should
-        // be positive, hence all the fabs() below
-        // (w1 and h1 will be negative when mirroring, w0 and h0 will probably never be negative)
-        gdouble A = -fabs(w0*h0) + fabs(r0)*(fabs(w0) + fabs(h0));
-        gdouble B = -(fabs(w1) + fabs(h1)) * r0*r0;
-        gdouble C = fabs(w1 * h1 * r0*r0);
+        // matrix. We're trying to solve this equation:
+        // r1 = r0 * sqrt (((w1-r0)/(w0-r0))*((h1-r0)/(h0-r0)))
+        // The operant of the sqrt() must be positive, which is ensured by the fabs() a few lines above
+        gdouble A = -w0*h0 + r0*(w0 + h0);
+        gdouble B = -(w1 + h1) * r0*r0;
+        gdouble C = w1 * h1 * r0*r0;
         if (B*B - 4*A*C > 0) {
-            gdouble r1 = (-B - sqrt (B*B - 4*A*C))/(2*A);
+            gdouble r1 = fabs((-B - sqrt(B*B - 4*A*C))/(2*A));
             //gdouble r2 = (-B + sqrt (B*B - 4*A*C))/(2*A);
             //std::cout << "r0" << r0 << " r1" << r1 << " r2" << r2 << "\n";
             //
-            // I think r1 will always be positive if r0 is (mathematical proof?)
-            // but if w1 becomes negative, then the scale will be wrong if we just do  
+            // If w1 < 0 then the scale will be wrong if we just do
             // gdouble scale_x = (w1 - r1)/(w0 - r0);
-            // gdouble scale_y = (h1 - r1)/(h0 - r0);
-            // So let's do it like this: Calculate the absolute scale
-            gdouble scale_x = (fabs(w1) - fabs(r1))/(fabs(w0) - fabs(r0));
-            gdouble scale_y = (fabs(h1) - fabs(r1))/(fabs(h0) - fabs(r0));
-            scale *= NR::scale(flip_x*scale_x, flip_y*scale_y);
-            unbudge *= NR::translate (-flip_x * 0.5 * (fabs(r0) * scale_x - fabs(r1)), -flip_y * 0.5 * (fabs(r0) * scale_y - fabs(r1)));
+            // Here we also need the absolute values of w0, w1, h0, h1, and r1
+            gdouble scale_x = (w1 - r1)/(w0 - r0);
+            gdouble scale_y = (h1 - r1)/(h0 - r0);
+            scale *= NR::scale(flip_x * scale_x, flip_y * scale_y);
+            unbudge *= NR::translate (-flip_x * 0.5 * (r0 * scale_x - r1), -flip_y * 0.5 * (r0 * scale_y - r1));
         } else {
             scale *= direct;
         }
@@ -159,13 +158,45 @@ get_scale_transform_with_stroke (NR::Rect &bbox_param, gdouble strokewidth, bool
             scale *= direct;
         } else {// nonscaling strokewidth
             scale *= direct_constant_r;
-            unbudge *= NR::translate (flip_x * 0.5 * fabs(r0) * (1 - ratio_x), flip_y * 0.5 * fabs(r0) * (1 - ratio_y));
+            unbudge *= NR::translate (flip_x * 0.5 * r0 * (1 - ratio_x), flip_y * 0.5 * r0 * (1 - ratio_y));
         }
     }
 
     return (p2o * scale * unbudge * o2n);
 }
 
+NR::Rect
+get_visual_bbox (NR::Maybe<NR::Rect> const &initial_geom_bbox, NR::Matrix const &abs_affine, gdouble const initial_strokewidth, bool const transform_stroke)
+{
+    
+    g_assert(initial_geom_bbox);
+    
+    // Find the new geometric bounding box; Do this by transforming each corner of
+    // the initial geometric bounding box individually and fitting a new boundingbox
+    // around the transformerd corners  
+    NR::Point const p0 = initial_geom_bbox->corner(0) * abs_affine;    
+    NR::Rect new_geom_bbox = NR::Rect(p0, p0);
+    for (unsigned i = 1 ; i < 4 ; i++) {
+        new_geom_bbox.expandTo(initial_geom_bbox->corner(i) * abs_affine);
+    }
+
+    NR::Rect new_visual_bbox = new_geom_bbox; 
+    if (initial_strokewidth > 0 && initial_strokewidth < NR_HUGE) {
+        if (transform_stroke) {
+            // scale stroke by: sqrt (((w1-r0)/(w0-r0))*((h1-r0)/(h0-r0))) (for visual bboxes, see get_scale_transform_with_stroke)
+            // equals scaling by: sqrt ((w1/w0)*(h1/h0)) for geometrical bboxes            
+            // equals scaling by: sqrt (area1/area0) for geometrical bboxes
+            gdouble const new_strokewidth = initial_strokewidth * sqrt (new_geom_bbox.area() / initial_geom_bbox->area());
+            new_visual_bbox.growBy(0.5 * new_strokewidth);        
+        } else {
+            // Do not transform the stroke
+            new_visual_bbox.growBy(0.5 * initial_strokewidth);   
+        }
+    }
+    
+    return new_visual_bbox;
+}
+
 /*
   Local Variables:
   mode:c++