# Set to false for a short run
big_tests = false

verbose(true)

source(File.dirname(__FILE__) + "/drctest.gds", "TOPTOP")


0.1.um.to_s == "0.1" || raise("unexpected value")
0.1.micron.to_s == "0.1" || raise("unexpected value")
0.1.um2.to_s == "0.1" || raise("unexpected value")
0.1.mm2.to_s == "100000.0" || raise("unexpected value")
120.dbu.to_s == "0.12" || raise("unexpected value")
(0.1.um + 120.dbu).to_s == "0.22" || raise("unexpected value")
0.1.mm.to_s == "100.0" || raise("unexpected value")
1e-6.m.to_s == "1.0" || raise("unexpected value")
1.um.to_s == "1.0" || raise("unexpected value")
1.micron.to_s == "1.0" || raise("unexpected value")
1.um2.to_s == "1.0" || raise("unexpected value")
1.mm2.to_s == "1000000.0" || raise("unexpected value")
(1.um + 120.dbu).to_s == "1.12" || raise("unexpected value")
1.mm.to_s == "1000.0" || raise("unexpected value")

def run_testsuite(dm, ic, tiled = false)

  lb = 100
  
  a = input(RBA::LayerInfo::new(1, 0))
  b = input(2)
  c = input(3)
  x = input(10)
  y = input(11)

  h = {}
  layers.each { |l| h[l] = true }
  h[RBA::LayerInfo::new(1, 0)] || raise("missing layer 1/0 in layers list")
  
  c.is_merged? == false || raise("unexpected value")
  
  puts "--- general #{lb}"
  
  l1 = a&b
  l1.output(lb, dm)
  l1.is_empty? && raise("must not be empty")
  
  a.and(b).xor(l1).is_empty? || raise("xor not empty")
  tiled || (a.and(b).is_merged? == true || raise("unexpected value"))
  
  a.xor(b).output(RBA::LayerInfo::new(lb + 1, dm))
  a.xor(b).xor(a ^ b).is_empty? || raise("xor not empty")
  a.not(b).output(lb + 2, dm)
  a.not(b).xor(a - b).is_empty? || raise("xor not empty")
  a.or(b).output(lb + 3, dm)
  a.or(b).xor(a | b).is_empty? || raise("xor not empty")
  
  a.join(b).output(lb + 4, dm)
  a.join(b).xor(a + b).is_empty? || raise("xor not empty")

  if !tiled
    a.join(b).data.count == 16 * ic || raise("unexpected shape count")
  end
  
  c.raw.is_clean? == false || raise("unexpected value")
  # c.raw switched the semantics
  c.is_clean? == false || raise("unexpected value")
  (c.area / ic).to_s == "10.0" || raise("unexpected value")
  (c.edges.length / ic).to_s == "18.0" || raise("unexpected value")
  (c.perimeter / ic).to_s == "18.0" || raise("unexpected value")
  c.clean
  c.is_clean? == true || raise("unexpected value")
  (c.area / ic).to_s == "9.0" || raise("unexpected value")
  (c.edges.length / ic).to_s == "14.0" || raise("unexpected value")
  (c.perimeter / ic).to_s == "14.0" || raise("unexpected value")
  (c.dup.area / ic).to_s == "9.0" || raise("unexpected value")
  (c.raw.clean.area / ic).to_s == "9.0" || raise("unexpected value")
  (c.area / ic).to_s == "9.0" || raise("unexpected value")
  (c.raw.area / ic).to_s == "10.0" || raise("unexpected value")
  c.clean

  if ic == 1    
    c.bbox.to_s == "(-4,-5;0,-2)" || raise("unexpected value")
  end
  
  lb += 10
  puts "--- centers, start_segments, end_segments #{lb}"
  
  a.edges.centers(0.1).extended_out(0.01).output(lb, dm)
  a.edges.centers(0, 0.5).extended_out(0.01).output(lb + 1, dm)
  a.edges.centers(0.5, 0.5).extended_out(0.01).output(lb + 2, dm)
  a.edges.start_segments(0.1).extended_out(0.01).output(lb + 3, dm)
  a.edges.start_segments(0, 0.5).extended_out(0.01).output(lb + 4, dm)
  a.edges.start_segments(0.5, 0.5).extended_out(0.01).output(lb + 5, dm)
  a.edges.end_segments(0.1).extended_out(0.01).output(lb + 6, dm)
  a.edges.end_segments(0, 0.5).extended_out(0.01).output(lb + 7, dm)
  a.edges.end_segments(0.5, 0.5).extended_out(0.01).output(lb + 8, dm)
  
  a.polygons? == true || raise("unexpected value")
  a.edges? == false || raise("unexpected value")
  a.edge_pairs? == false || raise("unexpected value")
  a.edges.edge_pairs? == false || raise("unexpected value")
  a.edges.polygons? == false || raise("unexpected value")
  a.edges.edges? == true || raise("unexpected value")
  a.space(0.5).edge_pairs? == true || raise("unexpected value")
  a.space(0.5).polygons? == false || raise("unexpected value")
  a.space(0.5).edges? == false || raise("unexpected value")
  
  lb += 10
  puts "--- extended #{lb}"
  
  a.edges.with_length(0.6).extended(begin: 0.1, end: 0.15, out: 0.1, in: -0.05).output(lb, dm)
  a.edges.with_length(0.6).extended(out: 0.1, in: -0.05).output(lb + 1, dm)
  a.edges.with_length(0..0.7).extended(out: 0.1, in: -0.05).output(lb + 2, dm)
  a.edges.with_length(0..0.7).extended(out: 0.1, in: -0.05, joined: true).output(lb + 3, dm)
  a.edges.with_length(0..0.7).extended(begin: 0.1, end: 0.15, out: 0.1, in: 0.05, joined: true).output(lb + 4, dm)
  a.edges.with_length(0..0.7).extended(0.1, 0.15, 0.1, 0.05).output(lb + 5, dm)
  a.edges.with_length(0..0.7).extended(0.1, 0.15, 0.1, 0.05, true).output(lb + 6, dm)
  a.edges.with_length(0..0.7).extended_out(0.1).output(lb + 7, dm)
  a.edges.with_length(0..0.7).extended_in(0.1).output(lb + 8, dm)
  
  lb += 10
  puts "--- extents #{lb}"
  
  b.extents.output(lb, dm)
  
  lb += 10
  puts "--- first_edges, second_edges, edges #{lb}"
  
  a.enclosing(b, 0.5, whole_edges).first_edges.extended_out(0.01).output(lb, dm)
  a.enclosing(b, 0.5, whole_edges).second_edges.extended_out(0.01).output(lb + 1, dm)
  a.enclosing(b, 0.5, whole_edges).edges.extended_out(0.01).output(lb + 2, dm)
  
  lb += 10
  puts "--- hulls, holes #{lb}"
  
  c.xor(b).holes.output(lb, dm)
  c.xor(b).hulls.output(lb + 1, dm)
  c.xor(b).xor(c.xor(b).hulls - c.xor(b).holes).is_empty? || raise("xor not empty")
  
  lb += 10
  puts "--- interacting, in, not_in #{lb}"
  
  b.interacting(a).output(lb, dm)
  b.interacting(a).in(b).output(lb + 1, dm)
  (b|a).not_in(b).output(lb + 2, dm)
  x.in(b).output(lb + 3, dm)
  b.sized(0.1).in(b).is_empty? == true || raise("unexpected value")
  
  lb += 10
  puts "--- inside, outside, overlapping, interacting #{lb}"
  
  b.inside(c).output(lb, dm)
  b.outside(c).output(lb + 1, dm)
  b.overlapping(c).output(lb + 2, dm)
  b.interacting(c).output(lb + 3, dm)
  bdup = b.dup
  bdup.select_inside(c)
  bdup.xor(b.inside(c)).output(lb + 4, dm)
  bdup = b.dup
  bdup.select_outside(c)
  bdup.xor(b.outside(c)).output(lb + 5, dm)
  bdup = b.dup
  bdup.select_overlapping(c)
  bdup.xor(b.overlapping(c)).output(lb + 6, dm)
  bdup = b.dup
  bdup.select_interacting(c)
  bdup.xor(b.interacting(c)).output(lb + 7, dm)
  
  lb += 10
  puts "--- merge #{lb}"
  
  c.raw
  if !tiled
    c.merged.is_merged? == true || raise("unexpected value")
  end
  c.merged.output(lb, dm)
  c.merged(2).output(lb + 1, dm)
  c.merged(3).output(lb + 2, dm)
  c.clean
  cdup = c.dup
  cdup.merge.xor(c.merged).output(lb + 3, dm)
  cdup.xor(c.merged).output(lb + 4, dm)
  cdup = c.dup
  cdup.merge(2).xor(c.merged(2)).output(lb + 5, dm)
  cdup.xor(c.merged(2)).output(lb + 6, dm)
  
  lb += 10
  puts "--- move, transform, rotate, scale #{lb}"
  
  c.moved(0.2, -0.1).output(lb, dm)
  cdup = c.dup
  cdup.move(0.2, -0.1)
  cdup.xor(c.moved(0.2, -0.1)).is_empty? == true || raise("xor not empty")
  t = RBA::DCplxTrans::new(0.5)
  c.transformed(t).output(lb + 1, dm)
  cdup = c.dup
  cdup.transform(t)
  cdup.xor(c.transformed(t)).is_empty? == true || raise("xor not empty")
  c.scaled(0.5).output(lb + 2, dm)
  cdup = c.dup
  cdup.scale(0.5)
  cdup.xor(c.scaled(0.5)).is_empty? == true || raise("xor not empty")
  c.rotated(45).output(lb + 3, dm)
  cdup = c.dup
  cdup.rotate(45)
  cdup.xor(c.rotated(45)).is_empty? == true || raise("xor not empty")
  
  lb += 10
  puts "--- rectangles, rectilinear #{lb}"
  
  a.rectangles.output(lb, dm)
  a.non_rectangles.output(lb + 1, dm)
  a.rectilinear.output(lb + 2, dm)
  a.non_rectilinear.output(lb + 3, dm)
  c.raw
  c.non_rectangles.output(lb + 4, dm)
  c.clean
  c.rectangles.output(lb + 5, dm)
  
  lb += 10
  puts "--- enclosing #{lb}"
  
  a.enclosing(b, 0.5).polygons.output(lb, dm)
  a.enclosing(b, 0.5).polygons.xor(a.enc(b, 0.5).polygons).is_empty? == true || raise("xor not empty")
  a.enclosing(b, 0.5, euclidian).polygons.output(lb + 1, dm)
  a.enclosing(b, 0.5, square).polygons.output(lb + 2, dm)
  a.enclosing(b, 0.5, projection).polygons.output(lb + 3, dm)
  a.enclosing(b, 0.5, euclidian, whole_edges).polygons.output(lb + 4, dm)
  a.enclosing(b, 0.5, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  a.enclosing(b, 0.5, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  a.enclosing(b, 0.5, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- enclosing (edges) #{lb}"
  
  ae = a.edges
  be = b.edges
  ae.enclosing(be, 0.5).polygons.output(lb, dm)
  ae.enclosing(be, 0.5).polygons.xor(ae.enc(be, 0.5).polygons).is_empty? == true || raise("xor not empty")
  ae.enclosing(be, 0.5, euclidian).polygons.output(lb + 1, dm)
  ae.enclosing(be, 0.5, square).polygons.output(lb + 2, dm)
  ae.enclosing(be, 0.5, projection).polygons.output(lb + 3, dm)
  ae.enclosing(be, 0.5, euclidian, whole_edges).polygons.output(lb + 4, dm)
  ae.enclosing(be, 0.5, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  ae.enclosing(be, 0.5, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  ae.enclosing(be, 0.5, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- isolated #{lb}"
  
  b.isolated(0.4).polygons.output(lb, dm)
  b.isolated(0.4).polygons.xor(b.iso(0.4).polygons).is_empty? == true || raise("xor not empty")
  b.isolated(0.4, euclidian).polygons.output(lb + 1, dm)
  b.isolated(0.4, square).polygons.output(lb + 2, dm)
  b.isolated(0.4, projection).polygons.output(lb + 3, dm)
  b.isolated(0.4, euclidian, whole_edges).polygons.output(lb + 4, dm)
  b.isolated(0.4, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  b.isolated(0.4, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  b.isolated(0.4, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- notch #{lb}"
  
  b.notch(0.4).polygons.output(lb, dm)
  b.notch(0.4, euclidian).polygons.output(lb + 1, dm)
  b.notch(0.4, square).polygons.output(lb + 2, dm)
  b.notch(0.4, projection).polygons.output(lb + 3, dm)
  b.notch(0.4, euclidian, whole_edges).polygons.output(lb + 4, dm)
  b.notch(0.4, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  b.notch(0.4, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  b.notch(0.4, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- overlap #{lb}"
  
  a.overlap(b, 0.5).polygons.output(lb, dm)
  a.overlap(b, 0.5, euclidian).polygons.output(lb + 1, dm)
  a.overlap(b, 0.5, square).polygons.output(lb + 2, dm)
  a.overlap(b, 0.5, projection).polygons.output(lb + 3, dm)
  a.overlap(b, 0.5, euclidian, whole_edges).polygons.output(lb + 4, dm)
  a.overlap(b, 0.5, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  a.overlap(b, 0.5, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  a.overlap(b, 0.5, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- overlap (edges) #{lb}"
  
  ae = a.edges
  be = b.edges
  ae.overlap(be, 0.5).polygons.output(lb, dm)
  ae.overlap(be, 0.5, euclidian).polygons.output(lb + 1, dm)
  ae.overlap(be, 0.5, square).polygons.output(lb + 2, dm)
  ae.overlap(be, 0.5, projection).polygons.output(lb + 3, dm)
  ae.overlap(be, 0.5, euclidian, whole_edges).polygons.output(lb + 4, dm)
  ae.overlap(be, 0.5, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  ae.overlap(be, 0.5, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  ae.overlap(be, 0.5, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- separation #{lb}"
  
  a.separation(b, 0.4).polygons.output(lb, dm)
  a.separation(b, 0.4).polygons.xor(a.sep(b, 0.4).polygons).is_empty? == true || raise("xor not empty")
  a.separation(b, 0.4, euclidian).polygons.output(lb + 1, dm)
  a.separation(b, 0.4, square).polygons.output(lb + 2, dm)
  a.separation(b, 0.4, projection).polygons.output(lb + 3, dm)
  a.separation(b, 0.4, euclidian, whole_edges).polygons.output(lb + 4, dm)
  a.separation(b, 0.4, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  a.separation(b, 0.4, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  a.separation(b, 0.4, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- sized #{lb}"
  
  x.sized(0.1).output(lb, dm)
  x.sized(0.1, 0).output(lb + 1, dm)
  x.sized(0.1, diamond_limit).output(lb + 2, dm)
  x.sized(0.1, octagon_limit).output(lb + 3, dm)
  x.sized(0.1, square_limit).output(lb + 4, dm)
  x.sized(0.1, acute_limit).output(lb + 5, dm)
  x.sized(0.1, no_limit).output(lb + 6, dm)
  xdup = x.dup
  xdup.size(0.1, 0)
  xdup.xor(x.sized(0.1, 0)).is_empty? == true || raise("xor not empty")
  a.sized(4.0, no_limit).output(lb + 7, dm)
  
  lb += 10
  puts "--- space #{lb}"
  
  b.space(0.4).polygons.output(lb, dm)
  b.space(0.4, euclidian).polygons.output(lb + 1, dm)
  b.space(0.4, square).polygons.output(lb + 2, dm)
  b.space(0.4, projection).polygons.output(lb + 3, dm)
  b.space(0.4, euclidian, whole_edges).polygons.output(lb + 4, dm)
  b.space(0.4, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  b.space(0.4, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  b.space(0.4, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- space (edges) #{lb}"
  
  be = b.edges
  be.space(0.4).polygons.output(lb, dm)
  be.space(0.4, euclidian).polygons.output(lb + 1, dm)
  be.space(0.4, square).polygons.output(lb + 2, dm)
  be.space(0.4, projection).polygons.output(lb + 3, dm)
  be.space(0.4, euclidian, whole_edges).polygons.output(lb + 4, dm)
  be.space(0.4, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  be.space(0.4, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  be.space(0.4, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- width #{lb}"
  
  b.width(0.4).polygons.output(lb, dm)
  b.width(0.4, euclidian).polygons.output(lb + 1, dm)
  b.width(0.4, square).polygons.output(lb + 2, dm)
  b.width(0.4, projection).polygons.output(lb + 3, dm)
  b.width(0.4, euclidian, whole_edges).polygons.output(lb + 4, dm)
  b.width(0.4, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  b.width(0.4, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  b.width(0.4, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- width (edges) #{lb}"
  
  be = b.edges
  be.width(0.4).polygons.output(lb, dm)
  be.width(0.4, euclidian).polygons.output(lb + 1, dm)
  be.width(0.4, square).polygons.output(lb + 2, dm)
  be.width(0.4, projection).polygons.output(lb + 3, dm)
  be.width(0.4, euclidian, whole_edges).polygons.output(lb + 4, dm)
  be.width(0.4, euclidian, projection_limits(0.4, nil)).polygons.output(lb + 5, dm)
  be.width(0.4, euclidian, projection_limits(nil, 0.4)).polygons.output(lb + 6, dm)
  be.width(0.4, euclidian, projection_limits(0..0.4)).polygons.output(lb + 7, dm)
  
  lb += 10
  puts "--- with_angle #{lb}"
  
  a.edges.with_angle(90.degree).extended_out(0.01).output(lb, dm)
  a.edges.with_angle(45 .. 90.1).extended_out(0.01).output(lb + 1, dm)
  a.edges.with_angle(0, 45.1).extended_out(0.01).output(lb + 2, dm)
  a.edges.without_angle(90.degree).extended_out(0.01).output(lb + 3, dm)
  a.edges.without_angle(45 .. 90.1).extended_out(0.01).output(lb + 4, dm)
  a.edges.without_angle(0, 45.1).extended_out(0.01).output(lb + 5, dm)
  a.with_angle(135).polygons.output(lb + 6, dm)
  
  lb += 10
  puts "--- with_area #{lb}"
  
  b.with_area(0.49.um2).output(lb, dm)
  b.with_area(0 .. 0.5).output(lb + 1, dm)
  b.with_area(0.5, nil).output(lb + 2, dm)
  b.without_area(0.49.um2).output(lb + 3, dm)
  b.without_area(nil, 0.5).output(lb + 4, dm)
  b.without_area(0.5, nil).output(lb + 5, dm)
  
  lb += 10
  puts "--- with_perimeter #{lb}"
  
  b.with_perimeter(2.8.um).output(lb, dm)
  b.with_perimeter(0 .. 3.um).output(lb + 1, dm)
  b.with_perimeter(3.um, nil).output(lb + 2, dm)
  b.without_perimeter(2.8.um).output(lb + 3, dm)
  b.without_perimeter(nil, 3.um).output(lb + 4, dm)
  b.without_perimeter(3.um, nil).output(lb + 5, dm)
  
  lb += 10
  puts "--- with_bbox_height #{lb}"
  
  b.with_bbox_height(0.7.um).output(lb, dm)
  b.with_bbox_height(0 .. 0.7.um+1.dbu).output(lb + 1, dm)
  b.with_bbox_height(0.7.um, nil).output(lb + 2, dm)
  b.without_bbox_height(0.7.um).output(lb + 3, dm)
  b.without_bbox_height(nil, 0.7.um+1.dbu).output(lb + 4, dm)
  b.without_bbox_height(0.7.um, nil).output(lb + 5, dm)
  
  lb += 10
  puts "--- with_bbox_width #{lb}"
  
  b.with_bbox_width(0.7.um).output(lb, dm)
  b.with_bbox_width(0 .. 0.7.um+1.dbu).output(lb + 1, dm)
  b.with_bbox_width(0.7.um, nil).output(lb + 2, dm)
  b.without_bbox_width(0.7.um).output(lb + 3, dm)
  b.without_bbox_width(nil, 0.7.um+1.dbu).output(lb + 4, dm)
  b.without_bbox_width(0.7.um, nil).output(lb + 5, dm)
  
  lb += 10
  puts "--- with_bbox_min #{lb}"
  
  b.with_bbox_min(0.7.um).output(lb, dm)
  b.with_bbox_min(0 .. 0.7.um+1.dbu).output(lb + 1, dm)
  b.with_bbox_min(0.7.um, nil).output(lb + 2, dm)
  b.without_bbox_min(0.7.um).output(lb + 3, dm)
  b.without_bbox_min(nil, 0.7.um+1.dbu).output(lb + 4, dm)
  b.without_bbox_min(0.7.um, nil).output(lb + 5, dm)
  
  lb += 10
  puts "--- with_bbox_max #{lb}"
  
  b.with_bbox_max(0.7.um).output(lb, dm)
  b.with_bbox_max(0 .. 0.7.um+1.dbu).output(lb + 1, dm)
  b.with_bbox_max(0.7.um, nil).output(lb + 2, dm)
  b.without_bbox_max(0.7.um).output(lb + 3, dm)
  b.without_bbox_max(nil, 0.7.um+1.dbu).output(lb + 4, dm)
  b.without_bbox_max(0.7.um, nil).output(lb + 5, dm)
  
  lb += 10
  puts "--- with_length #{lb}"
  
  a.edges.with_length(0.6.um).extended_out(0.01).output(lb, dm)
  a.edges.with_length(0 .. 0.6.um+1.dbu).extended_out(0.01).output(lb + 1, dm)
  a.edges.with_length(0.6, nil).extended_out(0.01).output(lb + 2, dm)
  a.edges.without_length(0.6.um).extended_out(0.01).output(lb + 3, dm)
  a.edges.without_length(0 .. 0.6.um+1.dbu).extended_out(0.01).output(lb + 4, dm)
  a.edges.without_length(0.6, nil).extended_out(0.01).output(lb + 5, dm)
  
  lb += 10
  puts "--- ongrid, snap #{lb}"
  
  a.ongrid(0.2).polygons(1.dbu).output(lb, dm)
  a.ongrid(0, 0.1).polygons(1.dbu).output(lb + 1, dm)
  a.snapped(0.2).output(lb + 2, dm)
  a.snapped(0, 0.1).output(lb + 3, dm)
  adup = a.dup
  adup.snap(0.2)
  adup.xor(a.snapped(0.2)).output(lb + 4, dm)
  
  lb += 10
  puts "--- odd polygon check #{lb}"
  
  y.odd_polygons.output(lb, dm)
  
end

cell("TOP")

x = polygon_layer
x.is_empty? == true || raise("unexpected value")
x.is_box? == false || raise("unexpected value")
x.insert(box(4.0, 0, 4.7, 0.7))
x.is_empty? == false || raise("unexpected value")
x.is_box? == true || raise("unexpected value")
x.insert(polygon([ p(0, 0), p(2.0, 0), p(1.0, 1.0) ]))
x.insert(polygon([ p(0, -5.0), p(2.0, -5.0), p(1.0, -6.0) ]))
x.insert(path([ p(0, -2), p(2.0, -2) ], 0.2))
x.is_box? == false || raise("unexpected value")
x.output(10, 0)

y = polygon_layer
y.insert(polygon([ p(0, 0), p(0, 1.0), p(2.0, 1.0), p(2.0, 2.0), p(1.0, 2.0), p(1.0, 0) ]))
y.output(11, 0)

a1 = input(1)
b1 = input(2)
c1 = input(3)
  
puts "=== Simple tests ==="

cell("TOP")

c1.rounded_corners(0.5, 0.5, 16).output(1010, 0)
c1.smoothed(1.5).output(1011, 0)

a1.texts.output(1020, 0)
a1.texts("A*").output(1021, 0)
a1.texts(text("A*")).output(1022, 0)
a1.texts(pattern("A*")).output(1023, 0)
a1.texts(pattern("A*"), as_dots).extended(0.05, 0.05, 0.05, 0.05).output(1024, 0)
a1.texts(pattern("A*"), as_boxes).output(1025, 0)

a1.middle.sized(0.05).output(1030, 0)
a1.middle(as_dots).extended(0.05, 0.05, 0.05, 0.05).output(1031, 0)
a1.middle(as_boxes).sized(0.05).output(1032, 0)
a1.extent_refs(0.5, 0.5).sized(0.05).output(1040, 0)
a1.extent_refs(:center).sized(0.05).output(1040, 1)
a1.extent_refs(:c).sized(0.05).output(1040, 2)
a1.extent_refs(:bottom).sized(0.05).output(1041, 0)
a1.extent_refs(:b).sized(0.05).output(1041, 1)
a1.extent_refs(:b, as_edges).extended(0.05, 0.05, 0.05, 0.05).output(1041, 2)
a1.extent_refs(:top).sized(0.05).output(1042, 0)
a1.extent_refs(:t).sized(0.05).output(1042, 1)
a1.extent_refs(:left).sized(0.05).output(1043, 0)
a1.extent_refs(:l).sized(0.05).output(1043, 1)
a1.extent_refs(:right).sized(0.05).output(1044, 0)
a1.extent_refs(:r).sized(0.05).output(1044, 1)
a1.extent_refs(:bottom_left).sized(0.05).output(1045, 0)
a1.extent_refs(:bl).sized(0.05).output(1045, 1)
a1.extent_refs(:bottom_center).sized(0.05).output(1046, 0)
a1.extent_refs(:bc).sized(0.05).output(1046, 1)
a1.extent_refs(:bottom_right).sized(0.05).output(1047, 0)
a1.extent_refs(:br).sized(0.05).output(1047, 1)
a1.extent_refs(:top_left).sized(0.05).output(1048, 0)
a1.extent_refs(:tl).sized(0.05).output(1048, 1)
a1.extent_refs(:top_center).sized(0.05).output(1049, 0)
a1.extent_refs(:tc).sized(0.05).output(1049, 1)
a1.extent_refs(:top_right).sized(0.05).output(1050, 0)
a1.extent_refs(:tr).sized(0.05).output(1050, 1)
a1.extent_refs(:left_center).sized(0.05).output(1051, 0)
a1.extent_refs(:lc).sized(0.05).output(1051, 1)
a1.extent_refs(:right_center).sized(0.05).output(1052, 0)
a1.extent_refs(:rc).sized(0.05).output(1052, 1)
a1.extent_refs(0.25, 0.5, 0.5, 0.75).output(1053, 0)

a1.corners.sized(0.05).output(1060, 0)
a1.corners(-90.0, as_boxes).sized(0.05).output(1061, 0)
a1.corners(-90.0, as_dots).extended(0.05, 0.05, 0.05, 0.05).output(1062, 0)

a1.select { |p| p.bbox.width < 0.8 }.output(1100, 0)
a1.collect { |p| p.is_box? && p.bbox.enlarged(0.1, 0.1) }.output(1101, 0)
a1.collect_to_region { |p| p.is_box? && p.bbox.enlarged(0.1, 0.1) }.output(1102, 0)
a1.collect_to_edges { |p| p.is_box? && p.bbox.enlarged(0.1, 0.1) }.output(1103, 0)
a1.collect { |p| p.is_box? && p.bbox.transformed(RBA::VCplxTrans::new(1000.0)).enlarged(120, 120) }.output(1104, 0)
a1.collect { |p| p.is_box? && [ p.bbox.transformed(RBA::VCplxTrans::new(1000.0)).enlarged(150, 150), p.bbox.transformed(RBA::VCplxTrans::new(1000.0)).enlarged(120, 120) ] }.output(1105, 0)
lx = polygon_layer
a1.each { |p| p.is_box? && lx.insert(p) }
lx.output(1106, 0)
a1.collect { |p| p.is_box? && RBA::Region::new(p.bbox.transformed(RBA::VCplxTrans::new(1000.0))).sized(120) }.output(1107, 0)
a1.collect { |p| p.is_box? && RBA::Polygon::new(p.bbox.transformed(RBA::VCplxTrans::new(1000.0))) }.output(1108, 0)
a1.collect { |p| p.is_box? && RBA::DPolygon::new(p.bbox) }.output(1109, 0)
a1.collect { |p| p.is_box? && RBA::SimplePolygon::new(p.bbox.transformed(RBA::VCplxTrans::new(1000.0))) }.output(1110, 0)
a1.collect { |p| p.is_box? && RBA::DSimplePolygon::new(p.bbox) }.output(1111, 0)
a1.collect_to_edges { |p| p.is_box? && p.bbox.transformed(RBA::VCplxTrans::new(1000.0)).enlarged(120, 120) }.output(1112, 0)
a1.collect_to_edges { |p| p.is_box? && [ p.bbox.transformed(RBA::VCplxTrans::new(1000.0)).enlarged(150, 150), p.bbox.transformed(RBA::VCplxTrans::new(1000.0)).enlarged(120, 120) ] }.output(1113, 0)
a1.collect_to_edges { |p| p.is_box? && RBA::Region::new(p.bbox.transformed(RBA::VCplxTrans::new(1000.0))).sized(120) }.output(1114, 0)
a1.collect_to_edges { |p| p.is_box? && RBA::Polygon::new(p.bbox.transformed(RBA::VCplxTrans::new(1000.0))) }.output(1115, 0)
a1.collect_to_edges { |p| p.is_box? && RBA::DPolygon::new(p.bbox) }.output(1116, 0)
a1.collect_to_edges { |p| p.is_box? && RBA::SimplePolygon::new(p.bbox.transformed(RBA::VCplxTrans::new(1000.0))) }.output(1117, 0)
a1.collect_to_edges { |p| p.is_box? && RBA::DSimplePolygon::new(p.bbox) }.output(1118, 0)

a1.edges.select { |p| p.length < 0.8 }.output(1120, 0)
a1.edges.collect { |p| p.length < 0.8 && p.transformed(RBA::VCplxTrans::new(1000.0)) }.output(1121, 0)
a1.edges.collect_to_region { |p| p.length < 0.8 && p.bbox.enlarged(0.1, 0.1) }.output(1122, 0)
a1.edges.collect_to_region { |p| p.length < 0.8 && p.bbox.transformed(RBA::VCplxTrans::new(1000.0)).enlarged(100, 100) }.output(1123, 0)

# edge pair collect
a1.width(1.5).collect { |p| p.transformed(RBA::VCplxTrans::new(1000.0)) }.output(1120, 0)
a1.width(1.5).collect_to_edge_pairs { |p| p.transformed(RBA::VCplxTrans::new(1000.0)) }.output(1121, 0)

puts "=== Single-cell testsuite ==="

run_testsuite(0, 1)

output_cell("NEW")

puts "=== Foreign-cell test ==="

a1.and(b1).output(1000, 0)

cell("TOPTOP")

a1.and(b1).output(1001, 0)

if big_tests

  puts "=== Multi-cell testsuite ==="

  run_testsuite(1, 900)

  puts "=== Multi-cell testsuite (tiled) ==="

  tiles(4.1, 5.2)
  # must be that big ...
  tile_borders(1.um, 1.um)
  threads(4)

  run_testsuite(2, 900, true)

  puts "=== Multi-cell testsuite (pseudo-tiled) ==="

  tiles(10000.0, 10000.0)
  tile_borders(0, 0)
  threads(4)

  run_testsuite(3, 900, true)

  flat
  threads(1)

  puts "=== XOR ==="

  (100..900).each do |lo|
    a = input(lo, 0)
    b = input(lo, 1)
    if !a.is_empty? || !b.is_empty?
      a.xor(b).output(lo, 10)
    end
    a = input(lo, 1)
    b = input(lo, 2)
    if !a.is_empty? || !b.is_empty?
      a.xor(b).output(lo, 11)
    end
    a = input(lo, 1)
    b = input(lo, 3)
    if !a.is_empty? || !b.is_empty?
      a.xor(b).output(lo, 12)
    end
  end

end

source.layout.write("drctest_out.gds")