ariths-gen/chr2c.py

#!/usr/bin/python3
# A simple script to convert the input CGP chromosome to the corresponding output in C code
# Script originally written by Vojtech Mrazek
# Modifications to suite the appropriate needs of ArithsGen made by Jan Klhufek
import sys
import re
import numpy as np
import argparse


# Parse all nodes present in input CGP
def parse_node(n):
    return list(map(int, re.match(r"\[(\d+)\](\d+),(\d+),(\d+)", n).groups()))


# Recursively detect wires required to get the result wire with id 'id' and activate them (will be generated)
def activate(id, act, c_in, cdata_dict):
    if act[id]:
        return 0

    act[id] = True

    conf = cdata_dict[id-(c_in+2)]
    if conf:
        _, a, b, f = conf
        activate(a, act, c_in, cdata_dict)
        if f != 0 and f != 1:  # Identity, NOT
            activate(b, act, c_in, cdata_dict)


# Declare new output logic gate's wire if it has not been declared yet (and then assign it the result of corresponding logical function)
# Or return previously declared wire (which is used as an input into corresponding logical function)
def get_sig(name, trans, create=False):
    if create is True and int(name) not in list(trans.keys()):
        trans[name] = "sig_%d" % name
    return trans[name]


def parse_chromosome(chromosome, signed=False, function=None):
    
    # Load input CGP
    f = open(chromosome, "r")
    chrom = "\n".join([x for x in f])
    f.close()

    output = ["// " + chromosome]

    # Splitting input representation into three parts
    preamble, data, outputs = re.match(r"\{(.*)\}\((.*)\)\(([^\)]+)\)", chrom).groups()
    # Parsing CGP preamble parameters
    c_in, c_out, c_rows, c_cols, c_ni, c_no, c_lback = list(map(int, preamble.split(",")))
    # Parsing all present nodes
    cdata = list(map(parse_node, data.split(")(")))

    # Parse position and corresponding information regarding each chromosome
    cdata_dict = [None for i in range(0, c_cols * c_rows)]

    for id, chromosome in enumerate(cdata):
        cid, a, b, f = chromosome
        cdata_dict[id] = (id + 2 + c_in, a, b, f)

    # Reserve position for all wires present in the genotype
    # Later used to detect if they need to be activated for obtaining the result
    act = [False for i in range(0, (2 + c_in + c_cols * c_rows))]
    # Reserved constant wires indexes
    act[0] = True
    act[1] = True

    # Obtaining all output wires (MSB to LSB)
    outs = list(map(int, outputs.split(",")))
    # Recursively activate all (node) wires needed to obtain the corresponding result wire
    for o in outs:
        activate(o, act, c_in, cdata_dict)

    # Export converted input declarations into C code function header / body start
    output.append("#include <stdint.h>")

    if not function:
        function = f"circuit{c_in/2:.0f}"
    
    if signed:
        dtype = "int64_t"
    else:
        dtype = "uint64_t"

    output.append(f"{dtype} {function}({dtype} a, {dtype} b) {{")
    output.append("  int wa[%d];" % int(c_in/2))
    output.append("  int wb[%d];" % int(c_in/2))
    output.append(f"  {dtype} y = 0;")

    # Export converted input assignments into C code function body
    trans = {}
    trans[0] = "0x00"
    trans[1] = "0x01"
    for i in range(0, int(c_in/2)):
        trans[i+2] = "wa[%d]" % i
        output.append("  wa[%d] = (a >> %d) & 0x01;" % (i, i))
        trans[i + (c_in/2)+2] = "wb[%d]" % i
        output.append("  wb[%d] = (b >> %d) & 0x01;" % (i, i))

    # Lists representing individual nodes (lines), output wires (lines_end)
    lines = []
    lines_end = []

    # Export converted nodes into C code function body
    for dat, a in zip(cdata_dict, act):
        id, a, b, f = dat

        if f == 0:
            lines.append("  uint8_t " + get_sig(id, trans, True) + " = " + get_sig(a, trans) + ";")
        elif f == 1:
            lines.append("  uint8_t " + get_sig(id, trans, True) + " = !" + get_sig(a, trans) + ";")
        elif f == 2:
            lines.append("  uint8_t " + get_sig(id, trans, True) + " = " + get_sig(a, trans) + " & " + get_sig(b, trans) + ";")
        elif f == 3:
            lines.append("  uint8_t " + get_sig(id, trans, True) + " = " + get_sig(a, trans) + " | " + get_sig(b, trans) + ";")
        elif f == 4:
            lines.append("  uint8_t " + get_sig(id, trans, True) + " = " + get_sig(a, trans) + " ^ " + get_sig(b, trans) + ";")
        elif f == 5:
            lines.append("  uint8_t " + get_sig(id, trans, True) + " = ~(" + get_sig(a, trans) + " & " + get_sig(b, trans) + ") & 0x01;")
        elif f == 6:
            lines.append("  uint8_t " + get_sig(id, trans, True) + " = ~(" + get_sig(a, trans) + " | " + get_sig(b, trans) + ") & 0x01;")
        elif f == 7:
            lines.append("  uint8_t " + get_sig(id, trans, True) + " = ~(" + get_sig(a, trans) + " ^ " + get_sig(b, trans) + ") & 0x01;")

    # Export converted outputs into C code function body
    for i in range(0, c_out):
        if outs[i] in trans:
            lines_end.append("  y |=  (%s & 0x01ull) << %d; // default output" % (trans[outs[i]], i))
        else:
            assert False

    # Ensure proper sign extension if it is needed
    if signed is True:
        last_out_wire = trans[outs[-1]] # used for sign extension
        for i in range(c_out, 64):
            lines_end.append("  y |=  (%s & 0x01ull) << %d; // default output" % (last_out_wire, i))

    # Print final result return in C code function body and close it
    lines_end.append("  return y;")
    output += lines
    output += lines_end;
    output.append("}")

    return "\n".join(output)

# Parse and convert the input CGP chromosome to output C code representation
if __name__ == "__main__":
    argparse.ArgumentParser()
    parser = argparse.ArgumentParser()

    parser.add_argument("chromosome", help="Path to CGP file")
    parser.add_argument("--signed", action="store_true")
    parser.add_argument("--function", default=None)
    args = parser.parse_args()
    print(parse_chromosome(**vars(args)))