{"id":294874,"date":"2023-07-26T09:51:35","date_gmt":"2023-07-26T09:51:35","guid":{"rendered":"\/forum\/forums\/topic\/integrate-integrate2-command-not-working-as-intended-not-integrating\/"},"modified":"2023-07-26T09:51:35","modified_gmt":"2023-07-26T09:51:35","slug":"integrate-integrate2-command-not-working-as-intended-not-integrating","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/integrate-integrate2-command-not-working-as-intended-not-integrating\/","title":{"rendered":"Integrate\/Integrate2 command not working as intended\/not integrating"},"content":{"rendered":"

Hello, I am trying to calculate the field confinement factor of my SOI design with the formula in the image below (skipped the denominator while calculating, also know that this formulsa is not suitable for SOI but just trying to get an approx for anyone concerned):<\/p>\n

\"\"<\/p>\n

While trying to calculate the integration, I have used the integrate2 command to take the integral of the Ey_string which is equal tothe mode1 Ey:<\/p>\n

simname = “mode1”;<\/p>\n

Ey_string = getdata(simname, ‘Ey’); #and the usage of the integrate2 function was as follows:<\/p>\n

y_vector = [(-0.1925 * 10^-6):(-0.1025* 10^-6)]; # the y limits of the slot, 90nm<\/p>\n

EY_integ = integrate2(Ey_string,1,y_vector);<\/p>\n

After doing these steps, my expectation was to get a scalar number to see the confinement factor in a percentile fashion. Checking the Ey_string it is the same as the Ey of the mode1 itself of FDE, after that when I check the EY_integ, it is the same as Ey_string. Nothing changed, produced a matrix instead of a scalar number.<\/p>\n

I would appreciate any help that would broaden my understanding of LSF and the integrate command. Sorry for not using the code block for the code parts in the post it’s my first time here, couldn’t make it work somehow. The full script:<\/p>\n

simname = “mode1”;<\/p>\n

Ez_string = getdata(simname, ‘Ez’);
Ez_string = real(Ez_string * Ez_string);<\/p>\n

Ey_string = getdata(simname, ‘Ey’);
Ey_string = real(Ey_string * Ey_string);<\/p>\n

y_vector = [(-0.1925 * 10^-6):(-0.1025* 10^-6)];
z_vector = [(1.03 * 10^-6): (1.25* 10^-6)];<\/p>\n

EY_integ = integrate2(Ey_string,1,y_vector);<\/p>\n

EZ_integ = integrate2(Ez_string,1,z_vector);<\/p>\n

result_e = struct;
result_e.Name = “E”;
E_sum = EY_integ + EZ_integ;<\/p>\n

print(E_sum);<\/p>\n","protected":false},"template":"","class_list":["post-294874","topic","type-topic","status-closed","hentry","topic-tag-FDE-1","topic-tag-lsf","topic-tag-lumerical","topic-tag-photonics","topic-tag-script"],"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