diff --git a/ElectromagneticScattering/scattererTmatrix.pro b/ElectromagneticScattering/scattererTmatrix.pro index d288e8294aec8931aa6a6c06422bba06ab2b4d40..3e62be26e38ebfd3abe401a19bf6fc824de60373 100644 --- a/ElectromagneticScattering/scattererTmatrix.pro +++ b/ElectromagneticScattering/scattererTmatrix.pro @@ -204,8 +204,8 @@ PostProcessing { { Name Mnm_source~{pe} ; Value { Local { [ Mnm_source~{pe}[] ]; In All_domains; Jacobian JVol; } } } For po In {1:p_max} - { Name feM~{pe}~{po} ; Value { Integral { [ (1/r_pml_in)^2 *normalize_fenm_Z~{po}[]* {u}*Conj[Znm~{po}[]] ]; In SurfInt; Integration Int_1 ; Jacobian JSur; } } } - { Name fhM~{pe}~{po} ; Value { Integral { [ (1/r_pml_in)^2 *normalize_fhnm_X~{po}[]* {u}*Conj[Xnm~{po}[]] ]; In SurfInt; Integration Int_1 ; Jacobian JSur; } } } + { Name feM~{pe}~{po} ; Value { Integral { [ 1/(r_pml_in^2*epsilonr_Out[]) * normalize_fenm_Z~{po}[]* {u}*Conj[Znm~{po}[]] ]; In SurfInt; Integration Int_1 ; Jacobian JSur; } } } + { Name fhM~{pe}~{po} ; Value { Integral { [ 1/(r_pml_in^2*epsilonr_Out[]) * normalize_fhnm_X~{po}[]* {u}*Conj[Xnm~{po}[]] ]; In SurfInt; Integration Int_1 ; Jacobian JSur; } } } EndFor } }