quantum.cc

Go to the documentation of this file.

/* Copyright (C) 2013
   Oliver Lemke <olemke@core-dump.info>
 
   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, 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, write to the Free Software
   Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,
   USA. */
 
#include "quantum.h"
 
#include <stdexcept>
#include "mystring.h"
 
 
bool QuantumNumbers::Compare(const QuantumNumbers& qn) const
{
    const QuantumContainer& qnumbers2 = qn.GetNumbers();
 
    bool match = true;
 
    QuantumContainer::const_iterator qnr1it = mqnumbers.begin();
    QuantumContainer::const_iterator qnr2it;
 
    // Loop over all quantum numbers in mqnumbers and use their keys
    // to find the corresponding numbers in qnumbers2
    while (match && qnr1it != mqnumbers.end())
    {
        qnr2it = qnumbers2.find(qnr1it->first);
 
        // If one of the two numbers is undefined, it is considered as
        // a match.
        if (qnr2it != qnumbers2.end()
            && !(qnr1it->second).isUndefined()
            && !(qnr2it->second).isUndefined()
            && qnr1it->second != qnr2it->second)
            match = false;
 
        qnr1it++;
    }
 
    return match;
}
 
std::istream& operator>>(std::istream& is, QuantumNumbers& qn)
{
    String name;
    Rational r;
 
    is >> name >> r;
 
    // Define a helper macro to save some typing.
#define INPUT_QUANTUM(ID) \
    if (name == #ID) qn.Set(QN_ ## ID, r)
 
    INPUT_QUANTUM(J);
    else INPUT_QUANTUM(N);
    else INPUT_QUANTUM(S);
    else INPUT_QUANTUM(F);
    else INPUT_QUANTUM(Omega);
    else INPUT_QUANTUM(K1);
    else INPUT_QUANTUM(K2);
    else INPUT_QUANTUM(v1);
    else INPUT_QUANTUM(v2);
    else INPUT_QUANTUM(v3);
    else
    {
        std::ostringstream os;
        os << "Unknown quantum number: " << name << " (" << r << ").";
        throw std::runtime_error(os.str());
    }
 
#undef INPUT_QUANTUM
 
    return is;
}
 
std::ostream& operator<<(std::ostream& os, const QuantumNumbers& qn)
{
    bool first = true;
    // Define a helper macro to save some typing.
#define OUTPUT_QUANTUM(ID) \
    if (!qn[QN_ ## ID].isUndefined()) \
      { if (!first) os << " "; first = false; os << #ID << " " << qn[QN_ ## ID]; }
 
    OUTPUT_QUANTUM(J);
    OUTPUT_QUANTUM(N);
    OUTPUT_QUANTUM(S);
    OUTPUT_QUANTUM(F);
    OUTPUT_QUANTUM(Omega);
    OUTPUT_QUANTUM(K1);
    OUTPUT_QUANTUM(K2);
    OUTPUT_QUANTUM(v1);
    OUTPUT_QUANTUM(v2);
    OUTPUT_QUANTUM(v3);
 
#undef OUTPUT_QUANTUM
 
    return os;
}
 
std::ostream& operator<<(std::ostream& os, const QuantumNumberRecord& qr)
{
    os << "Upper: " << qr.Upper() << " ";
    os << "Lower: " << qr.Lower();
 
    return os;
}