![](\"\/forum\/wp-content\/uploads\/sites\/2\/2023\/06\/25-06-2023-1687676782-mceclip0.png\")
<\/p>\n\u8001\u5e08\u60a8\u597d\uff0c\u6211\u6700\u8fd1\u770b\u4e86\u5b98\u7f51\u6848\u4f8b\uff081x2_MMI\u548c1x4_CWDM_2D\uff09\uff0c\u4e8e\u662f\u6211\u5c31\u60f3\u5c1d\u8bd5\u4e00\u4e0b1×4\u529f\u5206\u5668\u4ef6\uff0c\u4f46\u662f\u5728\u4f18\u5316\u8fed\u4ee3\u7684\u8fc7\u7a0b\u4e2d\u9047\u5230\u4e86\u65e0\u6cd5\u89e3\u91ca\u7684\u95ee\u9898\u3002
\u4f18\u5316\u811a\u672c\u4ee3\u7801\u4e3a\uff1a<\/p>\n
def runSim(params, eps_bg, eps_wg, x_pos, y_pos, z_pos, filter_R, min_feature_size, working_dir, beta = 1):<\/p>\n
geometry = TopologyOptimization3DLayered(params=params, eps_min=eps_bg, eps_max=eps_wg, x=x_pos, y=y_pos, z=z_pos, filter_R=filter_R, min_feature_size=min_feature_size, beta=beta)<\/p>\n
######## DEFINE FIGURE OF MERIT FOR EACH OUTPUT WAVEGUIDE ########<\/p>\n
fom1 = ModeMatch(monitor_name = ‘fom_1’, mode_number = ‘Fundamental TE mode’, direction = ‘Forward’, norm_p = 2, target_fom=0.25)<\/p>\n
fom2 = ModeMatch(monitor_name = ‘fom_2’, mode_number = ‘Fundamental TE mode’, direction = ‘Forward’, norm_p = 2, target_fom=0.25)<\/p>\n
fom3 = ModeMatch(monitor_name = ‘fom_3’, mode_number = ‘Fundamental TE mode’, direction = ‘Forward’, norm_p = 2, target_fom=0.25)<\/p>\n
fom4 = ModeMatch(monitor_name = ‘fom_4’, mode_number = ‘Fundamental TE mode’, direction = ‘Forward’, norm_p = 2, target_fom=0.25)<\/p>\n
######## DEFINE OPTIMIZATION ALGORITHM ########<\/p>\n
optimizer = ScipyOptimizers(max_iter=80, method=’L-BFGS-B’, scaling_factor=1, pgtol=1e-6, ftol=1e-5, scale_initial_gradient_to=0.25)<\/p>\n
######## DEFINE SETUP SCRIPT AND INDIVIDUAL OPTIMIZERS ########<\/p>\n
script = os.path.join(os.path.dirname(__file__), ‘1x4_test_material.fsp’)<\/p>\n
size_x = max(x_pos)-min(x_pos)<\/p>\n
script = script.replace(‘opt_size_x=10e-6′,’opt_size_x={:1.6g}’.format(size_x))<\/p>\n
size_y = max(y_pos)-min(y_pos)<\/p>\n
script = script.replace(‘opt_size_y=16e-6′,’opt_size_y={:1.6g}’.format(size_y))<\/p>\n
wavelengths = Wavelengths(start = 1550e-9, stop = 1550e-9, points = 1)<\/p>\n
opt1 = Optimization(base_script=script, wavelengths = wavelengths, fom=fom1, geometry=geometry, optimizer=optimizer, use_deps=False, hide_fdtd_cad=False, plot_history=False, store_all_simulations=False)<\/p>\n
opt2 = Optimization(base_script=script, wavelengths = wavelengths, fom=fom2, geometry=geometry, optimizer=optimizer, use_deps=False, hide_fdtd_cad=False, plot_history=False, store_all_simulations=False)<\/p>\n
opt3 = Optimization(base_script=script, wavelengths = wavelengths, fom=fom3, geometry=geometry, optimizer=optimizer, use_deps=False, hide_fdtd_cad=False, plot_history=False, store_all_simulations=False)<\/p>\n
opt4 = Optimization(base_script=script, wavelengths = wavelengths, fom=fom4, geometry=geometry, optimizer=optimizer, use_deps=False, hide_fdtd_cad=False, plot_history=False, store_all_simulations=False)<\/p>\n
######## PUT EVERYTHING TOGETHER AND RUN ########<\/p>\n
opt = opt1+opt2+opt3+opt4<\/p>\n
opt.continuation_max_iter = 20<\/p>\n
opt.run(working_dir = working_dir)<\/p>\n
\u5728\u8fed\u4ee3\u8fc7\u7a0b\u4e2d\u7684FOM\u5982\u4e0b\u56fe\u6240\u793a\uff1a<\/p>\n
<\/p>\n
\u60f3\u95ee\u4e00\u4e0b\u8001\u5e08\u4e2d\u95f4\u8fd9\u4e24\u4e2a\u7aef\u53e3\u529f\u7387\u8fd9\u4e48\u5927\u662f\u6b63\u5e38\u7684\u5417\uff1f
\u5f80\u4e0b\u7ee7\u7eed\u8fdb\u884c\u4e86\u51e0\u6b21\u8fed\u4ee3\uff0c\u53d1\u73b04\u4e2a\u7aef\u53e3\u7684\u529f\u7387\u53d8\u5316\u4e0d\u5927\uff0c\u5f97\u5230\u7684\u7ed3\u679c\u6ca1\u6709\u5411P=0.25\u65b9\u5411\u8d70\u3002<\/p>\n
\u5e0c\u671b\u8001\u5e08\u5e2e\u5fd9\u89e3\u7b54\u4e00\u4e0b\u7591\u60d1\uff0c\u975e\u5e38\u611f\u8c22\u3002<\/p>\n
<\/p>\n","protected":false},"template":"","class_list":["post-290237","topic","type-topic","status-closed","hentry"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t