Dear All,
Can you check attached simple .ctl file and give me advice?
Problem: Plane wave pulse Ez (wave in Y direction) in _empty_ 3D space
(periodic boundaries, PML on Y dir).
When I set _BOTH_ symmetries:
(set! symmetries (list
(make mirror-sym (direction X) (phase 1))
(make mirror-sym (direction Z) (phase -1))
))
the fields go to infinity. BUT when I comment out one symmetry or all then
the results seems to be OK.
Thanks in advance,
Serge
begin 666 sph-3D-symXZtest.ctl
M.R!2969L96-T:6]N+51R86YS;6ES<VEO;B!F;W(@<W1R=6-T=7)E(&EN(#-D
M+@H[(&$],3!N;3TQ92TX;0H[(&=E;VUE=')Y('!A<F%M971E<G,Z"BAD969I
M;F4M<&%R86T@1%]S<&AE<F4@(#$I(" @(" ).R!D:6%M971E<B!O9B!T:&4@
M<W!H97)E"BAD969I;F4M<&%R86T@9'!M;" @,BD)"3L@=&AI8VMN97-S(&]F
M('!M;"!L87EE<G,@:6X@62!D:7)E8W1I;VX**&1E9FEN92UP87)A;2!097)I
M;V0@(#,@*0D).R!P97)I;V0@;V8@<W!H97)E<R!I;B!U;FET<R!O9B!$"CL@
M9&]M86EN.@HH9&5F:6YE+7!A<F%M('-X("@J(%!E<FEO9"!$7W-P:&5R92DI
M"3L@<VEZ92!I;B!8(&1I<F5C=&EO;@HH9&5F:6YE+7!A<F%M('-Z("@J(%!E
M<FEO9"!$7W-P:&5R92DI"3L@<VEZ92!I;B!:(&1I<F5C=&EO;@HH9&5F:6YE
M+7!A<F%M('-Y("@K("@J(#,@1%]S<&AE<F4I(&1P;6P@9'!M;"DI.R!S:7IE
M(&EN(%D@9&ER96-T:6]N( H[('-O=7)C92!P87)A;65T97)S.@HH9&5F:6YE
M+7!A<F%M(&9C96X@,"XR,"D)"3L@<'5L<V4@8V5N=&5R(&9R97%U96YC>2 *
M*&1E9FEN92UP87)A;2!D9B @(# N,3 I( D).R!P=6QS92!W:61T:" H:6X@
M9G)E<75E;F-Y*2 @"BAD969I;F4M<&%R86T@;F9R97$@,C I"0D[(&YU;6)E
M<B!O9B!F<F5Q=65N8VEE<R!A="!W:&EC:"!T;R!C;VUP=71E(&9L=7@@<W1E
M< H[(&UA=&5R:6%L('!A<F%M971E<G,Z"BAD969I;F4M<&%R86T@3C$@,BXP
M*0D)"3L@;65D,2 ]($XQ*DXQ(&5M8BYM961I=6T@?B!:;D\**'-E="UP87)A
M;2$@<F5S;VQU=&EO;B R,"D)"3L@/CT@."IS<7)T*$5P<U]M87@I*F9M87@*
M*&1E9FEN92UP87)A;2!R871I;R Q*0D)"3L@>2\H>"QZ*2!R871I;R!O9B!E
M;&QI<'-O:60@<&%R=&EC;&4@( H[(&QO9VEC86P@<&%R86UE=&5R<SH**&1E
M9FEN92UP87)A;2!E;7!T>3\@("!T<G5E*2 )"3L@:68@=')U92P@<V]U<F-E
M(&-A<FEB<F%T:6]N(')U;@H[(&-A;&,N(&5L;&EP<V]I9"!S:7IE.@HH9&5F
M:6YE('-P:"UY("@J($1?<W!H97)E("@D97AP="!R871I;R H+R R(#,I*2DI
M"BAD969I;F4@<W!H+7@@*"\@<W!H+7D@<F%T:6\I*0HH9&5F:6YE('-P:"UZ
M('-P:"UX*0HH<')I;G0@(G@](B!S<&@M>" B('D](B!S<&@M>2 B('H](B!S
M<&@M>B B('9O;#TB("@J('-P:"UX('-P:"UY('-P:"UZ*2 B/3TB("@D97AP
M="!$7W-P:&5R92 S*2 B/UQN(BD*.R!D969I;F4@;6%T97)I86QS.@HH9&5F
M:6YE(&UE9#$@"2AM86ME(&UE9&EU;2 H97!S:6QO;B H*B!.,2!.,2DI*2D*
M"CL@3&%T=&EC92 F($)O=6YD87)Y(&-O;F1I=&EO;G,Z"BAS970A(&=E;VUE
M=')Y+6QA='1I8V4@*&UA:V4@;&%T=&EC92 H<VEZ92!S>"!S>2!S>BDI*0D[
M(#-$(&-A<V4**'-E="$@:RUP;VEN=" H=F5C=&]R,R P(# @,"DI"BAS970A
M(&5N<W5R92UP97)I;V1I8VET>2!T<G5E*0HH<V5T(2!P;6PM;&%Y97)S("AL
M:7-T( H)"2AM86ME('!M;" H=&AI8VMN97-S(&1P;6PI(" H9&ER96-T:6]N
M(%DI("AS:61E($AI9V@I*0H)"2AM86ME('!M;" H=&AI8VMN97-S(&1P;6PI
M(" H9&ER96-T:6]N(%DI("AS:61E($QO=R I*0HI*0H[($=E;VUE=')Y.@HH
M<V5T(2!G96]M971R>0H)*&EF(&5M<'1Y/PH)"2AL:7-T"@D)"2AM86ME(&)L
M;V-K("AC96YT97(@," P(# I("AS:7IE(&EN9FEN:71Y(&EN9FEN:71Y(&EN
M9FEN:71Y*2 H;6%T97)I86P@;65D,2DI*0D[(&UE9#$@9FEL;"!A;&P*"0DH
M;&ES= H)"0DH;6%K92!B;&]C:R H8V5N=&5R(# @," P*2 H<VEZ92!I;F9I
M;FET>2!I;F9I;FET>2!I;F9I;FET>2D@*&UA=&5R:6%L(&UE9#$I*0D[(&UE
M9#$@9FEL;"!A;&P*"0D)*&UA:V4@96QL:7!S;VED("AC96YT97(@," P(# I
M(" @(" H<VEZ92!S<&@M>"!S<&@M>2!S<&@M>BD@*&UA=&5R:6%L(&UE=&%L
M*2D)"3L@<W!H97)E('!E<F9E8W0M;65T86P@"@D)*0H)*0HI"CL@4V]U<F-E
M(&EN(&UE9#$Z"BAS970A('-O=7)C97,@*&QI<W0*(" @(" @(" @(" @(" @
M*&UA:V4@<V]U<F-E"CL@(" @(" @(" @(" @(" @("AS<F,@*&UA:V4@8V]N
M=&EN=6]U<RUS<F,@*&9R97%U96YC>2!F8V5N*2 I*0H@(" @(" @(" @(" @
M(" @("AS<F,@*&UA:V4@9V%U<W-I86XM<W)C("AF<F5Q=65N8WD@9F-E;BD@
M*&9W:61T:"!D9BDI*0H@(" @(" @(" @(" @(" @("AC;VUP;VYE;G0@17HI
M("AC96YT97(@," H*R H+R!S>2 M,BD@9'!M;" H+R R(')E<V]L=71I;VXI
M*2 P*2 H<VEZ92!S>" P('-Z*2DI*0H**'-E="$@<WEM;65T<FEE<R H;&ES
M=" )"B @(" @(" H;6%K92!M:7)R;W(M<WEM("AD:7)E8W1I;VX@6"D@*'!H
M87-E(#$I*0H@(" @(" @*&UA:V4@;6ER<F]R+7-Y;2 H9&ER96-T:6]N(%HI
M("AP:&%S92 M,2DI"BDI"@H[(&9L=7AE<R F(&9L=7@@<F5G:6]N<PHH9&5F
M:6YE(')[email protected]!R969L96-T960@9FQU> H)*&%D9"UF;'5X(&9C96X@9&8@
M;F9R97$*"0DH;6%K92!F;'5X+7)E9VEO;B H8V5N=&5R(# @*"L@*"\@<WD@
M+3(I(&1P;6P@*"\@-"!R97-O;'5T:6]N*2D@,"D@*'-I>F4@<W@@,"!S>BD@
M*&1I<F5C=&EO;B!9*2DI*0HH9&5F:6YE('1R86YS(#L@=')A;G-M:71T960@
M9FQU> H)*&%D9"UF;'5X(&9C96X@9&8@;F9R97$*"0DH;6%K92!F;'5X+7)E
M9VEO;B H8V5N=&5R(# @*"T@*"\@<WD@(#(I(&1P;6P@*"\@-"!R97-O;'5T
M:6]N*2D@,"D@*'-I>F4@<W@@,"!S>BD@*&1I<F5C=&EO;B!9*2DI*0H*"CL@
M+2TM+2!R=6X@+2TM+2TM+2TM+2TM+2TM+2TM+2TM+2TM+2TM+2TM+2TM+2TM
M+2TM+2TM+2TM+2TM+2TM+2TM+2TM+2TM+2TM+2TM+0HH<V5T(2!P<F]G<F5S
M<RUI;G1E<G9A;" Q,"D)"3L@;W5T<'5T('!R;V=R97-S)R!M<V=S(&EN=&5R
M=F%L"BAU<V4M;W5T<'5T+61I<F5C=&]R>2D*"BAI9B H;F]T(&5M<'1Y/RD@
M(" @(" @(#L@9F]R(&YO<FUA;"!R=6XL(&QO860@;F5G871E9"!F:65L9',@
M=&\@<W5B=')A8W0@:6YC:61E;G0@9G)O;2!R969L+B!F:65L9',*"2AL:7-T
M( H)"2AL;V%D+6UI;G5S+69L=7@@(G)E9FPM9FQU>"(@<F5F;"DI*0H**')U
M;BUS;W5R8V5S*R *"2AS=&]P+7=H96XM9FEE;&1S+61E8V%Y960@,3 @17H@
M"@D)*'9E8W1O<C,@," H+R!S<&@M>2 M,2XY*2 P*2 Q92TS*0H)"2AA="UB
M96=I;FYI;F<@;W5T<'5T+65P<VEL;VXI"@D)*'1O+6%P<&5N9&5D(")E>BUS
M;&EC92(@"@D)"2AA="UE=F5R>2 Q("AI;BUV;VQU;64@*'9O;'5M92 H8V5N
M=&5R(# @," P*2 H<VEZ92!S>"!S>2 P*2D@;W5T<'5T+65F:65L9"UZ*2DI
M"BD*"BAI9B!E;7!T>3\@"0D).R!F<F]M(&YO<FUA;&EZ871I;VX@<G5N+"!S
M879E(&9L=7@@9FEE;&1S(&9O<B!R969L+B!P;&%N90H)*'-A=F4M9FQU>" B
M<F5F;"UF;'5X(B!R969L*2D*"0H[('!R:6YT(&9L=7AE<SH**&1I<W!L87DM
H9FQU>&5S(')E9FP@=')A;G,@*0D)"3L@<')I;G0@9FQU>&5S"@H*"@``
`
end
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
