|
DifferentialGeometry[DGbasis] - calculate the dual basis to a given basis of vectors or 1-forms
|
|
Calling Sequence
|
|
DGbasis(S, option)
|
|
|
Parameters
|
|
|
S
|
-
|
a list of vectors, forms or tensors
|
|
option
|
-
|
the keyword argument method = "real"
|
|
|
|
|
|
Description
|
|
|
•
|
Let S = [S_1, S_2, ..., S_k] be a list of vectors, matrices, differential forms, or tensors. Then DGbasis(S) returns a sublist B = [S_i_1, S_i_2, ..., S_i_r] of S such that the elements of B define a basis for the subspace spanned by the elements of S. Thus the elements of B are linearly independent and span(S) = span(B).
|
|
•
|
With the keyword argument method = "real", the set S is viewed as a vector space of vectors, forms or tensors over the real numbers and the basis is calculated using real number coefficients instead of general expression (function) coefficients.
|
|
•
|
This command is part of the DifferentialGeometry package, and so can be used in the form DGbasis(...) only after executing the command with(DifferentialGeometry). It can always be used in the long form DifferentialGeometry:-DGbasis.
|
|
|
|
Examples
|
|
|
>
|
|
Initialize a 4-dimensional manifold M with coordinates [x, y, z, w].
|
>
|
![DGsetup([x, y, z, w], M)](/support/helpjp/helpview.aspx?si=5573/file01426/math89.png)
|
Example 1.
Find a basis B1 for the span of the vectors in S1.
|
>
|
![S1 := evalDG([D_x, D_x+D_y, D_y, `&mult`(0, D_x), D_y-D_z, D_z, D_x+D_y+D_w])](/support/helpjp/helpview.aspx?si=5573/file01426/math103.png)
|
|
![[_DG([["vector", M, []], [[[1], 1]]]), _DG([["vector", M, []], [[[1], 1], [[2], 1]]]), _DG([["vector", M, []], [[[2], 1]]]), _DG([["vector", M, []], [[[1], 0]]]), _DG([["vector", M, []], [[[2], 1], [[3], -1]]]), _DG([["vector", M, []], [[[3], 1]]]), _DG([["vector", M, []], [[[1], 1], [[2], 1], [[4], 1]]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math106.png)
| (1) |
|
>
|
|
|
![[_DG([["vector", M, []], [[[1], 1]]]), _DG([["vector", M, []], [[[1], 1], [[2], 1]]]), _DG([["vector", M, []], [[[2], 1], [[3], -1]]]), _DG([["vector", M, []], [[[1], 1], [[2], 1], [[4], 1]]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math113.png)
| (2) |
Use the GetComponents command to check that each vector in S is a unique linear combination of the vectors in B. (GetComponents returns the coefficients of the vectors in S as linear combination of the vectors in B. If one of the vectors in S is not a linear combination of the vectors in B, an empty list [] is returned).
|
>
|
|
|
![[[1, 0, 0, 0], [0, 1, 0, 0], [-1, 1, 0, 0], [0, 0, 0, 0], [0, 0, 1, 0], [-1, 1, -1, 0], [0, 0, 0, 1]]](/support/helpjp/helpview.aspx?si=5573/file01426/math142.png)
| (3) |
The basis B1 is not the simplest basis for S1. Another choice of basis for S1 can be found using the command CanonicalBasis.
|
>
|
![Tools:-CanonicalBasis(B1, [D_x, D_y, D_z, D_w])](/support/helpjp/helpview.aspx?si=5573/file01426/math157.png)
|
|
![[_DG([["vector", M, []], [[[1], 1]]]), _DG([["vector", M, []], [[[2], 1]]]), _DG([["vector", M, []], [[[3], 1]]]), _DG([["vector", M, []], [[[4], 1]]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math160.png)
| (4) |
Example 2.
Find a basis for the span of the tensors in S2.
|
>
|
![S2 := evalDG([`&t`(dx, dx), `&t`(dy, dx), `&t`(dx, dx)+`&t`(dy, dx), `&t`(dy, dy)])](/support/helpjp/helpview.aspx?si=5573/file01426/math172.png)
|
|
![[_DG([["tensor", M, [["cov_bas", "cov_bas"], []]], [[[1, 1], 1]]]), _DG([["tensor", M, [["cov_bas", "cov_bas"], []]], [[[2, 1], 1]]]), _DG([["tensor", M, [["cov_bas", "cov_bas"], []]], [[[1, 1], 1], [[2, 1], 1]]]), _DG([["tensor", M, [["cov_bas", "cov_bas"], []]], [[[2, 2], 1]]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math175.png)
| (5) |
|
>
|
|
|
![[_DG([["tensor", M, [["cov_bas", "cov_bas"], []]], [[[1, 1], 1]]]), _DG([["tensor", M, [["cov_bas", "cov_bas"], []]], [[[2, 1], 1]]]), _DG([["tensor", M, [["cov_bas", "cov_bas"], []]], [[[2, 2], 1]]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math182.png)
| (6) |
Example 3.
DGbasis also accepts a list of Matrices or Vectors.
|
>
|
![S3 := [Matrix([[0, 1], [0, 0]]), Matrix([[1, 0], [0, 1]]), Matrix([[1, 1], [0, 1]]), Matrix([[1, 0], [0, -1]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math191.png)
|
|
![[Matrix(%id = 18446790986273155062), Matrix(%id = 18446790986273155182), Matrix(%id = 18446790986273155302), Matrix(%id = 18446790986273155422)]](/support/helpjp/helpview.aspx?si=5573/file01426/math194.png)
| (7) |
|
>
|
|
|
![[Matrix(%id = 18446790986273155062), Matrix(%id = 18446790986273155182), Matrix(%id = 18446790986273155422)]](/support/helpjp/helpview.aspx?si=5573/file01426/math203.png)
| (8) |
Example 4.
|
>
|
![S4 := evalDG([D_x, x*D_x, (3*x-4)*D_x])](/support/helpjp/helpview.aspx?si=5573/file01426/math209.png)
|
|
![[_DG([["vector", M, []], [[[1], 1]]]), _DG([["vector", M, []], [[[1], x]]]), _DG([["vector", M, []], [[[1], 3*x-4]]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math212.png)
| (9) |
Here is the basis for S4 using function coefficients.
|
>
|
|
|
![[_DG([["vector", M, []], [[[1], 1]]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math223.png)
| (10) |
Here is the basis for S4 using real number coefficients.
|
>
|
|
|
![[_DG([["vector", M, []], [[[1], 1]]]), _DG([["vector", M, []], [[[1], x]]])]](/support/helpjp/helpview.aspx?si=5573/file01426/math234.png)
| (11) |
|
M >
|
|
|
|
Download Help Document
Was this information helpful?
