HeatStructure.cpp

Go to the documentation of this file.

// HeatStructure.cpp: implementation of the HeatStructure class.
//
 
#include <fstream>
#include <boost/serialization/export.hpp>
 
#include "HeatStructure.h"
#include <boost/archive/binary_oarchive.hpp>
#include <boost/archive/binary_iarchive.hpp>
BOOST_CLASS_EXPORT_IMPLEMENT(HeatStructure)
 
 
using namespace std;
 
// Construction/Destruction
 
HeatStructure::HeatStructure()
{
    number_of_mesh_points = 0;
    unit_of_mesh_point_temperature = "NOT_SPECIFIED";
    THElement::setTypeOfCard(Relap5ComponentType::HeatStructure);
}
 
HeatStructure::~HeatStructure()
{
}
void HeatStructure::setMeshPointTemperature(int n, double x){
    int i, n_gap;
    if(mesh_point_temperature.size() >= (n + 1)){//n starts at 0 (zero)
        mesh_point_temperature[n] = x;
    }
    else{
        n_gap = n + 1 - mesh_point_temperature.size();
        for(i = 0; i < n_gap; i++){
            mesh_point_temperature.push_back(-1.0);
        }
        mesh_point_temperature[n] = x;
    }
}
void HeatStructure::setNumberOfMeshPoints(int n){
    number_of_mesh_points = n;
}
double HeatStructure::getMeshPointTemperature(int n){
    if(mesh_point_temperature.size() >= (n + 1)){
        return mesh_point_temperature[n];
    }
    else{
        ofstream fout("err messages.txt", ios::app);
        fout << r5pars_.timesno << "  HeatStructure::getMeshPointTemperature the n (" << n << ") exceeds the size of mesh_point_temperature array\n";
        fout.close();
        return -1.0;
    }
}
int HeatStructure::getNumberOfMeshPoints(){
    return number_of_mesh_points; 
}
 
 
void HeatStructure::setUnits(bool  si_unit){
    this->setUnitOfMeshPointTemperature(si_unit);
};
void HeatStructure::setUnitOfMeshPointTemperature(bool si_unit){
    if (si_unit){
        unit_of_mesh_point_temperature = "K";
    }
    else{
        unit_of_mesh_point_temperature = "F";
    }
}
string HeatStructure::getUnitOfMeshPointTemperature(){return unit_of_mesh_point_temperature;}
void HeatStructure::readInfoFromFile(ifstream &fin){
    THElement::readInfoFromFile(fin);
    fin >> this->number_of_mesh_points;
}
void HeatStructure::saveHeatStructureToFile(ofstream &fout){
    int i, num;
    num = mesh_point_temperature.size();
    fout << number_of_mesh_points << ' '
    << unit_of_mesh_point_temperature << ' '
    << num;
    for(i = 0; i < num; i++){
        fout << ' ' << mesh_point_temperature.at(i);
    }
    fout << '\n';
    this->saveTHElementToFile(fout);
}
void HeatStructure::retrieveHeatStructureFromFile(ifstream &fin){
    int i, num;
    double value;
    fin >>  number_of_mesh_points >> unit_of_mesh_point_temperature >> num;
    for(i = 0; i < num; i++){
        fin >> value;
        mesh_point_temperature.push_back(value);
    }
    this->retrieveTHElementFromFile(fin);
}
void HeatStructure::copyHeatStructure(HeatStructure &_hs){
    int i, num;
    this->copyTHElement(_hs);
    number_of_mesh_points = _hs.getNumberOfMeshPoints();
    unit_of_mesh_point_temperature = _hs.getUnitOfMeshPointTemperature();
    num = _hs.mesh_point_temperature.size();
    for(i = 0; i < num; i++){
        this->mesh_point_temperature.push_back(_hs.mesh_point_temperature.at(i));
    }
}