Construction of a vector frame from a family of vector fields

[egourgoulhon]

This ticket modifies DifferentiableManifold.vector_frame() to allow for constructing a vector frame from a spanning family of linearly independent vector fields:

sage: M = Manifold(2, 'M')
sage: X.<x,y> = M.chart()
sage: e0 = M.vector_field(1+x^2, 1+y^2)
sage: e1 = M.vector_field(2, -x*y)
sage: e = M.vector_frame('e', (e0, e1)); e
Vector frame (M, (e_0,e_1))
sage: e[0].display()
e_0 = (x^2 + 1) d/dx + (y^2 + 1) d/dy
sage: e[1].display()
e_1 = 2 d/dx - x*y d/dy
sage: (e[0], e[1]) == (e0, e1)
True

Previously, the only way to introduce the vector frame e was to first introduce the automorphism relating the frame (d/dx, d/dy) to (e0, e1) and to pass this automorphism to VectorFrame.new_frame():

sage: aut = M.automorphism_field()
sage: aut[:] = [[e0[0], e1[0]], [e0[1], e1[1]]]
sage: e = X.frame().new_frame(aut, 'e')

Implementation details: such functionality already existed for bases of finite rank free modules; the relevant code is extracted from the method FiniteRankFreeModule.basis() and put into the new method FreeModuleBasis._init_from_family(), in order to be used in DifferentiableManifold.vector_frame() as well.

CC: @tscrim

Component: geometry

Keywords: vector_frame

Author: Eric Gourgoulhon

Branch/Commit: 81e2f60

Reviewer: Michael Jung

Issue created by migration from https://trac.sagemath.org/ticket/28716

[egourgoulhon]

New commits:

013fb8b

Add construction of a vector frame from a family of vector fields

[egourgoulhon]

Commit: 013fb8b

[egourgoulhon]

Branch: public/manifolds/vector_frame_from_family-28716

[egourgoulhon]

Description changed:

--- 
+++ 
@@ -14,7 +14,7 @@
 sage: (e[0], e[1]) == (e0, e1)
 True
 ```
-Previously, the only way to introduce the vector frame `e` was to first introduce an automorphism relating the frame `(d/dx, d/dy)` to `(e0, e1)` and to pass this automorphism to `VectorFrame.new_frame()`:
+Previously, the only way to introduce the vector frame `e` was to first introduce the automorphism relating the frame `(d/dx, d/dy)` to `(e0, e1)` and to pass this automorphism to `VectorFrame.new_frame()`:
 
 ```
 sage: aut = M.automorphism_field()

[tscrim]

comment:4

Do we need the optional parameter? Basically, can we just use the fact that a tuple/list is being given and then assume it is suppose to be a family of vector fields? If it has to be a keyword, I would change from_family to the more descriptive from_vector_fields.

[DeRhamSource]

comment:5

This is a great idea and would be very useful for vector bundles, too. Sometimes I got really annoyed by this detour. Would you mind to adapt your code, if working, for vector bundles as-well?

By the way: We should combine vector bundles and the previous implementations really really soon (in this case inherit vector frames from local frames) otherwise things could get extremly messy.

Unfortunately, I am quite busy working at my master thesis right now. I can almost feel the deadline touching my skin. I promise to work on that as soon as I've gained some time back.

Even though I don't have the time now, I've opened the corresponding ticket #28718, just to keep this task in mind.

[egourgoulhon]

comment:6

Replying to @tscrim:

Thanks for your prompt feedback.

Do we need the optional parameter? Basically, can we just use the fact that a tuple/list is being given and then assume it is suppose to be a family of vector fields?

Good idea, this is much more user-friendly! I am on it...

[egourgoulhon]

comment:7

Replying to @DeRhamSource:

This is a great idea and would be very useful for vector bundles, too. Sometimes I got really annoyed by this detour.

Yes, this should have been done sooner...

Would you mind to adapt your code, if working, for vector bundles as-well?

OK, I'll try to do this (see below).

By the way: We should combine vector bundles and the previous implementations really really soon (in this case inherit vector frames from local frames) otherwise things could get extremly messy.

Yes, I agree. Note however that this ticket does not touch the class VectorFrame, only the user interface DifferentiableManifold.vector_frame(). I'll perform a similar change to the interfaces TopologicalVectorBundle.local_frame() and TensorBundle.local_frame().

Unfortunately, I am quite busy working at my master thesis right now. I can almost feel the deadline touching my skin.

Good luck with your master thesis!

I promise to work on that as soon as I've gained some time back.

Even though I don't have the time now, I've opened the corresponding ticket #28718, just to keep this task in mind.

Thanks.

[sagetrac-git]

Changed commit from 013fb8b to d0ef4d7