Demo Readonlymemory error for Triangulate

joyal2012_hybrid.jl

using PSSFSS

# From M.-A. Joyal and J.-J. Laurin, "Analysis and Design of Thin Circular Polarizers Based on
# Meander Lines", IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 60, NO. 6, JUNE 2012,
# pp. 3007-3011, Fig. 7 (hybrid design)

meand = meander(ntri=600, units=mil, a=160, b=201, h=100, w1=23, w2=23, rot=45)
rect = rectstrip(Nx=5, Ny=20, units=mil, Px=160, Py=201, Lx=31, Ly=139, rot=45)
pyralux = Layer(ϵᵣ=3.4, tanδ=0.002, width=5mil)
spacer = Layer(ϵᵣ=1.1, tanδ=0.01, width=74mil) # λ/8 in Rohacell 71 HF at 19 GHz

strata = [
      Layer()
      meand
      Layer(width=0.685mil)
      meand
      pyralux
      spacer
      rect
      Layer(width=0.685mil)
      rect
      pyralux
      spacer
      meand
      Layer(width=0.685mil)
      meand
      pyralux
      Layer()
      ]

flist = 16:23
steering = (θ=0, ϕ=0)
analyze(strata, flist, steering, logfile=devnull, resultfile=devnull)
nothing

li2009.jl

using PSSFSS

# From Figure 5 of Li, Jiao, and Zhao: "A Novel Microstrip Rectangular-Patch/Ring-
# Combination Reflectarray Element and Its Application", IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS, 
# VOL. 8, 2009, pp. 1119-1112.


L = 17 # Period of square lattice (mm)
k1 = 0.5
k2 = 0.125
FGHz = 10.0
L2s = LinRange(2.0, 16.5, 6)
ntris = round.(Int, LinRange(400,1000,6))
sheets = []
for (L2,ntri) in zip(L2s, ntris)
    println("L2 = $L2")
    L1 = k1 * L2
    w = k2 * L2
    a = [-1.0, (L2-2w)/√2]
    b = [L1/√2, L2/√2]
    sheet = polyring(sides=4, s1=[L, 0], s2=[0, L], ntri=ntri, orient=45, a=a, b=b, units=mm)
    push!(sheets,sheet)
    mystrata = [Layer(width=-7mm)
              sheet
              Layer(ϵᵣ=2.65, width=1mm)
              Layer(ϵᵣ=1.07, width=6mm)
              pecsheet() # PEC Wall
              Layer()]
    results = analyze(mystrata, FGHz, (ϕ=0, θ=0), logfile=devnull, resultfile=devnull,showprogress=false)
end

showerror.jl

for n in 1:20
    println("Trial ", n)
    include("joyal2012_hybrid.jl")
    include("li2009.jl")
    println("\n\n")
end