#!/usr/bin/env bash

# @(#) s1	Demonstrate arithmetic combination of column data, awk, dm.
# For dm utility in |stat, see:
# http://hcibib.org/perlman/stat/

# Utility functions: print-as-echo, print-line-with-visual-space, debug.
# export PATH="/usr/local/bin:/usr/bin:/bin"
LC_ALL=C ; LANG=C ; export LC_ALL LANG
pe() { for _i;do printf "%s" "$_i";done; printf "\n"; }
pl() { pe;pe "-----" ;pe "$*"; }
db() { ( printf " db, ";for _i;do printf "%s" "$_i";done;printf "\n" ) >&2 ; }
db() { : ; }
C=$HOME/bin/context && [ -f $C ] && $C paste awk dm

pl " Input data files data1 data2, combined into data3:"
paste data1 data2 > data3
head data1 data2 data3

pl " Results, awk:"
awk '{print $1,$2*$3}' data3

pl " Results, |stat dm:"
dm s1 'x2*x3' < data3

exit 0

$ ./s1

Environment: LC_ALL = C, LANG = C
(Versions displayed with local utility "version")
OS, ker|rel, machine: Linux, 2.6.26-2-amd64, x86_64
Distribution        : Debian 5.0.8 (lenny, workstation) 
bash GNU bash 3.2.39
paste (GNU coreutils) 6.10
awk GNU Awk 3.1.5
dm - ( local: ~/executable/dm, 2009-11-09 )

-----
 Input data files data1 data2, combined into data3:
==> data1 <==
AAAAA	3.465830E-12
BBBBB	4.263280E-08
CCCCC	1.113320E-17
DDDDD	0.000000E+00

==> data2 <==
3.932350E-12 
1.194380E-07 
4.901480E-17 
0.000000E+00 

==> data3 <==
AAAAA	3.465830E-12	3.932350E-12 
BBBBB	4.263280E-08	1.194380E-07 
CCCCC	1.113320E-17	4.901480E-17 
DDDDD	0.000000E+00	0.000000E+00 

-----
 Results, awk:
AAAAA 1.36289e-23
BBBBB 5.09198e-15
CCCCC 5.45692e-34
DDDDD 0

-----
 Results, |stat dm:
AAAAA	1.36289e-23
BBBBB	5.09198e-15
CCCCC	5.45692e-34
DDDDD	0