#!python
#
#    JKS - Measurement database system
#    Copyright (C) 2013-2024  Christoph Lehner (christoph.lehner@ur.de, https://github.com/lehner/jks)
#
#    This program is free software; you can redistribute it and/or modify
#    it under the terms of the GNU General Public License as published by
#    the Free Software Foundation; either version 2 of the License, or
#    (at your option) any later version.
#
#    This program is distributed in the hope that it will be useful,
#    but WITHOUT ANY WARRANTY; without even the implied warranty of
#    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
#    GNU General Public License for more details.
#
#    You should have received a copy of the GNU General Public License along
#    with this program; if not, write to the Free Software Foundation, Inc.,
#    51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
#
import jks, sys, os, math, glob, numpy as np
if len(sys.argv) != 5:
    print("%s out.jks A.jks B.jks op" % sys.argv[0])
    print("")
    print("Info:")
    print("- each joint tag of A and B is processed")
    print("- numpy is defined as np")
    print("- in op: r['A'] is the element from A, r['B'] is the element from B")
    sys.exit(0)

fno=sys.argv[1]
fnA=sys.argv[2]
fnB=sys.argv[3]
oop=sys.argv[4]

res=jks.resamples(fno)
inA=jks.resamples(fnA)
inB=jks.resamples(fnB)

jks2=jks.resamples()

for tag in inA.keys():
    if tag in inB.keys():
        jkA=inA.get(tag)
        jkB=inB.get(tag)

        jks2.add("A",jkA)
        jks2.add("B",jkB)
        op=eval("lambda r: %s" % oop)
        jk=jks2.apply(op)
        res.add(tag,jk)

res.save(fno)

