more reorg

Author: ksargsyan

apps/awkies/getrange.x

@@ -0,0 +1,48 @@
+#!/bin/bash
+#=====================================================================================
+
+# Given a Nxd matrix of samples, compute dx2 ranges
+# Example: getrange.x samples.dat > ranges.dat
+
+shopt -s expand_aliases
+alias awk="awk -v OFMT='%.15e'"
+
+if [ $# -lt 1 ]; then
+    echo "Number of arguments can not be less than 1"
+    echo "Syntax: $0 <filename> [<cushion_fraction>]"
+    exit
+elif [ $# -gt 2 ]; then
+    echo "Number of arguments can not be greater than 2"
+    echo "Syntax: $0 <filename> [<cushion_fraction>]"
+    exit
+elif [ $# -eq 1 ]; then
+    fr=0.0
+elif [ $# -eq 2 ]; then
+    fr=$2
+fi
+
+filename=$1
+
+DIM=`awk 'END{print NF}' $filename`
+
+for (( COL=1; COL<=$DIM ; COL++ )); do
+
+awk 'BEGIN {
+valmin=1e+100;
+valmax=-1e+100;
+line=1;
+}
+{
+if( $(col) < valmin )
+{
+  valmin=$(col);
+}
+if( $(col) > valmax )
+{
+  valmax=$(col);
+}
+}
+END{
+print valmin-fr*(valmax-valmin), valmax+fr*(valmax-valmin)
+}' col=$COL $filename
+done

apps/awkies/scale.x

@@ -0,0 +1,35 @@
+#!/bin/bash -e
+# Scales matrix data to or from a given parameter domain to [-1,1]^d
+# scale.x <input> <to or from> <domain> <output>
+
+shopt -s expand_aliases
+alias awk="awk -v OFMT='%.15e'"
+
+IN_FILE=$1
+TO_FROM=$2
+DOM_FILE=$3
+OUT_FILE=$4
+
+N=`awk 'END{print NR}' $IN_FILE`
+D=`awk 'END{print NR}' $DOM_FILE`
+
+DD=`awk 'END{print NF}' $IN_FILE`
+
+## check that DD=D
+
+echo "" > $OUT_FILE
+for (( i=1; i<=$D ; i++ )); do
+A=`awk 'NR==i{print $1}' i=$i $DOM_FILE`
+B=`awk 'NR==i{print $2}' i=$i $DOM_FILE`
+
+if [ "$TO_FROM" = "from" ]; then
+    awk '{print -1.+2.*($i-a)/(b-a)}' i=$i a=$A b=$B $IN_FILE > tmp
+elif [ "$TO_FROM" = "to" ]; then
+    awk '{print (a+b)/2.+$i*(b-a)/2.}' i=$i a=$A b=$B $IN_FILE > tmp
+else
+    echo "Second argument has to be 'to' or 'from'"
+    exit
+fi
+paste $OUT_FILE tmp > tmpp; mv tmpp $OUT_FILE
+
+done

apps/awkies/transpose.x

@@ -0,0 +1,26 @@
+#!/bin/bash -e
+# Transpose a matrix: assumes all lines have same number
+# of fields
+
+# Example: transpose_file.x file_in > file_out
+
+shopt -s expand_aliases
+alias awk="awk -v OFMT='%.15e'"
+
+exec awk '
+NR == 1 {
+	n = NF
+	for (i = 1; i <= NF; i++)
+		row[i] = $i
+	next
+}
+{
+	if (NF > n)
+		n = NF
+	for (i = 1; i <= NF; i++)
+		row[i] = row[i] " " $i
+}
+END {
+	for (i = 1; i <= n; i++)
+		print row[i]
+}' ${1+"$@"}

apps/pc_fit.py

@@ -113,7 +113,7 @@
     for isam in range(0, nsam, nevery):
         f = plt.figure(figsize=(12,4))
         plt.plot(xc, y[isam,:], 'bo-', ms=8, label='Model')
-        plt.plot(xc, ypred[isam,:], 'go-', ms=8, label='PC apprx.')
+        plt.plot(xc, ypred[isam,:], 'go-', ms=8, label='PC Apprx.')
         plt.title(f'Sample #{isam+1}')
         plt.xlabel(xlabel)
         plt.ylabel(ylabel)

apps/plot_ens.py

@@ -17,6 +17,8 @@
 args = parser.parse_args()
 
 ydata = np.loadtxt(args.ydata)
+if len(ydata.shape)==1:
+    ydata = ydata[:, np.newaxis]
 
 
 nsam, nout = ydata.shape

apps/uqpc/model.x

@@ -0,0 +1,26 @@
+#!/bin/bash
+
+INPUT_FILE=$1 
+OUTPUT_FILE=$2
+
+# Run the black-box model,
+# (a) Test function with 3 inputs and 5 outputs
+awk '{print exp($1), $3*log($1**2)+$1+($2**2)*(1-exp(-$3**2)), $1+$3**2, $2+$3, $1*$3}' $1 > $2
+
+# (b) Entry-wise exponential, same number of inputs and outputs
+#awk '{for (i=1; i<=NF; i++) { printf("%lg ",exp($i)) }; printf("\n")}' $1 > $2
+
+
+# # (c) A three parameter time-dependent model, i.e. as many outputs as there are time grid points
+# # Warning: this will generate a lot of figures after the whole workflow is done.
+# NGRID=50
+# echo -n"" > $2
+# for ((tgrid=1;tgrid<=$NGRID;tgrid++)); do
+#     awk '{print $3*$1**2+$2*exp((tgr/ngr)*$3)}' tgr=$tgrid ngr=$NGRID $1 > tmp
+#     paste $2 tmp > tmpp; mv tmpp $2; rm -f tmp
+# done
+# # Output names are time grid points
+# rm -rf outnames.txt
+# for ((i=1;i<=$NGRID;i++)); do
+#     echo "$i" >> outnames.txt
+# done