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

Commit

Permalink
Trac #30804: improved orientation example
Browse files Browse the repository at this point in the history
  • Loading branch information
mjungmath committed Nov 3, 2020
1 parent dfb6045 commit caa9d5f
Show file tree
Hide file tree
Showing 3 changed files with 20 additions and 4 deletions.
2 changes: 0 additions & 2 deletions src/doc/en/reference/manifolds/diff_manifold.rst
Original file line number Diff line number Diff line change
Expand Up @@ -10,8 +10,6 @@ Differentiable Manifolds

sage/manifolds/differentiable/examples/real_line

sage/manifolds/differentiable/examples/sphere

diff_scalarfield

diff_map
Expand Down
2 changes: 2 additions & 0 deletions src/doc/en/reference/manifolds/euclidean_space.rst
Original file line number Diff line number Diff line change
Expand Up @@ -8,4 +8,6 @@ Euclidean Spaces and Vector Calculus

sage/manifolds/differentiable/examples/euclidean

sage/manifolds/differentiable/examples/sphere

sage/manifolds/operators
20 changes: 18 additions & 2 deletions src/sage/manifolds/differentiable/examples/sphere.py
Original file line number Diff line number Diff line change
Expand Up @@ -131,15 +131,31 @@
on S^1(c)-{NP}: y1 |--> (x, y) = (2*y1/(y1^2 + 1) + 1, (y1^2 - 1)/(y1^2 + 1))
on S^1(c)-{SP}: yp1 |--> (x, y) = (2*yp1/(yp1^2 + 1) + 1, -(yp1^2 - 1)/(yp1^2 + 1))
Each sphere is orientable::
Since the sphere is a hypersurface, we can get a normal vector field by using
``normal``::
sage: n = S1c.normal(); n
Vector field n along the 1-sphere S^1(c) of radius 1 smoothly embedded in
the Euclidean plane E^2 centered at the Point c with values on the
Euclidean plane E^2
sage: n.display()
n = -cos(chi) e_x - sin(chi) e_y
Notice that this is just *one* normal field with arbitrary direction,
in this particular case `n` points inwards whereas `-n` points outwards.
However, the vector field `n` is indeed non-vanishing and hence the sphere
admits an orientation (as all spheres do)::
sage: orient = S1c.orientation(); orient
[Coordinate frame (S^1(c)-{SP}, (d/dyp1)), Vector frame (S^1(c)-{NP}, (f_1))]
sage: f = orient[1]
sage: f[1].display()
f_1 = -d/dy1
And hence admits a Volume form::
Notice that the orientation is chosen is such a way that `(\iota_*(f_1), -n)`
is oriented in the ambient Euclidean space, i.e. the last entry is the normal
vector field pointing outwards. Henceforth, the manifold admits
a volume form::
sage: h = S1c.induced_metric()
sage: h.display()
Expand Down

0 comments on commit caa9d5f

Please sign in to comment.