12345678910111213141516171819202122232425262728293031323334 |
- //---------------------------------------------------------------------
- //---------------------------------------------------------------------
- //
- // work_lhs.h
- //
- //---------------------------------------------------------------------
- //---------------------------------------------------------------------
- #ifndef __WORK_LHS_H
- #define __WORK_LHS_H
- #define fjac(m,n,i) fjac[(m-1)+5*((n-1)+5*(i+2))]
- #define njac(m,n,i) njac[(m-1)+5*((n-1)+5*(i+2))]
- #define lhsa(m,n,i) lhsa[(m-1)+5*((n-1)+5*(i+1))]
- #define lhsb(m,n,i) lhsb[(m-1)+5*((n-1)+5*(i+1))]
- #ifdef G_MAIN
- double fjac[5*5*(MAX_CELL_DIM+4)],
- njac[5*5*(MAX_CELL_DIM+4)],
- lhsa[5*5*(MAX_CELL_DIM+2)],
- lhsb[5*5*(MAX_CELL_DIM+2)],
- tmp1, tmp2, tmp3;
- // common /work_lhs/ fjac, njac, lhsa, lhsb, tmp1, tmp2, tmp3;
- #else
- extern double fjac[5*5*(MAX_CELL_DIM+4)],
- njac[5*5*(MAX_CELL_DIM+4)],
- lhsa[5*5*(MAX_CELL_DIM+2)],
- lhsb[5*5*(MAX_CELL_DIM+2)],
- tmp1, tmp2, tmp3;
- #endif /*G_MAIN*/
- #ifdef _OPENMP
- #pragma omp threadprivate (fjac, njac, lhsa, lhsb, tmp1, tmp2, tmp3)
- #endif
- #endif
|