#include<vector>
#include<iostream>
#include<cmath>
#include<algorithm>

#include "characteristic.hh"
#include "tetrahedron.hh"

double characteristic(const std::vector<float>& p1, const std::vector<float>& p2, const std::vector<float>& p3, const std::vector<float>& p4) {
	double elongation = 0.;
	double planarity = 0.;
	double characteristic = 0.;
	compute_tetrahedron_(p1, p2, p3, p4, elongation, planarity, characteristic);
	return characteristic;
}

double characteristic(const std::vector<double>& p1, const std::vector<double>& p2, const std::vector<double>& p3, const std::vector<double>& p4) {
	double elongation = 0.;
	double planarity = 0.;
	double characteristic = 0.;
	compute_tetrahedron_(p1, p2, p3, p4, elongation, planarity, characteristic);
	return characteristic;
}

double characteristic(const std::vector<long double>& p1, const std::vector<long double>& p2, const std::vector<long double>& p3, const std::vector<long double>& p4) {
	double elongation = 0.;
	double planarity = 0.;
	double characteristic = 0.;
	compute_tetrahedron_(p1, p2, p3, p4, elongation, planarity, characteristic);
	return characteristic;
}

double characteristic(const std::vector<int>& p1, const std::vector<int>& p2, const std::vector<int>& p3, const std::vector<int>& p4) {
	double elongation = 0.;
	double planarity = 0.;
	double characteristic = 0.;
	compute_tetrahedron_(p1, p2, p3, p4, elongation, planarity, characteristic);
	return characteristic;
}

double characteristic(const std::vector<short>& p1, const std::vector<short>& p2, const std::vector<short>& p3, const std::vector<short>& p4) {
	double elongation = 0.;
	double planarity = 0.;
	double characteristic = 0.;
	compute_tetrahedron_(p1, p2, p3, p4, elongation, planarity, characteristic);
	return characteristic;
}