c_interface_base.h
68bb631292e1master

c_interface_base.h
View Options


			#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+Nj] = source[i+Nj];
			}
			}
			}

			};

			#endif