From 3b9a6f54d08a4ef8f46b414a3cde4dfc3166014a Mon Sep 17 00:00:00 2001
From: Guillaume Demesy <guillaume.demesy@fresnel.fr>
Date: Mon, 29 Jul 2024 21:18:46 +0200
Subject: [PATCH] new test case updates
---
ElectromagneticScattering/scattererTmatrix.pro | 2 ++
ElectromagneticScattering/scattererTmatrix_data.geo | 2 +-
2 files changed, 3 insertions(+), 1 deletion(-)
diff --git a/ElectromagneticScattering/scattererTmatrix.pro b/ElectromagneticScattering/scattererTmatrix.pro
index 3e62be2..8d570c1 100644
--- a/ElectromagneticScattering/scattererTmatrix.pro
+++ b/ElectromagneticScattering/scattererTmatrix.pro
@@ -199,6 +199,7 @@ PostProcessing {
For pe In {1:p_max}
{ Name VPWM_postpro~{pe}; NameOfFormulation VPWMN_helmholtz_vector; NameOfSystem T;
Quantity {
+ // { Name epsr_xx; Value { Local { [ CompXX[epsilonr[]] ]; In All_domains; Jacobian JVol; } } }
{ Name E_scat ; Value { Local { [{u}]; In All_domains; Jacobian JVol; } } }
{ Name H_scat ; Value { Local { [siwt*I[]/(mur[]*mu0*omega0)*{Curl u}]; In All_domains; Jacobian JVol; } } }
{ Name Mnm_source~{pe} ; Value { Local { [ Mnm_source~{pe}[] ]; In All_domains; Jacobian JVol; } } }
@@ -227,6 +228,7 @@ PostOperation {
For pe In {1:p_max}
{Name VPWM_postop~{pe}; NameOfPostProcessing VPWM_postpro~{pe} ;
Operation {
+ // Print [ epsr_xx , OnElementsOf All_domains, File "epsr_xx.pos"];
If (flag_plotcuts==1)
Print [ E_scat , OnGrid
{(r_pml_in-1*nm)*Sin[$B]*Cos[$C], (r_pml_in-1*nm)*Sin[$B]*Sin[$C], (r_pml_in-1*nm)*Cos[$B]}
diff --git a/ElectromagneticScattering/scattererTmatrix_data.geo b/ElectromagneticScattering/scattererTmatrix_data.geo
index 159bb89..b601341 100644
--- a/ElectromagneticScattering/scattererTmatrix_data.geo
+++ b/ElectromagneticScattering/scattererTmatrix_data.geo
@@ -132,7 +132,7 @@ DefineConstant[
space2pml = {lambda_bg/2 , Name StrCat[pp3 , "0space around scatterer [nm]"] , Highlight Str[colorppNO] , Closed 1,Min (lambda_bg/10.), Max (3*lambda_bg)},
pml_size = {lambda_bg , Name StrCat[pp3 , "1PML thickness [nm]"] , Highlight Str[colorppWA] , Closed 0, Min (lambda_bg/10), Max (3*lambda_bg)},
paramaille = {5. , Name StrCat[pp3 , "2mesh size"], Highlight Str[colorppWA] , Closed 0},
- refine_scat = {1. , Name StrCat[pp3 , "3scatterer mesh refinement"], Highlight Str[colorppWA] , Closed 0}
+ refine_scat = {2. , Name StrCat[pp3 , "3scatterer mesh refinement"], Highlight Str[colorppWA] , Closed 0}
is_FEM_o2 = {1 , Name StrCat[pp3 , "4Interpolation order "] , Choices {0="order 1",1="order 2"}, Closed 0}
];
--
GitLab