struct.h

Go to the documentation of this file.

 
#define RESET   "\033[0m"
#define BLACK   "\033[30m"      /* Black */
#define RED     "\033[31m"      /* Red */
#define GREEN   "\033[32m"      /* Green */
#define YELLOW  "\033[33m"      /* Yellow */
#define BLUE    "\033[34m"      /* Blue */
#define MAGENTA "\033[35m"      /* Magenta */
#define CYAN    "\033[36m"      /* Cyan */
#define WHITE   "\033[37m"      /* White */
#define BOLDBLACK   "\033[1m"      /* Bold Black */
#define BOLDRED     "\033[1m\033[31m"      /* Bold Red */
#define BOLDGREEN   "\033[1m\033[32m"      /* Bold Green */
#define BOLDYELLOW  "\033[1m\033[33m"      /* Bold Yellow */
#define BOLDBLUE    "\033[1m\033[34m"      /* Bold Blue */
#define BOLDMAGENTA "\033[1m\033[35m"      /* Bold Magenta */
#define BOLDCYAN    "\033[1m\033[36m"      /* Bold Cyan */
#define BOLDWHITE   "\033[1m\033[37m"      /* Bold White */
 
 
/*-------------------------------------*/
/* Structures that store the variables */
/*-------------------------------------*/
 
 
 
typedef struct {
  PetscMPIInt size;
  PetscMPIInt rank;
  int         xs;
  int         ys;
  int         zs;
  int         xm;
  int         ym;
  int         zm;
  int         gxs;
  int         gys;
  int         gzs;
  int         gxm;
  int         gym;
  int         gzm;
} Proc;
 
typedef struct {
  DM        da;
  int       nx;
  int       ny;
  int       nz;
  double    dx;
  double    dy;
  double    dz;
  double    xmin;
  double    xmax;
  double    ymin;
  double    ymax;
  double    zmin;
  double    zmax;
  PetscBool xperiodic;
  PetscBool yperiodic;
  PetscBool zperiodic;
  int       ip;
  int       jp;
  int       kp;
  double    epsilon;
} Grid;
 
typedef struct {
  int       ni;
  int       kk,carac;
  int       ischeme;
  int       cscheme;
  int       fluid; /* 1 or 2 for viscous implit */
  double    mu1,rho1,mu2,rho2,dt,dtelast,lambda,cfl,time,timeac,timeac2,g,volume,st;
  double    Cc;
  double    Cv;
  double     Cg;
  double    Cs;
  double    E,rhoe,mue;
  double    hsurf;
  PetscBool force_exa;
  PetscBool force_ext1;
  PetscBool force_ext2;
  PetscBool bi_fluid;
  PetscBool redist;
  int       test;
} Simu;
 
typedef struct {
  int  nb; 
  int  nbe;
  int  ib; 
  int  *no;
  int  *ns; 
  int  *ne; 
  int  nexp,nexp2,nexp3;
  int  *choice;
  int  *file_format;
  int  test,*nst;
  char OBJ_FILE[30][PETSC_MAX_PATH_LEN];
  
  double *xtail,*ytail,*ztail,*xhead,*yhead,*zhead,xfood,yfood,zfood,thetafood;  
  double *St,*St_swim; /* the surface of the body (in 2D) */
 
  struct {
    double **u,**v;
    double ***Ub;
    double ***Vb;
    double ***Wb;
    double Uin;
    double Vin;
    double Win;
    double Vref;
    double *ubar;
    double *vbar;
    double *wbar;
    double *xbar;
    double *ybar;
    double *zbar;
    double *cp;
    double **Omega;
    double **Omegat;
    double **Theta;
    double **Theta_def;
    double **Thetai;
    double rpm,theta;
    double a,f;
  }mouv;
 
  struct {
    int  *is, *js, *ks, *ie, *je, *ke;
    int  *is1,*js1,*ks1,*ie1,*je1,*ke1;
    int  *is2,*js2,*ks2,*ie2,*je2,*ke2;
    int  *ist,*jst,*kst,*iet,*jet,*ket;
    double *Cx;
    double *Cy;
    double *Cz;
    double *Cxt;
    double *Cyt;
    double *Czt;
    double ***Cp,***Cpt;
    double ***Cs,***Cst;
    double ***C_useful,***Ct_useful,***C_useless,***Ct_useless,***C_other,***Ct_other,***C_total,***Ct_total,*Cn,Mn_exa,Mn_ext,*eta,*etas;
    double *Mx;
    double *My;
    double *Mz;
    double *Mxt;
    double *Myt;
    double *Mzt;
    double *Power,*Powert,*Power_useful,*Powert_useful,*Power_useless,*Powert_useless,*Power_other,*Powert_other;
    double *W_total,*W_useful,*W_useless,*W_other,**P_total,**P_useful,**P_useless,**P_other;
    double P,Pm,Pt,Pc,Pct,Pv,Pvt,Pb,Pbt,Pext,Pextt,Ec,Ect,Ecm1,Ecm1t;
    double **Fx;
    double **Fy;
    double **Fz;
    double **Fxt;
    double **Fyt;
    double **Fzt;
    double ***Fxx,***Fxxt,***Fyy,***Fyyt,***Fzz,***Fzzt;
    double ***dist,***distt,***distm1;
    int ***distk,***distkt;
  }coef;
 
  struct {
    int  nbarres[100];
    int  *nb_normal;
    int *nb_triangle;
    int *nb_quad;
    int *nb_vertex;
    int *nb_texture;
    int ***ind_vertex_q,***ind_vertex_t,***ind_texture_q,***ind_texture_t,***ind_normal_q,***ind_normal_t;
    double **vx,**vy,**vz,**nx,**ny,**nz,**tx,**ty,**tz;
    double *xi;
    double *yi;
    double *zi;
    double ***x;
    double ***y;
    double ***z;
    double ***xm1;
    double ***ym1;
    double ***zm1;
    double ***xs,***ys,***zs;
    double ***x0,***y0,***z0,***x0m1,***y0m1,***z0m1;
    double ***x1,***y1,***z1;
    double **xb;
    double **yb;
    double **zb;
    double **xb_m1,**yb_m1,**zb_m1,**ub,**vb,**wb;
    double **xbs,**ybs,**zbs;
    
    double ***Jinv;
    
    double ***Jxyz;
    
    double ***Jxyzm1;
    double ***n;
    double *length;
    double *V;
    double *Vt;
    double *Vs;
    double *S;
    double *St;
    double *rho;
    double *Amax;
    double *Amin;
    double *xg;
    double *yg;
    double *zg;
    double *xgd,*ygd,*zgd,*xgs,*ygs,*zgs,*xg1,*yg1,*zg1,*xg2,*yg2,*zg2;
    double *xmin,*xmax,*ymin,*ymax,*zmin,*zmax;
    
    double **theta;
    double **theta_m1;
    double **theta_m2;
    double re,ri,rc,yt,ym,yc,yl,r1,r2,m1,m2,rh,h;
 
  }geom;
 
  struct {
    int  n[100]; // nombre de barres
    double **x,**y,**xm1,**ym1,**xm2,**ym2;  //
    double **a,**b,**e;
    double **xc,**yc,**l;
    double **Wx,**Wy,**Wxt,**Wyt;
    double **m,**S,**In;
    double **phi,**phim1,**phim2,**dphi,**d2phi,**curv,**Aibphi;
    double **ce,**se,**cem1,**cem2,**sem1,**sem2,**d2ce,**d2se;
    double *E,*mu;
    double ***A,***A2,***Ap,***Hc,***Hs,***Gc,***Gs,***Jn,***Aphi;
    double ***GcAp,***GcApHc,***GsAp,***GsApHs,***GsApHc,***GcApHs;
    double **Mm,**Mmt,**bphi;
    double distance,dist;
    double *Cx,*Cy;
  }elast;
 
  //double *regul2;
 
} Obs;
 
typedef struct {
  Vec U;
  Vec lU;
  Vec Uf;
  Vec lUf;
  Vec V;
  Vec lV;
  Vec Vf;
  Vec lVf;
  Vec W;
  Vec lW;
  Vec Wf;
  Vec lWf;
  Vec lP;
  Vec ldPdy;
  Vec lLS; 
  Vec lLSind;
  Vec PHI;
  Vec lPHI; 
  Vec lOx; 
  Vec lOy;
  Vec lOz; 
  Vec lOn; 
  Vec lQ; 
  Vec lPSI;
  Vec lD;
  Vec lH;
  Vec DIV;
  Vec DU;
  Vec lDU;
  Vec Yx;
  Vec lYx;
  Vec Yy;
  Vec lYy;
  Vec Yz;
  Vec lYz;
  Vec lnx;
  Vec lny;
  Vec lnz;
  Vec a,la;
  Vec an,lan;
  Vec as,las;
} Vectors;
 
typedef struct {
  Vectors n;
  Vectors np1;
  Vectors nm1;
  Vectors m;
  Vectors e;
  Vectors t;
  Vectors i[10],k[10],l[10],im1[10],it[10],el[10],elm1[10],elt[10],mat[6][3][3];
} Fields;
 
typedef struct {
  int ***LSind;
  double ***U;
  double ***V;
  double ***W;
  double ***P;
  double ***dPdy;
  double ***Uf;
  double ***Vf;
  double ***Wf;
  double ***PHI;
  double ***LS;
  double ***Ox;
  double ***Oy;
  double ***Oz;
  double ***On;
  double ***Q;
  double ***PSI;
  double ***D;
  double ***H;
  double ***Yx;
  double ***Yy;
  double ***Yz;
  double ***nx;
  double ***ny;
  double ***nz;
  double ***a,***as,***an;
  double ***DU;
  double      ***DIV;
} Pointers;
 
typedef struct {
  Pointers n;
  Pointers np1;
  Pointers nm1;
  Pointers m;
  Pointers e;
  Pointers t,k[10],l[10],i[10],im1[10],it[10],cl,el[10],elm1[10],elt[10],mat[6][3][3];
} Arrays;
 
 
 
typedef struct {
  // Matrices
  Mat Ap;
  Mat Au;   
  // Right Hand Sides
  struct {
    Vec Bp;
    Vec Bu;
    Vec Bv;
    Vec Bw;   
  }vec;
  struct {
    double ***Bp;
    double ***Bu;
    double ***Bv;
    double ***Bw;
  }arr;
  // KSP Solvers
  KSP kspp;
  KSP kspu;
  int itp;
  int itu;
  int itv;
  int itw;
  MatNullSpace   nullsp;
} Linear;
 
 
 
 
/*-------------------------------------*/