wyj/blackhole-ray-trace
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

fix: use rmax instead of bmax

Browse Source
This commit is contained in:
Yingjie Wang
2025-12-03 15:58:08 -05:00
parent 455f15d309
commit 0c491d64ce
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/render.c
View File

@@ -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 h = (system->H)*(system->w)/(system->W);
double tanpsi2 = (system->w)*(system->w)+h*h; double tanpsi2 = (system->w)*(system->w)+h*h;
double cotpsi2 = 1/tanpsi2; 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 *x=NULL;
double *y = NULL; double *y = NULL;
int size; 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_interp_accel *acc = gsl_interp_accel_alloc();
gsl_spline *spline = gsl_spline_alloc(gsl_interp_steffen, size); gsl_spline *spline = gsl_spline_alloc(gsl_interp_steffen, size);