From bbbebf1fdf3033d8da6bcbef56a56b56fcf6eb01 Mon Sep 17 00:00:00 2001
From: Yingjie Wang <phywyj@gmail.com>
Date: Sun, 7 Dec 2025 18:53:12 -0500
Subject: [PATCH] update: rerun MC integral

---
 src/render.c | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/src/render.c b/src/render.c
index b1e95fc..3f6d85b 100644
--- a/src/render.c
+++ b/src/render.c
@@ -89,6 +89,11 @@ int MC_pixel_render(const System *system, int i, int j, const Spline_data spline
     gsl_monte_miser_init(miser_state);
     gsl_monte_miser_integrate(&F, xl, xu, 2, 100, r, miser_state, &color, &err);
     //printf("%g +- %g ", color, err);
+    if (err > pixel_render_err) {
+      int ncalls = 100*(err/pixel_render_err)*(err/pixel_render_err);
+      if (ncalls > pixel_render_max) ncalls = pixel_render_max;
+      gsl_monte_miser_integrate(&F, xl, xu, 2, ncalls, r, miser_state, &color, &err);
+    }
     rgb[c] = color;
   }
   //printf(")\n");