io.cpp

Go to the documentation of this file.

#include "io.hpp"
 
#include <filesystem>
#include <fstream>
#include <iomanip>
#include <stdexcept>
 
namespace fs = std::filesystem;
 
/**
 * @brief Save a temperature profile (x, T) to a CSV file.
 *
 * This function automatically creates the directory structure if missing.
 * Example output path: "result/dufort_frankel/t_0.10.csv"
 *
 * @param path  Relative or absolute path to the file.
 * @param x     Spatial coordinates [cm].
 * @param T     Temperature values [°C].
 */
void save_profile_csv(const std::string& path,
                      const std::vector<double>& x,
                      const std::vector<double>& T)
{
    // Validate input sizes
    if (x.size() != T.size())
        throw std::runtime_error("save_profile_csv: vector size mismatch");
 
    // --- Ensure the directory exists (Resolved Conflict) ---
    // We use the safer check from the 'dev' branch.
    fs::path fullPath(path);
    if (fullPath.has_parent_path())
    {
        fs::create_directories(fullPath.parent_path());
    }
 
    // --- Open file for writing ---
    std::ofstream file(fullPath);
    if (!file)
        throw std::runtime_error("save_profile_csv: cannot open file " + fullPath.string());
 
    // --- Write header and data ---
    file << "x_cm,T_C\n";
    file << std::fixed << std::setprecision(6);
    for (std::size_t i = 0; i < x.size(); ++i)
        file << x[i] << "," << T[i] << "\n";
}