Run benchmark script on 6 spiral approaches

main
mmcwilliams 3 years ago
parent 98467a4f6b
commit 5f592661ba
  1. 1
      .gitignore
  2. 3
      README.md
  3. 54
      scripts/benchmark.sh
  4. 67
      spiral/spiral_1.scad
  5. 47
      spiral/spiral_2.scad
  6. 29
      spiral/spiral_3.scad
  7. 35
      spiral/spiral_4.scad
  8. 24
      spiral/spiral_5.scad
  9. 81
      spiral/spiral_6.scad

1
.gitignore vendored

@ -1 +1,2 @@
*.DS_Store
benchmark

@ -83,8 +83,9 @@ PLA is not recommended but this doesn't mean you can't get an acceptable results
### Dependencies
* [OpenSCAD](https://www.openscad.org/downloads.html) + [cli](https://en.wikibooks.org/wiki/OpenSCAD_User_Manual/Using_OpenSCAD_in_a_command_line_environment) ([Mac](https://en.wikibooks.org/wiki/OpenSCAD_User_Manual/Using_OpenSCAD_in_a_command_line_environment#MacOS_notes)) ([Windows](https://en.wikibooks.org/wiki/OpenSCAD_User_Manual/Using_OpenSCAD_in_a_command_line_environment#Windows_notes))
* Bash
* [OpenSCAD](https://www.openscad.org/downloads.html) + [cli](https://en.wikibooks.org/wiki/OpenSCAD_User_Manual/Using_OpenSCAD_in_a_command_line_environment) ([Mac](https://en.wikibooks.org/wiki/OpenSCAD_User_Manual/Using_OpenSCAD_in_a_command_line_environment#MacOS_notes)) ([Windows](https://en.wikibooks.org/wiki/OpenSCAD_User_Manual/Using_OpenSCAD_in_a_command_line_environment#Windows_notes))
* [ADMesh](https://github.com/admesh/admesh)
### Build Scripts

@ -0,0 +1,54 @@
echo "Benchmarking GNAL spiral generation"
VERSION=`bash ./scripts/version.sh`
CPU=`bash ./scripts/cpu.sh`
DATE=`date -u +%s` #timestamp in seconds
NOTES="./notes/benchmark.csv"
ROTATIONS=( 1 2 5 10 20 )
DIAMETERS=( 47 225 300 )
COMPLETE=( 40 60 )
echo "version,cpu,date,source,diameter,rotations,$fn,size,time" > $NOTES
for SPIRAL in ./spiral/*.scad
do
FILENAME=`basename ${SPIRAL}`
echo "Rendering ${SPIRAL}..."
for ROT in "${ROTATIONS[@]}"
do
:
for D in "${DIAMETERS[@]}"
do
:
echo "Rendering ${SPIRAL} for ${ROT} rotations @ ${D}mm"
TMP=`mktemp`
start=`date +%s`
openscad -o "${TMP}.stl" -D N=${ROT} -D D=${D} -D FN=100 "${SPIRAL}"
end=`date +%s`
runtime=$((end-start))
size=`wc -c < "${TMP}.stl"`
size=`echo $size | xargs`
line="${VERSION},${CPU},${DATE},${FILENAME},${D},${ROT},100,$size,$runtime"
echo $line
echo $line >> $NOTES
rm "${TMP}.stl"
done
done
for C in "${COMPLETE[@]}"
do
:
echo "Rendering complete ${SPIRAL} with ${C} rotations"
TMP=`mktemp`
start=`date +%s`
echo -o "${TMP}.stl" -D N=${C} -D D=47 "${SPIRAL}"
end=`date +%s`
runtime=$((end-start))
echo "${TMP}.stl"
#rm "${TMP}.stl"
done
done

@ -0,0 +1,67 @@
D=300; // start diameter
N=1; // number of spirals
FN=100;
$fn=FN;
H = 2.1;
/**
* Render all spiral facets for as many rotations as supplied
*/
module spiral (START_D = 50, SPIRALS = 39) {
//STOP_D = 100;
SPACING = 0.86;//1.34;
TOP_T = 0.86; //thickness
BOTTOM_T = 1.4;
START_R = START_D / 2;
union () {
for (i = [0 : $fn]) {
rotate ([0, 0, i * (360 / $fn)]) {
for (x = [0: (SPIRALS - 1)]) {
spiral_facet(i, x, START_R, SPACING, BOTTOM_T, TOP_T, H);
}
}
}
}
}
/**
* Generates a single face of the spiral, in this case a trapazoidal
* shape. Issues are (1) performance (maybe use of hull()) and (2) all
* facet lenths are the same, despite the diameter. This means that
* there are excessive numbers of facets for the smaller spirals to
* compensate for the number of facets needed for the outer spiral.
*/
module spiral_facet (i, x, START_R, SPACING, BOTTOM_T, TOP_T, H) {
STEP_SIZE = ((SPACING + BOTTOM_T) / $fn);
STEP_OFFSET = i * STEP_SIZE;
SPIRAL_START_OFFSET = (x * (SPACING + BOTTOM_T));
ACTUAL_R = START_R + SPIRAL_START_OFFSET + STEP_OFFSET;
L = 2 * (ACTUAL_R * tan((360 / $fn) / 2));
ANGLE = -atan( STEP_SIZE / (L / 2) ) / 2;
OFFSET = START_R - (BOTTOM_T / 2) + SPIRAL_START_OFFSET + STEP_OFFSET;
translate ([OFFSET, 0, - H / 2]) {
rotate ([0, 0, ANGLE]) {
//replace hull for quick render?
//test spiral lib
hull () {
translate([0, 0, H])
cube([TOP_T, L, 0.1], center=true);
cube([BOTTOM_T, L, 0.1], center=true);
}
}
}
}
spiral(D, N);

@ -0,0 +1,47 @@
D=300; // start diameter
N=1; // number of spirals
FN=100;
$fn=FN;
module spiral_facet (i, SPIRAL, START_R = 48, SPACING = 2, FACET_SIZE = 1, FN = 360) {
BOTTOM_T = 1.2;
TOP_T = .3;
H = 2.1;
STEP_SIZE = SPACING / FN;
STEP_OFFSET = (SPIRAL * SPACING) + (i * STEP_SIZE);
ROT = i* (360 / FN);
ANGLE = 0;
OFFSET = START_R + STEP_OFFSET;
rotate([0, 0, ROT]) translate ([OFFSET, 0, - H / 2]) {
rotate ([0, 0, ANGLE]) {
hull () {
translate([(BOTTOM_T - TOP_T) / 2, 0, H])
cube([TOP_T, FACET_SIZE, 0.1], center=true);
cube([BOTTOM_T, FACET_SIZE, 0.1], center=true);
}
//cube([BOTTOM_T, FACET_SIZE, H], center = true);
}
}
}
//https://www.youtube.com/watch?v=lsSeRxpoIaY
//https://www.mathportal.org/calculators/plane-geometry-calculators/right-triangle-calculator.php
module spiral (SPIRALS = 40, START_D = 48) {
SPACING = 2.075;
FACET_SIZE = 1;
for (SPIRAL = [0 : SPIRALS - 1]) {
//C = PI * R^2
C = PI * pow(( (START_D / 2) + (SPIRAL * SPACING) ) / 2, 2);
FN = ceil( (C / FACET_SIZE) * .37 );
for (i = [0 : FN - 1]) {
spiral_facet(i, SPIRAL, START_D / 2, SPACING, FACET_SIZE, FN);
}
}
}
spiral(N, D);

@ -0,0 +1,29 @@
D=300; // start diameter
N=1; // number of spirals
FN=100;
$fn=FN;
include <../libraries/path_extrude.scad>;
module spiral (count = 40, start_d = 48, spacing = 2.075) {
//$fn = 80;
bottom = 1.2;
top = .3;
top_offset = (bottom - top);
h = 2.2;
angle_i = 360 / $fn;
increment = spacing / $fn;
facetProfile = [[0, 0], [top_offset, -h], [bottom, -h], [bottom, 0]];
spiralPath = [ for(t = [0 : $fn * count])
[((start_d / 2) + (t * increment)) * cos(t * angle_i), ((start_d / 2) + (t * increment)) * sin(t * angle_i), 0] ];
path_extrude(exShape=facetProfile, exPath=spiralPath);
}
spiral(N, D);

@ -0,0 +1,35 @@
D=300; // start diameter
N=1; // number of spirals
FN=100;
$fn=FN;
include <../libraries/path_extrude.scad>;
//Distinction from v3: use of for loop and union to join rotations
module spiral (count = 40, start_d = 48, spacing = 2.075) {
facet_size = 30;
bottom = 1.2;
top = .4;
top_offset = (bottom - top);
h = 2.2;
od = start_d + (spacing * 2 * count);
echo("SPIRAL LENGTH", PI * count * (od + start_d + 1) / 2);
echo("OUTER D", od);
facetProfile = [[0, 0], [top_offset, -h], [bottom, -h], [bottom, 0]];
union () {
for (s = [0 : count - 1]) {
d = start_d + (s * spacing * 2);
c = PI * pow(d / 2, 2);
$fn = ceil( c / facet_size );
angle_i = 360 / $fn;
increment = spacing / $fn;
spiralPath = [ for(t = [0 : $fn + 1]) [((d / 2) + (t * increment)) * cos(t * angle_i), ((d / 2) + (t * increment)) * sin(t * angle_i), 0] ];
path_extrude(exShape=facetProfile, exPath=spiralPath);
}
}
}
spiral(N, D);

@ -0,0 +1,24 @@
D=300; // start diameter
N=1; // number of spirals
FN=100;
$fn=FN;
include <../libraries/path_extrude.scad>;
//Distinction from v4 compressed spiralPath generation into single line
module spiral (count = 40, start_d = 48, spacing = 2.095) {
facet_size = 10;
bottom = 1.2;
top = .3;
top_offset = (bottom - top);
h = 2.2;
facetProfile = [[0, 0], [top_offset, -h], [bottom, -h], [bottom, 0]];
spiralPath = [ for (s = [0 : count - 1]) for(t = [0 : ceil( (PI * pow((start_d + (s * spacing * 2)) / 2, 2)) / facet_size ) - 1]) [(((start_d + (s * spacing * 2)) / 2) + (t * (spacing / ceil( (PI * pow((start_d + (s * spacing * 2)) / 2, 2)) / facet_size )))) * cos(t * (360 / ceil( (PI * pow((start_d + (s * spacing * 2)) / 2, 2)) / facet_size ))), (((start_d + (s * spacing * 2)) / 2) + (t * (spacing / ceil( (PI * pow((start_d + (s * spacing * 2)) / 2, 2)) / facet_size )))) * sin(t * (360 / ceil( (PI * pow((start_d + (s * spacing * 2)) / 2, 2)) / facet_size ))), 0] ];
path_extrude(exShape=facetProfile, exPath=spiralPath);
}
spiral(N, D);

@ -0,0 +1,81 @@
D=300; // start diameter
N=1; // number of spirals
FN=100;
$fn=FN;
R = (D + (N * 2.095)) / 2;
module ShapeToExtrude ()
{
bottom = 1.2;
top = .3;
top_offset = (bottom - top);
h = 2.2;
// Build in +x space. The outside edge of this shape must follow the extrusion path, or there will be open seams..
polygon ( points= [
[0, 0],
[top_offset, h],
[bottom, h],
[bottom, 0]
]);
}
module InwardSpiral (StepSize, Steps, StartRadius, Pitch, ShapeX) {
for (i=[0 : Steps - 1]) {
// This could be made more computationally efficient
// by collapsing intermediate values and by doing only
// essential calculations inside the loop, but for now
// let's just leave it easy to read.
ThisTheta = StepSize * i;
NextTheta = StepSize * (i + 1);
ThisRadius = StartRadius - i * (Pitch * (StepSize / 360));
// Spiral step approximated by arc of radius ThisRadius,
// passing through the start and end points calculated here.
NextRadius = StartRadius - (i + 1) * (Pitch * (StepSize / 360));
ThisX = ThisRadius * sin(ThisTheta);
ThisY = ThisRadius * cos(ThisTheta);
NextX = NextRadius * sin(NextTheta);
NextY = NextRadius * cos(NextTheta);
DeltaX = NextX - ThisX;
DeltaY = NextY - ThisY;
SlopeToNext = DeltaY / DeltaX;
BisecSlope = -1 / SlopeToNext;
ThisXYToBisector = sqrt(DeltaX * DeltaX + DeltaY * DeltaY) / 2;
BisectX = ThisX + (DeltaX / 2);
BisectY = ThisY + (DeltaY / 2);
BisectToCentre = sqrt(pow(ThisRadius, 2) - pow(ThisXYToBisector, 2));
AbsXComponent = sqrt(pow(BisectToCentre, 2) / ( 1 + pow(BisecSlope, 2)));
XComponent = NextY < ThisY ? AbsXComponent: -AbsXComponent;
YComponent = XComponent * BisecSlope;
CentreX = BisectX - XComponent;
CentreY = BisectY - YComponent;
ExtrudeAngle = -2 * acos(BisectToCentre / ThisRadius);
ArcOrientation = NextY < ThisY ? atan(BisecSlope) - (ExtrudeAngle / 2) : -180 + atan(BisecSlope) -(ExtrudeAngle / 2);
translate([CentreX, CentreY, 0]) {
rotate ([0, 0, ArcOrientation]) {
rotate_extrude (angle=ExtrudeAngle, $fn=300) translate ([ThisRadius - ShapeX, 0, 0])ShapeToExtrude();
}
}
}
}
module spiral () {
SPIRALROTATION=N * 360;
SPIRALSTEPS=180;
INITIALRADIUS= 300 / 2 ;
PITCH=2.095;
XMAXSHAPE=2;
InwardSpiral (
SPIRALROTATION/SPIRALSTEPS,
SPIRALSTEPS,
INITIALRADIUS,
PITCH,
XMAXSHAPE);
}
spiral();
Loading…
Cancel
Save