cs_cdo_connect.h

Go to the documentation of this file.

#ifndef __CS_CDO_CONNECT_H__
#define __CS_CDO_CONNECT_H__
 
/*============================================================================
 * Manage connectivity (Topological features of the mesh)
 *============================================================================*/
 
/*
  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 "fvm_defs.h"
 
#include "cs_base.h"
#include "cs_flag.h"
#include "cs_mesh.h"
#include "cs_mesh_adjacencies.h"
#include "cs_range_set.h"
 
/*----------------------------------------------------------------------------*/
 
BEGIN_C_DECLS
 
/*============================================================================
 * Macro definitions
 *============================================================================*/
 
/* Main categories to consider for high-level structures
   Remark: scalar-valued HHO-P1 and vector-valued CDO-Fb shares the same
   structures
*/
#define CS_CDO_CONNECT_VTX_SCAL   0 /* Vb or VCb scalar-valued eq. */
#define CS_CDO_CONNECT_VTX_VECT   1 /* Vb or VCb vector-valued eq. */
#define CS_CDO_CONNECT_FACE_SP0   2 /* Fb or HHO-P0 scalar-valued eq. */
#define CS_CDO_CONNECT_FACE_VP0   3 /* Fb vector-valued eq. */
#define CS_CDO_CONNECT_FACE_SP1   3 /* HHO-P1 scalar-valued */
#define CS_CDO_CONNECT_FACE_SP2   4 /* HHO-P2 scalar-valued eq. */
#define CS_CDO_CONNECT_FACE_VHP0  3 /* HHO-P0 vector-valued eq. */
#define CS_CDO_CONNECT_FACE_VHP1  5 /* HHO-P1 vector-valued eq. */
#define CS_CDO_CONNECT_FACE_VHP2  6 /* HHO-P2 vector-valued eq. */
 
#define CS_CDO_CONNECT_N_CASES    7
 
/* Additional macros */
#define CS_TRIANGLE_CASE          3 /* Number of vertices in a triangle */
 
/*============================================================================
 * Type definitions
 *============================================================================*/
 
/* Connectivity structure */
typedef struct {
 
  cs_lnum_t          n_vertices;
  cs_lnum_t          n_edges;
  cs_lnum_t          n_faces[3];  /* 0: all, 1: border, 2: interior */
  cs_lnum_t          n_cells;
 
  /* Edge-related members */
  cs_adjacency_t    *e2v;         /* edge --> vertices connectivity */
 
  /* Face-related members */
  cs_adjacency_t    *f2c;         /* face --> cells connectivity */
  cs_adjacency_t    *f2e;         /* face --> edges connectivity */
  cs_adjacency_t    *bf2v;        /* border face --> vertices connectivity
                                     (map from cs_mesh_t) */
  cs_adjacency_t    *if2v;        /* interior face --> vertices connectivity
                                     (map from cs_mesh_t) */
 
  /* Cell-related members */
  fvm_element_t     *cell_type;   /* type of cell */
  cs_flag_t         *cell_flag;   /* Flag (Border) */
  cs_adjacency_t    *c2f;         /* cell --> faces connectivity */
  cs_adjacency_t    *c2e;         /* cell --> edges connectivity */
  cs_adjacency_t    *c2v;         /* cell --> vertices connectivity */
 
  /* Delta of ids between the min./max. values of entities related to a cell
     Useful to store compactly the link between mesh ids and cell mesh ids
     needed during the cell mesh definition */
  cs_lnum_t  e_max_cell_range;
  cs_lnum_t  v_max_cell_range;
 
  /* Max. connectitivy size for cells */
  int  n_max_vbyc;   /* max. number of vertices in a cell */
  int  n_max_ebyc;   /* max. number of edges in a cell */
  int  n_max_fbyc;   /* max. number of faces in a cell */
  int  n_max_vbyf;   /* max. number of vertices in a face */
  int  n_max_v2fc;   /* max. number of faces connected to a vertex in a cell */
  int  n_max_v2ec;   /* max. number of edges connected to a vertex in a cell */
 
  /* Adjacency related to most of the linear systems */
  cs_adjacency_t       *v2v;    /* vertex to vertices through cells */
  cs_adjacency_t       *f2f;    /* face to faces through cells */
 
  /* Structures to handle parallelism/assembler */
  cs_range_set_t       *range_sets[CS_CDO_CONNECT_N_CASES];
  cs_interface_set_t   *interfaces[CS_CDO_CONNECT_N_CASES];
 
} cs_cdo_connect_t;
 
/*============================================================================
 * Global variables
 *============================================================================*/
 
/*============================================================================
 * Static Inline Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
static inline void
cs_connect_get_next_3_vertices(const cs_lnum_t   *f2e_ids,
                               const cs_lnum_t   *e2v_ids,
                               const cs_lnum_t    start_idx,
                               cs_lnum_t         *v0,
                               cs_lnum_t         *v1,
                               cs_lnum_t         *v2)
{
  const cs_lnum_t _2e0  = 2*f2e_ids[start_idx],
                  _2e1  = 2*f2e_ids[start_idx+1];
  const cs_lnum_t tmp = e2v_ids[_2e1];
 
  *v0 = e2v_ids[_2e0];
  *v1 = e2v_ids[_2e0+1];
  *v2 = ((tmp != *v0) && (tmp != *v1)) ? tmp : e2v_ids[_2e1+1];
}
 
/*============================================================================
 * Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_cdo_connect_t *
cs_cdo_connect_init(cs_mesh_t      *mesh,
                    cs_flag_t       vb_scheme_flag,
                    cs_flag_t       vcb_scheme_flag,
                    cs_flag_t       fb_scheme_flag,
                    cs_flag_t       hho_scheme_flag);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_cdo_connect_t *
cs_cdo_connect_free(cs_cdo_connect_t   *connect);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_cdo_connect_summary(const cs_cdo_connect_t  *connect);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_cdo_connect_dump(const cs_cdo_connect_t  *connect);
 
/*----------------------------------------------------------------------------*/
 
END_C_DECLS
 
#endif /* __CS_CDO_CONNECT_H__ */