src/coordinates/coordinates.cpp
//----------------------------------------------------------------------------------------
// FaceXArea functions: compute area of face with normal in X-dir as vector
FaceXArea functions: compute area of face with normal in X-dir as vector
FaceXArea functions : compute area of face with normal in X - dir as vector
void Coordinates::Face1Area(const int k, const int j, const int il, const int iu,
  AthenaArray<Real> &area)
{
#pragma simd
  for (int i=il; i<=iu; ++i){
    Real& area_i = area(i);
    area_i = dx2f(j)*dx3f(k);
  }
  return;
}
void Coordinates :: Face1Area ( const int k , const int j , const int il , const int iu , AthenaArray < Real > & area ) { #pragma  simd for ( int i = il ; i <= iu ; ++ i ) { Real & area_i = area ( i ) ; area_i = dx2f ( j ) * dx3f ( k ) ; } return ; }
FaceXArea functions: compute area of face with normal in X-dir as vector
FaceXArea functions : compute area of face with normal in X - dir as vector
// area1 = dy dz