interpolateHypercube.h

Go to the documentation of this file.

/*
  @copyright Russell Standish 2020
  @author Russell Standish
  This file is part of Civita.

  Civita 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.

  Civita 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 Civita.  If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef CIVITA_INTERPOLATE_HYPERCUBE_H
#define CIVITA_INTERPOLATE_HYPERCUBE_H
#include "tensorOp.h"
using namespace std;

namespace civita
{
  class InterpolateHC: public ITensor
  {
    TensorPtr arg;
    Hypercube interimHC;
    std::vector<std::size_t> strides; 
    std::vector<std::size_t> rotation; 
    
    std::vector<std::pair<XVector, std::vector<std::size_t>>> sortedArgHC;
    void sortAndAdd(const XVector&);

    //
    vector<std::size_t> splitAndRotate(std::size_t) const;
    
    struct WeightedIndex
    {
      std::size_t index;
      double weight;
      WeightedIndex(std::size_t index,double weight): index(index), weight(weight) {}
    };

    struct WeightedIndexVector: public vector<WeightedIndex>
    {
      bool argIsOnDestHypercube=false;
    };
  
    vector<WeightedIndexVector> weightedIndices;

    WeightedIndexVector bodyCentredNeighbourhood(std::size_t idx) const;

  public:
    void setArgument(const TensorPtr& a, const ITensor::Args& args={"",0}) override;
    double operator[](std::size_t) const override;
    Timestamp timestamp() const override {return arg? arg->timestamp(): Timestamp();}
  };
  
}

#endif