#ifndef BTL_C_INTERFACE_BASE_H
#define BTL_C_INTERFACE_BASE_H

#include "utilities.h"
#include <vector>

template<class real> class c_interface_base
{

public:

  typedef real                      real_type;
  typedef std::vector<real>         stl_vector;
  typedef std::vector<stl_vector >  stl_matrix;

  typedef real* gene_matrix;
  typedef real* gene_vector;

  static void free_matrix(gene_matrix & A, int /*N*/){
    delete[] A;
  }

  static void free_vector(gene_vector & B){
    delete[] B;
  }

  static inline void matrix_from_stl(gene_matrix & A, stl_matrix & A_stl){
    int N = A_stl.size();
    A = new real[N*N];
    for (int j=0;j<N;j++)
      for (int i=0;i<N;i++)
        A[i+N*j] = A_stl[j][i];
  }

  static inline void vector_from_stl(gene_vector & B, stl_vector & B_stl){
    int N = B_stl.size();
    B = new real[N];
    for (int i=0;i<N;i++)
      B[i] = B_stl[i];
  }

  static inline void vector_to_stl(gene_vector & B, stl_vector & B_stl){
    int N = B_stl.size();
    for (int i=0;i<N;i++)
      B_stl[i] = B[i];
  }

  static inline void matrix_to_stl(gene_matrix & A, stl_matrix & A_stl){
    int N = A_stl.size();
    for (int j=0;j<N;j++){
      A_stl[j].resize(N);
      for (int i=0;i<N;i++)
        A_stl[j][i] = A[i+N*j];
    }
  }

  static inline void copy_vector(const gene_vector & source, gene_vector & cible, int N){
    for (int i=0;i<N;i++)
      cible[i]=source[i];
  }

  static inline void copy_matrix(const gene_matrix & source, gene_matrix & cible, int N){
    for (int j=0;j<N;j++){
      for (int i=0;i<N;i++){
        cible[i+N*j] = source[i+N*j];
      }
    }
  }

};

#endif