diffalg
Rosenfeld_Groebner
compute a characteristic decomposition of the radical differential ideal generated by a finite set of differential polynomials
Calling Sequence
Parameters
Description
Examples
Rosenfeld_Groebner (S, H, R, J)
S
-
list or set of differential polynomials of R
H
(optional) list or a set of differential polynomials of R
R
differential polynomial ring
J
(optional) radical differential ideal
Important: The diffalg package has been deprecated. Use the superseding package DifferentialAlgebra instead.
For an informal presentation, see the diffalg overview.
Rosenfeld_Groebner computes a characteristic decomposition of the radical differential ideal P = {S}:(H)^infinity.
If the parameter H is omitted, Rosenfeld_Groebner computes a characteristic decomposition of the radical differential ideal P={S} generated by the differential polynomials of S.
R is a differential polynomial ring constructed with the differential_ring command.
The output of Rosenfeld_Groebner depends on the ranking defined on R.
Rosenfeld_Groebner returns a list of characterizable differential ideals.
The empty list denotes the unit ideal (meaning that there is no solution).
Each characterizable differential ideal is stored in a table. Only the name of the table (characterizable) is printed on the screen. To access their defining characteristic sets you can use the commands rewrite_rules, equations, and inequations.
If the fourth parameter J is present, it is assumed to be another representation of P with respect to another ranking. It is used to spare some splittings. It can be used to speed up the computation, for example, if there is a natural ranking to compute the representation of P.
The command with(diffalg,Rosenfeld_Groebner) allows the use of the abbreviated form of this command.
with⁡diffalg:
The first example illustrates how the Rosenfeld_Groebner command splits a system of differential equations into a system representing the general solution and systems representing the singular solutions.
R≔differential_ring⁡ranking=z,y,derivations=x,notation=diff
R≔ODE_ring
eq1≔−y⁡x+x⁢diff⁡y⁡x,x+diff⁡y⁡x,x2+diff⁡z⁡x,x
eq1≔−y⁡x+x⁢ⅆⅆxy⁡x+ⅆⅆxy⁡x2+ⅆⅆxz⁡x
eq2≔−z⁡x+x⁢diff⁡z⁡x,x+diff⁡y⁡x,x⁢diff⁡z⁡x,x
eq2≔−z⁡x+x⁢ⅆⅆxz⁡x+ⅆⅆxy⁡x⁢ⅆⅆxz⁡x
P≔Rosenfeld_Groebner⁡eq1,eq2,R
P≔characterizable,characterizable,characterizable
equations⁡P1,inequations⁡P1
−y⁡x+x⁢ⅆⅆxy⁡x+ⅆⅆxy⁡x2+ⅆⅆxz⁡x,ⅆⅆxy⁡x3+2⁢ⅆⅆxy⁡x2⁢x−ⅆⅆxy⁡x⁢y⁡x+ⅆⅆxy⁡x⁢x2+z⁡x−y⁡x⁢x,x2+4⁢x⁢ⅆⅆxy⁡x+3⁢ⅆⅆxy⁡x2−y⁡x
equations⁡P2,inequations⁡P2
6⁢ⅆⅆxy⁡x⁢y⁡x+2⁢ⅆⅆxy⁡x⁢x2−9⁢z⁡x+7⁢y⁡x⁢x+2⁢x3,27⁢z⁡x2−18⁢z⁡x⁢y⁡x⁢x−4⁢z⁡x⁢x3−4⁢y⁡x3−y⁡x2⁢x2,3⁢y⁡x+x2,−2⁢x3−9⁢y⁡x⁢x+27⁢z⁡x
equations⁡P3,inequations⁡P3
27⁢z⁡x+x3,3⁢y⁡x+x2,
To obtain the characterizable differential ideal representing the general solution alone, we can proceed as follows.
G≔Rosenfeld_Groebner⁡eq1,eq2,x2+4⁢x⁢diff⁡y⁡x,x−y⁡x+3⁢diff⁡y⁡x,x2,R
G≔characterizable
It is sometimes the case that the radical differential ideal P generated by S is prime. This can be proved by exhibiting a ranking for which the characteristic decomposition of P consists of only one orthonomic characterizable differential ideal.
Before computing a representation of P with respect to the ranking of R, it may be useful to proceed as follows. Search for a ranking for which the characteristic decomposition is as described above. Assign J this computed characteristic decomposition. Then call Rosenfeld_Groebner with J as fourth parameter.
With such a fourth parameter, whatever the ranking of R is, the computed representation of P consists of only one characterizable differential ideal.
If J consists of a single non-orthonomic component or has more than one characterizable component, Rosenfeld_Groebner uses the information to avoid unnecessary splittings.
The example below illustrates this behavior for Euler's equations for an incompressible fluid in two dimensions.
p1≔v1t+v1⁢v1x+v2⁢v1y+px:
p2≔v2t+v1⁢v2x+v2⁢v2y+py:
p3≔v1x+v2y:
R≔differential_ring⁡derivations=x,y,t,ranking=v1,v2,p
R≔PDE_ring
Rosenfeld_Groebner⁡p1,p2,p3,R
characterizable,characterizable,characterizable
S≔differential_ring⁡derivations=t,x,y,ranking=lexp,v1,v2
S≔PDE_ring
J≔Rosenfeld_Groebner⁡p1,p2,p3,S
J≔characterizable
is_orthonomic⁡J
true
Rosenfeld_Groebner⁡p1,p2,p3,R,J
characterizable
See Also
diffalg(deprecated)
diffalg(deprecated)/differential_algebra
diffalg(deprecated)/Rosenfeld_Groebner_options
diffalg(deprecated)[differential_ring]
diffalg(deprecated)[equations]
diffalg(deprecated)[inequations]
diffalg(deprecated)[is_orthonomic]
diffalg(deprecated)[rewrite_rules]
Download Help Document