#!/bin/bash 

# 14-08-19 20:02:41 
# this script launches a batch of gencase runs that draw floating objects
# and reads out the moments of inertia out of the inertia tensor
# (diagonal elements); it does not inspect inertia products off the diagonal
# Usage note: before starting select the paths, ids and data sets in the section user-defined parts


# 1 USER-DEFINED PARTS

## functions
  
fail () { 
	echo Execution aborted. 
	#    read -n1 -r -p "Press any key to continue..." key 
	exit 1 
}

## functions to extract a single diagonal element of inertia tensor from xml file

Ixx() {
	cat ${dirout2}/${name}.xml | grep "<floating mkbound=\"${mkbound}\"" -A7 | grep v11 | cut -d\" -f2 
}

Iyy() {
	cat ${dirout2}/${name}.xml | grep "<floating mkbound=\"${mkbound}\"" -A7 | grep v22 | cut -d\" -f4  
}

Izz() {
	cat ${dirout2}/${name}.xml | grep "<floating mkbound=\"${mkbound}\"" -A7 | grep v33 | cut -d\" -f6 
}

## data sets

### a sample of 30 random dp's 
### they lead to a given of L/Dp's uniformly distributed in the range 8-96

sample="0.0667467342 \
0.0424280136 \
0.0375202416 \
0.0326268291 \
0.0321952850 \
0.0302021326 \
0.0293327468 \
0.0289502166 \
0.0274384254 \
0.0267748444 \
0.0234328231 \
0.0227003922 \
0.0215810239 \
0.0167618759 \
0.0160500455 \
0.0160066984 \
0.0158548490 \
0.0155654362 \
0.0153154617 \
0.0152776139 \
0.0151462866 \
0.0147402413 \
0.0140387424 \
0.0139238217 \
0.0136520394 \
0.0130019603 \
0.0119289974 \
0.0117426064 \
0.0110519980 \
0.0106245671	
"

### dp's that lead to integer L/dp between 8 and 96

integers=" \
0.1250000000 \
0.1111111111 \
0.0833333333 \
0.0625000000 \
0.0555555556 \
0.0416666665 \
0.0312500000 \
0.0277777778 \
0.0208333333 \
0.0178571429 \
0.0156250000 \
0.0138888889 \
0.0104166667 \
"

### select the data set you want to use ()

given=${integers}
given=${sample}
given=0.1200000000  # just one value to check it runs

## select simulation paths, object ids, and parameters

export name=InertiaBug01converge
export dirbin=~/Solvers/sph/dsph/gitlab44/bin/linux/
export mkbound1=10
export mkbound2=30
export mkbound3=31

### select the gencase version if ready for use

export gencase="${dirbin}/GenCase4_linux64-40027" # 11 may 2016
export gencase="${dirbin}/GenCase4_linux64-40109" # 21 june 2019
export gencase="${dirbin}/GenCase4_MkWord_linux64" # 4.0.092 23 march 2019
export gencase="${dirbin}/GenCase4_linux64" # 4.0.092 23 march 2019

# 2 DONE FOR YOU BY THE COMPUTER IN SOLITUDE

export LD_LIBRARY_PATH=${LD_LIBRARY_PATH}:${dirbin}

## remove and remake directory

export dirout=${name}
export dirout2=${dirout}/gencase
if [ -e ${dirout2} ]; then rm -r ${dirout2}; fi
mkdir -p ${dirout2}

## clear old output files

rm $0*.out

## loops through the list of givens

for dx in $given; do

	echo
	echo particle spacing $dx

	${gencase} -dp:$dx ${name}_Def ${dirout2}/${name} #-save:vtkbound #-createdirs:1
	if [ $? -ne 0 ] ; then fail; fi
	echo

	### read out information on floating object (extended)

	cat ${dirout2}/${name}.xml | grep "<inertia" -A3 -B4 | sed s/' '*\</\</g 
	echo
	
	### read out diagonal moments of inertia only

	#### first object
	mkbound=${mkbound1}
	a1=$(Ixx)
	b1=$(Iyy)
	c1=$(Izz)
#	Ixx && Iyy && Izz
	echo $dx $a1 $b1 $c1 2>&1 | tee -a $0-${mkbound}.out
	echo

	#### second object
	mkbound=${mkbound2}
	a2=$(Ixx)
	b2=$(Iyy)
	c2=$(Izz)
#	Ixx && Iyy && Izz
	echo $dx $a2 $b2 $c2 2>&1 | tee -a $0-${mkbound}.out
	echo

	#### third object
	mkbound=${mkbound3}
	a3=$(Ixx)
	b3=$(Iyy)
	c3=$(Izz)
#	Ixx && Iyy && Izz
	echo $dx $a3 $b3 $c3 2>&1 | tee -a $0-${mkbound}.out
	echo

	sleep 0.5

done

## shows GenCase version and close

echo
echo version: $( ${gencase} -h | grep GenCase | grep v )

echo "End of script $0"