Tensor[SymmetricProductsOfKillingTensors] - form all possible symmetric tensors of a given rank (linearly independent over the real numbers) from a list of symmetric tensors
Calling Sequences
SymmetricProducstOfKillingTensors(K, p, ptlist )
Parameters
K - a list of lists of Killing tensors K ={K1, K2, ...], where Kr is a list of Killing tensors of rank r on a manifold M
p - a positive integer.
ptlist - (optional) a list of points on the manifold M
Description
Examples
This program first computes all partitions p1,p2,...of the integer p, that is, all non-decreasing sequences p1,p2,...with p1+p2+⋅⋅⋅=p. Then, for each partition, all possible rank p symmetric tensors of the form T1⊙T2⊙T3⊙... ,where the tensor Ti is taken from the list K, are generated. From the totality of tensors so obtained a maximal set of linearly independent tensors (over the real numbers) is selected. Each symmetric tensor in the returned list is a Killing tensor if each of the Kiare.
The independent tensors are generated by a call to the DifferentialGeometry command DGbasis. For tensors with coefficients which are not rational functions, the DGbasis program may work faster using a Wronskian approach which requires the specification of a list of points on the manifold M.
This command is part of the DifferentialGeometry:-Tensor package, and so can be used in the form SymmetricProductsOfKillingTensors(...) only after executing the commands with(DifferentialGeometry), with(Tensor) in that order. It can always be used in the long form DifferentialGeometry:-Tensor:-SymmetricProductsOfKillingTensors.
with⁡DifferentialGeometry:with⁡Tensor:
Example 1.
DGsetup⁡x,y,M
frame name: M
g≔evalDG⁡dx&sdy
g≔dx2⁢dy+dy2⁢dx
Calculate the rank 1 Killing tensors.
K1≔KillingTensors⁡g,1
K1≔−y⁢dx+x⁢dy,dy,dx
Calculate all the rank 2 and rank 3 Killing tensors generated by K1.
S2≔SymmetricProductsOfKillingTensors⁡K1,2
S2≔y2⁢dx⁢dx−y⁢x⁢dx⁢dy−y⁢x⁢dy⁢dx+x2⁢dy⁢dy,−y2⁢dx⁢dy−y2⁢dy⁢dx+x⁢dy⁢dy,−y⁢dx⁢dx+x2⁢dx⁢dy+x2⁢dy⁢dx,dy⁢dy,dx2⁢dy+dy2⁢dx,dx⁢dx
nops⁡S2
6
S3≔SymmetricProductsOfKillingTensors⁡K1,3
S3≔−y3⁢dx⁢dx⁢dx+y2⁢x⁢dx⁢dx⁢dy+y2⁢x⁢dx⁢dy⁢dx−y⁢x2⁢dx⁢dy⁢dy+y2⁢x⁢dy⁢dx⁢dx−y⁢x2⁢dy⁢dx⁢dy−y⁢x2⁢dy⁢dy⁢dx+x3⁢dy⁢dy⁢dy,y23⁢dx⁢dx⁢dy+y23⁢dx⁢dy⁢dx−2⁢y⁢x3⁢dx⁢dy⁢dy+y23⁢dy⁢dx⁢dx−2⁢y⁢x3⁢dy⁢dx⁢dy−2⁢y⁢x3⁢dy⁢dy⁢dx+x2⁢dy⁢dy⁢dy,y2⁢dx⁢dx⁢dx−2⁢y⁢x3⁢dx⁢dx⁢dy−2⁢y⁢x3⁢dx⁢dy⁢dx+x23⁢dx⁢dy⁢dy−2⁢y⁢x3⁢dy⁢dx⁢dx+x23⁢dy⁢dx⁢dy+x23⁢dy⁢dy⁢dx,−y3⁢dx⁢dy⁢dy−y3⁢dy⁢dx⁢dy−y3⁢dy⁢dy⁢dx+x⁢dy⁢dy⁢dy,−y3⁢dx⁢dx⁢dy−y3⁢dx⁢dy⁢dx+x3⁢dx⁢dy⁢dy−y3⁢dy⁢dx⁢dx+x3⁢dy⁢dx⁢dy+x3⁢dy⁢dy⁢dx,−y⁢dx⁢dx⁢dx+x3⁢dx⁢dx⁢dy+x3⁢dx⁢dy⁢dx+x3⁢dy⁢dx⁢dx,dy⁢dy⁢dy,dx3⁢dy⁢dy+dy3⁢dx⁢dy+dy3⁢dy⁢dx,dx3⁢dx⁢dy+dx3⁢dy⁢dx+dy3⁢dx⁢dx,dx⁢dx⁢dx
nops⁡S3
10
Check that all the tensors in S2 are indeed Killing tensors using the CheckKillingTensor command.
map2⁡CheckKillingTensor,g,S2
0⁢dx⁢dx⁢dx,0⁢dx⁢dx⁢dx,0⁢dx⁢dx⁢dx,0⁢dx⁢dx⁢dx,0⁢dx⁢dx⁢dx,0⁢dx⁢dx⁢dx
Calculate all the rank 3 Killing tensors generated by S3 and the metric g. We get the list S3 again because the metric is already a constant linear combination of the rank 1 Killing tensors.
S3a≔SymmetricProductsOfKillingTensors⁡K1,g,3:
See Also
DifferentialGeometry
Tensor
CheckKillingTensor
DGbasis
GenerateTensors
GetComponents
KillingTensors
SymmetrizeIndices
Download Help Document