ComputeIntersections.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <vtkType.h>
#include <vtkUnstructuredGrid.h>
 
#include <Eigen/Core>
#include <range/v3/view/enumerate.hpp>
#include <vector>
 
#include "ComputeNaturalCoordsResult.h"
 
namespace AU = ApplicationUtils;

struct IntersectionResult
{
    Eigen::Vector3d point;
    double t;  // parameter along the line from p0 to p1, i.e. p = p0 + t*(p1 -
               // p0) used for sorting
};

std::vector<std::vector<IntersectionResult>> getOrderedAnchorCoords(
    vtkUnstructuredGrid* grid,
    Eigen::MatrixX3d const& realcoords,
    Eigen::VectorXd const& free_fraction,
    double const tol);

AU::ComputeNaturalCoordsResult setPhysicalPropertiesForIntersectionPoints(
    std::vector<std::vector<IntersectionResult>> const& anchor_coords,
    AU::ComputeNaturalCoordsResult const& original_anchor_data);