1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
|
# Copyright (c) 1997-2024
# Ewgenij Gawrilow, Michael Joswig, and the polymake team
# Technische Universität Berlin, Germany
# https://polymake.org
#
# This program 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 2, or (at your option) any
# later version: http://www.gnu.org/licenses/gpl.txt.
#
# This program 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.
#-------------------------------------------------------------------------------
object Polytope {
# @category Visualization
# Take a 3-polytope and write ASCII STL output.
#
# @param String filename
# @example > dodecahedron()->write_stl("/tmp/dodecahedron.stl");
user_method write_stl(String) : VERTICES, FACETS, VIF_CYCLIC_NORMAL {
my ($self,$filename) = @_;
my $n_facets=$self->FACETS->rows();
my $vertices=new Matrix<Float>($self->VERTICES->minor(All,~[0]));
open STL, ">$filename" or die "can't write to $filename";
print STL << ".";
solid Default
.
for (my $f=0; $f<$n_facets; ++$f) {
# facet normals pointing outwards
my $facet_normal=new Vector<Float>(-$self->FACETS->[$f]->slice(range_from(1)));
# each facet needs to be triangulated
my $vif=$self->VIF_CYCLIC_NORMAL->[$f];
my $k=$vif->size();
# the first vertex is contained in each triangle
my $first=$vertices->[$vif->[0]];
for (my $i=2; $i<$k; ++$i) {
my $second=$vertices->[$vif->[$i-1]];
my $third=$vertices->[$vif->[$i]];
print STL << ".";
facet normal $facet_normal
outer loop
vertex $first
vertex $second
vertex $third
endloop
endfacet
.
}
}
print STL << ".";
endsolid Default
.
close STL;
}
precondition : CONE_DIM { $this->CONE_DIM==4 }
precondition : CONE_AMBIENT_DIM { $this->CONE_AMBIENT_DIM==4 }
precondition : BOUNDED;
}
# Local Variables:
# mode: perl
# cperl-indent-level:3
# indent-tabs-mode:nil
# End:
|
