Code

optimize cairo paths by not drawing segments that are outside of view (for stroke...
[inkscape.git] / src / display / nr-filter.cpp
index 026aa4ee00ece0243af2ee8286b6ddedfa0fbee2..07c2649289b128fa9b9657ab3cd171ce4b390251 100644 (file)
  */
 
 #include <glib.h>
+#include <cmath>
 
 #include "display/nr-filter.h"
 #include "display/nr-filter-primitive.h"
 #include "display/nr-filter-gaussian.h"
 #include "display/nr-filter-slot.h"
 #include "display/nr-filter-types.h"
+#include "display/pixblock-scaler.h"
+#include "display/pixblock-transform.h"
 
 #include "display/nr-arena-item.h"
 #include "libnr/nr-pixblock.h"
 #include "libnr/nr-blit.h"
+#include "libnr/nr-matrix.h"
+#include "libnr/nr-scale.h"
 #include "svg/svg-length.h"
 #include "sp-filter-units.h"
+#if defined (SOLARIS_2_8)
+#include "round.h"
+using Inkscape::round;
+#endif 
 
 //#include "display/nr-arena-shape.h"
 
+__attribute__ ((const))
+inline static int _max4(const double a, const double b,
+                        const double c, const double d) {
+    double ret = a;
+    if (b > ret) ret = b;
+    if (c > ret) ret = c;
+    if (d > ret) ret = d;
+    return (int)round(ret);
+}
+
+__attribute__ ((const))
+inline static int _min4(const double a, const double b,
+                        const double c, const double d) {
+    double ret = a;
+    if (b < ret) ret = b;
+    if (c < ret) ret = c;
+    if (d < ret) ret = d;
+    return (int)round(ret);
+}
+
 namespace NR {
 
 Filter::Filter()
 {
-    _primitive_count = 1;
+    _primitive_count = 0;
     _primitive_table_size = 1;
     _primitive = new FilterPrimitive*[1];
-    _primitive[0] = new FilterGaussian;
+       _primitive[0] = NULL;
+    //_primitive_count = 1;
+    //_primitive[0] = new FilterGaussian;
     _common_init();
 }
 
@@ -51,13 +82,19 @@ Filter::Filter(int n)
 
 void Filter::_common_init() {
     _slot_count = 1;
+    // Having "not set" here as value means the output of last filter
+    // primitive will be used as output of this filter
     _output_slot = NR_FILTER_SLOT_NOT_SET;
 
+    // These are the default values for filter region,
+    // as specified in SVG standard
+    // NB: SVGLength.set takes prescaled percent values: -.10 means -10%
     _region_x.set(SVGLength::PERCENT, -.10, 0);
     _region_y.set(SVGLength::PERCENT, -.10, 0);
     _region_width.set(SVGLength::PERCENT, 1.20, 0);
     _region_height.set(SVGLength::PERCENT, 1.20, 0);
 
+    // Filter resolution, negative value here stands for "automatic"
     _x_pixels = -1.0;
     _y_pixels = -1.0;
 
@@ -74,32 +111,141 @@ Filter::~Filter()
 
 int Filter::render(NRArenaItem const *item, NRPixBlock *pb)
 {
-    FilterSlot slot(_slot_count);
+    if(!_primitive[0]) { // if there are no primitives, do nothing
+       return 0; 
+    }
+
+    Matrix trans = *item->ctm;
+    Matrix paraller_trans = trans;
+    bool notparaller = false;
+    FilterSlot slot(_slot_count, item);
     NRPixBlock *in = new NRPixBlock;
-    nr_pixblock_setup_fast(in, pb->mode,
-                           pb->area.x0, pb->area.y0,
-                           pb->area.x1, pb->area.y1, true);
-    nr_blit_pixblock_pixblock(in, pb);
+
+    // If filter effects region is not paraller to viewport,
+    // we must first undo the rotation / shear.
+    // It will be redone after filtering.
+    // If there is only rotation and uniform scaling (zoom), let's skip this,
+    // as it will not make a difference with gaussian blur.
+    if ((fabs(trans[1]) > 1e-6 || fabs(trans[2]) > 1e-6) &&
+        !(fabs(trans[0] - trans[3]) < 1e-6 && fabs(trans[1] + trans[2]) < 1e-6)) {
+        notparaller = true;
+
+        // TODO: if filter resolution is specified, scaling should be set
+        // according to that
+        double scaling_factor = sqrt(trans.expansionX() * trans.expansionX() +
+                                     trans.expansionY() * trans.expansionY());
+        scale scaling(scaling_factor, scaling_factor);
+        scale scaling_inv(1.0 / scaling_factor, 1.0 / scaling_factor);
+        trans *= scaling_inv;
+        paraller_trans.set_identity();
+        paraller_trans *= scaling;
+
+        Matrix itrans = trans.inverse();
+        int x0 = pb->area.x0;
+        int y0 = pb->area.y0;
+        int x1 = pb->area.x1;
+        int y1 = pb->area.y1;
+        int min_x = _min4(itrans[0] * x0 + itrans[2] * y0 + itrans[4],
+                          itrans[0] * x0 + itrans[2] * y1 + itrans[4],
+                          itrans[0] * x1 + itrans[2] * y0 + itrans[4],
+                          itrans[0] * x1 + itrans[2] * y1 + itrans[4]);
+        int max_x = _max4(itrans[0] * x0 + itrans[2] * y0 + itrans[4],
+                          itrans[0] * x0 + itrans[2] * y1 + itrans[4],
+                          itrans[0] * x1 + itrans[2] * y0 + itrans[4],
+                          itrans[0] * x1 + itrans[2] * y1 + itrans[4]);
+        int min_y = _min4(itrans[1] * x0 + itrans[3] * y0 + itrans[5],
+                          itrans[1] * x0 + itrans[3] * y1 + itrans[5],
+                          itrans[1] * x1 + itrans[3] * y0 + itrans[5],
+                          itrans[1] * x1 + itrans[3] * y1 + itrans[5]);
+        int max_y = _max4(itrans[1] * x0 + itrans[3] * y0 + itrans[5],
+                          itrans[1] * x0 + itrans[3] * y1 + itrans[5],
+                          itrans[1] * x1 + itrans[3] * y0 + itrans[5],
+                          itrans[1] * x1 + itrans[3] * y1 + itrans[5]);
+        
+        nr_pixblock_setup_fast(in, pb->mode,
+                               min_x, min_y,
+                               max_x, max_y, true);
+        if (in->size != NR_PIXBLOCK_SIZE_TINY && in->data.px == NULL) // memory allocation failed
+            return 0;
+        transform_nearest(in, pb, itrans);
+    } else if (_x_pixels >= 0) {
+        // If filter resolution is not set to automatic, we should
+        // scale the input image to correct resolution
+        /* If filter resolution is zero, the object should not be rendered */
+        if (_x_pixels == 0 || _y_pixels == 0) {
+            int size = (pb->area.x1 - pb->area.x0)
+                * (pb->area.y1 - pb->area.y0)
+                * NR_PIXBLOCK_BPP(pb);
+            memset(NR_PIXBLOCK_PX(pb), 0, size);
+            return 0;
+        }
+        // Resolution is specified as pixel length of our internal buffer.
+        // Though, we might not be rendering the whole object at time,
+        // so we need to calculate the correct pixel size
+        int x_len = (int)round(((pb->area.x1 - pb->area.x0) * _x_pixels) / (item->bbox.x1 - item->bbox.x0));
+        if (x_len < 1) x_len = 1;
+        // If y-resolution is also set, count y-area in the same way as x-area
+        // Otherwise, make y-area so, that aspect ratio of input pixblock and
+        // internal pixblock are the same.
+        int y_len;
+        if (_y_pixels > 0) {
+            y_len = (int)round(((pb->area.y1 - pb->area.y0) * _y_pixels) / (item->bbox.y1 - item->bbox.y0));
+        } else {
+            y_len = (int)round((x_len * (pb->area.y1 - pb->area.y0)) / (double)(pb->area.x1 - pb->area.x0));
+        }
+        if (y_len < 1) y_len = 1;
+        nr_pixblock_setup_fast(in, pb->mode, 0, 0, x_len, y_len, true);
+        if (in->size != NR_PIXBLOCK_SIZE_TINY && in->data.px == NULL) // memory allocation failed
+            return 0;
+        scale_bicubic(in, pb);
+        scale res_scaling(x_len / (double)(pb->area.x1 - pb->area.x0),
+                          y_len / (double)(pb->area.y1 - pb->area.y0));
+        paraller_trans *= res_scaling;
+    } else {
+        // If filter resolution is automatic, just make copy of input image
+        nr_pixblock_setup_fast(in, pb->mode,
+                               pb->area.x0, pb->area.y0,
+                               pb->area.x1, pb->area.y1, true);
+        if (in->size != NR_PIXBLOCK_SIZE_TINY && in->data.px == NULL) // memory allocation failed
+            return 0;
+        nr_blit_pixblock_pixblock(in, pb);
+    }
     slot.set(NR_FILTER_SOURCEGRAPHIC, in);
     in = NULL; // in is now handled by FilterSlot, we should not touch it
 
-    _primitive[0]->render(slot, *item->ctm);
-
+    // TODO: loop through ALL the primitives and render them one at a time
+    _primitive[0]->render(slot, paraller_trans);
     NRPixBlock *out = slot.get(_output_slot);
 
+    // Clear the pixblock, where the output will be put
+    // -> the original image does not show through
     int size = (pb->area.x1 - pb->area.x0)
         * (pb->area.y1 - pb->area.y0)
         * NR_PIXBLOCK_BPP(pb);
     memset(NR_PIXBLOCK_PX(pb), 0, size);
 
-    nr_blit_pixblock_pixblock(pb, out);
+    if (notparaller) {
+        transform_nearest(pb, out, trans);
+    } else if (_x_pixels < 0) {
+        // If the filter resolution is automatic, just copy our final image
+        // to output pixblock, otherwise use bicubic scaling
+        nr_blit_pixblock_pixblock(pb, out);
+    } else {
+        scale_bicubic(pb, out);
+    }
 
+    // Take note of the amount of used image slots
+    // -> next time this filter is rendered, we can reserve enough slots
+    // immediately
     _slot_count = slot.get_slot_count();
     return 0;
 }
 
 int Filter::get_enlarge(Matrix const &m)
 {
+    // Just sum the enlargement factor of all filter elements.
+    // TODO: this both sucks and blows for filters like feOffset
+    // -> ditch this method and design a better one...
     int enlarge = 0;
     for ( int i = 0 ; i < _primitive_count ; i++ ) {
         if(_primitive[i]) enlarge += _primitive[i]->get_enlarge(m);
@@ -149,11 +295,17 @@ void Filter::bbox_enlarge(NRRectL &bbox)
     }
 }
 
+/* Constructor table holds pointers to static methods returning filter
+ * primitives. This table is indexed with FilterPrimitiveType, so that
+ * for example method in _constructor[NR_FILTER_GAUSSIANBLUR]
+ * returns a filter object of type NR::FilterGaussian.
+ */
 typedef FilterPrimitive*(*FilterConstructor)();
 static FilterConstructor _constructor[NR_FILTER_ENDPRIMITIVETYPE];
 
 void Filter::_create_constructor_table()
 {
+    // Constructor table won't change in run-time, so no need to recreate
     static bool created = false;
     if(created) return;
 
@@ -174,8 +326,13 @@ void Filter::_create_constructor_table()
     _constructor[NR_FILTER_SPECULARLIGHTING] = NULL;
     _constructor[NR_FILTER_TILE] = NULL;
     _constructor[NR_FILTER_TURBULENCE] = NULL;
+    created = true;
 }
 
+/** Helper method for enlarging table of filter primitives. When new
+ * primitives are added, but we have no space for them, this function
+ * makes some more space.
+ */
 void Filter::_enlarge_primitive_table() {
     FilterPrimitive **new_tbl = new FilterPrimitive*[_primitive_table_size * 2];
     for (int i = 0 ; i < _primitive_count ; i++) {
@@ -205,6 +362,7 @@ FilterPrimitive *Filter::add_primitive(FilterPrimitiveType type)
     }
 
     _primitive[_primitive_count] = created;
+    _primitive_count++;
     return created;
 }
 
@@ -246,6 +404,27 @@ void Filter::clear_primitives()
     _primitive_count = 0;
 }
 
+void Filter::set_x(SVGLength &length)
+{ 
+  if (length._set)
+      _region_x = length;
+}
+void Filter::set_y(SVGLength &length)
+{
+  if (length._set)
+      _region_y = length;
+}
+void Filter::set_width(SVGLength &length)
+{
+  if (length._set)
+      _region_width = length;
+}
+void Filter::set_height(SVGLength &length)
+{ 
+  if (length._set)
+      _region_height = length;
+}
+
 } /* namespace NR */
 
 /*