#!/bin/csh
#
# This script generates the AR data sets.  Note that the input MR files
# must be sorted by octant as a preprocessing measure.  See the 
# script FromASCIItoBinarySortCubeScript for more details
#
# Robert S Laramee
# Mon 19 Mar '01
# FromMRtoARscript
#
# Recall the AR data sets do *not* apply to wavelet data or the "combo"
set BASE = "../data.selective"
#
# Remove output from previous run of this script.
#
echo 'compiling rlaramee/FromMRtoARcube.java...'
javac -deprecation rlaramee/FromMRtoARcube.java
#
echo 'removing previous AR files...'
@ level = 6
@ resolution = 2
while ($resolution <= 128)

    rm -vf $BASE/level${level}/headRes${resolution}AR*.bin
    @ level = $level - 1
    @ resolution = $resolution * 2

end
#
echo 'generating AR files...'
echo 'AR THRESHOLD: 1%...'
@ level = 6
@ resolution = 2
while ($resolution <= 128)

	java -Xmx1024m rlaramee/FromMRtoARcube \
	$BASE/level${level}/headRes${resolution}cube.bin \
	$resolution $resolution $resolution 0.01 \
	$BASE/level${level}/headRes${resolution}AR01.bin

	@ level = $level - 1
	@ resolution = $resolution * 2

end
echo 'AR THRESHOLD: 2%...'
@ level = 6
@ resolution = 2
while ($resolution <= 128)

	java -Xmx1024m rlaramee/FromMRtoARcube \
	$BASE/level${level}/headRes${resolution}cube.bin \
	$resolution $resolution $resolution 0.02 \
	$BASE/level${level}/headRes${resolution}AR02.bin

	@ level = $level - 1
	@ resolution = $resolution * 2

end
echo 'AR THRESHOLD: 5%...'
@ level = 6
@ resolution = 2
while ($resolution <= 128)

	java -Xmx1024m rlaramee/FromMRtoARcube \
	$BASE/level${level}/headRes${resolution}cube.bin \
	$resolution $resolution $resolution 0.05 \
	$BASE/level${level}/headRes${resolution}AR05.bin

	@ level = $level - 1
	@ resolution = $resolution * 2

end
echo 'AR THRESHOLD: 10%...'
@ level = 6
@ resolution = 2
while ($resolution <= 128)

	java -Xmx1024m rlaramee/FromMRtoARcube \
	$BASE/level${level}/headRes${resolution}cube.bin \
	$resolution $resolution $resolution 0.10 \
	$BASE/level${level}/headRes${resolution}AR10.bin

	@ level = $level - 1
	@ resolution = $resolution * 2

end
echo 'AR THRESHOLD: 15%...'
@ level = 6
@ resolution = 2
while ($resolution <= 128)

	java -Xmx1024m rlaramee/FromMRtoARcube \
	$BASE/level${level}/headRes${resolution}cube.bin \
	$resolution $resolution $resolution 0.15 \
	$BASE/level${level}/headRes${resolution}AR15.bin

	@ level = $level - 1
	@ resolution = $resolution * 2

end
echo 'AR THRESHOLD: 20%...'
@ level = 6
@ resolution = 2
while ($resolution <= 128)

	java -Xmx1024m rlaramee/FromMRtoARcube \
	$BASE/level${level}/headRes${resolution}cube.bin \
	$resolution $resolution $resolution 0.20 \
	$BASE/level${level}/headRes${resolution}AR20.bin

	@ level = $level - 1
	@ resolution = $resolution * 2

end
echo 'AR THRESHOLD: 25%...'
@ level = 6
@ resolution = 2
while ($resolution <= 128)

	java -Xmx1024m rlaramee/FromMRtoARcube \
	$BASE/level${level}/headRes${resolution}cube.bin \
	$resolution $resolution $resolution 0.25 \
	$BASE/level${level}/headRes${resolution}AR25.bin

	@ level = $level - 1
	@ resolution = $resolution * 2

end
#
echo 'done generating AR files 2^3 - 128^3 cubes for 1-25% thresholds'
echo 'FromMRtoARscript finished'