Compute-energy-in-objects returns more than 100%
"Jack Holloway, MIT" <[email protected]> Sat, 25 Apr 2020 09:44:05 -0400
| Newsgroups | gmane.comp.science.photonic-bands |
|---|---|
| Message-ID | <[email protected]> |
--===============0105352211721366262==
Content-Type: multipart/signed;
boundary="Apple-Mail=_612C5818-E44A-4CCE-9988-3C42C04351C5";
protocol="application/pkcs7-signature";
micalg=sha-256
--Apple-Mail=_612C5818-E44A-4CCE-9988-3C42C04351C5
Content-Type: multipart/alternative;
boundary="Apple-Mail=_F1E0CC94-F145-4019-9342-66EBF355D45C"
--Apple-Mail=_F1E0CC94-F145-4019-9342-66EBF355D45C
Content-Transfer-Encoding: quoted-printable
Content-Type: text/plain;
charset=utf-8
Folks,
I just started using MPB to study various rectangular dielectric =
waveguide geometries. I=E2=80=99m interested in looking at the portion =
of the field energy in a rectangular core and a rectangular cladding. =
When I use compute-energy-in-objects, it curiously returns more than =
100%.=20
I=E2=80=99m using:
mpbi 1.7.0, Copyright (C) 1999-2012, MIT
Using libctl 4.2.0 and Guile 3.0.2.
On an Ubuntu installation.
When I=E2=80=99m looking at just the first mode:
(get-dfield 1) ; get d-field for first band
(compute-field-energy)
(print
"Gamma_guide: "
(compute-energy-in-objects (make block (size infinity w h) (center 0 0 =
0) (material air)))
"\n")
(print
"Gamma_clad: "
(compute-energy-in-objects
(make block (size infinity (+ w (* 2 tClad)) (+ h (* 2 tClad)))
(center 0 0 0) (material air)))
"\n")
(print
"Gamma_cell: "
(compute-energy-in-objects (make block (size infinity sc-z sc-y) =
(center 0 0 0) (material air)))
"\n")
I received the below output:
D-energy-components:, 102, 1, 0.0919373, 0.907724, 0.000338274
$35 =3D (0.9999999999997949 3.843388742030928e-9 0.09193731941226466 =
0.9077244029679912 9.523365114586835e-11 3.382689149641591e-4 =
4.766081590605888e-9)
Gamma_guide: 1.085048755617456
Gamma_clad: 1.0923722948863863
Gamma_cell: 1.0923723008598665
What am I missing?
-JW=
--Apple-Mail=_F1E0CC94-F145-4019-9342-66EBF355D45C
Content-Transfer-Encoding: quoted-printable
Content-Type: text/html;
charset=utf-8
<html><head><meta http-equiv=3D"Content-Type" content=3D"text/html; =
charset=3Dutf-8"></head><body style=3D"word-wrap: break-word; =
-webkit-nbsp-mode: space; line-break: after-white-space;" class=3D""><span=
style=3D"caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" =
class=3D"">Folks,</span><div class=3D"" style=3D"caret-color: rgb(0, 0, =
0); color: rgb(0, 0, 0);"><br class=3D""></div><div class=3D"" =
style=3D"caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">I just started =
using MPB to study various rectangular dielectric waveguide geometries. =
I=E2=80=99m interested in looking at the portion of the field energy in =
a rectangular core and a rectangular cladding. When I use =
compute-energy-in-objects, it curiously returns more than =
100%. </div><div class=3D"" style=3D"caret-color: rgb(0, 0, 0); =
color: rgb(0, 0, 0);"><br class=3D""></div><div class=3D"" =
style=3D"caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">I=E2=80=99m =
using:</div><div class=3D"" style=3D"caret-color: rgb(0, 0, 0); color: =
rgb(0, 0, 0);"> <span class=3D"" style=3D"color: rgb(16, 16, 16); =
font-family: Monaco; font-size: 10px; background-color: rgba(255, 255, =
255, 0.85);">mpbi 1.7.0, Copyright (C) 1999-2012, MIT</span><div =
class=3D"" style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: no-common-ligatures;">Using libctl =
4.2.0 and Guile 3.0.2.</span></div><div class=3D""><br =
class=3D"webkit-block-placeholder"></div><div class=3D""><div class=3D"" =
style=3D"font-variant-ligatures: normal; font-variant-east-asian: =
normal; font-variant-position: normal; line-height: normal; word-wrap: =
break-word; -webkit-nbsp-mode: space; line-break: =
after-white-space;"><div class=3D"">On an Ubuntu installation.</div><div =
class=3D""><br class=3D""></div><div class=3D"">When I=E2=80=99m looking =
at just the first mode:</div><div class=3D""><br class=3D""></div><div =
class=3D""><div class=3D"" style=3D"margin: 0px; font-stretch: normal; =
font-size: 10px; line-height: normal; font-family: Monaco; color: =
rgb(16, 16, 16); background-color: rgba(255, 255, 255, 0.85);"><span =
class=3D"" style=3D"font-variant-ligatures: =
no-common-ligatures;">(get-dfield 1) ; get d-field for first =
band</span></div><div class=3D"" style=3D"margin: 0px; font-stretch: =
normal; font-size: 10px; line-height: normal; font-family: Monaco; =
color: rgb(16, 16, 16); background-color: rgba(255, 255, 255, =
0.85);"><span class=3D"" style=3D"font-variant-ligatures: =
no-common-ligatures;">(compute-field-energy)</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;">(print</span></div><div class=3D"" style=3D"margin: =
0px; font-stretch: normal; font-size: 10px; line-height: normal; =
font-family: Monaco; color: rgb(16, 16, 16); background-color: rgba(255, =
255, 255, 0.85);"><span class=3D"" style=3D"font-variant-ligatures: =
no-common-ligatures;"> "Gamma_guide: "</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;"> (compute-energy-in-objects (make block (size =
infinity w h) (center 0 0 0) (material air)))</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;"> "\n")</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;">(print</span></div><div class=3D"" style=3D"margin: =
0px; font-stretch: normal; font-size: 10px; line-height: normal; =
font-family: Monaco; color: rgb(16, 16, 16); background-color: rgba(255, =
255, 255, 0.85);"><span class=3D"" style=3D"font-variant-ligatures: =
no-common-ligatures;"> "Gamma_clad: "</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;"> (compute-energy-in-objects</span></div><div =
class=3D"" style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: no-common-ligatures;"> =
(make block (size infinity (+ w (* 2 tClad)) (+ h (* 2 =
tClad)))</span></div><div class=3D"" style=3D"margin: 0px; font-stretch: =
normal; font-size: 10px; line-height: normal; font-family: Monaco; =
color: rgb(16, 16, 16); background-color: rgba(255, 255, 255, =
0.85);"><span class=3D"" style=3D"font-variant-ligatures: =
no-common-ligatures;"> =
(center 0 0 0) (material air)))</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;"> "\n")</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;">(print</span></div><div class=3D"" style=3D"margin: =
0px; font-stretch: normal; font-size: 10px; line-height: normal; =
font-family: Monaco; color: rgb(16, 16, 16); background-color: rgba(255, =
255, 255, 0.85);"><span class=3D"" style=3D"font-variant-ligatures: =
no-common-ligatures;"> "Gamma_cell: "</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;"> (compute-energy-in-objects (make block (size =
infinity sc-z sc-y) (center 0 0 0) (material air)))</span></div><div =
class=3D"" style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;"> "\n")</span></div></div><div class=3D""><br =
class=3D""></div><div class=3D"">I received the below output:</div><div =
class=3D""><br class=3D""></div><div class=3D""><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: =
no-common-ligatures;">D-energy-components:, 102, 1, 0.0919373, 0.907724, =
0.000338274</span></div><div class=3D"" style=3D"margin: 0px; =
font-stretch: normal; font-size: 10px; line-height: normal; font-family: =
Monaco; color: rgb(16, 16, 16); background-color: rgba(255, 255, 255, =
0.85);"><span class=3D"" style=3D"font-variant-ligatures: =
no-common-ligatures;">$35 =3D (0.9999999999997949 3.843388742030928e-9 =
0.09193731941226466 0.9077244029679912 9.523365114586835e-11 =
3.382689149641591e-4 4.766081590605888e-9)</span></div><div class=3D"" =
style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: no-common-ligatures;">Gamma_guide: =
1.085048755617456</span></div><div class=3D"" style=3D"margin: 0px; =
font-stretch: normal; font-size: 10px; line-height: normal; font-family: =
Monaco; color: rgb(16, 16, 16); background-color: rgba(255, 255, 255, =
0.85);"><span class=3D"" style=3D"font-variant-ligatures: =
no-common-ligatures;">Gamma_clad: 1.0923722948863863</span></div><div =
class=3D"" style=3D"margin: 0px; font-stretch: normal; font-size: 10px; =
line-height: normal; font-family: Monaco; color: rgb(16, 16, 16); =
background-color: rgba(255, 255, 255, 0.85);"><span class=3D"" =
style=3D"font-variant-ligatures: no-common-ligatures;">Gamma_cell: =
1.0923723008598665</span></div></div></div></div></div><div class=3D"" =
style=3D"caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><br =
class=3D"Apple-interchange-newline"></div><div class=3D"" =
style=3D"caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">What am I =
missing?</div><div class=3D"" style=3D"caret-color: rgb(0, 0, 0); color: =
rgb(0, 0, 0);"><br class=3D""></div><div class=3D"" style=3D"caret-color: =
rgb(0, 0, 0); color: rgb(0, 0, 0);">-JW</div></body></html>=
--Apple-Mail=_F1E0CC94-F145-4019-9342-66EBF355D45C--
--Apple-Mail=_612C5818-E44A-4CCE-9988-3C42C04351C5
Content-Disposition: attachment;
filename=smime.p7s
Content-Type: application/pkcs7-signature;
name=smime.p7s
Content-Transfer-Encoding: base64
MIAGCSqGSIb3DQEHAqCAMIACAQExDzANBglghkgBZQMEAgEFADCABgkqhkiG9w0BBwEAAKCCA8Qw
ggPAMIIDKaADAgECAhEA27Z1/WjjOHpmBWXD2rMGHjANBgkqhkiG9w0BAQsFADBsMQswCQYDVQQG
EwJVUzEWMBQGA1UECBMNTWFzc2FjaHVzZXR0czEuMCwGA1UEChMlTWFzc2FjaHVzZXR0cyBJbnN0
aXR1dGUgb2YgVGVjaG5vbG9neTEVMBMGA1UECxMMQ2xpZW50IENBIHYxMB4XDTE5MDgzMDIwMTkx
NFoXDTIwMDczMDIwMTkxNFowgacxCzAJBgNVBAYTAlVTMRYwFAYDVQQIEw1NYXNzYWNodXNldHRz
MS4wLAYDVQQKEyVNYXNzYWNodXNldHRzIEluc3RpdHV0ZSBvZiBUZWNobm9sb2d5MRUwEwYDVQQL
EwxDbGllbnQgQ0EgdjExGDAWBgNVBAMTD0phY2sgVyBIb2xsb3dheTEfMB0GCSqGSIb3DQEJARYQ
aG9sbG93YXlATUlULkVEVTCCASIwDQYJKoZIhvcNAQEBBQADggEPADCCAQoCggEBAMPLnm4aEmCZ
10AWBw3UC8qHqBo0ocTB2bUohT+O6Rvzr7GS0WrT/UmJhSeRe4238KWdjchzXkYUEDUkEed4oz2a
9Fkg64Kr5UXkTzpAlxpcRNhF/BYpIFuNPs+/jJc/3TpjDbMurkI2hq1fNR++1VxxZ7zNkvU4Qhzl
XJlohDDasVes/FJklpKZIyNl8AOXjkqKZz8OGlhGG5M81If0ly8WHAO5kqq5Y3thY2IJ4qtE1mns
qFcAE+cBABZvX4wnwn8PrkhEk/DZ5lzv6udljSXSrlVITTPElIc8vggUcjbt3fJD5NnBrqgb1lmF
wICdpy0/936BoWpFZFIHGeEWYHECAwEAAaOBoTCBnjAJBgNVHRMEAjAAMBEGCWCGSAGG+EIBAQQE
AwIFoDAdBgNVHSUEFjAUBggrBgEFBQcDBAYIKwYBBQUHAwIwCwYDVR0PBAQDAgXgMB0GA1UdDgQW
BBQEmw95hOg6t8wJcQ/7QdSq/YcygzAzBgNVHR8ELDAqMCigJqAkhiJodHRwOi8vY2EubWl0LmVk
dS9jYS9taXRjbGllbnQuY3JsMA0GCSqGSIb3DQEBCwUAA4GBAAGnr5+jfs9sfzfYWu+q+JRLteWC
F91IytY/tB0bhpWH2b48qb6H3HX/NAzEUlpTImE69IrYn+vzpTm0X6T+VQKfrZbS7rptRRtvXaQG
6hpzi49m+2VoPBz2z/bwtXVc3+FTsBjY2h3hLQ2JH3ks/qifNM1uVX/zEtMmuiaMXZHCMYIDRjCC
A0ICAQEwgYEwbDELMAkGA1UEBhMCVVMxFjAUBgNVBAgTDU1hc3NhY2h1c2V0dHMxLjAsBgNVBAoT
JU1hc3NhY2h1c2V0dHMgSW5zdGl0dXRlIG9mIFRlY2hub2xvZ3kxFTATBgNVBAsTDENsaWVudCBD
QSB2MQIRANu2df1o4zh6ZgVlw9qzBh4wDQYJYIZIAWUDBAIBBQCgggGVMBgGCSqGSIb3DQEJAzEL
BgkqhkiG9w0BBwEwHAYJKoZIhvcNAQkFMQ8XDTIwMDQyNTEzNDQwNVowLwYJKoZIhvcNAQkEMSIE
ID+VmG1Z91sGqNYVUbyi9y0BPrgcudnf+1O0+HsE4P67MIGSBgkrBgEEAYI3EAQxgYQwgYEwbDEL
MAkGA1UEBhMCVVMxFjAUBgNVBAgTDU1hc3NhY2h1c2V0dHMxLjAsBgNVBAoTJU1hc3NhY2h1c2V0
dHMgSW5zdGl0dXRlIG9mIFRlY2hub2xvZ3kxFTATBgNVBAsTDENsaWVudCBDQSB2MQIRANu2df1o
4zh6ZgVlw9qzBh4wgZQGCyqGSIb3DQEJEAILMYGEoIGBMGwxCzAJBgNVBAYTAlVTMRYwFAYDVQQI
Ew1NYXNzYWNodXNldHRzMS4wLAYDVQQKEyVNYXNzYWNodXNldHRzIEluc3RpdHV0ZSBvZiBUZWNo
bm9sb2d5MRUwEwYDVQQLEwxDbGllbnQgQ0EgdjECEQDbtnX9aOM4emYFZcPaswYeMA0GCSqGSIb3
DQEBAQUABIIBAJ8h7BDY7bDyvFvw27RCpxWBEonY9biaTqek6HXn6SvmTJ2ljUUoT+b8OSXW8dWR
bNE6PRLxxFQE2WxVhFJ2WjOoiD/zsLBYKbmfXo/Zils8ZKXmPHh47DyRUXY9TppcOXpIWvzrbT7a
QV5TFMMb5z/S66+spY2fdN6/W7owQlhIXuI1vFZT/5sH5qUQgYrbik9x9fHHglHR5f8UhgWslAOT
jQhp3DLYjdhp3T+Q5VVCWXDyQBE+yb9gZvJRJHjuSq1zuKE5a8dwQq8luJCGZMpyI3UJ8QXZKv5r
uboqf7FhoSTOh5xTsH/DYWsLyIdktp1+Y2BT4dHYu+F0l43vw20AAAAAAAA=
--Apple-Mail=_612C5818-E44A-4CCE-9988-3C42C04351C5--
--===============0105352211721366262==
Content-Type: text/plain; charset="utf-8"
MIME-Version: 1.0
Content-Transfer-Encoding: base64
Content-Disposition: inline
X19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX18KbXBiLWRpc2N1
c3MgbWFpbGluZyBsaXN0Cm1wYi1kaXNjdXNzQGFiLWluaXRpby5taXQuZWR1Cmh0dHA6Ly9hYi1p
bml0aW8ubWl0LmVkdS9jZ2ktYmluL21haWxtYW4vbGlzdGluZm8vbXBiLWRpc2N1c3M=
--===============0105352211721366262==--