Made corrections to the check_density.py script so that partial
areas on the right and top sides of the design that don't fill out
the full 70um x 70um of a check tile are prorated when computing
the total fill density for the area.
diff --git a/scripts/check_density.py b/scripts/check_density.py
index 76f62a3..129a848 100755
--- a/scripts/check_density.py
+++ b/scripts/check_density.py
@@ -153,6 +153,12 @@
print('puts stdout "YTILES: $ytiles"', file=ofile)
print('', file=ofile)
+ # Need to know what fraction of a full tile is the last row and column
+ print('set xfrac [expr {($xtiles * $stepsizex - $fullwidth + 0.0) / $stepsizex}]', file=ofile)
+ print('set yfrac [expr {($ytiles * $stepsizey - $fullheight + 0.0) / $stepsizey}]', file=ofile)
+ print('puts stdout "XFRAC: $xfrac"', file=ofile)
+ print('puts stdout "YFRAC: $yfrac"', file=ofile)
+
print('cif ostyle density', file=ofile)
# Process density at steps. For efficiency, this is done in 70x70 um
@@ -229,6 +235,10 @@
met3fill = []
met4fill = []
met5fill = []
+ xtiles = 0
+ ytiles = 0
+ xfrac = 0.0
+ yfrac = 0.0
while mproc:
status = mproc.poll()
@@ -270,6 +280,10 @@
xtiles = int(dpair[1].strip())
elif layer == 'YTILES':
ytiles = int(dpair[1].strip())
+ elif layer == 'XFRAC':
+ xfrac = float(dpair[1].strip())
+ elif layer == 'YFRAC':
+ yfrac = float(dpair[1].strip())
for line in errlines.splitlines():
print(line)
@@ -302,16 +316,32 @@
print('')
print('Density results (total tiles = ' + str(total_tiles) + '):')
+ # Full areas are 10 x 10 tiles = 100. But the right and top sides are
+ # not full tiles, so the full area must be prorated.
+
+ sideadjust = 90.0 + (10.0 * xfrac)
+ topadjust = 90.0 + (10.0 * yfrac)
+ corneradjust = 81.0 + (9.0 * xfrac) + (9.0 * yfrac) + (xfrac * yfrac)
+
print('')
print('FOM Density:')
for y in range(0, ytiles - 9):
+ if y == ytiles - 10:
+ atotal = topadjust
+ else:
+ atotal = 100.0
for x in range(0, xtiles - 9):
+ if x == xtiles - 10:
+ if y == ytiles - 10:
+ atotal = corneradjust
+ else:
+ atotal = sideadjust
fomaccum = 0
for w in range(y, y + 10):
base = xtiles * w + x
fomaccum += sum(fomfill[base : base + 10])
- print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(fomaccum / 100.0))
+ print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(fomaccum / atotal))
if fomaccum < 33.0:
print('***Error: FOM Density < 33%')
elif fomaccum > 57.0:
@@ -320,24 +350,42 @@
print('')
print('POLY Density:')
for y in range(0, ytiles - 9):
+ if y == ytiles - 10:
+ atotal = topadjust
+ else:
+ atotal = 100.0
for x in range(0, xtiles - 9):
+ if x == xtiles - 10:
+ if y == ytiles - 10:
+ atotal = corneradjust
+ else:
+ atotal = sideadjust
polyaccum = 0
for w in range(y, y + 10):
base = xtiles * w + x
polyaccum += sum(polyfill[base : base + 10])
- print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(polyaccum / 100.0))
+ print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(polyaccum / atotal))
print('')
print('LI Density:')
for y in range(0, ytiles - 9):
+ if y == ytiles - 10:
+ atotal = topadjust
+ else:
+ atotal = 100.0
for x in range(0, xtiles - 9):
+ if x == xtiles - 10:
+ if y == ytiles - 10:
+ atotal = corneradjust
+ else:
+ atotal = sideadjust
liaccum = 0
for w in range(y, y + 10):
base = xtiles * w + x
liaccum += sum(lifill[base : base + 10])
- print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(liaccum / 100.0))
+ print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(liaccum / atotal))
if liaccum < 35.0:
print('***Error: LI Density < 35%')
elif liaccum > 70.0:
@@ -346,13 +394,22 @@
print('')
print('MET1 Density:')
for y in range(0, ytiles - 9):
+ if y == ytiles - 10:
+ atotal = topadjust
+ else:
+ atotal = 100.0
for x in range(0, xtiles - 9):
+ if x == xtiles - 10:
+ if y == ytiles - 10:
+ atotal = corneradjust
+ else:
+ atotal = sideadjust
met1accum = 0
for w in range(y, y + 10):
base = xtiles * w + x
met1accum += sum(met1fill[base : base + 10])
- print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met1accum / 100.0))
+ print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met1accum / atotal))
if met1accum < 35.0:
print('***Error: MET1 Density < 35%')
elif met1accum > 70.0:
@@ -361,13 +418,22 @@
print('')
print('MET2 Density:')
for y in range(0, ytiles - 9):
+ if y == ytiles - 10:
+ atotal = topadjust
+ else:
+ atotal = 100.0
for x in range(0, xtiles - 9):
+ if x == xtiles - 10:
+ if y == ytiles - 10:
+ atotal = corneradjust
+ else:
+ atotal = sideadjust
met2accum = 0
for w in range(y, y + 10):
base = xtiles * w + x
met2accum += sum(met2fill[base : base + 10])
- print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met2accum / 100.0))
+ print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met2accum / atotal))
if met2accum < 35.0:
print('***Error: MET2 Density < 35%')
elif met2accum > 70.0:
@@ -376,13 +442,22 @@
print('')
print('MET3 Density:')
for y in range(0, ytiles - 9):
+ if y == ytiles - 10:
+ atotal = topadjust
+ else:
+ atotal = 100.0
for x in range(0, xtiles - 9):
+ if x == xtiles - 10:
+ if y == ytiles - 10:
+ atotal = corneradjust
+ else:
+ atotal = sideadjust
met3accum = 0
for w in range(y, y + 10):
base = xtiles * w + x
met3accum += sum(met3fill[base : base + 10])
- print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met3accum / 100.0))
+ print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met3accum / atotal))
if met3accum < 35.0:
print('***Error: MET3 Density < 35%')
elif met3accum > 70.0:
@@ -391,13 +466,22 @@
print('')
print('MET4 Density:')
for y in range(0, ytiles - 9):
+ if y == ytiles - 10:
+ atotal = topadjust
+ else:
+ atotal = 100.0
for x in range(0, xtiles - 9):
+ if x == xtiles - 10:
+ if y == ytiles - 10:
+ atotal = corneradjust
+ else:
+ atotal = sideadjust
met4accum = 0
for w in range(y, y + 10):
base = xtiles * w + x
met4accum += sum(met4fill[base : base + 10])
- print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met4accum / 100.0))
+ print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met4accum / atotal))
if met4accum < 35.0:
print('***Error: MET4 Density < 35%')
elif met4accum > 70.0:
@@ -406,13 +490,22 @@
print('')
print('MET5 Density:')
for y in range(0, ytiles - 9):
+ if y == ytiles - 10:
+ atotal = topadjust
+ else:
+ atotal = 100.0
for x in range(0, xtiles - 9):
+ if x == xtiles - 10:
+ if y == ytiles - 10:
+ atotal = corneradjust
+ else:
+ atotal = sideadjust
met5accum = 0
for w in range(y, y + 10):
base = xtiles * w + x
met5accum += sum(met5fill[base : base + 10])
- print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met5accum / 100.0))
+ print('Tile (' + str(x) + ', ' + str(y) + '): ' + str(met5accum / atotal))
if met5accum < 45.0:
print('***Error: MET5 Density < 45%')
elif met5accum > 86.0:
