a_unit_system.cxx

Go to the documentation of this file.

/*
         Copyright 2019 Pierre SMARS (smars@yuntech.edu.tw)
         This program is free software: you can redistribute it and/or modify
         it under the terms of the GNU General Public License as published by
         the Free Software Foundation, either version 2 of the License, or
         (at your option) any later version.
 
         This program is distributed in the hope that it will be useful,
         but WITHOUT ANY WARRANTY; without even the implied warranty of
         MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
         GNU General Public License for more details.
 
         You should have received a copy of the GNU General Public License
         along with this program.  If not, see <http://www.gnu.org/licenses/>.
         */
#include "a_unit_system.h"
#include <sstream>
#include <iomanip>
#include <math.h>
#include <fstream>
 
//---------------------------------------------------------------------------
const std::string a_unit_system::help()
{
        std::ostringstream o;
        o << "***************" << std::endl;
        o << "a_unit_system:" << std::endl;
        o << "***************" << std::endl;
        o << "This is a 'a_unit_system' class" << std::endl;
        o << std::endl;
        o << "Create a unit system:" << std::endl;
        o << "a_unit_system s" << std::endl;
        o << std::endl;
        o << "Commands:" << std::endl;
        o << "name _name: set name of system of units" << std::endl;
        o << "name: get name of system of units" << std::endl;
        o << "decimals1 _val: set number of decimals in 1D" << std::endl;
        o << "decimals1: get number of decimals in 1D" << std::endl;
        o << "decimals2 _val: set number of decimals in 2D" << std::endl;
        o << "decimals2: get number of decimals in 2D" << std::endl;
        o << "addlevel1 _unit: add a unit to 1D unit set" << std::endl;
        o << "addlevel2 _unit: add a unit to 2D unit set" << std::endl;
        o << "level1: get pointer to a unit of the 1D unit set" << std::endl;
        o << "level2: get pointer to a unit of the 2D unit set" << std::endl;
        o << "dim1 _val: format a 1D measure" << std::endl;
        o << "dim1_det _val: format a 1D measure (+ system name)" << std::endl;
        o << "dim2 _val: format a 2D measure" << std::endl;
        o << "dim2_det _val: format a 2D measure (+ system name)" << std::endl;
        o << "dim1 _point: format a point" << std::endl;
        o << "print: set parameters" << std::endl;
        o << "read _file: read parameters from file" << std::endl;
        o << "write _file: write parameters to file" << std::endl;
        o << std::endl;
        o << a_unit_base::help();
        return o.str();
}
//***************************************************************************
//---------------------------------------------------------------------------
void a_unit_system::read(std::istream &in)
{
        std::getline(in,name_);
        int n;
        in >> n;
        for (int i=0; i<n; i++)
        {
                auto u = new a_unit;
                in >> u;
                this->addlevel1(*u);
        }
        in >> decimals1_;
        in >> n;
        for (int i=0; i<n; i++)
        {
                auto u = new a_unit;
                in >> u;
                this->addlevel2(*u);
        }
        in >> decimals2_;
}
//---------------------------------------------------------------------------
void a_unit_system::write(std::ostream &o) const
{
        o << this->name() << std::endl;
        int n1=this->levels1();
        o << n1 << std::endl;
        for (int i=0; i<n1; i++)
                o << level1(i) << std::endl;
        o << decimals1_ << std::endl;
        int n2=this->levels2();
        o << n2 << std::endl;
        for (int i=0; i<n2; i++)
                o << level2(i) << std::endl;
        o << decimals2_ << std::endl;
}
//---------------------------------------------------------------------------
void a_unit_system::reset()
{
        for (std::deque<a_unit *>::iterator u = levels1_.begin(); u < levels1_.end(); u++)
                delete *u;              
        levels1_.clear();
        for (std::deque<a_unit *>::iterator u = levels2_.begin(); u < levels2_.end(); u++)
                delete *u;              
        levels2_.clear();
}
//***************************************************************************
//---------------------------------------------------------------------------
a_unit* a_unit_system::level1(int n) const
{
        return levels1_[n];
}
//---------------------------------------------------------------------------
a_unit* a_unit_system::level2(int n) const
{
        return levels2_[n];
}
//---------------------------------------------------------------------------
a_unit_system::a_unit_system(const a_unit_system& u)
{
        name_ = u.name();
        decimals1_ = u.decimals1();
        for (int tel = 0; tel < u.levels1(); tel++)
        {
                auto temp = new a_unit(*u.level1(tel));
                levels1_.push_back(temp);
        }
        decimals2_ = u.decimals2();
        for (int tel = 0; tel < u.levels2(); tel++)
        {
                auto temp = new a_unit(*u.level2(tel));
                levels2_.push_back(temp);
        }
}
//---------------------------------------------------------------------------
a_unit_system::~a_unit_system()
{
        this->reset();
}
//---------------------------------------------------------------------------
void a_unit_system::addlevel1(a_unit& u)
{
        double r = u.ratio();
        for (std::deque<a_unit *>::iterator it=levels1_.begin(); it<levels1_.end(); it++)
        {
                if ((*it)->ratio()>r)
                {
                        levels1_.insert(it,&u);
                        return;
                }
        }
        levels1_.push_back(&u);
}
//---------------------------------------------------------------------------
void a_unit_system::addlevel2(a_unit& u)
{
        double r = u.ratio();
        for (std::deque<a_unit *>::iterator it=levels2_.begin(); it<levels2_.end(); it++)
        {
                if ((*it)->ratio()>r)
                {
                        levels2_.insert(it,&u);
                        return;
                }
        }
        levels2_.push_back(&u);
}
//---------------------------------------------------------------------------
std::string a_unit_system::dim1(const double v)
{
        double rem = fabs(v);
        std::ostringstream s;
        if (v<0) {s << "-";}
        int num = levels1_.size();
        for (int tel = 0; tel < num; tel++)
        {
                int v0;
                double v1;
                a_unit * unit = levels1_[tel];
                v1 = rem*unit->ratio();
                v0 = int(v1);
                rem -= v0/unit->ratio();
                if ((v0 != 0)||(tel == num-1)) 
                {
                        if ((tel == num-1)&&(decimals1_ != 0))
                        {
                                s << std::fixed << std::setprecision(decimals1_);
                                s << v0+rem*unit->ratio();
                        }
                        else
                                s << v0;        
                        s << unit->short_name() << " ";
                }
        }
        return s.str();
}
//---------------------------------------------------------------------------
std::string a_unit_system::dim1(const a_point& p)
{
        std::ostringstream    s;
        s << this->dim1(p.x()) << " " << this->dim1(p.y()) << " " << this->dim1(p.z()); 
        return s.str();
}
//---------------------------------------------------------------------------
std::string a_unit_system::dim1_det(const double v)
{
        std::string s0 = dim1(v);
        std::ostringstream s;
        s << "[" << this->name() << "]";
        return (s0+s.str());
}
//---------------------------------------------------------------------------
std::string a_unit_system::dim2(const double v)
{
        double rem = fabs(v);
        std::ostringstream s;
        if (v<0) {s << "-";}
        int num = levels2_.size();
        for (int tel = 0; tel < num; tel++)
        {
                int v0;
                double v1;
                a_unit * unit = levels2_[tel];
                v1 = rem*unit->ratio();
                v0 = int(v1);
                rem -= v0/unit->ratio();
                if ((v0 != 0)||(tel == num-1)) 
                {
                        if ((tel == num-1)&&(decimals2_ != 0))
                        {
                                s << std::fixed << std::setprecision(decimals2_);
                                s << v0+rem*unit->ratio();
                        }
                        else
                                s << v0;        
                        s << unit->short_name() << " ";
                }
        }
        return s.str();
}
//---------------------------------------------------------------------------
std::string a_unit_system::dim2_det(const double v)
{
        std::string s0 = dim2(v);
        std::ostringstream s;
        s << "[" << this->name() << "]";
        return (s0+s.str());
}