&graph
1 D Distribution|Orbit
  -> D   rd__o dist_ss pd1_or start_dsor intpar
  -> SS  rd__ss dist_ss pd1_ss start_dss conpar
2 SS Steady state
  -> SS  rd__ss ss_ss pd1_ss start_ssss conpar
3 BP Branching point
8 LP Limit point
&compclass !the way curves of the given type are computed
D   1 !trajectories
SS  2 !implicit functions
BP  2
LP  2
&menu
Initial state{
 D
 SS
}
Curve{
 D|Orbit
 SS
}
&deftypes
! CurveType PointType NewCurveType
      D        D         D
      SS       SS        SS
&stdlinalg
&&$(Paths.KERNEL)stdalg$(Extensions.DESC)
&utilities
_util(CheckParameters)
=(ErrMessage)
=(ActivateTestFuns)
=(IndProcFunc)
&functions
xxx
RHS
 _Rhs+
 int $(Compiler.RHSATTR) _Rhs(double *_u, double *_u_x, double *_u_xx, double *_s, double *_p, double *_t, double *_f) {~
  @Phase {#define $Phase (*Phase(_u+&Phase))~}
  @Phase {#define $Phase_@Space {$Space} (*Phase(_u_x+&Phase))~}
  @Phase {#define $Phase_@Space {$Space_$Space} (*Phase(_u_xx+&Phase))~}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Space {#define $Space _s[0]~}
  @Phase {#define $Phase_@Time {$Time} (*Phase(_f+&Phase))~}
  /*
 return 0;~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Space {#undef $Space~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Phase {#undef $Phase_@Space {$Space_$Space}~}
  @Phase {#undef $Phase_@Time {$Time}~}
  /*
  @Phase { $Phase_@Time {$Time}}
  @Phase { $Phase}
  @Phase { $Phase_@Space {$Space}}
  @Phase { $Phase_@Space {$Space_$Space}}
  @Space { $Space}
  @Par   { $Par}
  /*
Derivatives{1st order
1st order derivatives of RHS
 _Rhs_Der1
 typedef double _d1T[COMPONENTSDIM][3*COMPONENTSDIM+PARAMETERSDIM+SPACEDIM];~static _d1T _d1;~double * $(Compiler.RHSATTR) _Rhs_Der1(double *_1u, double *_2u_x, double *_3u_xx, double *_4s, double *_5p, double *_t) {~
  @Phase {#define $Phase (*Phase(_1u+&Phase))~}
  @Phase {#define $Phase_@Space {$Space} (*Phase(_2u_x+&Phase))~}
  @Phase {#define $Phase_@Space {$Space_$Space} (*Phase(_3u_xx+&Phase))~}
  @Par   {#define $Par (*Par(_5p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Space {#define $Space _4s[0]~}
  @Phase {#define $Phase_@Time {$Time} (*Phase*(_d1T*)(_d1+&Phase))~}
  #define der1(f,u) _d1[&f-_1u][_c(u)]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u : (#v[4]=='2' ? _2u_x-COMPONENTSDIM : (#v[4]=='3' ? _3u_xx-2*COMPONENTSDIM : (#v[4]=='4' ? _4s-3*COMPONENTSDIM : _5p-3*COMPONENTSDIM-SPACEDIM)))))~
  /*
 return _d1[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Space {#undef $Space~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Phase {#undef $Phase_@Space {$Space_$Space}~}
  @Phase {#undef $Phase_@Time {$Time}~}
  #undef der1~
  #undef _c~
  /*
 xxx
  /*
2nd order
2nd order derivatives of RHS
 _Rhs_Der2
 #define _n 3*COMPONENTSDIM+PARAMETERSDIM+SPACEDIM-1~typedef double _d2T[COMPONENTSDIM][dersub2(_n,_n)+1];~static _d2T _d2;~#undef _n~double * $(Compiler.RHSATTR) _Rhs_Der2(double *_1u, double *_2u_x, double *_3u_xx, double *_4s, double *_5p, double *_t) {~
  @Phase {#define $Phase (*Phase(_1u+&Phase))~}
  @Phase {#define $Phase_@Space {$Space} (*Phase(_2u_x+&Phase))~}
  @Phase {#define $Phase_@Space {$Space_$Space} (*Phase(_3u_xx+&Phase))~}
  @Par   {#define $Par (*Par(_5p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Space {#define $Space _4s[0]~}
  @Phase {#define $Phase_@Time {$Time} (*Phase*(_d2T*)(_d2+&Phase))~}
  #define der2(f,u,v) _d2[&f-_1u][dersub2(_c(u),_c(v))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u : (#v[4]=='2' ? _2u_x-COMPONENTSDIM : (#v[4]=='3' ? _3u_xx-2*COMPONENTSDIM : (#v[4]=='4' ? _4s-3*COMPONENTSDIM : _5p-3*COMPONENTSDIM-SPACEDIM)))))~
  /*
 return _d2[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Space {#undef $Space~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Phase {#undef $Phase_@Space {$Space_$Space}~}
  @Phase {#undef $Phase_@Time {$Time}~}
  #undef der2~
  #undef _c~
  /*
 xxx
  /*
3rd order
3rd order derivatives of RHS
 _Rhs_Der3
 #define _n 3*COMPONENTSDIM+PARAMETERSDIM+SPACEDIM-1~typedef double _d3T[COMPONENTSDIM][dersub3(_n,_n,_n)+1];~static _d3T _d3;~#undef _n~double * $(Compiler.RHSATTR) _Rhs_Der3(double *_1u, double *_2u_x, double *_3u_xx, double *_4s, double *_5p, double *_t) {~
  @Phase {#define $Phase (*Phase(_1u+&Phase))~}
  @Phase {#define $Phase_@Space {$Space} (*Phase(_2u_x+&Phase))~}
  @Phase {#define $Phase_@Space {$Space_$Space} (*Phase(_3u_xx+&Phase))~}
  @Par   {#define $Par (*Par(_5p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Space {#define $Space _4s[0]~}
  @Phase {#define $Phase_@Time {$Time} (*Phase*(_d3T*)(_d3+&Phase))~}
  #define der3(f,u,v,w) _d3[&f-_1u][dersub3(_c(u),_c(v),_c(w))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u : (#v[4]=='2' ? _2u_x-COMPONENTSDIM : (#v[4]=='3' ? _3u_xx-2*COMPONENTSDIM : (#v[4]=='4' ? _4s-3*COMPONENTSDIM : _5p-3*COMPONENTSDIM-SPACEDIM)))))~
  /*
 return _d3[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Space {#undef $Space~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Phase {#undef $Phase_@Space {$Space_$Space}~}
  @Phase {#undef $Phase_@Time {$Time}~}
  #undef der3~
  #undef _c~
  /*
 xxx
  /*
4th order
4th order derivatives of RHS
 _Rhs_Der4
 #define _n 3*COMPONENTSDIM+PARAMETERSDIM+SPACEDIM-1~typedef double _d4T[COMPONENTSDIM][dersub4(_n,_n,_n,_n)+1];~static _d4T _d4;~#undef _n~double * $(Compiler.RHSATTR) _Rhs_Der4(double *_1u, double *_2u_x, double *_3u_xx, double *_4s, double *_5p, double *_t) {~
  @Phase {#define $Phase (*Phase(_1u+&Phase))~}
  @Phase {#define $Phase_@Space {$Space} (*Phase(_2u_x+&Phase))~}
  @Phase {#define $Phase_@Space {$Space_$Space} (*Phase(_3u_xx+&Phase))~}
  @Par   {#define $Par (*Par(_5p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Space {#define $Space _4s[0]~}
  @Phase {#define $Phase_@Time {$Time} (*Phase*(_d4T*)(_d4+&Phase))~}
  #define der4(f,u,v,w,x) _d4[&f-_1u][dersub4(_c(u),_c(v),_c(w),_c(x))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u : (#v[4]=='2' ? _2u_x-COMPONENTSDIM : (#v[4]=='3' ? _3u_xx-2*COMPONENTSDIM : (#v[4]=='4' ? _4s-3*COMPONENTSDIM : _5p-3*COMPONENTSDIM-SPACEDIM)))))~
  /*
 return _d4[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Space {#undef $Space~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Phase {#undef $Phase_@Space {$Space_$Space}~}
  @Phase {#undef $Phase_@Time {$Time}~}
  #undef der4~
  #undef _c~
  /*
 xxx
  /*
5th order}
5th order derivatives of RHS
 _Rhs_Der5
 #define _n 3*COMPONENTSDIM+PARAMETERSDIM+SPACEDIM-1~typedef double _d5T[COMPONENTSDIM][dersub5(_n,_n,_n,_n,_n)+1];~static _d5T _d5;~#undef _n~double * $(Compiler.RHSATTR) _Rhs_Der5(double *_1u, double *_2u_x, double *_3u_xx, double *_4s, double *_5p, double *_t) {~
  @Phase {#define $Phase (*Phase(_1u+&Phase))~}
  @Phase {#define $Phase_@Space {$Space} (*Phase(_2u_x+&Phase))~}
  @Phase {#define $Phase_@Space {$Space_$Space} (*Phase(_3u_xx+&Phase))~}
  @Par   {#define $Par (*Par(_5p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Space {#define $Space _4s[0]~}
  @Phase {#define $Phase_@Time {$Time} (*Phase*(_d5T*)(_d5+&Phase))~}
  #define der5(f,u,v,w,x,y) _d5[&f-_1u][dersub5(_c(u),_c(v),_c(w),_c(x),_c(y))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u : (#v[4]=='2' ? _2u_x-COMPONENTSDIM : (#v[4]=='3' ? _3u_xx-2*COMPONENTSDIM : (#v[4]=='4' ? _4s-3*COMPONENTSDIM : _5p-3*COMPONENTSDIM-SPACEDIM)))))~
  /*
 return _d5[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Space {#undef $Space~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Phase {#undef $Phase_@Space {$Space_$Space}~}
  @Phase {#undef $Phase_@Time {$Time}~}
  #undef der5~
  #undef _c~
  /*
 xxx
  /*
Left end BC{BC0
Boundary conditions at the left end of the unit interval
 _Bc0+
 int $(Compiler.RHSATTR) _Bc0(double *_u0, double *_u0p, double *_p, double *_t, double *_f) {~
  @Phase {#define $Phase(s) (*Phase(_u##s+&Phase))~}
  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_u##s##p+&Phase))~}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (*(_f+&Phase))~}
  /*
 return 0;~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  /*
  @Phase&{ BC_&Phase}
  @Phase { $Phase(0)}
  @Phase { $Phase_@Space {$Space}(0)}
  @Par   { $Par}
  @Time  { $Time}
  /*
BC0'
1st order derivatives of boundary conditions at the left side
 _Bc0_Der1
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM~typedef double _Bc0_d1T[COMPONENTSDIM][_n];~#undef _n~static _Bc0_d1T _Bc0_d1;~double * $(Compiler.RHSATTR) _Bc0_Der1(double *_1u0, double *_2u0p, double *_3p, double *_t) {~
  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc0_d1T*)(_Bc0_d1+&Phase))~}
  #define der1(f,u) _Bc0_d1[f][_c(u)]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  /*
 return _Bc0_d1[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der1~
  #undef _c~
  /*
 xxx
  /*
BC0''
2nd order derivatives of boundary conditions at the left side
 _Bc0_Der2
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM-1~typedef double _Bc0_d2T[COMPONENTSDIM][dersub2(_n,_n)+1];~#undef _n~static _Bc0_d2T _Bc0_d2;~double * $(Compiler.RHSATTR) _Bc0_Der2(double *_1u0, double *_2u0p, double *_3p, double *_t) {~
  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc0_d2T*)(_Bc0_d2+&Phase))~}
  #define der2(f,u,v) _Bc0_d2[f][dersub2(_c(u),_c(v))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  /*
 return _Bc0_d2[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der2~
  #undef _c~
  /*
 xxx
  /*
BC0'''
3rd order derivatives of boundary conditions at the left side
 _Bc0_Der3
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM-1~typedef double _Bc0_d3T[COMPONENTSDIM][dersub3(_n,_n,_n)+1];~#undef _n~static _Bc0_d3T _Bc0_d3;~double * $(Compiler.RHSATTR) _Bc0_Der3(double *_1u0, double *_2u0p, double *_3p, double *_t) {~
  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc0_d3T*)(_Bc0_d3+&Phase))~}
  #define der3(f,u,v,w) _Bc0_d3[f][dersub3(_c(u),_c(v),_c(w))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  /*
 return _Bc0_d3[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der3~
  #undef _c~
  /*
 xxx
  /*
BC0''''
4th order derivatives of boundary conditions at the left side
 _Bc0_Der4
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM-1~typedef double _Bc0_d4T[COMPONENTSDIM][dersub4(_n,_n,_n,_n)+1];~#undef _n~static _Bc0_d4T _Bc0_d4;~double * $(Compiler.RHSATTR) _Bc0_Der4(double *_1u0, double *_2u0p, double *_3p, double *_t) {~
!!  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
!!  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Phase {#define _1u0_&Phase (*Phase(_1u0+&Phase))~}
  @Phase {#define _2u0p_&Phase (*Phase(_2u0p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc0_d4T*)(_Bc0_d4+&Phase))~}
  #define der4(f,u,v,w,x) _Bc0_d4[f][dersub4(_c(u),_c(v),_c(w),_c(x))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  @Phase {#define $Phase(s) "_1u0_&Phase" /*== don't touch me*/~}
  @Phase {#define $Phase_@Space {$Space}(s) "_2u0p_&Phase" /*== don't touch me*/~}
  /*
 return _Bc0_d4[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der4~
  #undef _c~
  /*
 xxx
  /*
BC0'''''}
5th order derivatives of boundary conditions at the left side
 _Bc0_Der5
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM-1~typedef double _Bc0_d5T[COMPONENTSDIM][dersub5(_n,_n,_n,_n,_n)+1];~#undef _n~static _Bc0_d5T _Bc0_d5;~double * $(Compiler.RHSATTR) _Bc0_Der5(double *_1u0, double *_2u0p, double *_3p, double *_t) {~
!!  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
!!  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Phase {#define _1u0_&Phase (*Phase(_1u0+&Phase))~}
  @Phase {#define _2u0p_&Phase (*Phase(_2u0p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc0_d5T*)(_Bc0_d5+&Phase))~}
  #define der5(f,u,v,w,x,y) _Bc0_d5[f][dersub5(_c(u),_c(v),_c(w),_c(x),_c(y))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  @Phase {#define $Phase(s) "_1u0_&Phase" /*== don't touch me*/~}
  @Phase {#define $Phase_@Space {$Space}(s) "_2u0p_&Phase" /*== don't touch me*/~}
  /*
 return _Bc0_d5[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der5~
  #undef _c~
  /*
 xxx
  /*
Right end BC{BC1
Boundary conditions at the right side of the unit interval
 _Bc1+
 int $(Compiler.RHSATTR) _Bc1(double *_u1, double *_u1p, double *_p, double *_t, double *_f) {~
  @Phase {#define $Phase(s) (*Phase(_u##s+&Phase))~}
  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_u##s##p+&Phase))~}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (*(_f+&Phase))~}
  /*
 return 0;~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  /*
  @Phase&{ BC_&Phase}
  @Phase { $Phase(1)}
  @Phase { $Phase_@Space {$Space}(1)}
  @Par   { $Par}
  @Time  { $Time}
  /*
BC1'
1st order derivatives of boundary conditions at the right side
 _Bc1_Der1
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM~typedef double _Bc1_d1T[COMPONENTSDIM][_n];~#undef _n~static _Bc1_d1T _Bc1_d1;~double * $(Compiler.RHSATTR) _Bc1_Der1(double *_1u1, double *_2u1p, double *_3p, double *_t) {~
  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc1_d1T*)(_Bc1_d1+&Phase))~}
  #define der1(f,u) _Bc1_d1[f][_c(u)]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  /*
 return _Bc1_d1[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der1~
  #undef _c~
  /*
 xxx
  /*
BC1''
2nd order derivatives of boundary conditions at the right side
 _Bc1_Der2
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM-1~typedef double _Bc1_d2T[COMPONENTSDIM][dersub2(_n,_n)+1];~#undef _n~static _Bc1_d2T _Bc1_d2;~double * $(Compiler.RHSATTR) _Bc1_Der2(double *_1u1, double *_2u1p, double *_3p, double *_t) {~
  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc1_d2T*)(_Bc1_d2+&Phase))~}
  #define der2(f,u,v) _Bc1_d2[f][dersub2(_c(u),_c(v))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  /*
 return _Bc1_d2[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der2~
  #undef _c~
  /*
 xxx
  /*
BC1'''
3rd order derivatives of boundary conditions at the right side
 _Bc1_Der3
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM-1~typedef double _Bc1_d3T[COMPONENTSDIM][dersub3(_n,_n,_n)+1];~#undef _n~static _Bc1_d3T _Bc1_d3;~double * $(Compiler.RHSATTR) _Bc1_Der3(double *_1u1, double *_2u1p, double *_3p, double *_t) {~
  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc1_d3T*)(_Bc1_d3+&Phase))~}
  #define der3(f,u,v,w) _Bc1_d3[f-][dersub3(_c(u),_c(v),_c(w))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  /*
 return _Bc1_d3[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der3~
  #undef _c~
  /*
 xxx
  /*
BC1''''
4th order derivatives of boundary conditions at the right side
 _Bc1_Der4
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM-1~typedef double _Bc1_d4T[COMPONENTSDIM][dersub4(_n,_n,_n,_n)+1];~#undef _n~static _Bc1_d4T _Bc1_d4;~double * $(Compiler.RHSATTR) _Bc1_Der4(double *_1u1, double *_2u1p, double *_3p, double *_t) {~
!!  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
!!  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Phase {#define _1u1_&Phase (*Phase(_1u1+&Phase))~}
  @Phase {#define _2u1p_&Phase (*Phase(_2u1p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc1_d4T*)(_Bc1_d4+&Phase))~}
  #define der4(f,u,v,w,x) _Bc1_d4[f][dersub4(_c(u),_c(v),_c(w),_c(x))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  @Phase {#define $Phase(s) "_1u1_&Phase" /*== don't touch me*/~}
  @Phase {#define $Phase_@Space {$Space}(s) "_2u1p_&Phase" /*== don't touch me*/~}
  /*
 return _Bc1_d4[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der4~
  #undef _c~
  /*
 xxx
  /*
BC1'''''}
5th order derivatives of boundary conditions at the right side
 _Bc1_Der5
 #define _n 2*COMPONENTSDIM+PARAMETERSDIM-1~typedef double _Bc1_d5T[COMPONENTSDIM][dersub5(_n,_n,_n,_n,_n)+1];~#undef _n~static _Bc1_d5T _Bc1_d5;~double * $(Compiler.RHSATTR) _Bc1_Der5(double *_1u1, double *_2u1p, double *_3p, double *_t) {~
!!  @Phase {#define $Phase(s) (*Phase(_1u##s+&Phase))~}
!!  @Phase {#define $Phase_@Space {$Space}(s) (*Phase(_2u##s##p+&Phase))~}
  @Phase {#define _1u1_&Phase (*Phase(_1u1+&Phase))~}
  @Phase {#define _2u1p_&Phase (*Phase(_2u1p+&Phase))~}
  @Par   {#define $Par (*Par(_3p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase&{#define BC_&Phase (**(_Bc1_d5T*)(_Bc1_d5+&Phase))~}
  #define der5(f,u,v,w,x,y) _Bc1_d5[f][dersub5(_c(u),_c(v),_c(w),_c(x),_c(y))]~
  #define _c(v) (&v-(#v[4]=='1' ? _1u0 : (#v[4]=='2' ? _2u0p-COMPONENTSDIM : _3p-2*COMPONENTSDIM)))~
  @Phase {#define $Phase(s) "_1u1_&Phase" /*== don't touch me*/~}
  @Phase {#define $Phase_@Space {$Space}(s) "_2u1p_&Phase" /*== don't touch me*/~}
  /*
 return _Bc1_d5[0];~}~
  @Phase {#undef $Phase~}
  @Phase {#undef $Phase_@Space {$Space}~}
  @Par   {#undef $Par~}
  @Time  {#undef $Time~}
  @Phase&{#undef BC_&Phase~}
  #undef der5~
  #undef _c~
  /*
 xxx
  /*
&derivatives
  1
&visiblenames
2 *Phase Components~ ~>= 1
3 Par   Parameters~ ~>= 0
1 Time  Time~ ~== 1
10 *Space Space coordinate~ ~== 1
0
4 Eigv Eigenvalues~@Phase& { Re_&Phase} @Phase& { Im_&Phase}~
6 TestEp- Test functions (EP)~Branching Fold Hopf~
7 TestH- Test functions (H)~Generalized_Hopf Bogdanov_Takens Zero_Hopf Double_Hopf~
8 TestLp- Test functions (LP)~Cusp Bogdanov_Takens Hopf~
9 TestLc- Test functions (LC)~CycleFold Period_doubling Neimark_Sacker~
5 UserFunc- User functions~ ~
11 StdFunc- Standard functions ~L2_Norm @Phase { Min_$Phase} @Phase { Max_$Phase}~
! Special classes (with 255>code>200) MUST follow those with code<200
254 SPECCL_AUS d1~ ~ ! ActiveUserSing
253 SPECCL_DIM d2~ ~ ! G-dim of point
&help
cont_pd1.hlp
&windows
!