#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