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
73
74
75
76
77
78
79
80
81
82
83
84
|
# 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.
#-------------------------------------------------------------------------------
CREDIT vinci
Vinci - Computing Volumes of Convex Polytopes.
Copyright by Benno Büeler, Andreas Enge, and Komei Fukuda
http://www.lix.polytechnique.fr/Labo/Andreas.Enge/Vinci.html
# path to the "vinci" executable
custom $vinci;
CONFIGURE {
find_program($vinci, "vinci");
}
object Polytope<Float> {
# @category Triangulation and volume
# Use the [[wiki:external_software#vinci]] program for volume computation.
label vinci
rule vinci.volume: VOLUME : CONE_AMBIENT_DIM, VERTICES, N_VERTICES, FACETS, N_FACETS {
my $dim = $this->CONE_AMBIENT_DIM;
# write the ext-file
my $tempname=new Tempfile;
open my $P, ">$tempname.ext"
or die "can't create temporary file $tempname.ext: $!";
print $P "Automatically generated by polymake.\n",
"begin\n", $this->N_VERTICES, " $dim real\n",
dense($this->VERTICES),
"end\n";
close $P;
# write the ine-file
open $P, ">$tempname.ine"
or die "can't create temporary file $tempname.ine: $!";
print $P "Automatically generated by polymake.\n",
"begin\n", $this->N_FACETS, " $dim real\n",
dense($this->FACETS),
"end\n";
close $P;
# call vinci
open $P, "$vinci $tempname |"
or die "couldn't run $vinci: $!";
local $_;
while (<$P>) {
if (/^Volume:\s+/) {
$this->VOLUME=0+$';
close $P;
return;
}
}
die "can't parse output from $vinci";
}
precondition : N_FACETS, CONE_AMBIENT_DIM {
$this->CONE_AMBIENT_DIM <= 11 && $this->N_FACETS <= 200
}
precondition : CONE_DIM, CONE_AMBIENT_DIM {
$this->CONE_DIM==$this->CONE_AMBIENT_DIM;
}
precondition : BOUNDED;
}
# Local Variables:
# mode: perl
# cperl-indent-level:3
# End:
|
