1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
|
function export_empire( CSX, FDTD, filename, options )
% export_empire( CSX, FDTD, filename, options )
%
% Exports the geometry defined in CSX to filename as an Empire python file.
%
% CSX: CSX-object created by InitCSX()
% FDTD: FDTD-object created by InitFDTD()
% filename: export filename (e.g. '/tmp/export.py')
% options (optional): struct
% options.ignore : cell array with property names to ignore
%
% CSXCAD matlab interface
% -----------------------
% author: Sebastian Held <sebastian.held@gmx.de>
% 17. Jun 2010: initial version
% 28. Sep 2010: rewritten using Empire python scripting
% See also InitCSX InitFDTD
if nargin < 3
options = [];
end
fid = fopen( filename, 'w' );
% write header
fprintf( fid, ['# Empire python script\n\n' ...
'# start Empire and select File/"Run Script" and select this file\n\n'] );
fprintf( fid, 'from scripting import *\n\n' );
export_empire_meshlines( fid, CSX );
export_empire_settings( fid, CSX, FDTD );
if isfield(CSX.Properties,'Material')
% process material
Material = CSX.Properties.Material;
options.Property = 'Material';
process_primitives( fid, Material, options );
end
if isfield(CSX.Properties,'Metal')
% process PEC
Metal = CSX.Properties.Metal;
options.Property = 'Metal';
process_primitives( fid, Metal, options );
end
fclose( fid );
% -----------------------------------------------------------------------------
function str = gym_vect( start,stop )
str = ['[' num2str(start(1)) ',' num2str(stop(1)) ',' num2str(start(2)) ',' num2str(stop(2)) ',' num2str(start(3)) ',' num2str(stop(3)) ']'];
function str = gym_vector( v )
str = ['[' num2str(v(1))];
for a=2:numel(v), str = [str ', ' num2str(v(a))]; end
str = [str ']'];
% -----------------------------------------------------------------------------
function primitive_box( fid, CSX_box, layertype )
%primitive_box( fid, CSX_box, layertype )
% layertype may be: 'dielectric Dielectric' or 'conductor Conductor'
global current_layer_prio
properties = '';
if current_layer_prio ~= CSX_box.ATTRIBUTE.Priority
properties = [properties ' prio ' num2str(CSX_box.ATTRIBUTE.Priority)];
end
if ~isempty(properties)
fprintf( fid, 'g.objecttext("%s")\n', [layertype properties] );
end
start = [CSX_box.P1.ATTRIBUTE.X CSX_box.P1.ATTRIBUTE.Y CSX_box.P1.ATTRIBUTE.Z];
stop = [CSX_box.P2.ATTRIBUTE.X CSX_box.P2.ATTRIBUTE.Y CSX_box.P2.ATTRIBUTE.Z];
fprintf( fid, 'g.layerextrude("z",%g,%g)\n', start(3), stop(3) );
fprintf( fid, 'g.box(%g,%g,%g,%g)\n', start(1), start(2), stop(1), stop(2) );
% -----------------------------------------------------------------------------
function primitive_cylinder( fid, CSX_cylinder, layertype )
%primitive_cylinder( fid, CSX_cylinder, layertype )
% layertype may be: 'dielectric Dielectric' or 'conductor Conductor'
global current_layer_prio
properties = '';
if current_layer_prio ~= CSX_cylinder.ATTRIBUTE.Priority
properties = [properties ' prio ' num2str(CSX_cylinder.ATTRIBUTE.Priority)];
end
if ~isempty(properties)
fprintf( fid, 'g.objecttext("%s")\n', [layertype properties] );
end
start = [CSX_cylinder.P1.ATTRIBUTE.X CSX_cylinder.P1.ATTRIBUTE.Y CSX_cylinder.P1.ATTRIBUTE.Z];
stop = [CSX_cylinder.P2.ATTRIBUTE.X CSX_cylinder.P2.ATTRIBUTE.Y CSX_cylinder.P2.ATTRIBUTE.Z];
radius = CSX_cylinder.ATTRIBUTE.Radius;
if start([1 2]) == stop([1 2])
fprintf( fid, 'g.layerextrude("z",%g,%g)\n', start(3), stop(3) );
fprintf( fid, 'g.circle(%g,%g,%g)\n', start(1), start(2), radius );
elseif start([1 3]) == stop([1 3])
fprintf( fid, 'g.layerextrude("y",%g,%g)\n', start(2), stop(2) );
fprintf( fid, 'g.circle(%g,%g,%g)\n', start(3), start(1), radius );
elseif start([2 3]) == stop([2 3])
fprintf( fid, 'g.layerextrude("x",%g,%g)\n', start(1), stop(1) );
fprintf( fid, 'g.circle(%g,%g,%g)\n', start(2), start(3), radius );
else
disp( ['cylinder coordinates: (' num2str(start) ') -> (' num2str(stop) ')'] );
error 'cylinder orientation not supported'
end
% -----------------------------------------------------------------------------
function str = primitive_wire( CSX_wire, options )
radius = CSX_wire.ATTRIBUTE.WireRadius;
str = ['sphere_sweep { linear_spline, ' num2str(numel(CSX_wire.Vertex))];
for a=1:numel(CSX_wire.Vertex)
% iterate over all vertices
v = [CSX_wire.Vertex{a}.ATTRIBUTE.X CSX_wire.Vertex{a}.ATTRIBUTE.Y CSX_wire.Vertex{a}.ATTRIBUTE.Z];
str = [str ', ' pov_vect(v) ', ' num2str(radius)];
end
str = [str ' ' options '}'];
% -----------------------------------------------------------------------------
function primitive_polygon( fid, CSX_polygon, layertype )
%primitive_polygon( fid, CSX_polygon, layertype )
% layertype may be: 'dielectric Dielectric' or 'conductor Conductor'
global current_layer_prio
properties = '';
if current_layer_prio ~= CSX_polygon.ATTRIBUTE.Priority
properties = [properties ' prio ' num2str(CSX_polygon.ATTRIBUTE.Priority)];
end
if ~isempty(properties)
fprintf( fid, 'g.objecttext("%s")\n', [layertype properties] );
end
Elevation = CSX_polygon.ATTRIBUTE.Elevation;
if (CSX_polygon.NormDir.ATTRIBUTE.X ~= 0)
NormDir = 'x';
elseif (CSX_polygon.NormDir.ATTRIBUTE.Y ~= 0)
NormDir = 'y';
elseif (CSX_polygon.NormDir.ATTRIBUTE.Z ~= 0)
NormDir = 'z';
end
coords = [];
for a=1:numel(CSX_polygon.Vertex)
% iterate over all vertices
coords = [coords CSX_polygon.Vertex{a}.ATTRIBUTE.X1 CSX_polygon.Vertex{a}.ATTRIBUTE.X2];
end
fprintf( fid, 'g.layerextrude("%s",%g,%g)\n', NormDir, Elevation, Elevation+1 );
fprintf( fid, 'g.poly(' );
for n=1:numel(coords)
if n > 1
fprintf( fid, ', %g', coords(n) );
else
fprintf( fid, '%g', coords(n) );
end
end
fprintf( fid, ')\n' );
% -----------------------------------------------------------------------------
function process_primitives( fid, prop, options )
global current_layer_prio
ignore = {};
if isfield(options,'ignore'), ignore = options.ignore; end
% iterate over all properties and extract the priority
prio = [];
for num=1:numel(prop)
if isfield(prop{num}.Primitives,'Box')
for a=1:numel(prop{num}.Primitives.Box)
% iterate over all boxes
prio(end+1) = prop{num}.Primitives.Box{a}.ATTRIBUTE.Priority;
end
end
if isfield(prop{num}.Primitives,'Cylinder')
for a=1:numel(prop{num}.Primitives.Cylinder)
% iterate over all cylinders
prio(end+1) = prop{num}.Primitives.Cylinder{a}.ATTRIBUTE.Priority;
end
end
if isfield(prop{num}.Primitives,'Polygon')
for a=1:numel(prop{num}.Primitives.Polygon)
% iterate over all polygons
prio(end+1) = prop{num}.Primitives.Polygon{a}.ATTRIBUTE.Priority;
end
end
end
% Now prio is a vector full of priorities. Extract the priority, which is
% most often used.
u_prios = unique(prio);
count = histc(prio,u_prios);
[~,idx] = max(count);
current_layer_prio = u_prios(idx);
% now export all properties/primitives
for num=1:numel(prop)
Name = prop{num}.ATTRIBUTE.Name;
if any( strcmp(Name,ignore) )
disp( ['omitting ' Name '...'] );
continue
end
disp( ['processing ' prop{num}.ATTRIBUTE.Name '...'] );
if strcmp(options.Property,'Metal')
layertype = 'conductor';
properties = [layertype ' sigma infinite sides auto rough 0 prio ' num2str(current_layer_prio) ' automodel auto'];
else
layertype = 'dielectric';
epsr = 1;
sigma = 0;
if isfield(prop{num},'PropertyX')
if isfield(prop{num}.PropertyX.ATTRIBUTE,'Epsilon'), epsr = prop{num}.PropertyX.ATTRIBUTE.Epsilon; end
if isfield(prop{num}.PropertyX.ATTRIBUTE,'Kappa'), sigma = prop{num}.PropertyX.ATTRIBUTE.Kappa; end
end
properties = [layertype ' epsr ' num2str(epsr) ' sigma ' num2str(sigma) ' prio ' num2str(current_layer_prio) ' tand 0 density 0'];
end
export_empire_layer( fid, Name, properties )
if isfield(prop{num}.Primitives,'Box')
for a=1:numel(prop{num}.Primitives.Box)
% iterate over all boxes
Box = prop{num}.Primitives.Box{a};
primitive_box( fid, Box, layertype );
end
end
if isfield(prop{num}.Primitives,'Cylinder')
for a=1:numel(prop{num}.Primitives.Cylinder)
% iterate over all cylinders
Cylinder = prop{num}.Primitives.Cylinder{a};
primitive_cylinder( fid, Cylinder, layertype );
end
end
if isfield(prop{num}.Primitives,'Wire')
for a=1:numel(prop{num}.Primitives.Wire)
% iterate over all wires
Wire = prop{num}.Primitives.Wire{a};
% str = primitive_wire( Wire, obj_modifier );
% fprintf( fid, '%s\n', str );
end
warning( 'CSXCAD:export_empire', 'primitive wire currently unsupported' );
end
if isfield(prop{num}.Primitives,'Polygon')
for a=1:numel(prop{num}.Primitives.Polygon)
% iterate over all polygons
Polygon = prop{num}.Primitives.Polygon{a};
primitive_polygon( fid, Polygon, layertype );
end
end
fprintf( fid, 'g.load()\n\n' );
end
function export_empire_layer( fid, name, layertext )
fprintf( fid, '# create a new layer\n' );
fprintf( fid, 'g = gmf()\n' );
fprintf( fid, 'g.layer("%s")\n', name );
fprintf( fid, 'g.layerextrude("z",0,0)\n' );
fprintf( fid, 'g.layertext("%s")\n', layertext );
function export_empire_meshlines(fid,CSX)
fprintf( fid, '# create the mesh\n' );
fprintf( fid, 'empty_disc()\n' );
fprintf( fid, 'g = gmf()\n' );
fprintf( fid, 'g.raw_write("DISC X %g")\n', CSX.RectilinearGrid.XLines );
fprintf( fid, 'g.raw_write("DISC Y %g")\n', CSX.RectilinearGrid.YLines );
fprintf( fid, 'g.raw_write("DISC Z %g")\n', CSX.RectilinearGrid.ZLines );
fprintf( fid, 'g.load()\n\n' );
function empire_BC = convert_BC( openEMS_BC )
if ischar(openEMS_BC)
if strcmp(openEMS_BC,'PEC')
empire_BC = 'electric';
elseif strcmp(openEMS_BC,'PMC')
empire_BC = 'magnetic';
elseif strncmp(openEMS_BC,'MUR',3) || strncmp(openEMS_BC,'PML',3)
empire_BC = 'pml 6';
else
empire_BC = 'UNKNOWN';
end
else
if openEMS_BC == 0
empire_BC = 'electric';
elseif openEMS_BC == 1
empire_BC = 'magnetic';
elseif openEMS_BC == 2 || openEMS_BC == 3
empire_BC = 'pml 6';
else
empire_BC = 'UNKNOWN';
end
end
function export_empire_settings(fid,CSX,FDTD)
fprintf( fid, '# settings\n' );
fprintf( fid, 'set_fs_parameter( ''dx'', %g )\n', CSX.RectilinearGrid.ATTRIBUTE.DeltaUnit );
fprintf( fid, 'set_fs_parameter( ''xmin'', ''%s'' )\n', convert_BC(FDTD.BoundaryCond.ATTRIBUTE.xmin) );
fprintf( fid, 'set_fs_parameter( ''xmax'', ''%s'' )\n', convert_BC(FDTD.BoundaryCond.ATTRIBUTE.xmax) );
fprintf( fid, 'set_fs_parameter( ''ymin'', ''%s'' )\n', convert_BC(FDTD.BoundaryCond.ATTRIBUTE.ymin) );
fprintf( fid, 'set_fs_parameter( ''ymax'', ''%s'' )\n', convert_BC(FDTD.BoundaryCond.ATTRIBUTE.ymax) );
fprintf( fid, 'set_fs_parameter( ''zmin'', ''%s'' )\n', convert_BC(FDTD.BoundaryCond.ATTRIBUTE.zmin) );
fprintf( fid, 'set_fs_parameter( ''zmax'', ''%s'' )\n', convert_BC(FDTD.BoundaryCond.ATTRIBUTE.zmax) );
fprintf( fid, 'set_fs_parameter( ''f_0'', %g )\n', FDTD.Excitation.ATTRIBUTE.f0 );
fprintf( fid, 'set_fs_parameter( ''f_c'', %g )\n', FDTD.Excitation.ATTRIBUTE.fc );
f_start = max(0,FDTD.Excitation.ATTRIBUTE.f0 - FDTD.Excitation.ATTRIBUTE.fc);
f_end = max(0,FDTD.Excitation.ATTRIBUTE.f0 + FDTD.Excitation.ATTRIBUTE.fc);
fprintf( fid, 'set_fs_parameter( ''F_START'', %g )\n', f_start );
fprintf( fid, 'set_fs_parameter( ''F_END'', %g )\n', f_end );
fprintf( fid, 'set_fs_parameter( ''F_MID'', %g )\n', (f_start+f_end)/2 );
|
