weiye
/
Droplet_LTORC
Watch 1
Star 0
Fork 0
Code Issues 0 Pull Requests 0 Releases 0 Wiki Activity
Droplet Lagrangian Transient One-dimensional Reacting Code Implementation of both liquid and gas phase governing equations.
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
4 Commits
1 Branch
1.7MB
Tree: 0aba332165
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from '0aba332165'
${ noResults }
Droplet_LTORC/macros.h

78 lines
2.7KB
Raw Blame History 

  1. #ifndef MACRO_DEF

  2. #define MACRO_DEF

  3. //#define SUNDIALS_DOUBLE_PRECISION 1

  4. //#define SUNDIALS_SINGLE_PRECISION 1

  5. 

  6. #define ZERO  RCONST(0.0)

  7. #define HALF  RCONST(0.5)

  8. #define ONE   RCONST(1.0)

  9. #define TWO   RCONST(2.0)

  10. #define THREE RCONST(3.0)

  11. #define FOUR  RCONST(4.0)

  12. #define TEN   RCONST(10.0)

  13. 

  14. /* In order to keep begin the index numbers from 1 instead of 0, we define

  15.  * macros here. Also, we define macros to ease referencing various variables in

  16.  * the sundials nvector.

  17.  */

  18. #define psi(i)   psidata[i-1]

  19. 

  20. #define T(i)    ydata[((i-1)*data->nvar)+data->nt]

  21. #define Y(i,k)  ydata[((i-1)*data->nvar)+data->ny+k-1]

  22. #define R(i)    ydata[((i-1)*data->nvar)+data->nr]

  23. #define P(i)    ydata[((i-1)*data->nvar)+data->np]

  24. #define Mdot(i) ydata[((i-1)*data->nvar)+data->nm]

  25. 

  26. #define Tdot(i)    ydotdata[((i-1)*data->nvar)+data->nt]

  27. #define Ydot(i,k)  ydotdata[((i-1)*data->nvar)+data->ny+k-1]

  28. #define Rdot(i)    ydotdata[((i-1)*data->nvar)+data->nr]

  29. #define Pdot(i)    ydotdata[((i-1)*data->nvar)+data->np]

  30. #define Mdotdot(i) ydotdata[((i-1)*data->nvar)+data->nm]

  31. 

  32. #define Tres(i)    resdata[((i-1)*data->nvar)+data->nt]

  33. #define Yres(i,k)  resdata[((i-1)*data->nvar)+data->ny+k-1]

  34. #define Rres(i)    resdata[((i-1)*data->nvar)+data->nr]

  35. #define Pres(i)    resdata[((i-1)*data->nvar)+data->np]

  36. #define Mdotres(i) resdata[((i-1)*data->nvar)+data->nm]

  37. 

  38. #define Tid(i)    iddata[((i-1)*data->nvar)+data->nt]

  39. #define Yid(i,k)  iddata[((i-1)*data->nvar)+data->ny+k-1]

  40. #define Rid(i)    iddata[((i-1)*data->nvar)+data->nr]

  41. #define Pid(i)    iddata[((i-1)*data->nvar)+data->np]

  42. #define Mdotid(i) iddata[((i-1)*data->nvar)+data->nm]

  43. 

  44. #define Yav(i) Yav[i-1]

  45. #define YAvg(i) YAvg[i-1]

  46. #define YVmhalf(i) YVmhalf[i-1]

  47. #define YVphalf(i) YVphalf[i-1]

  48. #define X(i) X[i-1]

  49. #define Xp(i) Xp[i-1]

  50. #define Xgradhalf(i) Xgradhalf[i-1]

  51. #define XLeft(i) XLeft[i-1]

  52. #define XRight(i) XRight[i-1]

  53. #define gradX(i) gradX[i-1]

  54. #define wdot(i) wdot[i-1]

  55. #define enthalpy(i) enthalpy[i-1]

  56. #define energy(i) energy[i-1]

  57. #define Cp(i) Cp[i-1]

  58. 

  59. #define atolT(i)   atolvdata[((i-1)*data->nvar)+data->nt]

  60. #define atolY(i,k)   atolvdata[((i-1)*data->nvar)+data->ny+k-1]

  61. #define atolR(i)   atolvdata[((i-1)*data->nvar)+data->nr]

  62. #define atolP(i)   atolvdata[((i-1)*data->nvar)+data->np]

  63. #define atolMdot(i)   atolvdata[((i-1)*data->nvar)+data->nm]

  64. 

  65. 

  66. #define constraintsY(i,k)   constraintsdata[((i-1)*data->nvar)+data->ny+k-1]

  67. #define constraintsR(i)   constraintsdata[((i-1)*data->nvar)+data->nr]

  68. 

  69. /*Following marcos are defined to calculate the characteristic time-scale */

  70. #define wdot_mole(i)    wdot_mole[i-1]

  71. #define wdot_mass(i)    wdot_mass[i-1]

  72. #define MW(i)         MW[i-1]

  73. #define time_scale(i,k)    time_scale[(i-1)*data->nsp+k-1]

  74. #define concentra(i)     concentra[i-1]

  75. 

  76. 

  77. #endif