Skip to content
This repository has been archived by the owner on Jan 30, 2023. It is now read-only.

Commit

Permalink
Use new format of affine hull data
Browse files Browse the repository at this point in the history
  • Loading branch information
Matthias Koeppe committed Apr 19, 2021
1 parent 6e686a4 commit 643392c
Showing 1 changed file with 4 additions and 6 deletions.
10 changes: 4 additions & 6 deletions src/sage/geometry/polyhedron/base.py
Original file line number Diff line number Diff line change
Expand Up @@ -10576,19 +10576,17 @@ def affine_hull_manifold(self, name=None, latex_name=None, start_index=0, ambien
CH = H.chart(names=names)

data = self.affine_hull_projection(return_all_data=True, **kwds)
projection_linear_map, projection_translation_vector = data['projection_map']
projection_matrix = projection_linear_map.matrix().transpose()
section_linear_map, section_translation_vector = data['section_map']
section_matrix = section_linear_map.matrix().transpose()
projection_matrix = data.projection_linear_map.matrix().transpose()
projection_translation_vector = data.projection_translation
section_matrix = data.section_linear_map.matrix().transpose()
section_translation_vector = data.section_translation

from sage.symbolic.ring import SR
# We use the slacks of the (linear independent) equations as the foliation parameters
foliation_parameters = vector(SR.var(f't{i}') for i in range(self.ambient_dim() - self.dim()))
normal_matrix = matrix(equation.A() for equation in self.equation_generator()).transpose()
slack_matrix = normal_matrix.pseudoinverse()

projection_linear_map.matrix().transpose().right_kernel_matrix().transpose()

phi = H.diff_map(ambient_space, {(CH, CE):
(section_matrix * vector(CH._xx) + section_translation_vector
+ normal_matrix * foliation_parameters).list()})
Expand Down

0 comments on commit 643392c

Please sign in to comment.