//:
// \file
// \author Isabel Restrepo
// \date 25-Jan-2011

#include "bvpl_taylor_basis_factory.h"
#include <iostream>
#include <map>
#include <utility>
#ifdef _MSC_VER
#  include "vcl_msvc_warnings.h"
#endif

//: Constructor from filename
bvpl_taylor_basis_factory::bvpl_taylor_basis_factory(std::string filename)
{
  //initialize variables
  filename_ = std::move(filename);
  angle_ = 0.0f;
  rotation_axis_ = canonical_rotation_axis_;
  parallel_axis_ = canonical_parallel_axis_;

  //this skernel is symmetric around main axis
  angular_resolution_=0;

  //create the default kernel
  create_canonical();
}

void bvpl_taylor_basis_factory::create_canonical()
{
  typedef bvpl_kernel_dispatch dispatch;

  std::ifstream ifs(filename_.c_str());

  //set the dimension of the 3-d bounding box containing the kernels
  if (!ifs.eof())
  {
    ifs >> min_point_;
    ifs >> max_point_;
  }
  std::cout << " Max point: " << max_point_ << " Min point: " << min_point_ << std::endl;

  while (!ifs.eof())
  {
    vgl_point_3d<float> this_loc;
    float weight;
    ifs >> this_loc;
    ifs >> weight;
    canonical_kernel_.emplace_back(this_loc, dispatch(weight));
  }

  //set the current kernel
  kernel_ = canonical_kernel_;
  factory_name_ = name();
}


/************************bvpl_taylor_basis_loader**************************************/

void bvpl_taylor_basis_loader::create_basis(std::map<std::string, bvpl_kernel_sptr> &taylor_basis)
{
  //get filenames, iterate through files reading the kernels
  std::vector<std::string> filenames;
  files(filenames);

  auto file_it =filenames.begin();

  for (; file_it != filenames.end(); file_it++)
  {
    std::string filename = path_ + '/' + *file_it + ".txt";
    std::cout << "Reading kernel file : " << filename << std::endl;
    bvpl_taylor_basis_factory factory(filename);
    bvpl_kernel_sptr kernel = new bvpl_kernel(factory.create());
    //kernel->print();
    taylor_basis.insert(std::pair<std::string, bvpl_kernel_sptr>(*file_it, kernel));
  }
}

//: Kernels needed for 2 degree approximation of 3D functions
void bvpl_taylor_basis_loader::files(std::vector<std::string> &filenames)
{
  if (degree_ == 2)
  {
    filenames.emplace_back("I0");
    filenames.emplace_back("Ix");
    filenames.emplace_back("Iy");
    filenames.emplace_back("Iz");
    filenames.emplace_back("Ixx");
    filenames.emplace_back("Iyy");
    filenames.emplace_back("Izz");
    filenames.emplace_back("Ixy");
    filenames.emplace_back("Ixz");
    filenames.emplace_back("Iyz");
  }
}