DirectConditionGenerator Class Reference

Detailed Description

Definition at line 13 of file DirectConditionGenerator.h.

#include <DirectConditionGenerator.h>

Public Member Functions
	DirectConditionGenerator ()=default
	~DirectConditionGenerator ()=default
const std::vector< std::pair< std::size_t, double > > &	directToSurfaceNodes (const MeshLib::Mesh &mesh, const std::string &filename)
const std::vector< std::pair< std::size_t, double > > &	directWithSurfaceIntegration (MeshLib::Mesh &mesh, const std::string &filename, double scaling)
int	writeToFile (const std::string &name) const

Private Attributes
std::vector< std::pair< std::size_t, double > >	_direct_values

Constructor & Destructor Documentation

DirectConditionGenerator()

DirectConditionGenerator::DirectConditionGenerator ( )

default

~DirectConditionGenerator()

DirectConditionGenerator::~DirectConditionGenerator ( )

default

Member Function Documentation

directToSurfaceNodes()

const std::vector< std::pair< std::size_t, double > > & DirectConditionGenerator::directToSurfaceNodes	(	const MeshLib::Mesh &	mesh,
		const std::string &	filename )

Definition at line 19 of file DirectConditionGenerator.cpp.

{
    if (_direct_values.empty())
    {
        GeoLib::Raster* raster(
            FileIO::AsciiRasterInterface::readRaster(filename));
        if (!raster)
        {
            ERR("Error in DirectConditionGenerator::directToSurfaceNodes() - "
                "could not load raster file.");
            return _direct_values;
        }
 
        Eigen::Vector3d const dir = -Eigen::Vector3d::UnitZ();
        const std::vector<MeshLib::Node*> surface_nodes(
            MeshToolsLib::MeshSurfaceExtraction::getSurfaceNodes(mesh, dir,
                                                                 90));
        const double no_data(raster->getHeader().no_data);
        _direct_values.reserve(surface_nodes.size());
        for (auto const* surface_node : surface_nodes)
        {
            double val(raster->getValueAtPoint(*surface_node));
            val = (val == no_data) ? 0 : val;
            _direct_values.emplace_back(surface_node->getID(), val);
        }
        delete raster;
 
        std::for_each(surface_nodes.begin(), surface_nodes.end(),
                      std::default_delete<MeshLib::Node>());
    }
    else
        ERR("Error in DirectConditionGenerator::directToSurfaceNodes() - Data "
            "vector contains outdated values.");
 
    return _direct_values;
}

References _direct_values, ERR(), GeoLib::Raster::getHeader(), MeshToolsLib::MeshSurfaceExtraction::getSurfaceNodes(), GeoLib::Raster::getValueAtPoint(), GeoLib::RasterHeader::no_data, and FileIO::AsciiRasterInterface::readRaster().

Referenced by CondFromRasterDialog::accept().

directWithSurfaceIntegration()

const std::vector< std::pair< std::size_t, double > > & DirectConditionGenerator::directWithSurfaceIntegration	(	MeshLib::Mesh &	mesh,
		const std::string &	filename,
		double	scaling )

Definition at line 58 of file DirectConditionGenerator.cpp.

{
    if (!_direct_values.empty())
    {
        ERR("Error in DirectConditionGenerator::directWithSurfaceIntegration()"
            "- Data vector contains outdated values...");
        return _direct_values;
    }
 
    std::unique_ptr<GeoLib::Raster> raster(
        FileIO::AsciiRasterInterface::readRaster(filename));
    if (!raster)
    {
        ERR("Error in DirectConditionGenerator::directWithSurfaceIntegration()"
            "- could not load raster file.");
        return _direct_values;
    }
 
    Eigen::Vector3d const dir({0.0, 0.0, -1.0});
    double const angle(90);
    auto prop_name = MeshLib::getBulkIDString(MeshLib::MeshItemType::Node);
    std::unique_ptr<MeshLib::Mesh> surface_mesh(
        MeshToolsLib::MeshSurfaceExtraction::getMeshSurface(mesh, dir, angle,
                                                            prop_name));
 
    std::vector<double> node_area_vec =
        MeshToolsLib::MeshSurfaceExtraction::getSurfaceAreaForNodes(
            *surface_mesh);
    const std::vector<MeshLib::Node*>& surface_nodes(surface_mesh->getNodes());
    const std::size_t nNodes(surface_mesh->getNumberOfNodes());
    const double no_data(raster->getHeader().no_data);
 
    MeshLib::PropertyVector<int> const* node_id_pv = nullptr;
    try
    {
        node_id_pv = surface_mesh->getProperties().getPropertyVector<int>(
            prop_name, MeshLib::MeshItemType::Node, 1);
    }
    catch (std::runtime_error const& e)
    {
        WARN("{:s}", e.what());
        return _direct_values;
    }
 
    _direct_values.reserve(nNodes);
    for (std::size_t i = 0; i < nNodes; ++i)
    {
        double val(raster->getValueAtPoint(*surface_nodes[i]));
        val = (val == no_data) ? 0 : ((val * node_area_vec[i]) / scaling);
        _direct_values.emplace_back((*node_id_pv)[i], val);
    }
 
    return _direct_values;
}

References _direct_values, ERR(), MeshLib::getBulkIDString(), MeshToolsLib::MeshSurfaceExtraction::getMeshSurface(), MeshToolsLib::MeshSurfaceExtraction::getSurfaceAreaForNodes(), MeshLib::Node, FileIO::AsciiRasterInterface::readRaster(), and WARN().

Referenced by CondFromRasterDialog::accept().

writeToFile()

int DirectConditionGenerator::writeToFile

(

const std::string &

name

)

const

Definition at line 114 of file DirectConditionGenerator.cpp.

{
    std::ofstream out(name.c_str(), std::ios::out);
 
    if (out)
    {
        for (const auto& direct_value : _direct_values)
        {
            out << direct_value.first << "\t" << direct_value.second << "\n";
        }
 
        out.close();
    }
    return 0;
}

References _direct_values.

Member Data Documentation

_direct_values

std::vector< std::pair<std::size_t,double> > DirectConditionGenerator::_direct_values

private

Definition at line 29 of file DirectConditionGenerator.h.

Referenced by directToSurfaceNodes(), directWithSurfaceIntegration(), and writeToFile().

---

The documentation for this class was generated from the following files:

• DirectConditionGenerator.h
• DirectConditionGenerator.cpp