cs_flag.h

Go to the documentation of this file.

#ifndef __CS_FLAG_H__
#define __CS_FLAG_H__

/*============================================================================
 * Manage the definition/setting of a computation
 *============================================================================*/

/*
  This file is part of code_saturne, a general-purpose CFD tool.

  Copyright (C) 1998-2022 EDF S.A.

  This program is free software; you can redistribute it and/or modify it under
  the terms of the GNU General Public License as published by the Free Software
  Foundation; either version 2 of the License, or (at your option) any later
  version.

  This program is distributed in the hope that it will be useful, but WITHOUT
  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
  details.

  You should have received a copy of the GNU General Public License along with
  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
  Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/

/*----------------------------------------------------------------------------*/

/*----------------------------------------------------------------------------
 *  Local headers
 *----------------------------------------------------------------------------*/

#include "cs_defs.h"

/*----------------------------------------------------------------------------*/

BEGIN_C_DECLS

/*============================================================================
 * Macro definitions
 *============================================================================*/

#define CS_FLAG_BOUNDARY_CELL_BY_FACE     (1 << 0)
#define CS_FLAG_BOUNDARY_CELL_BY_VERTEX   (1 << 1)
#define CS_FLAG_BOUNDARY_CELL_BY_EDGE     (1 << 2)
#define CS_FLAG_SOLID_CELL                (1 << 3)

#define CS_FLAG_CURRENT_TO_PREVIOUS       (1 << 0)
#define CS_FLAG_INITIALIZATION            (1 << 1)

#define CS_FLAG_SYS_MASS_MATRIX  (1 << 0) 
#define CS_FLAG_SYS_SYM          (1 << 1) 
#define CS_FLAG_SYS_TIME_DIAG    (1 << 2) 
#define CS_FLAG_SYS_REAC_DIAG    (1 << 3) 
#define CS_FLAG_SYS_SOURCES_HLOC (1 << 4) 
#define CS_FLAG_SYS_VECTOR       (1 << 5) 
#define CS_FLAG_STATE_UNIFORM      (1 << 0) 
#define CS_FLAG_STATE_CELLWISE     (1 << 1) 
#define CS_FLAG_STATE_FACEWISE     (1 << 2) 
#define CS_FLAG_STATE_STEADY       (1 << 3) 
#define CS_FLAG_STATE_POTENTIAL    (1 << 4) 
#define CS_FLAG_STATE_CIRCULATION  (1 << 5) 
#define CS_FLAG_STATE_FLUX         (1 << 6) 
#define CS_FLAG_STATE_DENSITY      (1 << 7) 
#define CS_FLAG_FULL_LOC  (1 <<  0) 
#define CS_FLAG_SCALAR    (1 <<  1) 
#define CS_FLAG_VECTOR    (1 <<  2) 
#define CS_FLAG_TENSOR    (1 <<  3) 
#define CS_FLAG_VERTEX    (1 <<  4) 
#define CS_FLAG_EDGE      (1 <<  5) 
#define CS_FLAG_FACE      (1 <<  6) 
#define CS_FLAG_CELL      (1 <<  7) 
#define CS_FLAG_PRIMAL    (1 <<  8) 
#define CS_FLAG_DUAL      (1 <<  9) 
#define CS_FLAG_BORDER    (1 << 10) 
#define CS_FLAG_BY_CELL   (1 << 11) 
#define CS_FLAG_BY_FACE   (1 << 12) 
#define CS_FLAG_SCHEME_SCALAR  (1 << 0) 
#define CS_FLAG_SCHEME_VECTOR  (1 << 1) 
#define CS_FLAG_SCHEME_NAVSTO  (1 << 2) 
#define CS_FLAG_SCHEME_POLY0   (1 << 3) 
#define CS_FLAG_SCHEME_POLY1   (1 << 4) 
#define CS_FLAG_SCHEME_POLY2   (1 << 5) 
/*============================================================================
 * Type definitions
 *============================================================================*/

typedef unsigned char       cs_mask_t;   /* Elementary flag */

/*============================================================================
 * Global variables
 *============================================================================*/

/* Default locations */
extern const cs_flag_t  cs_flag_primal_vtx;
extern const cs_flag_t  cs_flag_primal_edge;
extern const cs_flag_t  cs_flag_primal_face;
extern const cs_flag_t  cs_flag_primal_cell;
extern const cs_flag_t  cs_flag_vertex; /* equal to cs_flag_primal_vtx */
extern const cs_flag_t  cs_flag_cell;   /* equal to cs_flag_primal_cell */
extern const cs_flag_t  cs_flag_boundary_face;

extern const cs_flag_t  cs_flag_dual_vtx;
extern const cs_flag_t  cs_flag_dual_face;
extern const cs_flag_t  cs_flag_dual_cell;
extern const cs_flag_t  cs_flag_dual_cell_byc;
extern const cs_flag_t  cs_flag_dual_face_byc;

/* Part of dual cell closure belonging to a boundary primal face */

extern const cs_flag_t  cs_flag_dual_closure_byf;

/* According to the extended flag defined below one can identify which set of
 * quantities or connectivities have to be built on-the-fly and stored in a
 * local structure possibly owned by each thread and with a cellwise scope */

typedef unsigned int  cs_eflag_t;

/* Store predefined flags */

extern const cs_eflag_t  cs_flag_need_v;
extern const cs_eflag_t  cs_flag_need_e;
extern const cs_eflag_t  cs_flag_need_f;
extern const cs_eflag_t  cs_flag_need_fe;
extern const cs_eflag_t  cs_flag_need_ef;
extern const cs_eflag_t  cs_flag_need_peq;
extern const cs_eflag_t  cs_flag_need_dfq;
extern const cs_eflag_t  cs_flag_need_pfq;
extern const cs_eflag_t  cs_flag_need_deq;
extern const cs_eflag_t  cs_flag_need_pfc;

typedef enum {

  /* Compute simple and cellwise information for vertices */

  CS_FLAG_COMP_PV      = 1 << 0,   /* = 1 */

  /* Compute cellwise quantities for vertices */

  CS_FLAG_COMP_PVQ     = 1 << 1,   /* = 2 */

  /* Compute simple and cellwise information for edges */

  CS_FLAG_COMP_PE      = 1 << 2,   /* = 4 */

  /* Compute cellwise quantities for edges */

  CS_FLAG_COMP_PEQ     = 1 << 3,   /* = 8 */

  /* Compute cellwise quantities for dual faces (associated to edges) */

  CS_FLAG_COMP_DFQ     = 1 << 4,   /* = 16 */

  /* Compute simple and cellwise information for faces */

  CS_FLAG_COMP_PF      = 1 << 5,   /* = 32 */

  /* Compute cellwise quantities for faces */

  CS_FLAG_COMP_PFQ     = 1 << 6,   /* = 64 */

  /* Compute cellwise quantities for dual edges (associated to faces) */

  CS_FLAG_COMP_DEQ     = 1 << 7,   /* = 128 */

  /* Compute the cellwise connectivity edge to vertices */

  CS_FLAG_COMP_EV      = 1 << 8,   /* = 256 */

  /* Compute cellwise connectivity face to edges */

  CS_FLAG_COMP_FE      = 1 << 9,   /* = 512 */

  /* Compute cellwise quantities associated to the couple (face, edge) */

  CS_FLAG_COMP_FEQ     = 1 << 10,  /* = 1024 */

  /* Compute cellwise connectivity face to vertices */

  CS_FLAG_COMP_FV      = 1 << 11,  /* = 2048 */

  /* Compute cellwise connectivity edge to faces */

  CS_FLAG_COMP_EF      = 1 << 12,  /* = 4096 */

  /* Compute elemental portion of dual faces associated to the couple
     (edge, face) */

  CS_FLAG_COMP_SEF     = 1 << 13,  /* = 8192 */

  /* Compute cellwise quantities related to the height of the pyramid with
     basis spanned by a face and with apex the cell center */

  CS_FLAG_COMP_HFQ     = 1 << 14,  /* = 16384 */

  /* Compute cellwise orientation of oriented edges belonging to a face */

  CS_FLAG_COMP_FES     = 1 << 15,  /* = 32768 */

  /* Compute cellwise quantities related to the volume of the pyramid with
     basis spanned by a face and with apex the cell center */

  CS_FLAG_COMP_PFC     = 1 << 16,  /* = 65536 */

  /* Compute cellwise quantities related to the volume surrounding an edge */

  CS_FLAG_COMP_PEC     = 1 << 17,  /* = 131072 */

  /* Compute cellwise diameters */

  CS_FLAG_COMP_DIAM    = 1 << 18,  /* = 262144 */

} cs_flag_comp_bits_t;

typedef enum {

  CS_FLAG_LOCATION_PRIMAL_VTX = CS_FLAG_PRIMAL | CS_FLAG_VERTEX,
  CS_FLAG_LOCATION_PRIMAL_EDGE = CS_FLAG_PRIMAL | CS_FLAG_EDGE,
  CS_FLAG_LOCATION_PRIMAL_FACE = CS_FLAG_PRIMAL | CS_FLAG_FACE,
  CS_FLAG_LOCATION_PRIMAL_CELL = CS_FLAG_PRIMAL | CS_FLAG_CELL,
  CS_FLAG_LOCATION_DUAL_VTX  = CS_FLAG_DUAL | CS_FLAG_VERTEX,
  CS_FLAG_LOCATION_DUAL_EDGE = CS_FLAG_DUAL | CS_FLAG_EDGE,
  CS_FLAG_LOCATION_DUAL_FACE = CS_FLAG_DUAL | CS_FLAG_FACE,
  CS_FLAG_LOCATION_DUAL_CELL = CS_FLAG_DUAL | CS_FLAG_CELL,

  CS_FLAG_N_LOCATIONS

} cs_flag_location_t;

/*============================================================================
 * Public function prototypes
 *============================================================================*/

/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/

static inline bool
cs_flag_test(cs_flag_t    flag_to_check,
             cs_flag_t    reference)
{
  if ((flag_to_check & reference) == reference)
    return true;
  else
    return false;
}

/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/

static inline bool
cs_flag_at_least(cs_flag_t    flag_to_check,
                 int          n_masks,
                 cs_flag_t    masks[])
{
  for (int i = 0; i < n_masks; i++)
    if ((flag_to_check & masks[i]) == masks[i])
      return true;
  return false;
}

/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/

static inline bool
cs_eflag_test(cs_eflag_t    flag_to_check,
              cs_eflag_t    reference)
{
  if ((flag_to_check & reference) == reference)
    return true;
  else
    return false;
}

/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/

const char *
cs_flag_str_location(cs_flag_t  loc);

/*----------------------------------------------------------------------------*/

END_C_DECLS

#endif /* __CS_FLAG_H__ */