condition_number.cpp
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include "catch.hpp"

	#include <Eigen/Dense>
	#include "micpsolver/estimate_cond_number.hpp"

	TEST_CASE("Condition number") {

	SECTION("Lower matrix") {
	Eigen::MatrixXd mat(4,4);
	mat << 1, 0, 0 , 0,
	1, 1, 0, 0,
	1, 1, 1, 0,
	1, 1, 1, 1;
	double cond = specmicp::micpsolver::estimate_condition_number(
	mat.triangularView<Eigen::Lower>());
	CHECK(cond == 4);
	}

	SECTION("Upper matrix") {
	Eigen::MatrixXd mat(4,4);
	mat << 1, 0, 0 , 0,
	1, 1, 0, 0,
	1, 1, 1, 0,
	1, 1, 1, 1;
	double cond = specmicp::micpsolver::estimate_condition_number(
	mat.transpose().triangularView<Eigen::Upper>());
	CHECK(cond == 4);
	}

	// TODO : more tests
	}