Changeset 26183

Timestamp:

Nov 18, 2009, 11:01:03 AM (16 years ago)

Author:

bills

Message:

prevent division by zero

File:

• trunk/ippScripts/scripts/magic_tree.pl (modified) (1 diff)

trunk/ippScripts/scripts/magic_tree.pl

@@ -169 +169 @@
     # Coordinates on rotated celestial sphere
     my $phi = atan2($eta,$xi) + pi/2;
-    my $theta = atan(180 / pi / sqrt($xi**2 + $eta**2));
+    my $denominator = sqrt($xi**2 + $eta**2);
+    &my_die("denominator is zero!!", $magic_id, $PS_EXIT_PROG_ERROR) if $denominator == 0;
+
+    my $theta = atan(180 / pi / $denominator);
 
     # Coordinates on celestial sphere