From f0238cfb6997c4acfc2bd200de7295f3fa36968f Mon Sep 17 00:00:00 2001 From: Stanislaw Halik Date: Sun, 3 Mar 2019 21:09:10 +0100 Subject: don't index Eigen --- eigen/demos/mix_eigen_and_c/binary_library.cpp | 185 ------------------------- 1 file changed, 185 deletions(-) delete mode 100644 eigen/demos/mix_eigen_and_c/binary_library.cpp (limited to 'eigen/demos/mix_eigen_and_c/binary_library.cpp') diff --git a/eigen/demos/mix_eigen_and_c/binary_library.cpp b/eigen/demos/mix_eigen_and_c/binary_library.cpp deleted file mode 100644 index 15a2d03..0000000 --- a/eigen/demos/mix_eigen_and_c/binary_library.cpp +++ /dev/null @@ -1,185 +0,0 @@ -// This file is part of Eigen, a lightweight C++ template library -// for linear algebra. -// -// Copyright (C) 2009 Benoit Jacob -// -// This Source Code Form is subject to the terms of the Mozilla -// Public License v. 2.0. If a copy of the MPL was not distributed -// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. - -// This C++ file compiles to binary code that can be linked to by your C program, -// thanks to the extern "C" syntax used in the declarations in binary_library.h. - -#include "binary_library.h" - -#include - -using namespace Eigen; - -/************************* pointer conversion methods **********************************************/ - -////// class MatrixXd ////// - -inline MatrixXd& c_to_eigen(C_MatrixXd* ptr) -{ - return *reinterpret_cast(ptr); -} - -inline const MatrixXd& c_to_eigen(const C_MatrixXd* ptr) -{ - return *reinterpret_cast(ptr); -} - -inline C_MatrixXd* eigen_to_c(MatrixXd& ref) -{ - return reinterpret_cast(&ref); -} - -inline const C_MatrixXd* eigen_to_c(const MatrixXd& ref) -{ - return reinterpret_cast(&ref); -} - -////// class Map ////// - -inline Map& c_to_eigen(C_Map_MatrixXd* ptr) -{ - return *reinterpret_cast*>(ptr); -} - -inline const Map& c_to_eigen(const C_Map_MatrixXd* ptr) -{ - return *reinterpret_cast*>(ptr); -} - -inline C_Map_MatrixXd* eigen_to_c(Map& ref) -{ - return reinterpret_cast(&ref); -} - -inline const C_Map_MatrixXd* eigen_to_c(const Map& ref) -{ - return reinterpret_cast(&ref); -} - - -/************************* implementation of classes **********************************************/ - - -////// class MatrixXd ////// - - -C_MatrixXd* MatrixXd_new(int rows, int cols) -{ - return eigen_to_c(*new MatrixXd(rows,cols)); -} - -void MatrixXd_delete(C_MatrixXd *m) -{ - delete &c_to_eigen(m); -} - -double* MatrixXd_data(C_MatrixXd *m) -{ - return c_to_eigen(m).data(); -} - -void MatrixXd_set_zero(C_MatrixXd *m) -{ - c_to_eigen(m).setZero(); -} - -void MatrixXd_resize(C_MatrixXd *m, int rows, int cols) -{ - c_to_eigen(m).resize(rows,cols); -} - -void MatrixXd_copy(C_MatrixXd *dst, const C_MatrixXd *src) -{ - c_to_eigen(dst) = c_to_eigen(src); -} - -void MatrixXd_copy_map(C_MatrixXd *dst, const C_Map_MatrixXd *src) -{ - c_to_eigen(dst) = c_to_eigen(src); -} - -void MatrixXd_set_coeff(C_MatrixXd *m, int i, int j, double coeff) -{ - c_to_eigen(m)(i,j) = coeff; -} - -double MatrixXd_get_coeff(const C_MatrixXd *m, int i, int j) -{ - return c_to_eigen(m)(i,j); -} - -void MatrixXd_print(const C_MatrixXd *m) -{ - std::cout << c_to_eigen(m) << std::endl; -} - -void MatrixXd_multiply(const C_MatrixXd *m1, const C_MatrixXd *m2, C_MatrixXd *result) -{ - c_to_eigen(result) = c_to_eigen(m1) * c_to_eigen(m2); -} - -void MatrixXd_add(const C_MatrixXd *m1, const C_MatrixXd *m2, C_MatrixXd *result) -{ - c_to_eigen(result) = c_to_eigen(m1) + c_to_eigen(m2); -} - - - -////// class Map_MatrixXd ////// - - -C_Map_MatrixXd* Map_MatrixXd_new(double *array, int rows, int cols) -{ - return eigen_to_c(*new Map(array,rows,cols)); -} - -void Map_MatrixXd_delete(C_Map_MatrixXd *m) -{ - delete &c_to_eigen(m); -} - -void Map_MatrixXd_set_zero(C_Map_MatrixXd *m) -{ - c_to_eigen(m).setZero(); -} - -void Map_MatrixXd_copy(C_Map_MatrixXd *dst, const C_Map_MatrixXd *src) -{ - c_to_eigen(dst) = c_to_eigen(src); -} - -void Map_MatrixXd_copy_matrix(C_Map_MatrixXd *dst, const C_MatrixXd *src) -{ - c_to_eigen(dst) = c_to_eigen(src); -} - -void Map_MatrixXd_set_coeff(C_Map_MatrixXd *m, int i, int j, double coeff) -{ - c_to_eigen(m)(i,j) = coeff; -} - -double Map_MatrixXd_get_coeff(const C_Map_MatrixXd *m, int i, int j) -{ - return c_to_eigen(m)(i,j); -} - -void Map_MatrixXd_print(const C_Map_MatrixXd *m) -{ - std::cout << c_to_eigen(m) << std::endl; -} - -void Map_MatrixXd_multiply(const C_Map_MatrixXd *m1, const C_Map_MatrixXd *m2, C_Map_MatrixXd *result) -{ - c_to_eigen(result) = c_to_eigen(m1) * c_to_eigen(m2); -} - -void Map_MatrixXd_add(const C_Map_MatrixXd *m1, const C_Map_MatrixXd *m2, C_Map_MatrixXd *result) -{ - c_to_eigen(result) = c_to_eigen(m1) + c_to_eigen(m2); -} -- cgit v1.2.3