ROperator_Concat.hxx

Go to the documentation of this file.

#ifndef TMVA_SOFIE_ROPERATOR_Concat
 #define TMVA_SOFIE_ROPERATOR_Concat
 
 
 #include "TMVA/SOFIE_common.hxx"
 #include "TMVA/ROperator.hxx"
 #include "TMVA/RModel.hxx"
 
 #include <sstream>
 #include <algorithm>
 #include <iterator>
 #include <iomanip>
 #include <limits>
 
 namespace TMVA{
 namespace Experimental{
 namespace SOFIE{
 
     class ROperator_Concat final : public ROperator
     {
     private:
         int fAxis=0;
         int fnewAxis=0;
         std::vector<std::string> fInputs;
         std::string fOutput;
         std::vector<Dim>fOutputShape;
         std::vector<Dim> fOutputShapeData; // in case output is a shape tensor we store here the output shape value data (can be parametric)
         std::vector<std::vector<Dim>> fInputShapes;
 
     public:
 
         ROperator_Concat(){}
         ROperator_Concat(std::vector<std::string> inputs, int axis, int newAxis, std::string output):
         fAxis(axis), fnewAxis(newAxis), fOutput(UTILITY::Clean_name(output)) {
            fInputs.reserve(inputs.size());
            for (auto & name : inputs)
               fInputs.push_back(UTILITY::Clean_name(name));
 
         fInputTensorNames.resize(fInputs.size());
         std::transform(fInputs.begin(), fInputs.end(), fInputTensorNames.begin(),
                   [](const std::string& s) -> std::string_view { return s; });
         fOutputTensorNames = { fOutput };
         }
 
         std::vector<ETensorType> TypeInference(std::vector<ETensorType> input) override {
             return input;
         }
 
         // get shape of output given inputs. It is going to be called after initialized
         std::vector<std::vector<size_t>> ShapeInference(std::vector<std::vector<size_t>> inputs) override {
             std::vector<std::vector<size_t>> ret(1);
            // treat negative axis case
            if (fAxis<0) {
               fAxis = inputs[0].size()+fAxis;
            }
            if (fAxis < 0 || fAxis >= (int) inputs[0].size())
               throw std::runtime_error("TMVA SOFIE Concat Op - invalid axis value ");
 
            int concat_dim=0;
            // case of Concat (fNewAxis = 0) and not ConcatFromSequence
            if(fnewAxis == 0){
               for (size_t i = 0; i < inputs.size(); i++) {
                  if (i > 0 && inputs[i].size() != inputs[i - 1].size())
                     throw std::runtime_error("TMVA SOFIE Concat Op - input tensors have different shapes " +
                                              ConvertShapeToString(inputs[i]) + " and " + ConvertShapeToString(inputs[i - 1]));
                  for (size_t iaxis = 0; iaxis < inputs[i].size(); iaxis++) {
                     if ((int)iaxis == fAxis)
                        concat_dim += inputs[i][iaxis];
                     else if (i > 0 && inputs[i][iaxis] != inputs[i - 1][iaxis])
                        throw std::runtime_error("TMVA SOFIE Concat Op - input tensors have wrong shapes " +
                                                 ConvertShapeToString(inputs[i]) + " and " +
                                                 ConvertShapeToString(inputs[i - 1]));
                  }
               }
 
               // output shape
               ret[0] = inputs[0];
               ret[0][fAxis] = concat_dim;
            }
            std::vector<int> stack;
            // case ConCatFromSequence
            if(fnewAxis == 1){
               for(size_t i = 0; i < inputs.size(); i++) {
                  if (i > 0 && inputs[i].size() != inputs[i-1].size() )
                  throw std::runtime_error("TMVA SOFIE Concat Op - input tensors have different shapes " + fInputs[i] + " : " +
                     ConvertShapeToString(inputs[i]) + " and " + fInputs[i-1] + " : " + ConvertShapeToString(inputs[i-1]));
                  for (size_t iaxis = 0; iaxis < inputs[i].size(); iaxis++) {
                     if ((int) iaxis == fAxis)
                        stack.push_back(inputs[i][iaxis]);
                     else
                     if (i> 0 && inputs[i][iaxis] != inputs[i-1][iaxis])
                        throw std::runtime_error("TMVA SOFIE Concat Op - input tensors have wrong shapes " +
                        ConvertShapeToString(inputs[i]) + " and " + ConvertShapeToString(inputs[i-1]));
                  }
 
               }
               for(auto it:stack)
               ret[0].push_back(it);
            }
 
            return ret;
         }
 
         // get shape of output given inputs. It is going to be called after initialized
         std::vector<Dim> ShapeInference(const std::vector<std::vector<Dim>> & inputs, const RModel & model) {
            std::vector<Dim> ret(inputs[0].size());
            // treat negative axis case
            if (fAxis<0) {
               fAxis = inputs[0].size()+fAxis;
            }
            if (fAxis < 0 || fAxis >= (int) inputs[0].size())
               throw std::runtime_error("TMVA SOFIE Concat Op - invalid axis value ");
 
            Dim concat_dim;
            if(fnewAxis == 0){
               for (size_t i = 0; i < inputs.size(); i++) {
                  if (i > 0 && inputs[i].size() != inputs[i - 1].size())
                     throw std::runtime_error("TMVA SOFIE Concat Op - input tensors have different shapes " + fInputs[i] + " : " +
                                              ConvertDimShapeToString(inputs[i]) + " and " + fInputs[i-1] + " : " + ConvertDimShapeToString(inputs[i - 1]));
                  for (size_t iaxis = 0; iaxis < inputs[i].size(); iaxis++) {
                     if ((int)iaxis == fAxis) {
                        // support both integer and params shape for the concatenation axis
                        if (concat_dim.param.empty() && concat_dim.dim == 0)
                           concat_dim = inputs[i][iaxis];
                        else if (inputs[i][iaxis].isParam || concat_dim.isParam) {
                           concat_dim =
                              Dim{ concat_dim.GetVal() + std::string(" + ") + inputs[i][iaxis].GetVal(),
                                 static_cast<size_t>(-1)};
                        } else {
                           concat_dim = Dim { concat_dim.dim + inputs[i][iaxis].dim };
                        }
                     }
                     else if (i == 0) {
                        ret[iaxis] = inputs[i][iaxis];
                     }
                     else if ((!inputs[i][iaxis].isParam && !ret[iaxis].isParam) && (inputs[i][iaxis].dim != ret[iaxis].dim)) {
                        throw std::runtime_error("TMVA SOFIE Concat Op - input tensors have wrong shapes " +
                                                 ConvertDimShapeToString(inputs[i]) + " and " +
                                                 ConvertDimShapeToString(inputs[i - 1]));
                     }
                     else if (!inputs[i][iaxis].isParam && ret[iaxis].isParam){
                        // if shape is not parametric use it
                        ret[iaxis] = inputs[i][iaxis];
                     }
                     else if (inputs[i][iaxis].isParam && ret[iaxis].isParam) {
                        // check which parameter is first in RModel list
                        auto & dimNames = model.GetDimShapeNames();
                        auto p1 = std::find(dimNames.begin(), dimNames.end(), inputs[i][iaxis].param);
                        auto p2 = std::find(dimNames.begin(), dimNames.end(), ret[iaxis].param);
                        if (p1 < p2) ret[iaxis] = inputs[i][iaxis];
                     }
 
                  }
                  // add parenthesis in case is an expression
                  if (concat_dim.isParam && concat_dim.dim == static_cast<size_t>(-1))
                     concat_dim =  Dim{ std::string("(") + concat_dim.GetVal() +  std::string(")"), concat_dim.dim };
               }
 
               // output shape for concatenated axis
               ret[fAxis] = concat_dim;
 
            }
            // case of stacking (not supported yet)
            // here we need to check that input shapes are the same
            // for example for fAxis == 0
            // output shapes: [inputs.size(), inputs[0][0], inputs[0][1],....]
            if(fnewAxis == 1){
               throw std::runtime_error("TMVA SOFIE Concat Op - stacking (i.e. COncatFromSequence with new_axis=1) is not supported ");
            }
            return ret;
         }
 
         void Initialize(RModel& model) override {
            std::vector<std::vector<size_t>> inputIntShapes;
            for (auto &it : fInputs) {
               if (model.CheckIfTensorAlreadyExist(it) == false) {
                  throw std::runtime_error("TMVA SOFIE Concat Op Input Tensor " + it + " is not found in model");
               }
               fInputShapes.push_back(model.GetDimTensorShape(it));
               if (!model.IsDynamicTensor(it)) {
                  inputIntShapes.push_back(ConvertShapeToInt(fInputShapes.back()));
               }
            }
            if (inputIntShapes.size() == fInputs.size()) {
               // if all input shapes are static we can compute output shape at initialization time
               auto outputIntShape = ShapeInference(inputIntShapes)[0];
               fOutputShape = ConvertShapeToDim(outputIntShape);
               if (model.Verbose())
                  std::cout << "Initialize Concat operator with defined inputs shapes, "
                           << "output has shape " << ConvertShapeToString(outputIntShape) << std::endl;
 
            } else {
               // if at least one input shape is dynamic we need to compute output shape using the symbolic expression for the dimensions
               fOutputShape = ShapeInference(fInputShapes, model);
               if (model.Verbose())
                  std::cout << "Initialize Concat operator with dynamic inputs shapes, "
                           << "output has shape " << ConvertDimShapeToString(fOutputShape) << std::endl;
            }
 
            // check if concat has constant inputs , axis 0(concat contigous memory and type is integer)
            bool isOutputShape = false;
 
            // if (model.GetTensorType(fInputs[0]) == ETensorType::INT64 && fAxis == 0) {
            fIsOutputConstant = true;
            isOutputShape = true;
 
            for (auto &input : fInputs) {
               if (model.IsDynamicTensor(input)) {
                  fIsOutputConstant = false;
                  isOutputShape = false;
                  break;
               }
               if (!model.IsInitializedTensor(input)) {
                  if (model.IsShapeTensor(input)) {
                     // if it is a shape tensor we can have constant output if the shapes are defined)
                     auto shapeData = model.GetShapeTensorValues(input);
                     bool isShapeFullyDefined = ConvertShapeToInt(shapeData).size() == shapeData.size();
                     if (!isShapeFullyDefined) {
                        fIsOutputConstant = false;
                     } else {
                        // if shape is fully defined we can consider output as constant and we can compute the output
                        // shape at initialization time
                        fIsOutputConstant = fIsOutputConstant && true;
                     }
                     // inputs are then shape tensors and output is a shape tensor
                     isOutputShape = true;
                  } else {
                     // case of standard intermediate tensor
                     fIsOutputConstant = false;
                     isOutputShape = false;
                     break;
                  }
               } else {
                  fIsOutputConstant = fIsOutputConstant && true;
               }
            }
            //}
 
            if (fIsOutputConstant) {
               auto outputShape = ConvertShapeToInt(fOutputShape); // conversion must be possible
               std::vector<int64_t> outputData(ConvertShapeToLength(outputShape));
               size_t offset = 0;
               for (auto &input : fInputs) {
                  auto inputData = static_cast<int64_t *>(model.GetInitializedTensorData(input).get());
                  auto inputShape = model.GetTensorShape(input); // shape is not dynamic if it is constant
                  size_t inputLength = ConvertShapeToLength(inputShape);
                  std::copy(inputData, inputData + inputLength, outputData.begin() + offset);
                  offset += inputLength;
                  // the data of the input tensor don't need to be written in the generated code and data file
                  model.SetNotWritableInitializedTensor(input);
               }
               model.AddConstantTensor<int64_t>(fOutput, outputShape, outputData.data());
               if (model.Verbose()) {
                  std::cout << "output of Concat is a constant tensor " << ConvertShapeToString(outputShape) << " : "
                            << ConvertValuesToString(outputData) << " (constant)" << std::endl;
               }
            } else if (isOutputShape) {
               auto outputShape = ConvertShapeToInt(fOutputShape); // conversion must be possible
               if (outputShape.size() != 1)
                  throw std::runtime_error("TMVA SOFIE Concat Op - output shape for shape tensor must have rank 1");
               // output shape is a rank 1 tensor with size equal to the output rank
               std::vector<Dim> outputData(outputShape[0]);
               size_t offset = 0;
               for (auto &input : fInputs) {
                  std::vector<Dim> inputData;
                  auto inputShape = model.GetTensorShape(input);         // shape is not dynamic
                  size_t inputLength = ConvertShapeToLength(inputShape); // shape can be a scalar
                  if (model.IsShapeTensor(input)) {
                     inputData = model.GetShapeTensorValues(input);
                  } else if (model.IsInitializedTensor(input)) {
                     inputData.resize(inputLength);
                     auto intData = static_cast<int64_t *>(model.GetInitializedTensorData(input).get());
                     for (size_t i = 0; i < inputData.size(); i++)
                        inputData[i] = Dim{static_cast<size_t>(intData[i])};
                  } else {
                     // this should not happen
                     throw std::runtime_error("TMVA SOFIE Concat Operator- invalid tensor input " + input +
                                              " for shape output type");
                  }
                  std::copy(inputData.begin(), inputData.end(), outputData.begin() + offset);
                  offset += inputLength;
               }
               // add output tensor
               model.AddShapeTensor(fOutput, outputData, false); // cannot be a  scalar
               fOutputShapeData = outputData;
               if (model.Verbose()) {
                  std::cout << "output of Concat is a shape tensor " << ConvertShapeToString(outputShape) << " : "
                            << ConvertDimShapeToString(outputData) << " (shape)" << std::endl;
               }
               fIsOutputParamShape = true;
            }
            if (!fIsOutputConstant && !fIsOutputParamShape) {
               model.AddIntermediateTensor(fOutput, model.GetTensorType(fInputs[0]), fOutputShape);
               if (model.Verbose()) {
                  std::cout << "Concat ---> " << fOutput << " " <<  ConvertDimShapeToString(fOutputShape) << std::endl;
               }
            }
         }
 
         std::string Generate(std::string opName) override {
            opName = "op_" + opName;
            std::stringstream out;
            out<<"\n//--------- Concat " << opName << " --> " << fOutput << "  " << ConvertDimShapeToString(fOutputShape) << "\n";
 
            if (fIsOutputConstant) return out.str();
 
            if (fIsOutputParamShape) {
               // output is a shape tensor defined by the concatenation of the input shapes
               out << "// output is a shape tensor defined by the concatenation of the input shapes\n";
               for (int i = 0; i < static_cast<int>(fOutputShape
                  [0].dim); i++) {
                  out << SP << "tensor_" << fOutput << "[" << i << "] = " << fOutputShapeData[i] << ";\n";
               }
               return out.str();
            }
            // special case when memory is contiguous
            bool hasShapeOnes = true;
            for(int i = 0; i<fAxis; ++i){
               if(fInputShapes[0][i].dim !=1){
                  hasShapeOnes = false;
                  break;
               }
            }
            if (fAxis == 0 || hasShapeOnes) {
               std::string offset;
               for(size_t i=0; i<fInputs.size(); ++i) {
                  auto length = ConvertDimShapeToLength(fInputShapes[i]);
                  out << SP << "TMVA::Experimental::SOFIE::Copy(tensor_" << fOutput;
                  if (i > 0)
                     out << offset;
                  offset += " + " + length;
                  out << ", " << "tensor_" << fInputs[i] << ", " + length << ");\n";
               }
            }
            else {
 
               std::vector<Dim> outStride = UTILITY::ComputeStrideFromShape(fOutputShape);
               std::vector<std::vector<Dim>> inStrides(fInputs.size());
               int idx = 0;
               for ( auto &s : inStrides) {
                  s = UTILITY::ComputeStrideFromShape(fInputShapes[idx]);
                  idx++;
               }
               for (int i = 0; i < fAxis; ++i) {
                  // loop on dimensions
                  out << SP << "for (size_t i" << i << " = 0; i" << i << " < " << fOutputShape[i].GetVal() << "; ++i" << i <<") {\n";
               }
 
               out << SP << SP << SP << "int idxOut = ";
               for (int k = 0; k < fAxis; k++) {
                  if (k > 0) out << " + ";
                  out << outStride[k].GetVal() << "*i" << k;
               }
               out << ";\n";
 
               for (size_t j = 0; j < fInputs.size(); j++) {
                  if (j>0)
                  out << SP << SP << SP << "idxOut += " << inStrides[j-1][fAxis-1].GetVal() << ";\n";
                  out << SP << SP << SP << "int idxIn" << j <<" = ";
                  for (int k = 0; k < fAxis; k++) {
                     if (k > 0) out << " + ";
                     out << inStrides[j][k].GetVal() << "*i" << k;
                  }
                  out << ";\n";
                  out << SP << SP << SP << "for (size_t iC = 0; iC < " << inStrides[j][fAxis-1].GetVal() << "; ++iC) {\n";
                  out << SP << SP << SP << SP << "tensor_" << fOutput << "[idxOut+iC] = tensor_" << fInputs[j] << "[idxIn" << j << "+iC];\n";
                  out << SP << SP << SP << "}\n";
               // concatenate the axis values
               }
                for (int i = 0; i < fAxis; ++i) {
                    out << SP << "}\n";
                }
            }
 
            return out.str();
         }
     };
 }//SOFIE
 }//Experimental
 }//TMVA
 
 #endif //TMVA_SOFIE_ROPERATOR_CONCAT