eigen-2.0.1: cbits/eigen-proxy.h

#define EIGEN_MPL2_ONLY
#define EIGEN2_SUPPORT
#define EIGEN_NO_EIGEN2_DEPRECATED_WARNING
#define EIGEN_INITIALIZE_MATRICES_BY_ZERO
void eigen_assert_fail(const char* condition, const char* function, const char* file, int line);
#define eigen_assert(x) do {\
    if (!(x)) eigen_assert_fail(#x, __PRETTY_FUNCTION__, __FILE__, __LINE__);\
} while(false)
#include <Eigen/Core>

using namespace Eigen;

extern "C" {

const char* eigen_add(int, void*, int, int, const void*, int, int, const void*, int, int);
const char* eigen_sub(int, void*, int, int, const void*, int, int, const void*, int, int);
const char* eigen_mul(int, void*, int, int, const void*, int, int, const void*, int, int);
const char* eigen_norm(int, void*, const void*, int, int);
const char* eigen_squaredNorm(int, void*, const void*, int, int);
const char* eigen_blueNorm(int, void*, const void*, int, int);
const char* eigen_hypotNorm(int, void*, const void*, int, int);
const char* eigen_sum(int, void*, const void*, int, int);
const char* eigen_prod(int, void*, const void*, int, int);
const char* eigen_mean(int, void*, const void*, int, int);
const char* eigen_trace(int, void*, const void*, int, int);
const char* eigen_determinant(int, void*, const void*, int, int);
const char* eigen_diagonal(int, void*, int, int, const void*, int, int);
const char* eigen_transpose(int, void*, int, int, const void*, int, int);
const char* eigen_normalize(int, void*, int, int);
const char* eigen_random(int, void*, int, int);
const char* eigen_identity(int, void*, int, int);
const char* eigen_inverse(int, void*, int, int, const void*, int, int);
const char* eigen_conjugate(int, void*, int, int, const void*, int, int);
const char* eigen_adjoint(int, void*, int, int, const void*, int, int);
bool eigen_initParallel();
int eigen_getNbThreads();
void eigen_setNbThreads(int);

enum Decomposition {
	PartialPivLU,
	FullPivLU,
	HouseholderQR,
	ColPivHouseholderQR,
	FullPivHouseholderQR,
	LLT,
	LDLT,
	JacobiSVD,
	// SelfAdjointEigenSolver,
	// ComplexEigenSolver,
	// EigenSolver,
	// GeneralizedSelfAdjointEigenSolver
};

const char* eigen_rank(int, Decomposition d, int*, const void*, int, int);
const char* eigen_kernel(int, Decomposition d, void**, int*, int*, const void*, int, int);
const char* eigen_image(int, Decomposition d, void**, int*, int*, const void*, int, int);

const char* eigen_solve(int, Decomposition d, // Ax=b 
	void*, int, int, // x
	const void*, int, int, // A
	const void*, int, int); // b


const char* eigen_relativeError(int, void* e, 
	const void*, int, int, // x
	const void*, int, int, // A
	const void*, int, int); // b;


} // end extern "C"

#define GUARD_START try { assert(inited); do {
#define GUARD_END } while(false); return 0; } catch (const std::exception& ex) { return strdup(ex.what()); }