1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
// *****************************************************************************
/*!
  \file      src/IO/NetgenMeshWriter.cpp
  \copyright 2012-2015 J. Bakosi,
             2016-2018 Los Alamos National Security, LLC.,
             2019-2021 Triad National Security, LLC.,
             2022-2025 J. Bakosi
             All rights reserved. See the LICENSE file for details.
  \brief     Netgen mesh writer class definition
  \details   Netgen mesh writer class definition. Only supports tetrahedra.
*/
// *****************************************************************************

#include <iomanip><--- Include file:  not found. Please note: Cppcheck does not need standard library headers to get proper results.
#include <ostream><--- Include file:  not found. Please note: Cppcheck does not need standard library headers to get proper results.
#include <vector><--- Include file:  not found. Please note: Cppcheck does not need standard library headers to get proper results.
#include <algorithm><--- Include file:  not found. Please note: Cppcheck does not need standard library headers to get proper results.
#include <iterator><--- Include file:  not found. Please note: Cppcheck does not need standard library headers to get proper results.
#include <utility><--- Include file:  not found. Please note: Cppcheck does not need standard library headers to get proper results.
#include <cstddef><--- Include file:  not found. Please note: Cppcheck does not need standard library headers to get proper results.

#include "Types.hpp"
#include "UnsMesh.hpp"<--- Include file: "UnsMesh.hpp" not found.
#include "Exception.hpp"
#include "NetgenMeshWriter.hpp"

using tk::NetgenMeshWriter;

void
NetgenMeshWriter::writeMesh( const UnsMesh& mesh )
// *****************************************************************************
//  Public interface for writing Netgen mesh
//! \param[in] mesh Unstructured mesh object
// *****************************************************************************
{
  writeNodes( mesh );
  writeElements( mesh );
}

void
NetgenMeshWriter::writeNodes( const UnsMesh& mesh )
// *****************************************************************************
//  Write nodes
//! \param[in] mesh Unstructured mesh object
// *****************************************************************************
{
  // Write out number of nodes
  m_outFile << mesh.nnode() << std::endl;

  // Write node coordinates: x-coord y-coord z-coord
  m_outFile << std::setprecision(6) << std::fixed;
  for ( std::size_t i=0; i<mesh.nnode(); ++i ) {
    m_outFile << '\t' << mesh.x()[i]
              << '\t' << mesh.y()[i]
              << '\t' << mesh.z()[i] << std::endl;
  }
}

void
NetgenMeshWriter::writeElements( const UnsMesh& mesh )
// *****************************************************************************
//  Write elements, i.e., connectivity
//! \param[in] mesh Unstructured mesh object
// *****************************************************************************
{
  if (mesh.tetinpoel().empty()) return;

  // Make sure tetrahedron element connectivity starts with zero
  Assert( *std::minmax_element( begin(mesh.tetinpoel()),
                                end(mesh.tetinpoel()) ).first == 0,
          "tetrahedron node ids should start from zero" );

  // Get number of tetrahedra in mesh
  auto n = mesh.tetinpoel().size()/4;  

  // Write out number of tetrahedra
  m_outFile << n << std::endl;

  // Create empty tag vector
  std::vector< std::vector< int > > tg;
  tg.resize( n );
  for (auto& t : tg) t.push_back( 1 );

  // Write out tetrehadra element tags and connectivity
  for (std::size_t i=0; i<n; ++i) {
    // tag n[1-4]
    m_outFile << '\t' << tg[i][0]
              << '\t' << mesh.tetinpoel()[i*4+3]+1
              << '\t' << mesh.tetinpoel()[i*4+0]+1
              << '\t' << mesh.tetinpoel()[i*4+1]+1
              << '\t' << mesh.tetinpoel()[i*4+2]+1 << std::endl;
  }

  if (mesh.triinpoel().empty()) return;

  // Make sure triangle element connectivity starts with zero
  Assert( *std::minmax_element( begin(mesh.triinpoel()),
                                end(mesh.triinpoel()) ).first == 0,
          "triangle node ids should start from zero" );

  // Get number of triangles in mesh
  n = mesh.triinpoel().size()/3;

  // Write out number of triangles
  m_outFile << n << std::endl;

  // Create empty tag vector if there is no tag
  tg.clear();
  tg.resize( n );
  for (auto& t : tg) t.push_back( 1 );

  // Write out triangle element tags and connectivity
  for (std::size_t i=0; i<n; ++i) {
    // tag n[1-4]
    m_outFile << '\t' << tg[i][0]
              << '\t' << mesh.triinpoel()[i*3+0]+1
              << '\t' << mesh.triinpoel()[i*3+1]+1
              << '\t' << mesh.triinpoel()[i*3+2]+1 << std::endl;
  }
}