Important: The diffalg package has been deprecated. Use the superseding package DifferentialAlgebra instead.

A number of environment variables can be set to modify the computations made by Rosenfeld_Groebner.

_Env_diffalg_gen_order

The input differential polynomials to Rosenfeld_Groebner are treated in sequence. By default they are treated by increasing rank. Setting _Env_diffalg_gen_order := 'decreasingly' might change the computation time as well as the characteristic decomposition in the output.

To return to the default behavior, use _Env_diffalg_gen_order := 'increasingly'.

_Env_diffalg_char

Rosenfeld_Groebner proceeds in two fundamental steps. It first computes a decomposition into regular differential ideals. The second step consists of computing a decomposition into characterizable differential ideals of these regular differential ideals. The default output is  therefore a characteristic decomposition of the radical differential ideal generated by the input differential polynomials.

Though it is rarely  the bottleneck of the algorithm, you might want to compute only the decomposition into regular differential ideals. This can be achieved by setting _Env_diffalg_char := false. The default is true.

Regular differential ideals are represented in the same way as characterizable differential ideals. They are tables and appear as regular. The regular differential systems of equations and inequations defining them can be accessed by using the commands equations and inequations. Unlike characterizable differential ideals, the inequations can contain differential polynomials that are unrelated to the initials and separants of the equations.

_Env_diffalg_charpres_method

This option allows to choose between two different methods to compute a decomposition into characterizable differential ideals from the regular differential ideals.

By setting _Env_diffalg_charpres_method:='groebner', the final decomposition is processed using Groebner basis computations. By setting _Env_diffalg_charpres_method:='regchar', the final decomposition is processed using algebraic triangular technics. It might happen that one method succeeds whereas the other does not (see example below).

The default is _Env_diffalg_charpres_method:='groebner'.

_Env_diffalg_charpres_normalize

This option works together with the option _Env_diffalg_charpres_method.  If you use _Env_diffalg_charpres_method:='regchar', you can choose to compute normalized or non-normalized characterizable differential ideals. A normalized characterizable differential ideal satisfies the property : the initials of its equations are free of leaders. A non-normalized characterizable differential ideal does not satisfy this property. Setting _Env_diffalg_charpres_normalize:=true produces normalized outputs and _Env_diffalg_charpres_normalize:=false produces non-normalized ones.

The size of the output can be smaller by computing non-normalized characterizable differential ideals (see example below).

The default is _Env_diffalg_charpres_normalize:=true.

_Env_diffalg_primestudy

This environment variable affects the computation if Rosenfeld_Groebner is called with a fourth parameter that is a characteristic decomposition, with respect to another ranking, of the radical ideal to be treated. If this characteristic decomposition  consists of a single characterizable differential ideal that you know is prime, setting _Env_diffalg_primestudy := true takes advantage of this knowledge.

If the characterizable differential ideal passes the is_orthonomic test, it is not necessary to set _Env_diffalg_primestudy.

The default is _Env_diffalg_primestudy := false.

_Env_diffalg_split

Setting _Env_diffalg_split := false eliminates the branching on the initials and separants at each step. The computation seeks what is deduced to be the most general component. This saves computations. Nonetheless, the result might be misleading. The computation is indeed incomplete and does not provide a characteristic decomposition. The result may lead you to infer that there is no common zero to the input differential polynomials when there is one in fact.

$\mathrm{with}\left(\mathrm{diffalg}\right)\:$

$\mathrm{\_Env\_diffalg\_uses\_DifferentialAlgebra}≔\mathrm{false}$

$\mathrm{\_Env\_diffalg\_uses\_DifferentialAlgebra}≔\mathrm{false}$

$R≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[\left[v\,u\right]\,\left[p\,q\,w\,r\right]\right]\,\mathrm{derivations}=\left[y\,x\right]\right)\:$

$S≔\left[\left(u\left[x\right]-u\left[y\right]\right)v\left[x,x\right]-q\left[\right]\,\left(u\left[x\right]-u\left[y\right]\right)v\left[x\right]-p\left[\right]\,{u\left[x\right]}^2-w\left[\right]\,{u\left[y\right]}^2+w\left[\right]\right]$

$S≔\left[\left(u_x-u_y\right)v_{x,x}-q\left[\right]\,\left(u_x-u_y\right)v_x-p\left[\right]\,u_x^2-w\left[\right]\,u_y^2+w\left[\right]\right]$

$\mathrm{\_Env\_diffalg\_char}≔\mathrm{false}\:$$\mathrm{RegularDecomp}≔\mathrm{Rosenfeld\_Groebner}\left(S\,R\right)$

$\mathrm{RegularDecomp}≔\left[\mathrm{regular}\,\mathrm{regular}\,\mathrm{regular}\,\mathrm{regular}\,\mathrm{regular}\,\mathrm{regular}\,\mathrm{regular}\right]$

$\mathrm{\_Env\_diffalg\_char}≔\mathrm{true}\:$$\mathrm{CharacteristicDecomp}≔\mathrm{Rosenfeld\_Groebner}\left(S\,R\right)$

$\mathrm{CharacteristicDecomp}≔\left[\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\right]$

$R≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[y\,x\,p\right]\,\mathrm{derivations}=\left[t\right]\right)\:$

$S≔\left\{\left(1+{x\left[\right]}^2+2x\left[\right]\right){y\left[\right]}^2-\left(1+{x\left[\right]}^2-2x\left[\right]\right)\,\left(x\left[\right]+1\right)y\left[\right]+x\left[\right]-1\,{p\left[\right]}^2-1\,{x\left[\right]}^2-p\left[\right]\right\}$

$S≔\left\{\left(x\left[\right]+1\right)y\left[\right]+x\left[\right]-1\,\left({x\left[\right]}^2+2x\left[\right]+1\right){y\left[\right]}^2-{x\left[\right]}^2+2x\left[\right]-1\,{p\left[\right]}^2-1\,{x\left[\right]}^2-p\left[\right]\right\}$

$\mathrm{\_Env\_diffalg\_char}≔\mathrm{false}\:$$\mathrm{RegularDecomp}≔\mathrm{Rosenfeld\_Groebner}\left(S\,R\right)$

$\mathrm{RegularDecomp}≔\left[\mathrm{regular}\right]$

$\mathrm{\_Env\_diffalg\_char}≔\mathrm{true}\:$$\mathrm{CharacteristicDecomp}≔\mathrm{Rosenfeld\_Groebner}\left(S\,R\right)$

$\mathrm{CharacteristicDecomp}≔\left[\mathrm{characterizable}\,\mathrm{characterizable}\right]$

$\mathrm{rewrite\_rules}\left(\mathrm{RegularDecomp}\right)$

$\left[\left[y\left[\right]=-\frac{x\left[\right]-1}{x\left[\right]+1}\,{x\left[\right]}^2=p\left[\right]\,{p\left[\right]}^2=1\right]\right]$

$\mathrm{rewrite\_rules}\left(\mathrm{CharacteristicDecomp}\right)$

$\left[\left[y\left[\right]=-x\left[\right]+\frac{p\left[\right]}{2}+\frac{1}{2}\,x\left[\right]=1\,p\left[\right]=1\right]\,\left[y\left[\right]=-x\left[\right]+\frac{p\left[\right]}{2}+\frac{1}{2}\,{x\left[\right]}^2=p\left[\right]\,p\left[\right]=\mathrm{-1}\right]\right]$

$R≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[y\,x\right]\,\mathrm{derivations}=\left[t\right]\right)\:$

$S≔\left[{y\left[t\right]}^2+y\left[t\right]y\left[\right]-x\left[\right]y\left[t\right]+tx\left[\right]y\left[t\right]\,y\left[\right]-x\left[\right]+x\left[t\right]\right]\;$$H≔\left[x\left[t,t,t\right]-x\left[t\right]\right]$

$S≔\left[tx\left[\right]y_t-x\left[\right]y_t+y_{t}y\left[\right]+y_t^2\,y\left[\right]-x\left[\right]+x_t\right]$

$H≔\left[x_{t,t,t}-x_t\right]$

$\mathrm{\_Env\_diffalg\_char}≔\mathrm{false}\:$$\mathrm{RegularDecomp}≔\mathrm{Rosenfeld\_Groebner}\left(S\,H\,R\right)$

$\mathrm{RegularDecomp}≔\left[\mathrm{regular}\right]$

$\mathrm{\_Env\_diffalg\_char}≔\mathrm{true}\:$$\mathrm{CharacteristicDecomp}≔\mathrm{Rosenfeld\_Groebner}\left(S\,H\,R\right)$

$\mathrm{CharacteristicDecomp}≔\left[\mathrm{characterizable}\right]$

$\mathrm{rewrite\_rules}\left(\mathrm{RegularDecomp}\right)$

$\left[\left[y\left[\right]=x\left[\right]-x_t\,x_{t,t}^2=-tx\left[\right]x_t+tx\left[\right]x_{t,t}+x_t{x}_{t,t}\right]\right]$

$\mathrm{rewrite\_rules}\left(\mathrm{CharacteristicDecomp}\right)$

$\left[\left[y\left[\right]=x\left[\right]-x_t\,x_{t,t}=tx\left[\right]\right]\right]$

$\mathrm{\_Env\_diffalg\_charpres\_method}≔'\mathrm{regchar}'$

$\mathrm{\_Env\_diffalg\_charpres\_method}≔\mathrm{regchar}$

$R≔\mathrm{differential\_ring}\left(\mathrm{derivations}=\left[x\right]\,\mathrm{ranking}=\left[\left[\mathrm{u1}\,\mathrm{u2}\,\mathrm{u3}\,\mathrm{u4}\,\mathrm{u5}\right]\right]\right)$

$R≔\mathrm{ODE\_ring}$

$\mathrm{p1}≔\left(\mathrm{u2}\left[\right]\mathrm{u3}\left[\right]+\mathrm{u4}\left[\right]+\mathrm{u2}\left[\right]\right){\mathrm{u1}\left[\right]}^2+\left(\mathrm{u2}\left[\right]\mathrm{u3}\left[\right]+3\right)\mathrm{u1}\left[\right]+1$

$\mathrm{p1}≔\left(\mathrm{u2}\left[\right]\mathrm{u3}\left[\right]+\mathrm{u2}\left[\right]+\mathrm{u4}\left[\right]\right){\mathrm{u1}\left[\right]}^2+\left(\mathrm{u2}\left[\right]\mathrm{u3}\left[\right]+3\right)\mathrm{u1}\left[\right]+1$

$\mathrm{p2}≔\left(1+\mathrm{u4}\left[\right]+\mathrm{u5}\left[\right]\right){\mathrm{u3}\left[\right]}^2+\mathrm{u3}\left[\right]\left({\mathrm{u5}\left[\right]}^2+\mathrm{u4}\left[\right]+2\right)+\mathrm{u5}\left[\right]+2$

$\mathrm{p2}≔\left(1+\mathrm{u4}\left[\right]+\mathrm{u5}\left[\right]\right){\mathrm{u3}\left[\right]}^2+\mathrm{u3}\left[\right]\left({\mathrm{u5}\left[\right]}^2+\mathrm{u4}\left[\right]+2\right)+\mathrm{u5}\left[\right]+2$

$\mathrm{Rosenfeld\_Groebner}\left(\left[\mathrm{p1}\,\mathrm{p2}\right]\,R\right)$

$\left[\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\,\mathrm{characterizable}\right]$

$\mathrm{\_Env\_diffalg\_charpres\_method}≔'\mathrm{regchar}'$

$\mathrm{\_Env\_diffalg\_charpres\_method}≔\mathrm{regchar}$

$R≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[u\,v\right]\,\mathrm{derivations}=\left[x\,y\right]\right)$

$R≔\mathrm{PDE\_ring}$

$\mathrm{sys}≔\left[u\left[x,x\right]\left(v\left[y\right]+v\left[\right]+1\right)+u\left[y\right]\,v\left[\right]{v\left[y\right]}^5-2v\left[\right]\right]$

$\mathrm{sys}≔\left[u_{x,x}\left(v_y+v\left[\right]+1\right)+u_y\,v\left[\right]v_y^5-2v\left[\right]\right]$

$\mathrm{\_Env\_diffalg\_charpres\_normalize}≔\mathrm{true}$

$\mathrm{\_Env\_diffalg\_charpres\_normalize}≔\mathrm{true}$

$\mathrm{equations}\left(\mathrm{Rosenfeld\_Groebner}\left(\mathrm{sys}\,R\right)\right)$

$\left[\left[u_y+u_{x,x}\,v\left[\right]\right]\,\left[u_{x,x}{v\left[\right]}^5+u_y{v\left[\right]}^4-u_y{v\left[\right]}^3{v}_y+u_y{v\left[\right]}^2{v}_y^2-u_{y}v\left[\right]v_y^3+u_y{v}_y^4+5u_{x,x}{v\left[\right]}^4+4u_y{v\left[\right]}^3-3u_y{v\left[\right]}^2{v}_y+2u_{y}v\left[\right]v_y^2-u_y{v}_y^3+10u_{x,x}{v\left[\right]}^3+6u_y{v\left[\right]}^2-3u_{y}v\left[\right]v_y+u_y{v}_y^2+10u_{x,x}{v\left[\right]}^2+4u_{y}v\left[\right]-u_y{v}_y+5u_{x,x}v\left[\right]+u_y+3u_{x,x}\,v_y^5-2\right]\right]$

$\mathrm{\_Env\_diffalg\_charpres\_normalize}≔\mathrm{false}$

$\mathrm{\_Env\_diffalg\_charpres\_normalize}≔\mathrm{false}$

$\mathrm{equations}\left(\mathrm{Rosenfeld\_Groebner}\left(\mathrm{sys}\,R\right)\right)$

$\left[\left[u_y+u_{x,x}\,v\left[\right]\right]\,\left[u_{x,x}v\left[\right]+u_{x,x}{v}_y+u_y+u_{x,x}\,v_y^5-2\right]\right]$

$\mathrm{R0}≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[\left[y\,x\right]\,p\right]\,\mathrm{derivations}=\left[t\right]\right)\:$

$S≔\left\{x\left[\right]y\left[t\right]+x\left[t\right]+1\,{x\left[t\right]}^2+y\left[\right]+p\left[\right]\right\}\;$$H≔\left\{x\left[\right]\,x\left[t\right]\right\}$

$S≔\left\{x_t^2+p\left[\right]+y\left[\right]\,x\left[\right]y_t+x_t+1\right\}$

$H≔\left\{x\left[\right]\,x_t\right\}$

$\mathrm{C0}≔\mathrm{Rosenfeld\_Groebner}\left(S\,H\,\mathrm{R0}\right)$

$\mathrm{C0}≔\left[\mathrm{characterizable}\right]$

$\mathrm{equations}\left(\mathrm{C0}\left[1\right]\right)$

$\left[x\left[\right]y_t+x_t+1\,x_t^2+p\left[\right]+y\left[\right]\right]$

$R≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[x\,y\,p\right]\,\mathrm{derivations}=\left[t\right]\right)\:$

$\mathrm{\_Env\_diffalg\_prime\_study}≔\mathrm{true}\:$$\mathrm{C2}≔\mathrm{Rosenfeld\_Groebner}\left(S\,H\,R\,\mathrm{C0}\left[1\right]\right)$

$\mathrm{C2}≔\left[\mathrm{characterizable}\right]$

$R≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[\mathrm{\lambda }\,\left[x\,y\right]\right]\,\mathrm{derivations}=\left[t\right]\right)\:$

$\mathrm{Pendulum}≔\left[y\left[t,t\right]+\mathrm{\lambda }\left[\right]y\left[\right]+9\,x\left[t,t\right]+\mathrm{\lambda }\left[\right]x\left[\right]\,{x\left[\right]}^2+{y\left[\right]}^2-1\right]$

$\mathrm{Pendulum}≔\left[\mathrm{\lambda }\left[\right]y\left[\right]+y_{t,t}+9\,\mathrm{\lambda }\left[\right]x\left[\right]+x_{t,t}\,{x\left[\right]}^2+{y\left[\right]}^2-1\right]$

$\mathrm{\_Env\_diffalg\_split}≔\mathrm{false}\:$$\mathrm{Cg}≔\mathrm{Rosenfeld\_Groebner}\left(\mathrm{Pendulum}\,R\right)\;$$\mathrm{rewrite\_rules}\left(\mathrm{Cg}\right)$

$\mathrm{Cg}≔\left[\mathrm{characterizable}\right]$

$\left[\left[\mathrm{\lambda }\left[\right]=-\frac{9{y\left[\right]}^3+y_t^2-9y\left[\right]}{{y\left[\right]}^2-1}\,y_{t,t}=\frac{9{y\left[\right]}^4+y\left[\right]y_t^2-18{y\left[\right]}^2+9}{{y\left[\right]}^2-1}\,{x\left[\right]}^2=-{y\left[\right]}^2+1\right]\right]$

$\mathrm{\_Env\_diffalg\_split}≔\mathrm{true}\:$$C≔\mathrm{Rosenfeld\_Groebner}\left(\mathrm{Pendulum}\,R\right)\;$$\mathrm{rewrite\_rules}\left(C\right)$

$C≔\left[\mathrm{characterizable}\,\mathrm{characterizable}\right]$

$\left[\left[\mathrm{\lambda }\left[\right]=-\frac{9{y\left[\right]}^3+y_t^2-9y\left[\right]}{{y\left[\right]}^2-1}\,y_{t,t}=\frac{9{y\left[\right]}^4+y\left[\right]y_t^2-18{y\left[\right]}^2+9}{{y\left[\right]}^2-1}\,{x\left[\right]}^2=-{y\left[\right]}^2+1\right]\,\left[\mathrm{\lambda }\left[\right]=-\frac{9}{y\left[\right]}\,x\left[\right]=0\,{y\left[\right]}^2=1\right]\right]$

$R≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[y\,x\right]\,\mathrm{derivations}=\left[t\right]\right)\:$

$S≔\left\{\left(x\left[t\right]+1\right){y\left[t\right]}^2+x\left[t\right]-1\,{x\left[t\right]}^2-1\right\}$

$S≔\left\{\left(x_t+1\right)y_t^2+x_t-1\,x_t^2-1\right\}$

$\mathrm{\_Env\_diffalg\_split}≔\mathrm{false}\:$$\mathrm{Rosenfeld\_Groebner}\left(S\,R\right)$

$\mathrm{\_Env\_diffalg\_split}≔\mathrm{true}\:$$\mathrm{Rosenfeld\_Groebner}\left(S\,R\right)\;$$\mathrm{rewrite\_rules}\left(\right)$

$\left[\mathrm{characterizable}\right]$