--- ray/src/util/genambpos.pl	2014/04/24 23:15:42	2.1
+++ ray/src/util/genambpos.pl	2014/05/01 04:51:12	2.5
@@ -1,17 +1,19 @@
 #!/usr/bin/perl -w
-# RCSid $Id: genambpos.pl,v 2.1 2014/04/24 23:15:42 greg Exp $
+# RCSid $Id: genambpos.pl,v 2.5 2014/05/01 04:51:12 greg Exp $
 #
 # Visualize ambient positions and gradients
 #
 use strict;
 sub userror {
-	print STDERR "Usage: genambpos [-l lvl][-s scale][-p][-d] scene.amb > ambloc.rad\n";
+	print STDERR "Usage: genambpos [-l lvl][-w minwt][-r rad][-s sf][-p][-d] scene.amb > ambloc.rad\n";
 	exit 1;
 }
 my $lvlsel = -1;
 my $scale = 0.25;
 my $doposgrad = 0;
 my $dodirgrad = 0;
+my $minwt = 0.5001**6;
+my $fixedrad="";
 my $savedARGV = "genambpos @ARGV";
 # Get options
 while ($#ARGV >= 0) {
@@ -22,9 +24,15 @@ while ($#ARGV >= 0) {
 	} elsif ("$ARGV[0]" =~ /^-l/) {
 		$lvlsel = $ARGV[1];
 		shift @ARGV;
+	} elsif ("$ARGV[0]" =~ /^-w/) {
+		$minwt = $ARGV[1];
+		shift @ARGV;
 	} elsif ("$ARGV[0]" =~ /^-s/) {
 		$scale = $ARGV[1];
 		shift @ARGV;
+	} elsif ("$ARGV[0]" =~ /^-r/) {
+		$fixedrad = "-e psiz:$ARGV[1]";
+		shift @ARGV;
 	} elsif ("$ARGV[0]" =~ /^-./) {
 		userror();
 	} else {
@@ -34,11 +42,11 @@ while ($#ARGV >= 0) {
 }
 userror() if ($#ARGV != 0);
 my $cmd = "getinfo < $ARGV[0] " . 
-		q[| sed -n 's/^.* -aa \([.0-9][.0-9]*\) .*$/\1/p'];
+		q[| sed -n 's/^.* -aa \([.0-9][^ ]*\) .*$/\1/p'];
 my $ambacc=`$cmd`;
 die "Missing -aa setting in header\n" if (! $ambacc );
-$scale *= $ambacc**.25;
-my $outfmt = '
+$scale *= $ambacc;
+my $ambfmt = '
 void glow posglow
 0
 0
@@ -48,16 +56,21 @@ posglow sphere position${recno}
 0
 0
 4 ${px} ${py} ${pz} ${psiz}
+';
+my $posgradfmt = '
+void glow arrglow
+0
+0
+4 ${wt*agr} ${wt*agg} ${wt*agb} 0
 
-posglow cone pgarrow${recno}
+arrglow cone pgarrow${recno}
 0
 0
 8
 	${ cx0 }	${ cy0 }	${ cz0 }
 	${ cx1 }	${ cy1 }	${ cz1 }
 	${ cr0 }	0
-';
-my $posgradfmt = '
+
 void brightfunc pgpat
 2 posfunc ambpos.cal
 0
@@ -82,8 +95,7 @@ pgeval polygon pgellipse${recno}
 	${ px3 } ${ py3 } ${ pz3 }
 	${ px4 } ${ py4 } ${ pz4 }
 ';
-$outfmt .= $posgradfmt if ($doposgrad);
-my $dirgradfmt='
+$posgradfmt .= '
 void glow tipglow
 0
 0
@@ -93,7 +105,8 @@ tipglow sphere atip
 0
 0
 4 ${ cx1 } ${ cy1 } ${ cz1 } ${psiz/7}
-
+' if ($dodirgrad);
+my $dirgradfmt='
 void brightfunc dgpat
 2 dirfunc ambpos.cal
 0
@@ -108,7 +121,7 @@ dgval ring dgdisk${recno}a
 0
 0
 8
-	${ px+dgx*.0001 } ${ py+dgy*.0001 } ${ pz+dgz*.0001 }
+	${ px+dgx*.001 } ${ py+dgy*.001 } ${ pz+dgz*.001 }
 	${ dgx } ${ dgy } ${ dgz }
 	0	${ r0/2 }
 
@@ -116,12 +129,22 @@ dgval ring dgdisk${recno}b
 0
 0
 8
-	${ px-dgx*.0001 } ${ py-dgy*.0001 } ${ pz-dgz*.0001 }
+	${ px-dgx*.001 } ${ py-dgy*.001 } ${ pz-dgz*.001 }
 	${ -dgx } ${ -dgy } ${ -dgz }
 	0	${ r0/2 }
 ';
-$outfmt .= $dirgradfmt if ($dodirgrad);
 # Load & convert ambient values
 print "# Output produced by: $savedARGV\n";
-system "lookamb -h -d $ARGV[0] | rcalc -e 'LV:$lvlsel;SF:$scale' -f rambpos.cal -o '$outfmt'\n";
+system "lookamb -h -d $ARGV[0] | rcalc -e 'LV:$lvlsel;MW:$minwt;SF:$scale'" .
+		" -f rambpos.cal -e cond=acond $fixedrad -o '$ambfmt'";
+if ($doposgrad) {
+	system "lookamb -h -d $ARGV[0] " .
+		"| rcalc -e 'LV:$lvlsel;MW:$minwt;SF:$scale'" .
+		" -f rambpos.cal -e cond=pcond $fixedrad -o '$posgradfmt'";
+}
+if ($dodirgrad) {
+	system "lookamb -h -d $ARGV[0] " .
+		"| rcalc -e 'LV:$lvlsel;MW:$minwt;SF:$scale'" .
+		" -f rambpos.cal -e cond=dcond -o '$dirgradfmt'";
+}
 exit;