use <dotSCAD/src/loft.scad>;
include <BOSL2/std.scad> // for offset
stages =200;
stage_height = 1.25;
rad = 50;
f1 = 25;
f2 = 25;
phase1 = 0;
phase2 = 180;
height_depth=5;
depth1 = 20;
depth2 = 20;
thickness = 2;
bottom_thickness = 3;
myslices = 5;
angle_step=0.1;
// generate outer points
points_base = [for (i = [0:1:stages]) let(f = pow(sin(i/stages * 120),2) * 7 + 1, var = 1 , a=((sin((i/stages*360*f)%360) * 0.5 + 0.5) * (var * height_depth))) [for(j = [0:angle_step:360-angle_step]) [sin(j) * (rad+a+(pauw(sin(j *f1+phase1),0.5)*0.5+0.5)*depth1*i/stages+(pauw(sin(j *f2+phase2),0.5)*0.5+0.5)*depth2*(1-i/stages)), cos(j) *(rad+a+(pauw(sin(j *f1+phase1),0.5)*0.5+0.5)*depth1*i/stages+(pauw(sin(j *f2+phase2),0.5)*0.5+0.5)*depth2*(1-i/stages))]]];
// generate inner points with offset from outer points for hollowing out
points_base_e = [for (i = [0:1:stages]) offset(points_base[i],r=-thickness, closed=true, quality=1)];
// points - outside
points = [for (i = [0:1:stages-1]) [for (j = [0:1:(len(points_base[i])-1)]) concat(points_base[i][j][0],points_base[i][j][1],i * stage_height)]];
// points - inside
points_e = [for (i = [0:1:stages]) [for (j = [0:1:(len(points_base_e[i])-1)]) concat(points_base_e[i][j][0],points_base_e[i][j][1],i * stage_height)]];
//echo("points: ", points);
// separate points in packs for stepped extrusion
steps = [for (i = [0:1:stages-1]) points[i]];
steps_e = [for (i = [0:1:stages]) points_e[i]];
// extrusion with dotSCAD's loft function
difference() {
loft(steps,slices = myslices);
difference() {
loft(steps_e,slices = myslices);
cylinder(h = bottom_thickness, r = rad \* 2 + 20, center=true);
}
}
//======================== FUNCTIONS
function elbow(i,a,b,c) = lookup(i, [[0,a],[0.5,b],[1,c]]);
function elbow2(i,a,b,c,d) = lookup(i, [[0,a],[0.33333,b],[0.66666,c],[1,d]]);
function pauw(x,p) = sign(x) * pow(abs(x),p);
