vpr/src/base/stats.h

#pragma once
#include <vector>
#include <limits>
#include <algorithm>
#include "vpr_types.h"

void routing_stats(bool full_stats, enum e_route_type route_type,
		int num_rr_switch, t_segment_inf * segment_inf, int num_segment,
		float R_minW_nmos, float R_minW_pmos,
        float grid_logic_tile_area,
		enum e_directionality directionality, int wire_to_ipin_switch,
		bool timing_analysis_enabled,
		vtr::vector_map<ClusterNetId, float *> &net_delay
#ifdef ENABLE_CLASSIC_VPR_STA
        , t_slack * slacks, const t_timing_inf &timing_inf
#endif
        );

void print_wirelen_prob_dist();

void print_lambda();

void get_num_bends_and_length(ClusterNetId inet, int *bends, int *length, int *segments);

int count_netlist_clocks();

// template functions must be defined in header, or explicitely instantiated in definition file (defeats the point of template)
template <typename T>
double linear_regression_vector(const std::vector<T>& vals, size_t start_x = 0) {
	// returns slope; index is x, val is y
	size_t n {vals.size() - start_x};

	double x_avg {0}, y_avg {0};
	for (size_t x = start_x; x < vals.size(); ++x) {
		x_avg += x;
		y_avg += vals[x];
	} 
	x_avg /= (double) n;
	y_avg /= (double) n;

	double numerator = 0, denominator = 0;
	for (size_t x = start_x; x < vals.size(); ++x) {
		numerator += (x - x_avg) * (vals[x] - y_avg);
		denominator += (x - x_avg) * (x - x_avg);
	}

	if (denominator == 0) return std::numeric_limits<double>::max();
	return numerator / denominator;
}