Group { front = Region[27] ; sides = Region[25] ; back = Region[26] ; air = Region[30] ; Domain3D = Region[{air}]; Domain2D = Region[{front}]; Domain = Region[{Domain3D,Domain2D}] ; } Function { nu0 = 1 / (4.e-7 * Pi) ; ep0 = 8.854187817e-12 ; epsilon [Domain3D] = ep0 ; nu [Domain3D] = nu0 ; NxH[front ] = Vector[0,0,Sin[Pi*Y[]/50]] ;//w=50; m=1; kt = Pi*m/W; Freq = 4.e9 ; } Constraint { { Name e ;Type Assign ; Case { { Region sides ; Value 0.0 ; } } } } Jacobian { { Name JVol ; Case { { Region All ; Jacobian Vol ; } } } { Name JSur ; Case { { Region All ; Jacobian Sur ; } } } } Integration { { Name I1 ; Case { { Type Gauss ; Case { { GeoElement Point ; NumberOfPoints 1 ; } { GeoElement Line ; NumberOfPoints 3 ; } { GeoElement Triangle ; NumberOfPoints 4 ; } { GeoElement Quadrangle ; NumberOfPoints 4 ; } { GeoElement Tetrahedron ; NumberOfPoints 4 ; } { GeoElement Hexahedron ; NumberOfPoints 6 ; } { GeoElement Prism ; NumberOfPoints 9 ; } } } } } } /* -------------------------------------------------------------------------- MW_e_3D -------------------------------------------------------------------------- */ FunctionSpace { { Name Hcurl_e_3D; Type Form1; BasisFunction { { Name se; NameOfCoef ee; Function BF_Edge; Support Domain; Entity EdgesOf[All]; } } Constraint { { NameOfCoef ee; EntityType EdgesOf ; NameOfConstraint e; } // { NameOfCoef ee ;EntityType EdgesOfTreeIn ; EntitySubType StartingOn ; NameOfConstraint Gauge; } } } } /*FunctionSpace { // 1st order nodal elements { Name form0_1st; Type Form0 ; BasisFunction { { Name se ; NameOfCoef ae ; Function BF_Node; Support Domain3D ; Entity NodesOf[ All ] ; } } Constraint { { NameOfCoef ae ; EntityType NodesOf ; NameOfConstraint e ; } } } // 2nd order nodal elements { Name form0_2nd ; Type Form0 ; BasisFunction { { Name s1 ; NameOfCoef v1 ; Function BF_Node; Support Domain ; Entity NodesOf[ All ] ; } { Name s2; NameOfCoef v2; Function BF_Node_2E; Support Domain; Entity EdgesOf[All]; } } Constraint { { NameOfCoef v1 ; EntityType NodesOf ; NameOfConstraint e ; } { NameOfCoef v2 ; EntityType EdgesOf ; NameOfConstraint e ; } } }*/ Formulation { { Name MW_e_3D; Type FemEquation; Quantity { { Name e; Type Local; NameOfSpace Hcurl_e_3D; } } Equation { Galerkin { DtDt [ epsilon[] * Dof{e} , {e} ]; In Domain3D; Integration I1; Jacobian JVol; } Galerkin { [ nu[] * Dof{d e} , {d e} ]; In Domain3D; Integration I1; Jacobian JVol; } Galerkin { Dt [ - NxH[] , {e} ]; In Domain2D; Integration I1; Jacobian JSur; } } } } Resolution { { Name MW_e_3D; System { { Name A; NameOfFormulation MW_e_3D; Type ComplexValue; Frequency Freq; } } Operation { Generate [A]; Solve [A]; SaveSolution [A]; } } } PostProcessing { { Name MW_e_3D; NameOfFormulation MW_e_3D; NameOfSystem A; Quantity { { Name e; Value{ Local{ [ {e} ] ; In Domain3D;Jacobian JVol; } } } { Name e2; Value{ Local{ [ Norm[{e}] ] ; In Domain3D;Jacobian JVol; } } } } } } PostOperation { { Name e ; NameOfPostProcessing MW_e_3D ; Operation { Print[ e, OnElementsOf Domain3D , File "e.pos" ] ; Print[ e2, OnElementsOf Domain3D , File "e2.pos" ] ; } } }