diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2018-07-04 12:22:21 -0400 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2018-07-04 12:22:21 -0400 |
| commit | 9646e43e6016536911f0703dc7e3af7a0eae56db (patch) | |
| tree | 99c744bb1e589a68b66b8fd54860b0c189aef677 /solid_angle.c | |
| parent | 6b8e727b10ff7a1ae219987fd4ff40f359a0c655 (diff) | |
| download | sddm-9646e43e6016536911f0703dc7e3af7a0eae56db.tar.gz sddm-9646e43e6016536911f0703dc7e3af7a0eae56db.tar.bz2 sddm-9646e43e6016536911f0703dc7e3af7a0eae56db.zip | |
add a test of the get_solid_angle function
Diffstat (limited to 'solid_angle.c')
| -rw-r--r-- | solid_angle.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/solid_angle.c b/solid_angle.c index 5fc6033..cbc71f1 100644 --- a/solid_angle.c +++ b/solid_angle.c @@ -15,7 +15,7 @@ double get_solid_angle(double *pos, double *pmt, double *n, double r) dir[1] = pos[1] - pmt[1]; dir[2] = pos[2] - pmt[2]; - L = dir[0]*n[0] + dir[1]*n[1] + dir[2]*n[2]; + L = fabs(dir[0]*n[0] + dir[1]*n[1] + dir[2]*n[2]); R = sqrt(dir[0]*dir[0] + dir[1]*dir[1] + dir[2]*dir[2]); r0 = sqrt(R*R - L*L); |