use <experimental/foliage_scroll.scad>
use <polyline_join.scad>
use <stereographic_extrude.scad>
use <sweep.scad>
use <util/choose.scad>
use <util/reverse.scad>
// use <surface/sf_hull.scad>

$fn = 12;
width = 150;
height = 150;
max_spirals = 4; 
angle_step = 360 / $fn; 
min_radius = 15; 
init_radius = rands(min_radius * 1.75, min_radius * 2.25, 1)[0];

theta_step = 0.1;
rf_step = 0.1;

climbing_rose();

module climbing_rose() {
    thickness = 0.02;

    theta_from = PI * 1.75;
    theta_to = PI * 15;

    rf_to = 1;
    module rose(thickness, theta_from, theta_to, rf_to, rf_step) {
        function phi(theta) =
            (PI / 2) * exp(-theta / (8 * PI));
            
        function wave(theta) = 
             1 - 0.5 * pow(1.25 * pow(1- ((3.6 * theta) % (2 * PI)) / PI - (theta > 0 ? 0 : 2), 2) - 0.25, 2);

        function g(rf, theta) =
            let(angle = phi(theta) * 57.2958)
            1.95653 * pow(rf, 2) * pow(1.27689 * rf - 1, 2) * sin(angle);

        function r(rf,theta) = 
            let(angle = phi(theta) * 180 / PI)
            wave(theta) * (rf * sin(angle) + g(rf, theta) * cos(angle));

        sf = 1 + thickness;
        sections = [
            for(theta = [theta_from:theta_step:theta_to])
            let(
                path = [ 
                    for(rf = [0.1:rf_step:rf_to])
                    let(
                        r = r(rf,theta),
                        angle = theta * 57.2958, 
                        angle2 = phi(theta)* 57.2958,
                        x = r*sin(angle), 
                        y = r*cos(angle),
                        z = wave(theta)*(rf*cos(angle2)-g(rf,theta)*sin(angle2))
                    )
                    [x, y, z]
                ]
            )
            concat(
                [
                    for(p = path)
                    [p[0] * sf, p[1] * sf, p[2]] - [0, 0, thickness]
                ],
                reverse(path)
            )
        ];
        
        sweep(reverse(sections));

        /*
            // very slow, but a thicker rose is allowed
            sf_pts = [
                for(theta = [theta_from:theta_step:theta_to]) [
                    for(rf = [0.1:rf_step:rf_to])
                    let(
                        r = r(rf,theta),
                        angle = theta * 57.2958, 
                        angle2 = phi(theta)* 57.2958,
                        x = r*sin(angle), 
                        y = r*cos(angle),
                        z = wave(theta)*(rf*cos(angle2)-g(rf,theta)*sin(angle2))
                    )
                    [x, y, z]
                ]
            ];
            
            sf_hull(sf_pts, thickness, $fn = 3);

        */
    }

    module draw(spirals) {    
        r = spirals[0][0] / 5;
        for(i = [0:len(spirals) - 1]) {
            rr = r * ((len(spirals) - pow(i == 0 || i == 1 ? .5 : i, 0.25)) / len(spirals));
            path = spirals[i][1];
            
            polyline_join(path)
            rotate([45, 0, 0])
                sphere(rr, $fn = 4);
            
            center = spirals[i][2];
            leng = len(path);
            pathes = [
                [for(i = [leng / 8:leng / 6]) path[i]],
                [for(i = [leng / 4:leng / 2]) path[i]]
            ];
            for(i = [0:1]) {
                p = choose(pathes[i]);
                
                v = (p - center);
                a = atan2(v.y, v.x);
                
                translate([each p, 0] + [0, 0, rr * 0.55])
                rotate([0, rands(0, 1, 1)[0] > 0.5 ? 17 : -17, a])
                rotate(a)
                scale([rr * 6, rr * 6, rr * 5]) 
                union() {
                    rose(thickness, theta_from * rands(2.75, 3.25, 1)[0] * r / rr, rands(0.75, 1.25, 1)[0] * theta_to * r / rr, rf_to, rf_step);
                    translate([0, 0, .125])
                        sphere(.1);
                }
            }
        }
    }

    spirals = foliage_scroll([width, height], max_spirals, init_radius, min_radius);

    draw(spirals);
}