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
85
86
87
88
89
90
91
92
93
94
95
|
/* 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.
--------------------------------------------------------------------------------
*/
#pragma once
#include "polymake/client.h"
#include "polymake/IncidenceMatrix.h"
#include "polymake/linalg.h"
#include "polymake/Array.h"
#include "polymake/Matrix.h"
namespace polymake { namespace polytope {
/** Transform a coordinate section by applying a given transformation matrix.
@param p_out destination file
@param p_in source file
@param section_req read request for the section to be transformed
@param tau transformation matrix
*/
template <typename TransMatrix> inline
void transform_section(BigObject& p_out, BigObject& p_in,
const AnyString& section_req,
const GenericMatrix<TransMatrix>& tau)
{
Matrix<typename TransMatrix::element_type> M;
std::string name;
if (p_in.lookup_with_property_name(section_req, name, allow_conversion()) >> M) {
if (M.rows())
p_out.take(name) << M*tau;
else
p_out.take(name) << M;
}
}
/** Transform all coordinate sections according to the given transformation matrix.
@param p object to transform
@param tau transformation matrix for points
*/
template <typename Scalar, typename TransMatrix>
BigObject transform(BigObject p_in, const GenericMatrix<TransMatrix>& tau,
bool store_reverse_transformation = true)
{
const auto tau_inv = inv(tau);
BigObject p_out("Polytope", mlist<Scalar>());
transform_section(p_out, p_in, "VERTICES | POINTS", tau);
transform_section(p_out, p_in, "LINEALITY_SPACE | INPUT_LINEALITY", tau);
transform_section(p_out, p_in, "ZONOTOPE_INPUT_POINTS", tau);
transform_section(p_out, p_in, "FACETS | INEQUALITIES", T(tau_inv));
transform_section(p_out, p_in, "AFFINE_HULL | EQUATIONS", T(tau_inv));
/* FIXME
transform_section(p_out, p_in, "LINEAR_OBJECTIVE", T(tau_inv));
*/
IncidenceMatrix<> VIF;
if (p_in.lookup("VERTICES_IN_FACETS") >> VIF)
p_out.take("VERTICES_IN_FACETS") << VIF;
Array<std::string> labels;
if (p_in.lookup("VERTEX_LABELS") >> labels) p_out.take("VERTEX_LABELS") << labels;
if (p_in.lookup("FACET_LABELS") >> labels) p_out.take("FACET_LABELS") << labels;
if (store_reverse_transformation) {
Matrix<typename TransMatrix::element_type> tau_rev;
if (p_in.get_attachment("REVERSE_TRANSFORMATION") >> tau_rev)
tau_rev= tau_inv * tau_rev;
else
tau_rev= tau_inv;
p_out.attach("REVERSE_TRANSFORMATION") << tau_rev;
}
return p_out;
}
} } // end namespace polymake
// Local Variables:
// mode:C++
// c-basic-offset:3
// indent-tabs-mode:nil
// End:
|
