/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D E   L A   D Y N A M I Q U E   D E   V E R H U L S T   T R I D I M E N S I O N N E L L E  :         */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrk/verhulst.31$I' :                                                                                           */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, 20190213114145).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D E   L A   D Y N A M I Q U E   D E   V E R H U L S T   T R I D I M E N S I O N N E L L E  :         */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Definition :                                                                                                               */
/*                                                                                                                                   */
/*                    Elle est definie par le systeme                                                                                */
/*                  d'equations differentielles :                                                                                    */
/*                                                                                                                                   */
/*                                       dx                                                                                          */
/*                                      ---- = R  .x.(A  -x) + R  .y.(A  -y) + R  .z.(A  -z)                                         */
/*                                       dt     xx     xx       xy     xy       xz     xz                                            */
/*                                                                                                                                   */
/*                                       dy                                                                                          */
/*                                      ---- = R  .x.(A  -x) + R  .y.(A  -y) + R  .z.(A  -z)                                         */
/*                                       dt     yx     yx       yy     yy       yz     yz                                            */
/*                                                                                                                                   */
/*                                       dz                                                                                          */
/*                                      ---- = R  .x.(A  -x) + R  .y.(A  -y) + R  .z.(A  -z)                                         */
/*                                       dt     zx     zx       zy     zy       zz     zz                                            */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   CRxx                                                                                                                          \
                    FTROIS
#define   CRxy                                                                                                                          \
                    FZERO
#define   CRxz                                                                                                                          \
                    FZERO
#define   CAxx                                                                                                                          \
                    FU
#define   CAxy                                                                                                                          \
                    FU
#define   CAxz                                                                                                                          \
                    FU
DEFV(Local,DEFV(Float,INIT(Rxx,CRxx)));
DEFV(Local,DEFV(Float,INIT(Rxy,CRxy)));
DEFV(Local,DEFV(Float,INIT(Rxz,CRxz)));
DEFV(Local,DEFV(Float,INIT(Axx,CAxx)));
DEFV(Local,DEFV(Float,INIT(Axy,CAxy)));
DEFV(Local,DEFV(Float,INIT(Axz,CAxz)));
                                        /* Definition des coefficients de 'dx'.                                                      */

#define   CRyx                                                                                                                          \
                    FZERO
#define   CRyy                                                                                                                          \
                    FTROIS
#define   CRyz                                                                                                                          \
                    FZERO
#define   CAyx                                                                                                                          \
                    FU
#define   CAyy                                                                                                                          \
                    FU
#define   CAyz                                                                                                                          \
                    FU
DEFV(Local,DEFV(Float,INIT(Ryx,CRyx)));
DEFV(Local,DEFV(Float,INIT(Ryy,CRyy)));
DEFV(Local,DEFV(Float,INIT(Ryz,CRyz)));
DEFV(Local,DEFV(Float,INIT(Ayx,CAyx)));
DEFV(Local,DEFV(Float,INIT(Ayy,CAyy)));
DEFV(Local,DEFV(Float,INIT(Ayz,CAyz)));
                                        /* Definition des coefficients de 'dy'.                                                      */

#define   CRzx                                                                                                                          \
                    FZERO
#define   CRzy                                                                                                                          \
                    FZERO
#define   CRzz                                                                                                                          \
                    FTROIS
#define   CAzx                                                                                                                          \
                    FU
#define   CAzy                                                                                                                          \
                    FU
#define   CAzz                                                                                                                          \
                    FU
DEFV(Local,DEFV(Float,INIT(Rzx,CRzx)));
DEFV(Local,DEFV(Float,INIT(Rzy,CRzy)));
DEFV(Local,DEFV(Float,INIT(Rzz,CRzz)));
DEFV(Local,DEFV(Float,INIT(Azx,CAzx)));
DEFV(Local,DEFV(Float,INIT(Azy,CAzy)));
DEFV(Local,DEFV(Float,INIT(Azz,CAzz)));
                                        /* Definition des coefficients de 'dz'.                                                      */

#define   CX0                                                                                                                           \
                    FDU
#define   CY0                                                                                                                           \
                    FDU
#define   CZ0                                                                                                                           \
                    FDU

#define   DCT                                                                                                                           \
                    FU
DEFV(Local,DEFV(Float,INIT(dct,DCT)));
                                        /* Definition de 'dt'.                                                                       */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D E S   T R O I S   F O N C T I O N S   ' F '  :                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   Fx(cx,cy,cz,t)                                                                                                                \
                    ADD3(MUL3(Rxx,cx,SOUS(Axx,cx))                                                                                      \
                        ,MUL3(Rxy,cy,SOUS(Axy,cy))                                                                                      \
                        ,MUL3(Rxz,cz,SOUS(Axz,cz))                                                                                      \
                         )                                                                                                              \
                                        /* Definition de la fonction F (x,y,z).                                                      */ \
                                        /*                            x                                                              */
#define   Fy(cx,cy,cz,t)                                                                                                                \
                    ADD3(MUL3(Ryx,cx,SOUS(Ayx,cx))                                                                                      \
                        ,MUL3(Ryy,cy,SOUS(Ayy,cy))                                                                                      \
                        ,MUL3(Ryz,cz,SOUS(Ayz,cz))                                                                                      \
                         )                                                                                                              \
                                        /* Definition de la fonction F (x,y,z).                                                      */ \
                                        /*                            y                                                              */
#define   Fz(cx,cy,cz,t)                                                                                                                \
                    ADD3(MUL3(Rzx,cx,SOUS(Azx,cx))                                                                                      \
                        ,MUL3(Rzy,cy,SOUS(Azy,cy))                                                                                      \
                        ,MUL3(Rzz,cz,SOUS(Azz,cz))                                                                                      \
                         )                                                                                                              \
                                        /* Definition de la fonction F (x,y,z).                                                      */ \
                                        /*                            z                                                              */