From 0c491d64ce223cabc86e475705b1a8d179967bf3 Mon Sep 17 00:00:00 2001 From: Yingjie Wang Date: Wed, 3 Dec 2025 15:58:08 -0500 Subject: [PATCH] fix: use rmax instead of bmax --- src/render.c | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/render.c b/src/render.c index 18690f3..ebae401 100644 --- a/src/render.c +++ b/src/render.c @@ -29,13 +29,14 @@ int render(System *system, double **buffer, int pixel_render_max, double pixel_r double h = (system->H)*(system->w)/(system->W); double tanpsi2 = (system->w)*(system->w)+h*h; double cotpsi2 = 1/tanpsi2; - double bmax = R0/sqrt(f(R0)*(1+cotpsi2)); + //double bmax = R0/sqrt(f(R0)*(1+cotpsi2)); + double rmax = sqrt(tanpsi2); double *x=NULL; double *y = NULL; int size; - init(bmax, chi_rela_err, &size, &x, &y); + init(rmax, chi_rela_err, &size, &x, &y); gsl_interp_accel *acc = gsl_interp_accel_alloc(); gsl_spline *spline = gsl_spline_alloc(gsl_interp_steffen, size);