cs_param_cdo.h

Go to the documentation of this file.

#ifndef __CS_PARAM_CDO_H__
#define __CS_PARAM_CDO_H__
 
/*============================================================================
 * Manage the definition/setting of a computation
 *============================================================================*/
 
/*
  This file is part of Code_Saturne, a general-purpose CFD tool.
 
  Copyright (C) 1998-2019 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
 *============================================================================*/
 
/* Specifications for OpenMP loops */
#define CS_CDO_OMP_CHUNK_SIZE     128
#define CS_CDO_OMP_SCHEDULE       schedule(static, CS_CDO_OMP_CHUNK_SIZE)
#define CS_CDO_OMP_SYNC_SECTIONS  0 /* > 0 --> critical sections
                                       otherwise atomic sections */
 
/* Avoid issues with assert in some OpenMp contructs using gcc 9 */
#if defined(HAVE_OPENMP) && defined(__GNUC__)
  #if __GNUC__ == 9
    #define CS_CDO_OMP_ASSERT(e)
  #else
    #define CS_CDO_OMP_ASSERT(e)  assert(e)
  #endif
#else
  #define CS_CDO_OMP_ASSERT(e)  assert(e)
#endif
 
/* Size of the buffer used to collect global ids for rows and columns
   when assembling the values in the global matrix from the local cellwise
   matrices */
#define CS_CDO_ASSEMBLE_BUF_SIZE  99
 
/* The following limitation only results from an optimization in the size of
   the bit mask (can be changed if needed by changing the definition of
   the type cs_mask_t) */
#define CS_CDO_N_MAX_REACTIONS  8 // Max number of reaction terms in an equation
 
#define CS_ALL_FACES   0        /* All faces: interior + border */
#define CS_BND_FACES   1        /* Boundary faces */
#define CS_INT_FACES   2        /* Interior faces */
 
/* Number of DoFs on faces and cells according to the polynomial space */
#define CS_N_FACE_DOFS_0TH  1
#define CS_N_FACE_DOFS_1ST  3
#define CS_N_FACE_DOFS_2ND  6
 
#define CS_N_CELL_DOFS_0TH  1
#define CS_N_CELL_DOFS_1ST  4
#define CS_N_CELL_DOFS_2ND  10
 
/*============================================================================
 * Type definitions
 *============================================================================*/
 
/* OpenMP STRATEGY FOR THE ASSEMBLY STEP */
/* ===================================== */
 
typedef enum {
 
  CS_PARAM_ASSEMBLE_OMP_ATOMIC,
  CS_PARAM_ASSEMBLE_OMP_CRITICAL,
  CS_PARAM_ASSEMBLE_OMP_N_STRATEGIES
 
} cs_param_assemble_omp_strategy_t;
 
/* DISCRETE HODGE OPERATORS */
/* ======================== */
 
typedef enum {
 
  /* Hodge operator between dual spaces */
  CS_PARAM_HODGE_TYPE_VPCD, /* from primal vertices to dual cells */
  CS_PARAM_HODGE_TYPE_EPFD, /* from primal edges to dual faces */
  CS_PARAM_HODGE_TYPE_FPED, /* from primal faces to dual edges */
  CS_PARAM_HODGE_TYPE_EDFP, /* from dual edges to primal faces */
  CS_PARAM_HODGE_TYPE_CPVD, /* from primal cells to dual vertices */
 
  /* Hodge operator for hybrid spaces */
  CS_PARAM_HODGE_TYPE_FB,   /* primal face + primal cells */
  CS_PARAM_HODGE_TYPE_VC,   /* primal vertices + primal cells */
  CS_PARAM_N_HODGE_TYPES
 
} cs_param_hodge_type_t;
 
typedef enum {
 
  CS_PARAM_HODGE_ALGO_VORONOI, // Under othogonality condition gives a diag. op.
  CS_PARAM_HODGE_ALGO_WBS,     // WBS: Whitney Barycentric Subdivision
  CS_PARAM_HODGE_ALGO_COST,    // COST: COnsistency & STabilization splitting
  CS_PARAM_HODGE_ALGO_AUTO,    /* Switch between the previous algo. according to
                                  the type of cell and its property */
  CS_PARAM_N_HODGE_ALGOS
 
} cs_param_hodge_algo_t;
 
typedef struct {
 
  bool   is_unity; /* No associated property. Property is equalt to the unity */
  bool   is_iso;   /* Associated property is isotroopic ? */
  bool   inv_pty;  /* Definition based on the property or its inverse */
 
  cs_param_hodge_type_t   type;  /* type of discrete Hodge operator */
  cs_param_hodge_algo_t   algo;  /* type of algorithm used to build this op. */
  double                  coef;  /* Value of the stabilization parameter
                                    if the COST algo. is used, otherwise 0. */
 
} cs_param_hodge_t;
 
/*============================================================================
 * Global variables
 *============================================================================*/
 
/* Separation lines: header1, header2 (compatible with markdown), other */
extern const char h1_sep[80];
extern const char h2_sep[80];
extern const char sepline[80];
extern const char msepline[50];
 
/* Activation of the CDO/HHO module */
extern int  cs_param_cdo_mode;
 
/*============================================================================
 * Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
const char *
cs_param_hodge_get_algo_name(const cs_param_hodge_t   h_info);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
const char *
cs_param_hodge_get_type_name(const cs_param_hodge_t   h_info);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_param_hodge_log(const char               *prefix,
                   const cs_param_hodge_t    hp);
 
/*----------------------------------------------------------------------------*/
 
END_C_DECLS
 
#endif /* __CS_PARAM_CDO_H__ */