nnc stuff.py

# split

for i in range(107):
  a[i] = b[i] + c[i]

for i0 in range(10):
  for i1 in range(10):
    a[i0 * 10 + i1] = b[i0 * 10 + i1] + c[i0 * 10 + i1]
for i0 in range(7):
  a[100 + i0] = b[100 + i0] + c[100 + i0]
  
# rfactor

for i in range(16):
  a[0] += b[i]

a1 = [0 for _ in range(4)]
for i0 in range(4):
  for i1 in range(4):
    a1[i0] += b[i0 * 4 + i1]
for i0 in range(4):
  a[0] += a1[i0]
  
# compute at

for i in range(16):
  for j in range(16):
    a[i, j] = b[i, j]
for i in range(16):
  for j in range(16):
    for k in range(16):
      c[k] += a[i, j]

for i in range(16):
  for j in range(16):
    a[i, j] = b[i, j]
    for k in range(16):
      c[k] += a[i, j]