dependencies.hpp

/*
 * dependencies.hpp
 *
 * Created on: 10 Jul 2012
 * Author: Holger Schmitz
 * Email: holger@notjustphysics.com
 *
 * Copyright 2012 Holger Schmitz
 *
 * This file is part of Schnek.
 *
 * Schnek 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.
 *
 * Schnek 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 Schnek.  If not, see <http://www.gnu.org/licenses/>.
 *
 */

#ifndef SCHNEK_DEPENDENCIES_HPP_
#define SCHNEK_DEPENDENCIES_HPP_

#include "variables.hpp"
#include "blockparameters.hpp"

#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wdeprecated-declarations"

#include <boost/variant.hpp>

#pragma GCC diagnostic pop

#include <memory>
#include <map>
#include <set>
#include <list>

namespace schnek {

class DependencyUpdater;

class DependencyMap
{
  private:
    typedef std::set<long> DependencySet;
    struct VarInfo
    {
      pVariable v;
      DependencySet dependsOn;
      DependencySet modifies;
      int counter;
      VarInfo() {}
      VarInfo(pVariable v_, DependencySet dependsOn_, DependencySet modifies_)
        : v(v_), dependsOn(dependsOn_), modifies(modifies_), counter(0) {
        assert(v);
      }
    };

    typedef std::map<long, VarInfo> DepMap;

    typedef std::map<long, VarInfo*> RefDepMap;
    typedef std::shared_ptr<RefDepMap> pRefDepMap;

    typedef std::set<pVariable> VariableSet;
    typedef std::list<pVariable> VariableList;

    DepMap dependencies;
    pBlockVariables blockVars;

    pVariable dummyVar;

    friend class DependencyUpdater;

    void constructMapRecursive(const pBlockVariables vars);
    void constructMap(const pBlockVariables vars);
    void resetCounters();

    void makeUpdateList(const VariableSet &independentVars, const VariableSet &dependentVars, VariableList &updateList);
    pRefDepMap makeUpdatePredecessors(const VariableSet &independentVars, const VariableSet &dependentVars);
    pRefDepMap makeUpdateFollowers(const VariableSet &independentVars, pRefDepMap reverseDeps);
    void makeUpdateOrder(pRefDepMap deps, VariableList &updateList);

  public:
    DependencyMap(const pBlockVariables vars);
    void recreate() { constructMap(blockVars); }
    pBlockVariables getBlockVariables();
    void updateAll();

//    bool hasRoots(pVariable v, pParametersGroup roots);
};

typedef std::shared_ptr<DependencyMap> pDependencyMap;

class DependencyUpdater
{
  private:
    typedef std::set<pParameter> ParameterSet;
    typedef std::set<pVariable> VariableSet;
    typedef std::list<pVariable> VariableList;

    VariableList updateList;
    VariableSet independentVars;
    VariableSet dependentVars;
    ParameterSet dependentParameters;

    pDependencyMap dependencies;
    bool isValid;
  public:
    DependencyUpdater(pDependencyMap dependencies_);
    void addIndependent(pParameter v);
    void addDependent(pParameter v);
    void clearDependent();

    template<size_t rank, template<size_t> class CheckingPolicy>
    void addIndependentArray(Array<pParameter, rank, CheckingPolicy> v)
    { for (size_t i=0; i<rank; ++i) addIndependent(v[i]); }

    template<size_t rank, template<size_t> class CheckingPolicy>
    void addDependentArray(Array<pParameter, rank, CheckingPolicy> v)
    { for (size_t i=0; i<rank; ++i) addDependent(v[i]); }

    void update()
    {
      if (!isValid) {
        dependencies->makeUpdateList(independentVars, dependentVars, updateList);
        isValid = true;
      }
      for(pVariable v: updateList) v->evaluateExpression();
      for(pParameter p: dependentParameters) p->update();
    }
};

typedef std::shared_ptr<DependencyUpdater> pDependencyUpdater;

} // namespace


#endif // SCHNEK_DEPENDENCIES_HPP_