MeshToolsLib::MeshGenerator::AddFaultToVoxelGrid Namespace Reference

Functions
bool	isVoxelGrid (MeshLib::Mesh const &mesh)
bool	addFaultToVoxelGrid (MeshLib::Mesh *mesh, MeshLib::Mesh const *fault, int const fault_id)

Function Documentation

addFaultToVoxelGrid()

bool MeshToolsLib::MeshGenerator::AddFaultToVoxelGrid::addFaultToVoxelGrid	(	MeshLib::Mesh *	mesh,
		MeshLib::Mesh const *	fault,
		int const	fault_id )

Definition at line 190 of file AddFaultToVoxelGrid.cpp.

{
    if (mesh == nullptr)
    {
        ERR("Input mesh not found...");
        return false;
    }
    if (!isVoxelGrid(*mesh))
    {
        ERR("The input mesh is not a voxel grid. The input mesh must be "
            "a voxel grid (i.e. an equally sized axis "
            "aligned hexahedra mesh).");
        return false;
    }
 
    if (fault == nullptr)
    {
        ERR("Fault mesh not found...");
        return false;
    }
    if (fault->getDimension() != 2)
    {
        ERR("Fault needs to be a 2D mesh.");
        return false;
    }
 
    Eigen::Vector3d half_cell_size;
    {
        auto const n = *mesh->getElement(0)->getNode(0);
        auto const c = MeshLib::getCenterOfGravity(*mesh->getElement(0));
        half_cell_size[0] = std::abs(c[0] - n[0]);
        half_cell_size[1] = std::abs(c[1] - n[1]);
        half_cell_size[2] = std::abs(c[2] - n[2]);
    }
 
    markFaults(*mesh, *fault, fault_id, half_cell_size);
 
    return true;
}

References ERR(), MeshLib::getCenterOfGravity(), MeshLib::Mesh::getDimension(), MeshLib::Mesh::getElement(), MeshLib::Element::getNode(), and isVoxelGrid().

isVoxelGrid()

bool MeshToolsLib::MeshGenerator::AddFaultToVoxelGrid::isVoxelGrid

(

MeshLib::Mesh const &

mesh

)

Definition at line 159 of file AddFaultToVoxelGrid.cpp.

{
    auto const& elements = mesh.getElements();
    if (std::any_of(elements.cbegin(), elements.cend(),
                    [&](auto const& e) {
                        return (e->getGeomType() !=
                                MeshLib::MeshElemType::HEXAHEDRON);
                    }))
    {
        return false;
    }
    auto is_voxel = [](auto const& e)
    {
        auto const n = e->getNodes();
        return ((*n[0])[2] != (*n[1])[2] || (*n[1])[2] != (*n[2])[2] ||
                (*n[4])[2] != (*n[5])[2] || (*n[5])[2] != (*n[6])[2] ||
                (*n[1])[0] != (*n[2])[0] || (*n[2])[0] != (*n[5])[0] ||
                (*n[0])[0] != (*n[3])[0] || (*n[3])[0] != (*n[7])[0]);
    };
 
    if (std::any_of(elements.cbegin(), elements.cend(), is_voxel))
    {
        ERR("Input mesh needs to be voxel grid (i.e. equally sized axis "
            "aligned hexahedra).");
        return false;
    }
    return true;
}

References ERR(), and MeshLib::Mesh::getElements().

Referenced by AddFaultsToVoxelGridDialog::AddFaultsToVoxelGridDialog(), and addFaultToVoxelGrid().