Added erorr bars

Author: EduMenges

4/generate_graphs.py

@@ -3,25 +3,43 @@
 import numpy as np
 from scipy.optimize import curve_fit
 
+
 def model(x, a, b):
     return a * (x ** 2) + b * (x ** 2) * np.log(x)
 
 
 if __name__ == "__main__":
     df = pd.read_csv("results.csv")
 
-    x = df["n"]
-    y = df["time"]
-    plt.scatter(x, y)
+    grouped = df.groupby("n")["time"]
+
+    # Média (e desvio‑padrão) do tempo para cada tamanho
+    x_mean = grouped.mean().index.values.astype(float)
+    y_mean = grouped.mean().values
+    y_std = grouped.std().values
+
+    params, covariance = curve_fit(model, x_mean, y_mean)
+    a_fit, b_fit = params
+
+    y_fit = model(x_mean, a_fit, b_fit)
 
-    params, covariance = curve_fit(model, x, y)
+    ss_res = np.sum((y_mean - y_fit) ** 2)
+    ss_tot = np.sum((y_mean - np.mean(y_mean)) ** 2)
+    r_squared = 1 - ss_res / ss_tot
 
-    a_fit = params[0]
-    b_fit = params[1]
-    y_fit = model(x.unique(), a_fit, b_fit)
+    rmse = np.sqrt(np.mean((y_mean - y_fit) ** 2))
 
-    plt.plot(x.unique(), y_fit, label="Regressão não linear de quadrados mínimos", color="red")
+    plt.scatter(df["n"], df["time"], label="Observações individuais")
+
+    plt.errorbar(x_mean, y_mean, yerr=y_std, fmt="o", color="black",
+                 capsize=4, label="Média ± desvio‑padrão")
+
+    x_fit = np.linspace(x_mean.min(), x_mean.max(), 300)
+    y_fit = model(x_fit, a_fit, b_fit)
+    plt.plot(x_fit, y_fit, color="red", label=f"Ajuste (R²={r_squared:.4f})")
 
     plt.xlabel("Tamanho da matriz")
     plt.ylabel("Tempo (ms)")
+    plt.legend()
+
     plt.savefig("graphs/n_time.eps")

4/report.tex

@@ -9,6 +9,7 @@
 \usepackage{graphicx}
 \usepackage{fontspec}
 \usepackage{float}
+%\usepackage{listings}
 \setlength{\parskip}{10pt}
 \setmainfont{texgyreheros}[
     UprightFont = *-regular,
@@ -44,12 +45,13 @@ \section{Detalhes de implementação}
 \section{Ambiente de teste}
 
 A máquina de teste possui Linux Mint 22.1, processador Intel i5-12450H e 16GB de memória RAM.
-O compilador utilizado foi o GCC, versão 14.2, com todas as otimizações padrões habilitadas, e a implementação da biblioteca padrão do C++ é a libstdc++ versão 13.
+O compilador utilizado foi o GCC, versão 14.2, com todas as otimizações padrões habilitadas, e a implementação da biblioteca padrão do C++ é a libstdc++ versão 14.
 
 \section{Plano de teste}
 
 O gerador utilizado para gerar os grafos toma uma entrada $n$ e gera uma matriz $N \times N$ com valores no intervalo $[0, n \cdot n]$.
-Foram gerados e mensurados $10$ testes para cada $n$ pertencente ao intervalo $[1000, 18000]$, incrementando $1000$ a cada passo.
+
+Foram gerados $10$ testes para cada $n$ pertencente ao intervalo $[1000, 20000]$, incrementando $1000$ a cada passo. Para cada teste, mensurou-se o tempo de execução.
 
 \section{Resultados}
 

4/test.sh

@@ -0,0 +1,48 @@
+#!/usr/bin/env bash
+generator="/home/menges/repos/advanced_algorithms/cmake-build-gcc/bin/generator"
+target="/home/menges/repos/advanced_algorithms/cmake-build-gcc/bin/hungarian_algorithm"
+
+# Your commands here
+
+sizes=(
+    1000
+    2000
+    3000
+    4000
+    5000
+    6000
+    7000
+    8000
+    9000
+    10000
+    11000
+    12000
+    13000
+    14000
+    15000
+    16000
+    17000
+    18000
+    19000
+    20000
+)
+
+echo "n,time,result"
+
+for i in {1..10}; do
+        for size in "${sizes[@]}"; do
+        (
+            $generator $size >"tests/test_${size}"
+
+            start=$(date +%s%N)
+
+            result=$($target <"tests/test_${size}")
+
+            end=$(date +%s%N)
+            elapsed=$(((end - start) / 1000000))
+
+            echo "$size,$elapsed,$result"
+        ) &
+        done
+    wait
+done