Overview of the LAVF Object - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : Lie Symmetry Method : Commands for PDEs (and ODEs) : LieAlgebrasOfVectorFields : LAVF : LieAlgebrasOfVectorFields/LAVF

Overview of the LAVF Object

 

Description

LAVF Object Methods

Examples

Description

• 

The LAVF object is designed and created to represent a Lie algebras of vector field (LAVF). An LAVF is usually of the form before solving the symmetry determining system to get solution vector fields.

• 

There are more than 60 methods that are available for a LAVF object, types of methods including (i) basic and solution properties of a LAVF, (ii) Lie algebraic and geometrical properties of a LAVF, (iii) Lie algebraic and geometrical relationships among LAVFs, (iv) utility methods for manipulating a LAVF (or LAVFs). Some Maple existing builtins are extended for allowing LAVF object.

• 

A valid LAVF object is automatically in a rif-reduced form with respect to a ranking that is in total degree ranking. See methods IsRifReduced, IsTotalDegreeRanking and GetRanking of a LHPDE object.

• 

The LAVF object is allowed to be either of finite or infinite type. And the form of its determining system can be un-solved, or solved, or any forms in between (i.e. partially solved).

• 

All methods of the LAVF object become available only once a valid LAVF object is constructed successfully. To construct a LAVF object, see LieAlgebrasOfVectorFields[LAVF].

• 

The LAVF object is the main Maple object exported by the LieAlgebrasOfVectorFields package. See Overview of the LieAlgebrasOfVectorFields package for more detail.

• 

A LAVF object is mathematically represented by the formal vector field as a VectorField object and the determining system as a LHPDE object. Here the components of the formal vector field must be (part of) the dependent variables of the determining system.These two data attributes can be accessed via the GetVectorField and GetDeterminingSystem methods.

• 

After a LAVF object L is successfully constructed, each method in L can be accessed by either the short form method(L, arguments) or the long form L:-method(L, arguments).

LAVF Object Methods

• 

After a LAVF object is constructed, the following methods are available:

AdjointMatrix

AreCommuting

AreSame

AreSameSpace

Centraliser

Centre

CleanDependencies

Copy

DChange

DerivedAlgebra

DerivedSeries

GetDeterminingSystem

GetIDBasis

GetSpace

GetVectorField

Hypercentre

ImplicitForm

Intersection

InvariantCount

Invariants

IsAbelian

IsCommutative

IsFiniteType

IsFlat

IsIdeal

IsInvariant

IsLieAlgebra

IsNilpotent

IsPerfect

IsReductive

IsSemiSimple

IsSolvable

IsSubspace

IsTransitive

IsTrivial

KillingForm

KillingOrthogonal

KillingPolynomial

KillingRadical

LAVFSolve

LieProduct

LowerCentralSeries

NilpotentRadical

NilRadical

Normaliser

OrbitDimension

OrbitDistribution

ParametricDerivatives

SetIDBasis

SolutionDimension

SolvableRadical

StructureCoefficients

StructureConstants

Transporter

UpperCentralSeries

VectorSpaceSum

• 

The following Maple builtins functions are extended so that they work for a LAVF object: type, has, hastype, indets. See LAVF Object Overloaded Builtins for more detail.

Examples

with(LieAlgebrasOfVectorFields):

Typesetting:-Settings(userep=true):

Typesetting:-Suppress([xi(x,y),eta(x,y)]):

Example1: we first require to construct a vector field and a LHPDE object for representing the determining system for E(2).

V := VectorField([[xi(x,y),x], [eta(x,y),y]], space = [x,y]);

Vξx+ηy

(1)

E2Sys := LHPDE([diff(xi(x,y),y,y)=0, diff(eta(x,y),x)+diff(xi(x,y),y)=0, diff(eta(x,y),y)=0, diff(xi(x,y),x)=0], dep = [xi,eta]);

E2Sysξy,y=0,ηx+ξy=0,ηy=0,ξx=0,indep=x,y,dep=ξ,η

(2)

E2 := LAVF(V, E2Sys);

E2ξx+ηy&whereξy,y=0,ξx=0,ηx=ξy,ηy=0

(3)

E2 is a valid LAVF object and a collection of methods are available for E(2).

type(E2, 'LAVF');

true

(4)

exports(E2, 'static');

indets,has,hastype,type,GetVectorField,GetDeterminingSystem,ImplicitForm,SolutionDimension,IsFiniteType,IsTrivial,ParametricDerivatives,GetRanking,SetIDBasis,GetIDBasis,GetSpace,IsFlat,OrbitDistribution,OrbitDimension,InvariantCount,IsTransitive,Invariants,IsLieAlgebra,IsPerfect,DerivedAlgebra,IsSolvable,IsSoluble,DerivedSeries,SolvableRadical,SolubleRadical,Radical,IsNilpotent,Hypercentre,Hypercenter,NilRadical,Nilradical,LowerCentralSeries,UpperCentralSeries,IsAbelian,IsCommutative,Centre,Center,IsSemiSimple,IsReductive,NilpotentRadical,StructureConstants,StructureCoefficients,KillingRadical,KillingPolynomial,KillingForm,KillingOrthogonal,AdjointMatrix,AreCommuting,AreSame,AreSameSpace,Centraliser,Centralizer,Normaliser,CleanDependencies,Copy,DChange,dchange,Intersection,IsIdeal,IsInvariant,IsotropyRepresentation,IsSubspace,LAVFSolve,VectorSpaceSum,LieProduct,ProjectToSpace,Transporter,ModuleCopy,ModulePrint,ModuleApply

(5)

Basic properties of E2 can be obtained by:

X := GetVectorField(E2);

Xξx+ηy

(6)

S := GetDeterminingSystem(E2);

Sξy,y=0,ξx=0,ηx=ξy,ηy=0,indep=x,y,dep=ξ,η

(7)

Both X and S are Maple objects too, and have access to various methods too. For example,

type(S, 'LHPDE');

true

(8)

IsRifReduced(S);

true

(9)

IsTotalDegreeRanking(S);

true

(10)

GetRanking(S);

ξ,η

(11)

We can check algebraic properties of E2.

IsLieAlgebra(E2);

true

(12)

IsAbelian(E2);

false

(13)

IsSolvable(E2);

true

(14)

Fetch structure constants of E2 & display it.

C := StructureConstants(E2);

(15)

DisplayStructure(C, 'Z');

00Z200Z1Z2Z10

(16)

Some properties of E2 are also represented as LAVF objects (or lists of LAVF objects).

Centre(E2);

ξx+ηy&whereξ=0,η=0

(17)

SolvableRadical(E2);

ξx+ηy&whereξy,y=0,ξx=0,ηx=ξy,ηy=0

(18)

DerivedSeries(E2);

ξx+ηy&whereξy,y=0,ξx=0,ηx=ξy,ηy=0,ξx+ηy&whereξx=0,ηx=0,ξy=0,ηy=0,ξx+ηy&whereξ=0,η=0

(19)

A LAVF object can be converted into a PDO for vector fields.

VFPDO(E2);

X→yyXx,ⅆⅆxXx,xXy+yXx,ⅆⅆyXy

(20)

Solving E2 vector fields system to find vector fields.

LAVFSolve(E2);

_C1y+_C3x+_C1x+_C2y

(21)

Example2: now we consider three 3-d rotation vector fields that generating standard SO(3) action on R3(x,y,z). First we construct these vector fields.

Typesetting:-Suppress({xi,eta,zeta}(x,y,z)):

R1 := VectorField(x*D[y] - y*D[x], space = [x,y,z]);

R1yx+xy

(22)

R2 := VectorField(x*D[z] - z*D[x], space = [x,y,z]);

R2zx+xz

(23)

R3 := VectorField(y*D[z] - z*D[y], space = [x,y,z]);

R3zy+yz

(24)

X := VectorField(xi(x,y,z)*D[x] + eta(x,y,z)*D[y] + zeta(x,y,z)*D[z], space = [x,y,z]);

Xξx+ηy+ζz

(25)

Then a LAVF object as the vector fields system for SO(3) can be constructed by

SO3 := EliminationLAVF(X, [R1, R2, R3]);

SO3ξx+ηy+ζz&whereξ=ηyζzx,ηx=ζyz+ηx,ηy=0,ηz=ζy,ζy,y=0,ζx=ζyy+ζx,ζz=0

(26)

Some geometric properties of SO(3).

OrbitDimension(SO3);

2

(27)

IsTransitive(SO3);

false

(28)

InvariantCount(SO3);

1

(29)

Invariants may be found via integration.

Invariants(SO3);

x2+y2+z2

(30)

We can explore the geometric properties of SO3 in further detail by extracting its orbit distribution.

OD := OrbitDistribution(SO3);

ODyxx+y,zxx+z

(31)

OD is also a Maple object named Distribution, and it has access to various methods.

type(OD, 'Distribution');

true

(32)

exports(OD, 'static');

GetSpace,GetVectorFields,GetAnnihilator,AreSameSpace,IsTrivial,Dimension,Codimension,IsInvolutive,IsIntegrable,Integrals,IsSubspace,Intersection,VectorSpaceSum,IsInvariant,DerivedDistribution,CauchyDistribution,dchange,DChange,indets,has,hastype,type,convert

(33)

See Also

LieAlgebrasOfVectorFields (Package overview)

LAVF (Object overview)

LHPDE (Object overview)

VectorField (Object overview)

LieAlgebrasOfVectorFields[VectorField]

LieAlgebrasOfVectorFields[LHPDE]

LieAlgebrasOfVectorFields[LAVF]