From victor@cwi.nl Fri Jan 16 14:38:38 1998
From: victor@cwi.nl
Date: Fri, 16 Jan 1998 14:38:02 +0100 (MET)
To: Yu.A.KUZNETSOV@umpa.ens-lyon.fr
Subject: dde_g.con
Content-Length: 5759

&graph
1 O Point|Orbit
  -> O   dd__o point_or dde_or start_or integrpar
&compclass !the way curves of the given type are computed
O   1 !trajectories
!EP  2 !implicit functions
&menu
Initial point{
 O|Point
}
Curve{
 O|Orbit
}
&deftypes
! CurveType PointType NewCurveType
      O        O           O
&stdlinalg
&&$(Paths.KERNEL)stdalg$(Extensions.DESC)
&utilities
_util(CheckParameters)
=(ErrMessage)
=(ActivateTestFuns)
=(IndProcFunc)
&functions
xxx
RHS
 _Rhs+
 int $(Compiler.RHSATTR) _Rhs(double *_x, double *_p, double *_t, double *_f) {~
  @Phase {#define $Phase (*Phase(_x+&Phase))~}
  @Delay {@Phase {#define $Phase_$Delay (*Phase(_x+(@(%Phase)*(&Delay+1)+&Phase)))~}}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Time  {#define $Time _t[0]~}
  @Phase {#define $Phase_ (*Phase(_f+&Phase))~}
  /*
 return 0;~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Phase {#undef $Phase_~}
  @Delay {@Phase {#undef $Phase_$Delay~}}
  /*
  @Phase { $Phase_}
  @Phase { $Phase}
  @Delay {@Phase { $Phase_$Delay}}
  @Par   { $Par}
  @Time  { $Time}
  /*
Derivatives of RHS{1st order
1st order derivatives of RHS
 _Rhs_Der1
 #define _n COORDINATESDIM*(DELAYSDIM+1)+PARAMETERSDIM~typedef double _d1T[COORDINATESDIM][_n];~#undef _n~static _d1T _d1;~double * $(Compiler.RHSATTR) _Rhs_Der1(double *_x, double *_p, double *_t) {~
  @Time  {#define $Time _t[0]~}
  @Phase {#define $Phase (*Phase(_x+&Phase))~}
  @Delay {@Phase {#define $Phase_$Delay (*Phase(_x+(@(%Phase)*(&Delay+1)+&Phase)))~}}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Phase {#define $Phase_ (*Phase*(_d1T*)(_d1+&Phase))~}
  #define der1(f,u) _d1[&f-_x][_c(u)]~
  #define _c(v) (&v-(#v[4]=='p' ? _p-COORDINATESDIM*(DELAYSDIM+1) : _x))~
  /*
 return _d1[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Phase {#undef $Phase_~}
  #undef der1~
  #undef _c~
  /*
 xxx
  /*
2nd order
2nd order derivatives of RHS
 _Rhs_Der2
 #define _n COORDINATESDIM*(DELAYSDIM+1)+PARAMETERSDIM-1~typedef double _d2T[COORDINATESDIM][dersub2(_n,_n)+1];~static _d2T _d2;~#undef _n~double * $(Compiler.RHSATTR) _Rhs_Der2(double *_x, double *_p, double *_t) {~
  @Time  {#define $Time _t[0]~}
  @Phase {#define $Phase (*Phase(_x+&Phase))~}
  @Delay {@Phase {#define $Phase_$Delay (*Phase(_x+(@(%Phase)*(&Delay+1)+&Phase)))~}}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Phase {#define $Phase_ (*Phase*(_d2T*)(_d2+&Phase))~}
  #define der2(f,u,v) _d2[&f-_x][dersub2(_c(u),_c(v))]~
  #define _c(v) (&v-(#v[4]=='p' ? _p-COORDINATESDIM*(DELAYSDIM+1) : _x))~
  /*
 return _d2[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Phase {#undef $Phase_~}
  #undef der2~
  #undef _c~
  /*
 xxx
  /*
3rd order
3rd order derivatives of RHS
 _Rhs_Der3
 #define _n COORDINATESDIM*(DELAYSDIM+1)+PARAMETERSDIM-1~typedef double _d3T[COORDINATESDIM][dersub3(_n,_n,_n)+1];~static _d3T _d3;~#undef _n~double * $(Compiler.RHSATTR) _Rhs_Der3(double *_x, double *_p, double *_t) {~
  @Time  {#define $Time _t[0]~}
  @Phase {#define $Phase (*Phase(_x+&Phase))~}
  @Delay {@Phase {#define $Phase_$Delay (*Phase(_x+(@(%Phase)*(&Delay+1)+&Phase)))~}}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Phase {#define $Phase_ (*Phase*(_d3T*)(_d3+&Phase))~}
  #define der3(f,u,v,w) _d3[&f-_x][dersub3(_c(u),_c(v),_c(w))]~
  #define _c(v) (&v-(#v[4]=='p' ? _p-COORDINATESDIM*(DELAYSDIM+1) : _x))~
  /*
 return _d3[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Phase {#undef $Phase_~}
  #undef der3~
  #undef _c~
  /*
 xxx
  /*
4th order
4th order derivatives of RHS
 _Rhs_Der4
 #define _n COORDINATESDIM*(DELAYSDIM+1)+PARAMETERSDIM-1~typedef double _d4T[COORDINATESDIM][dersub4(_n,_n,_n,_n)+1];~static _d4T _d4;~#undef _n~double * $(Compiler.RHSATTR) _Rhs_Der4(double *_x, double *_p, double *_t) {~
  @Time  {#define $Time _t[0]~}
  @Phase {#define $Phase (*Phase(_x+&Phase))~}
  @Delay {@Phase {#define $Phase_$Delay (*Phase(_x+(@(%Phase)*(&Delay+1)+&Phase)))~}}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Phase {#define $Phase_ (*Phase*(_d4T*)(_d4+&Phase))~}
  #define der4(f,u,v,w,x) _d4[&f-_x][dersub4(_c(u),_c(v),_c(w),_c(x))]~
  #define _c(v) (&v-(#v[4]=='p' ? _p-COORDINATESDIM*(DELAYSDIM+1) : _x))~
  /*
 return _d4[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Phase {#undef $Phase_~}
  #undef der4~
  #undef _c~
  /*
 xxx
  /*
5th order}
5th order derivatives of RHS
 _Rhs_Der5
 #define _n COORDINATESDIM*(DELAYSDIM+1)+PARAMETERSDIM-1~typedef double _d5T[COORDINATESDIM][dersub5(_n,_n,_n,_n,_n)+1];~static _d5T _d5;~#undef _n~double * $(Compiler.RHSATTR) _Rhs_Der5(double *_x, double *_p, double *_t) {~
  @Time  {#define $Time _t[0]~}
  @Phase {#define $Phase (*Phase(_x+&Phase))~}
  @Delay {@Phase {#define $Phase_$Delay (*Phase(_x+(@(%Phase)*(&Delay+1)+&Phase)))~}}
  @Par   {#define $Par (*Par(_p+&Par))~}
  @Phase {#define $Phase_ (*Phase*(_d5T*)(_d5+&Phase))~}
  #define der5(f,u,v,w,x,y) _d5[&f-_x][dersub5(_c(u),_c(v),_c(w),_c(x),_c(y))]~
  #define _c(v) (&v-(#v[4]=='p' ? _p-COORDINATESDIM*(DELAYSDIM+1) : _x))~
  /*
 return _d5[0];~}~
  @Time  {#undef $Time~}
  @Phase {#undef $Phase~}
  @Par   {#undef $Par~}
  @Phase {#undef $Phase_~}
  #undef der5~
  #undef _c~
  /*
 xxx
  /*
&derivatives
  3
&visiblenames
2 Phase Coordinates~ ~> 0
6 Delay Delays~ ~>= 0
3 Par   Parameters~ ~>= 0
1 Time  Time~ ~== 1
0
4 Eigv- Eigenvalues~@Phase& { Re_&Phase} @Phase& { Im_&Phase}~
13 StdFunc- Standard functions ~L2_Norm @Phase { Min_$Phase} @Phase { Max_$Phase}~
5 UserFunc- User functions~ ~
! 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_dde.hlp
&windows
!