src/tetmeshbuilder.cpp

Go to the documentation of this file.

/***************************************************************************
                          tetmeshbuilder.h.cpp  -  description
                             -------------------
    begin                : Fri Dec 12 2003
    copyright            : (C) 2003 by Roman Geus
    email                : roman.geus@psi.ch
***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   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.                                   *
 *                                                                         *
 ***************************************************************************/

#include <new>
#include "tetmeshbuilder.h"

TetMeshBuilder::TetMeshBuilder() {
    _mesh = new mesh::TetMesh();
}

TetMeshBuilder::TetMeshBuilder(MPI_Comm comm) {
    _mesh = new mesh::ParallelTetMesh(comm);
}

TetMeshBuilder::~TetMeshBuilder() {
}

void TetMeshBuilder::init_coord(int nof_node) {
    _mesh->initPoint(nof_node);
}

void TetMeshBuilder::set_coord(int i, double x, double y, double z) {
    _mesh->insertPoint(i, x, y, z);
}

void TetMeshBuilder::finalize_coord() {
    _mesh->generateEdges();
    _mesh->generateFaces();
}

void TetMeshBuilder::init_tet(int nof_tet) {
    _mesh->initTet(nof_tet);
}

void TetMeshBuilder::set_tet(int t, int id0, int id1, int id2, int id3, int material) {
  _mesh->insertTet(t, id0, id1, id2, id3, material);
}

void TetMeshBuilder::finalize_tet() {
    _mesh->generateEdges();
    _mesh->generateFaces();
}

void TetMeshBuilder::init_bc(int nof_bc_face) {
//Dag: This is not used in the ParallelTetMesh class since it appearently has no function and the local no. of faces is not known at this time in the parallel reader.
}

void TetMeshBuilder::set_bc(int id0, int id1, int id2, int bc_id) {
    mesh::id_t face = _mesh->lookupFace(id0, id1, id2);
    if (bc_id == 0) {
        _mesh->setBoundaryOnFace(face);
    } else {
        assert((0 < bc_id) && (bc_id <= 3));
        _mesh->setSymmetryPlaneOnFace(face, bc_id - 1);
    }
}

void TetMeshBuilder::finalize_bc(int nof_sym) {
    _mesh->finalize_mesh(nof_sym);
}

---