SpinAmplitudeCache.h

Go to the documentation of this file.

/*  YAP - Yet another PWA toolkit
    Copyright 2015, Technische Universitaet Muenchen,
    Authors: Daniel Greenwald, Johannes Rauch

    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 3 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/>.
 */


#ifndef yap_SpinAmplitudeCache_h
#define yap_SpinAmplitudeCache_h

#include "fwd/Spin.h"

#include "SpinAmplitude.h"
#include "UnitSpinAmplitude.h"
#include "WeakPtrCache.h"

#include <memory>

namespace yap {

class Model;

class SpinAmplitudeCache :
    public WeakPtrCache<SpinAmplitude>
{
public:

    SpinAmplitudeCache(Model* model = nullptr) :
        WeakPtrCache(), Model_(model) {}

    bool equal(const std::shared_ptr<SpinAmplitude>& A, const std::shared_ptr<SpinAmplitude>& B) const override
    { return (A.get() == B.get()) or A->equalTo(*B); }

    using WeakPtrCache::find;

    weak_ptr_type find(unsigned two_J, const SpinVector& two_j, unsigned L, unsigned two_S) const;

    std::shared_ptr<SpinAmplitude> spinAmplitude(unsigned two_J, const SpinVector& two_j, unsigned L, unsigned two_S);

    bool consistent() const;

    friend class Model;

protected:

    void setModel(Model& model);

    std::shared_ptr<SpinAmplitude> unit(Model& m, unsigned two_J, const SpinVector& two_j, unsigned L, unsigned two_S) const
    { return std::shared_ptr<SpinAmplitude>(new UnitSpinAmplitude(m, two_J, two_j, L, two_S)); }

private:

    virtual std::shared_ptr<SpinAmplitude> create(Model& m, unsigned two_J, const SpinVector& two_j, unsigned L, unsigned two_S) const = 0;

    Model* Model_;

};

}

#endif