Revision fb6a558c8d9304637825079d884e437c0c0dc820 authored by Jaime Garibay on 21 May 2021, 14:25:30 UTC, committed by Jaime Garibay on 03 June 2021, 06:24:29 UTC
This feature adds the capability to model exchange reactions using the Phreeqc solver
1 parent d49bec2
SurfaceGrid.h
/*
* \date 2012-09-22
* \brief Declaration of the SurfaceGrid class.
*
* \file
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*/
#pragma once
#include <array>
#include <limits>
#include <optional>
#include <vector>
#include "AABB.h"
#include "Point.h"
namespace GeoLib {
// forward declarations
class Triangle;
class Surface;
class SurfaceGrid : public AABB {
public:
explicit SurfaceGrid(GeoLib::Surface const*const sfc);
bool isPointInSurface(MathLib::Point3d const & pnt,
double eps = std::numeric_limits<double>::epsilon()) const;
private:
void sortTrianglesInGridCells(GeoLib::Surface const* const sfc);
bool sortTriangleInGridCells(GeoLib::Triangle const*const triangle);
std::optional<std::array<std::size_t, 3>> getGridCellCoordinates(
MathLib::Point3d const& p) const;
std::array<double,3> _step_sizes{};
std::array<double,3> _inverse_step_sizes{};
std::array<std::size_t,3> _n_steps;
std::vector<std::vector<GeoLib::Triangle const*>> _triangles_in_grid_box;
};
} // end namespace GeoLib
Computing file changes ...