---

more/phys/skyrme.h

Go to the documentation of this file.

//  more/skyrme.h -- Skyrme parameters of the effective nuclear interaction
//  Copyright (C) 2001--2002  Petter Urkedal (petter.urkedal@matfys.lth.se)
//
//  This file 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 file 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
//
//  As a special exception, you may use this file as part of a free
//  software library without restriction.  Specifically, if other files
//  instantiate templates or use macros or inline functions from this
//  file, or you compile this file and link it with other files to
//  produce an executable, this file does not by itself cause the
//  resulting executable to be covered by the GNU General Public
//  License.  This exception does not however invalidate any other
//  reasons why the executable file might be covered by the GNU General
//  Public License.
//
//  $Id: skyrme.h,v 1.1 2002/05/30 18:01:40 petter_urkedal Exp $


#ifndef MORE_PHYS_SKYRME_H
#define MORE_PHYS_SKYRME_H

#include <more/bits/config.h>
#include <stdexcept>
#include <iosfwd>
#include <string>
#include <more/io/syncstream.h>
#include <list>

namespace more {
namespace phys {

  extern std::string skyrme_datadir;

  /** \class skyrme_interaction skyrme.h more/phys/skyrme.h
   **
   ** This class holds a Skyrme type parametrization of the nuclear
   ** interaction. It provides read and write methods, the former
   ** which may be applied to the supplied siii.skyrme, sgii.skyrme,
   ** etc. files. */
  template<typename T>
    struct skyrme_interaction
    {
        typedef T real_type;

        T get_t0() const { return t0; }
        T get_t1() const { return t1; }
        T get_t2() const { return t2; }
        T get_t3() const { return t3; }
        T get_x0() const { return x0; }
        T get_x1() const { return x1; }
        T get_x2() const { return x2; }
        T get_x3() const { return x3; }
        T get_alpha() const { return alpha; }

        T get_W0() const
        {
            if (!isospin_dependent_LS_p())
                return w;
            else
                throw std::logic_error("No W0 term in interaction.");
        }

        T get_W1() const
        {
            if (isospin_dependent_LS_p())
                return w;
            else
                throw std::logic_error("No W1 term in interaction.");
        }

        T get_W2() const
        {
            if (isospin_dependent_LS_p())
                return w2;
            else
                throw std::logic_error("No W2 term in interaction.");
        }

        bool isospin_dependent_LS_p() const { return m_force == 10; }

        void sync(more::io::syncstream&);

        /** Load interaction named \c name.  */
        void load(std::string name);

        /** A name suitable for printing. */
        std::string print_name() const { return m_name; }

        /** A name suitable for mangling into file names. */
        std::string path_comp_name() const { return m_fname; }

        friend std::istream&
        operator>>(std::istream& is, skyrme_interaction& V)
        { V.read(is); return is; }

        friend std::ostream&
        operator<<(std::ostream& os, const skyrme_interaction& V)
        { V.write(os); return os; }

      private:
        std::string m_name, m_fname;
        int m_force;
        real_type t0, x0, t1, x1, t2, x2, t3, x3, alpha, w, w2;

        void read(std::istream& sm);
        void write(std::ostream& sm) const;

        static void get_available();
        static bool s_done_init;
    };

  /** \class skyrme_name_list skyrme.h more/phys/skyrme.h
   **
   ** A list of the names of the skyrme interactions.  This is an
   ** empty structure which gives access to static data.  The default
   ** constructor makes sure the static data is constructed. */
  struct skyrme_name_list
  {
    private:
      typedef std::list<std::string> subcontainer;
    public:
      typedef subcontainer::const_iterator iterator;

      skyrme_name_list();

      static std::string name_to_path_comp(std::string name);

      iterator begin() const { return s_names.begin(); }
      iterator end() const { return s_names.end(); }

    private:
      static bool               s_done_init;
      static subcontainer       s_names;
  };

}}

#endif

---