NopSCADlib/tests/maths.scad

103 lines
2.5 KiB
OpenSCAD

//
// NopSCADlib Copyright Chris Palmer 2018
// nop.head@gmail.com
// hydraraptor.blogspot.com
//
// This file is part of NopSCADlib.
//
// NopSCADlib is free software: you can redistribute it and/or modify it under the terms of the
// GNU General Public License as published by the Free Software Foundation, either version 3 of
// the License, or (at your option) any later version.
//
// NopSCADlib is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
// without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
// See the GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License along with NopSCADlib.
// If not, see <https://www.gnu.org/licenses/>.
//
include <../utils/core/core.scad>
use <../utils/maths.scad>
use <../utils/annotation.scad>
tip = [0, 0, 20];
module shape() {
cylinder(d1 = 20, d2 = 0, h = tip.z);
}
t = [10, 20, 30];
r = [45, 20, 70];
s = [1, 0.5, 0.75];
module maths() {
//
// Translate, rotate and scale the shape
//
translate(t)
rotate(r)
scale(s)
shape();
//
// Apply the same transformations to the vector position of the tip
//
p = transform(tip, translate(t) * rotate(r) * scale(s));
//
// Place an arrow where the tip ends up
//
translate(p)
arrow();
//
// Unit vector pointing at p
//
u = unit(p);
//
// Point arrow in same direction
//
z = [0, 0, 1];
v = cross(u, z);
a = angle_between(u, z);
l = 20;
rotate(-a, v)
translate_z(l)
vflip()
arrow(l);
//
// Test Euler
//
assert(euler(rotate(r)) == r, "euler() failed");
//
// Circle intersect
//
r1 = 10;
c1 = [50, 0, 10];
r2 = 20;
c2 = [67, 0, 0];
p1 = circle_intersect(c1, r1, c2, r2);
p2 = circle_intersect(c2, r2, c1, r1);
rotate(90) {
color(grey(90))
translate(c1) rotate([90, 0, 0]) cylinder(r = r1, h = 4 * eps, center = true);
color(grey(80))
translate(c2) rotate([90, 0, 0]) cylinder(r = r2, h = eps, center = true);
color("red")
translate(p1) rotate([90, 0, 0]) cylinder(r = 0.1, h = 6 * eps, center = true);
color("blue")
translate(p2) rotate([90, 0, 0]) cylinder(r = 0.1, h = 6 * eps, center = true);
translate(p1) arrow();
translate(p2) vflip() arrow();
}
}
rotate(45)
maths();