unsaturated

condition_number.cpp
View Options

			#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
			}