// simnibs automatically generated file for getdp simulation // file created with getdp.py on 20140415_120207 // Andre Antunes WM = 1; WM_surface = 1001; GM = 2; GM_surface = 1002; CSF = 3; CSF_surface = 1003; Bone = 4; Bone_surface = 1004; Scalp = 5; Scalp_surface = 1005; fn_dadt = "almi5_surface.msh"; Group { WM = Region[WM]; WM_surface = Region[WM_surface]; GM = Region[GM]; GM_surface = Region[GM_surface]; CSF = Region[CSF]; CSF_surface = Region[CSF_surface]; Bone = Region[Bone]; Bone_surface = Region[Bone_surface]; Scalp = Region[Scalp]; Scalp_surface = Region[Scalp_surface]; Omega = Region[{WM,GM,CSF,Bone,Scalp}]; Gama = Region[{WM_surface, GM_surface, CSF_surface, Bone_surface, Scalp_surface}]; } Function { dadt[] = VectorField[XYZ[]]{0}; sigma[WM] = 0.126; sigma[WM_surface] = 0.126; sigma[GM] = 0.275; sigma[GM_surface] = 0.275; sigma[CSF] = 1.654; sigma[CSF_surface] = 1.654; sigma[Bone] = 0.01; sigma[Bone_surface] = 0.01; sigma[Scalp] = 0.465; sigma[Scalp_surface] = 0.465; } Jacobian { { Name Volume; Case {{ Region Omega; Jacobian Vol; }} } { Name Surface; Case {{ Region Gama; Jacobian Sur; }} } } Integration {{ Name GradGrad; Case {{ Type Gauss; Case { {GeoElement Triangle; NumberOfPoints 4;} {GeoElement Tetrahedron; NumberOfPoints 1;} } }} }} FunctionSpace {{ Name Hgrad_vf_Ele; Type Form0; BasisFunction {{ Name sn; NameOfCoef vn; Function BF_Node; Support Region[{Omega,Gama}]; Entity NodesOf[All]; }} }} Formulation {{ Name QS_Formulation; Type FemEquation; Quantity {{ Name v; Type Local; NameOfSpace Hgrad_vf_Ele; }} Equation { Galerkin { [sigma[] * Dof{Grad v}, {Grad v}]; In Omega; Jacobian Volume; Integration GradGrad; } Galerkin { [sigma[] * dadt[], {Grad v}]; In Omega; Jacobian Volume; Integration GradGrad; } Galerkin { [sigma[] * Dof{Grad v}, {Grad v}]; In Gama; Jacobian Surface; Integration GradGrad; } Galerkin { [sigma[] * dadt[], {Grad v}]; In Gama; Jacobian Surface; Integration GradGrad; } } }} Resolution {{ Name QS; System {{ Name QS_dAdt; NameOfFormulation QS_Formulation; }} Operation { GmshRead[fn_dadt]; Generate QS_dAdt; Solve QS_dAdt; SaveSolution QS_dAdt; } }} PostProcessing {{ Name QS_PostProcessing; NameOfFormulation QS_Formulation; NameOfSystem QS_dAdt; Quantity { {Name v; Value {Local { [ {v} ]; In Omega; Jacobian Volume;}}} {Name normE; Value {Local { [ Norm[-{Grad v}-dadt[]] ]; In Omega; Jacobian Volume;}}} } }} PostOperation V UsingPost QS_PostProcessing { Print [v, OnElementsOf Omega, Smoothing, File "almi5_surface_v.msh"]; Print [normE, OnElementsOf Omega, Smoothing, File "almi5_surface_normE.msh"]; }