From f08d661a5f959fc54cd3409c4c618324cdd31be4 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Wed, 4 Jul 2018 12:25:36 -0400 Subject: fix a bug in get_solid_angle() --- solid_angle.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/solid_angle.c b/solid_angle.c index cbc71f1..6959e2d 100644 --- a/solid_angle.c +++ b/solid_angle.c @@ -20,7 +20,7 @@ double get_solid_angle(double *pos, double *pmt, double *n, double r) r0 = sqrt(R*R - L*L); a = 4*r0*r/pow(r0+r,2); - Rmax = L*L + (r0+r)*(r0+r); + Rmax = sqrt(L*L + (r0+r)*(r0+r)); k = sqrt(4*r0*r/(L*L + pow(r0+r,2))); if (r0 <= r) { -- cgit