ERF_IBSEBMaterials.H File Reference

Material library of the building faces. More...

#include <AMReX_REAL.H>
#include <AMReX_Vector.H>
#include <string>

Include dependency graph for ERF_IBSEBMaterials.H:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	IBSEBMaterial

Functions
amrex::Vector< IBSEBMaterial >	read_ibseb_materials (const std::string &file)

Detailed Description

Material library of the building faces.

The CSV schema is the one of the ERF-SLUCM branch's material library, so the same file serves both models:

mat_id,name,albedo,emissivity,k_therm_W_per_mK,rho_cp_J_per_m3K,thickness_m,description

The I/O rank reads the file and broadcasts the entries as plain data; every rank then holds the same table. A face gets its material from its building (erf.ibseb.material_by_building) or the default (erf.ibseb.material_default); without a file every face uses the uniform properties of the inputs.

Function Documentation

read_ibseb_materials()

amrex::Vector<IBSEBMaterial> read_ibseb_materials

(

const std::string &

file

)

Read the material CSV on the I/O rank and broadcast it. Aborts on a missing file, a header that does not match the schema, a duplicate id, or a value outside its physical range.

{
    static_assert(std::is_trivially_copyable<IBSEBMaterial>::value, "IBSEBMaterial must be broadcastable");
    Vector<IBSEBMaterial> table;
    int n = 0;
    if (ParallelDescriptor::IOProcessor()) {
        std::ifstream in(file);
        if (!in.good()) { Abort("erf.ibseb.material_file: cannot open '" + file + "'"); }
        std::string line;
        if (!std::getline(in, line)) { Abort("erf.ibseb.material_file: empty file '" + file + "'"); }
        if (line.size() >= 3 && static_cast<unsigned char>(line[0]) == 0xEF) { line.erase(0, 3); }
        const std::string expected = "mat_id,name,albedo,emissivity,k_therm_W_per_mK,rho_cp_J_per_m3K,thickness_m,description";
        if (strip(line) != expected) {
            Abort("erf.ibseb.material_file: header must be '" + expected + "', got '" + line + "'");
        }
        std::set<int> seen;
        while (std::getline(in, line)) {
            if (strip(line).empty() || line[0] == '#') { continue; }
            std::stringstream ss(line);
            std::string f;
            IBSEBMaterial m;
            std::getline(ss, f, ','); m.mat_id = std::stoi(strip(f));
            std::getline(ss, f, ','); f = strip(f); std::strncpy(m.name, f.c_str(), sizeof(m.name) - 1); m.name[sizeof(m.name) - 1] = '\0';
            std::getline(ss, f, ','); m.albedo     = std::stod(strip(f));
            std::getline(ss, f, ','); m.emissivity = std::stod(strip(f));
            std::getline(ss, f, ','); m.k_therm    = std::stod(strip(f));
            std::getline(ss, f, ','); m.rho_cp     = std::stod(strip(f));
            std::getline(ss, f, ','); m.thickness  = std::stod(strip(f));
            if (!seen.insert(m.mat_id).second) { Abort("erf.ibseb.material_file: duplicate mat_id " + std::to_string(m.mat_id)); }
            if (m.albedo < 0.0 || m.albedo > 1.0 || m.emissivity <= 0.0 || m.emissivity > 1.0 ||
                m.k_therm <= 0.0 || m.rho_cp <= 0.0 || m.thickness <= 0.0) {
                Abort("erf.ibseb.material_file: material " + std::to_string(m.mat_id) + " has a value outside its range");
            }
            table.push_back(m);
        }
        n = static_cast<int>(table.size());
        if (n == 0) { Abort("erf.ibseb.material_file: no materials in '" + file + "'"); }
    }
    ParallelDescriptor::Bcast(&n, 1, ParallelDescriptor::IOProcessorNumber());
    table.resize(n);
    ParallelDescriptor::Bcast(reinterpret_cast<char*>(table.data()), n * sizeof(IBSEBMaterial),
                              ParallelDescriptor::IOProcessorNumber());
    return table;
}

Referenced by IBFaceSet::assign_materials().

Here is the caller graph for this function: