Refactored the gyroid infill.

xs/src/libslic3r/Fill/FillGyroid.cpp

@@ -9,149 +9,74 @@
 
 namespace Slic3r {
 
-Polyline FillGyroid::makeLineVert(double xPos, double yPos, double width, double height, double currentXBegin, double segmentSize, coord_t scaleFactor, 
+static inline Polyline make_wave_vertical(
-        double zCs, double zSn, bool flip, double decal){
+    double width, double height, double x0,
-    double maxSlope = abs(abs(zCs)-abs(zSn));
+    double segmentSize, double scaleFactor,
-    Polyline polyline;
+    double z_cos, double z_sin, bool flip)
-    polyline.points.push_back(Point(coord_t((std::max(std::min(currentXBegin, xPos+width),xPos) + decal) * scaleFactor), coord_t(yPos * scaleFactor)));
-    for(double y=yPos;y<yPos+height+segmentSize;y+=segmentSize){
-        if(y>yPos+height) y = yPos+height;
-        double ySn = sin(y +(zCs<0?3.14:0) + 3.14);
-        double yCs = cos(y +(zCs<0?3.14:0) + 3.14+(!flip?0:3.14));
-
-        double a = ySn;
-        double b = -zCs;
-        double res = zSn*yCs;
-        double r = sqrt(a*a + b*b);
-        double x = asin(a/r) + asin(res/r) +3.14;
-        x += currentXBegin;
-        
-        double ydeviation = 0.5*(flip?-1:1)*(zSn>0?-1:1)*decal*(1-maxSlope)*(res/r - a/r);
-        polyline.points.push_back(Point(coord_t((std::max(std::min(x, xPos+width),xPos)+decal-ydeviation/2) * scaleFactor), coord_t((y + ydeviation) * scaleFactor)));
-    }
-    
-    return polyline;
-}
-
-Polyline FillGyroid::makeLineHori(double xPos, double yPos, double width, double height, double currentYBegin, double segmentSize, coord_t scaleFactor, 
-        double zCs, double zSn, bool flip, double decal){
-    double maxSlope = abs(abs(zCs)-abs(zSn));
-    Polyline polyline;
-    polyline.points.push_back(Point(coord_t(xPos * scaleFactor), coord_t((std::max(std::min(currentYBegin, yPos+height),yPos)+decal) * scaleFactor)));
-    for(double x=xPos;x<xPos+width+segmentSize;x+=segmentSize){
-        if(x>xPos+width) x = xPos+width;
-        double xSn = sin(x +(zSn<0?3.14:0) +(flip?0:3.14));
-        double xCs = cos(x +(zSn<0?3.14:0) );
-        
-        double a = xCs;
-        double b = -zSn;
-        double res = zCs*xSn;
-        double r = sqrt(a*a + b*b);
-        double y = asin(a/r) + asin(res/r) +3.14/2;
-        y += currentYBegin;
-        
-        double xdeviation = 0.5*(flip?-1:1)*(zCs>0?-1:1)*decal*(1-maxSlope)*(res/r - a/r);
-        polyline.points.push_back(Point(coord_t((x + xdeviation) * scaleFactor), coord_t((std::max(std::min(y, yPos+height),yPos)+decal-xdeviation/2) * scaleFactor)));
-    }
-    
-    return polyline;
-}
-
-inline void FillGyroid::correctOrderAndAdd(const int num, Polyline poly, Polylines &array){
-    if(num%2==0){
-        Points temp(poly.points.rbegin(), poly.points.rend());
-        poly.points.assign(temp.begin(),temp.end());
-    }
-    array.push_back(poly);
-}
-
-// Generate a set of curves (array of array of 2d points) that describe a
-// horizontal slice of a truncated regular octahedron with a specified
-// grid square size.
-Polylines FillGyroid::makeGrid(coord_t gridZ, double density, double layer_width, size_t gridWidth, size_t gridHeight, size_t curveType)
 {
-    coord_t  scaleFactor = coord_t(scale_(layer_width) / density);
+    Polyline polyline;
-    Polylines result;
+    polyline.points.emplace_back(Point(coord_t(clamp(0., width, x0) * scaleFactor), 0));
-    Polyline *polyline2;
+    double phase_offset_sin = (z_cos < 0 ? M_PI : 0) + M_PI;
-    double segmentSize = density/2;
+    double phase_offset_cos = (z_cos < 0 ? M_PI : 0) + M_PI + (flip ? M_PI : 0.);
-    double decal = layer_width*density;
+    for (double y = 0.; y < height + segmentSize; y += segmentSize) {
-    double xPos = 0, yPos=0, width=gridWidth, height=gridHeight;
+        y = std::min(y, height);
-     //scale factor for 5% : 8 712 388
+        double a   = sin(y + phase_offset_sin);
-     // 1z = 10^-6 mm ?
+        double b   = - z_cos;
-    double z = gridZ/(1.0 * scaleFactor);
+        double res = z_sin * cos(y + phase_offset_cos);
-    double zSn = sin(z);
+        double r   = sqrt(sqr(a) + sqr(b));
-    double zCs = cos(z);
+        double x   = clamp(0., width, asin(a/r) + asin(res/r) + M_PI + x0);
-    
+        polyline.points.emplace_back(convert_to<Point>(Pointf(x, y) * scaleFactor));
-
-    int numLine = 0;
-    
-    if(abs(zSn)<=abs(zCs)){
-        //vertical
-        //begin to first one
-        int iter = 1;
-        double currentXBegin = xPos - PI/2;
-        currentXBegin = PI*(int)(currentXBegin/PI -1);
-        iter = (int)(currentXBegin/PI +1)%2;
-        bool flip = iter%2==1;
-        // bool needNewLine =false;
-        while(currentXBegin<xPos+width-PI/2){
-            
-            correctOrderAndAdd(numLine, makeLineVert(xPos, yPos, width, height, currentXBegin, segmentSize, scaleFactor, zCs, zSn, flip, 0), result);
-            numLine++;
-            
-            //then, return by the other side
-            iter++;
-            currentXBegin = currentXBegin + PI;
-            flip = iter%2==1;
-            
-            if(currentXBegin < xPos+width-PI/2){
-                
-                correctOrderAndAdd(numLine, makeLineVert(xPos, yPos, width, height, currentXBegin, segmentSize, scaleFactor, zCs, zSn, flip, 0), result);
-                numLine++;
-
-                // relance
-                iter++;
-                currentXBegin = currentXBegin + PI;
-                flip = iter%2==1;
-            }
-        }
-    }else{
-        //horizontal
-        
-
-        //begin to first one
-        int iter = 1;
-        //search first line output
-        double currentYBegin = yPos ;
-        currentYBegin = PI*(int)(currentYBegin/PI -0);
-        iter = (int)(currentYBegin/PI +1)%2;
-        
-        bool flip = iter%2==1;
-        
-        
-        while(currentYBegin < yPos+width){
-
-            correctOrderAndAdd(numLine, makeLineHori(xPos, yPos, width, height, currentYBegin, segmentSize, scaleFactor, zCs, zSn, flip, 0), result);
-            numLine++;
-        
-            //then, return by the other side
-            iter++;
-            currentYBegin = currentYBegin + PI;
-            flip = iter%2==1;
-            
-            if(currentYBegin<yPos+width){
-                
-                correctOrderAndAdd(numLine, makeLineHori(xPos, yPos, width, height, currentYBegin, segmentSize, scaleFactor, zCs, zSn, flip, 0), result);
-                numLine++;
-                
-                //relance
-                iter++;
-                currentYBegin = currentYBegin + PI;
-                flip = iter%2==1;
-            }
-        }
     }
+    if (flip)
+        std::reverse(polyline.points.begin(), polyline.points.end());
+    return polyline;
+}
 
+static inline Polyline make_wave_horizontal(
+    double width, double height, double y0, 
+    double segmentSize, double scaleFactor,
+    double z_cos, double z_sin, bool flip)
+{
+    Polyline polyline;
+    polyline.points.emplace_back(Point(0, coord_t(clamp(0., height, y0) * scaleFactor)));
+    double phase_offset_sin = (z_sin < 0 ? M_PI : 0) + (flip ? 0 : M_PI);
+    double phase_offset_cos = z_sin < 0 ? M_PI : 0.;
+    for (double x=0.; x < width + segmentSize; x += segmentSize) {
+        x = std::min(x, width);
+        double a   = cos(x + phase_offset_cos);
+        double b   = - z_sin;
+        double res = z_cos * sin(x + phase_offset_sin);
+        double r   = sqrt(sqr(a) + sqr(b));
+        double y   = clamp(0., height, asin(a/r) + asin(res/r) + 0.5 * M_PI + y0);
+        polyline.points.emplace_back(convert_to<Point>(Pointf(x, y) * scaleFactor));
+    }
+    if (flip)
+        std::reverse(polyline.points.begin(), polyline.points.end());
+    return polyline;
+}
+
+static Polylines make_gyroid_waves(double gridZ, double density, double layer_width, double width, double height)
+{
+    double scaleFactor = scale_(layer_width) / density;
+    double segmentSize = 0.5 * density;
+ //scale factor for 5% : 8 712 388
+ // 1z = 10^-6 mm ?
+    double z     = gridZ / scaleFactor;
+    double z_sin = sin(z);
+    double z_cos = cos(z);
+    Polylines result;
+    if (abs(z_sin) <= abs(z_cos)) {
+        // Vertical wave
+        double x0 = M_PI * (int)((- 0.5 * M_PI) / M_PI - 1.);
+        bool   flip          = ((int)(x0 / M_PI + 1.) & 1) != 0;
+        for (; x0 < width - 0.5 * M_PI; x0 += M_PI, flip = ! flip)
+            result.emplace_back(make_wave_vertical(width, height, x0, segmentSize, scaleFactor, z_cos, z_sin, flip));
+    } else {
+        // Horizontal wave
+        bool flip = true;
+        for (double y0 = 0.; y0 < width; y0 += M_PI, flip = !flip)
+            result.emplace_back(make_wave_horizontal(width, height, y0, segmentSize, scaleFactor, z_cos, z_sin, flip));
+    }
     return result;
 }
 
@@ -170,18 +95,16 @@ void FillGyroid::_fill_surface_single(
     bb.merge(_align_to_grid(bb.min, Point(2*M_PI*distance, 2*M_PI*distance)));
     
     // generate pattern
-    Polylines   polylines = makeGrid(
+    Polylines   polylines = make_gyroid_waves(
-        (coord_t)scale_(this->z),
+        scale_(this->z),
         params.density*this->scaling,
         this->spacing,
-        (size_t)(ceil(bb.size().x / distance) + 1),
+        ceil(bb.size().x / distance) + 1.,
-        (size_t)(ceil(bb.size().y / distance) + 1),
+        ceil(bb.size().y / distance) + 1.);
-        (size_t)(((this->layer_id/thickness_layers) % 2) + 1) );
     
     // move pattern in place
-    for (Polylines::iterator it = polylines.begin(); it != polylines.end(); ++ it)
+    for (Polyline &polyline : polylines)
-        it->translate(bb.min.x, bb.min.y);
+        polyline.translate(bb.min.x, bb.min.y);
-    
 
     // clip pattern to boundaries
     polylines = intersection_pl(polylines, (Polygons)expolygon);
@@ -198,18 +121,14 @@ void FillGyroid::_fill_surface_single(
             }
         }
         Polylines chained = PolylineCollection::chained_path_from(
-#if SLIC3R_CPPVER >= 11
             std::move(polylines), 
-#else
-            polylines,
-#endif
             PolylineCollection::leftmost_point(polylines), false); // reverse allowed
         bool first = true;
-        for (Polylines::iterator it_polyline = chained.begin(); it_polyline != chained.end(); ++ it_polyline) {
+        for (Polyline &polyline : chained) {
             if (! first) {
                 // Try to connect the lines.
                 Points &pts_end = polylines_out.back().points;
-                const Point &first_point = it_polyline->points.front();
+                const Point &first_point = polyline.points.front();
                 const Point &last_point = pts_end.back();
                 // TODO: we should also check that both points are on a fill_boundary to avoid 
                 // connecting paths on the boundaries of internal regions
@@ -217,17 +136,12 @@ void FillGyroid::_fill_surface_single(
                 if (first_point.distance_to(last_point) <= 5 * distance && 
                     expolygon_off.contains(Line(last_point, first_point))) {
                     // Append the polyline.
-                    pts_end.insert(pts_end.end(), it_polyline->points.begin(), it_polyline->points.end());
+                    pts_end.insert(pts_end.end(), polyline.points.begin(), polyline.points.end());
                     continue;
                 }
             }
             // The lines cannot be connected.
-#if SLIC3R_CPPVER >= 11
+            polylines_out.emplace_back(std::move(polyline));
-            polylines_out.push_back(std::move(*it_polyline));
-#else
-            polylines_out.push_back(Polyline());
-            std::swap(polylines_out.back(), *it_polyline);
-#endif
             first = false;
         }
     }

xs/src/libslic3r/Fill/FillGyroid.hpp

@@ -1,8 +1,6 @@
 #ifndef slic3r_FillGyroid_hpp_
 #define slic3r_FillGyroid_hpp_
 
-#include <map>
-
 #include "../libslic3r.h"
 
 #include "FillBase.hpp"
@@ -12,10 +10,8 @@ namespace Slic3r {
 class FillGyroid : public Fill
 {
 public:
-
+    FillGyroid() {}
-    FillGyroid(){ scaling = 1.75; }
+    virtual Fill* clone() const { return new FillGyroid(*this); }
-    virtual Fill* clone() const { return new FillGyroid(*this); };
-    virtual ~FillGyroid() {}
 
     // require bridge flow since most of this pattern hangs in air
     virtual bool use_bridge_flow() const { return true; }
@@ -23,8 +19,7 @@ public:
 protected:
     
     // mult of density, to have a good %of weight for each density parameter
-    float scaling;
+    float scaling = 1.75;
-
 
     virtual void _fill_surface_single(
         const FillParams                &params, 
@@ -32,22 +27,6 @@ protected:
         const std::pair<float, Point>   &direction, 
         ExPolygon                       &expolygon, 
         Polylines                       &polylines_out);
-    
-    // create the gyroid grid to clip.
-    Polylines makeGrid(coord_t gridZ, double density, double layer_width, size_t gridWidth, size_t gridHeight, size_t curveType);
-    //add line poly in reverse if needed into array
-    inline void correctOrderAndAdd(const int num, Polyline poly, Polylines &array);
-    //create a curved horinzontal line  (for each x, compute y)
-    Polyline makeLineHori(double xPos, double yPos, double width, double height, 
-        double currentYBegin, double segmentSize, coord_t scaleFactor, 
-        double zCs, double zSn, 
-        bool flip, double decal=0);
-    //create a curved vertival line (for each y, compute x)
-    Polyline makeLineVert(double xPos, double yPos, double width, double height, 
-        double currentXBegin, double segmentSize, coord_t scaleFactor, 
-        double zCs, double zSn, 
-        bool flip, double decal=0);
-
 };
 
 } // namespace Slic3r