Derivatives of RHS{1st order 1st order derivatives of RHS _Rhs_Der1 typedef double _d1T[COORDINATESDIM][COORDINATESDIM+PARAMETERSDIM];~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))~} @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 : _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+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))~} @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 : _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+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))~} @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 : _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+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))~} @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 : _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+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))~} @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 : _x))~ /* return _d5[0];~}~ @Time {#undef $Time~} @Phase {#undef $Phase~} @Par {#undef $Par~} @Phase {#undef $Phase_~} #undef der5~ #undef _c~ /* xxx /*