// simnibs automatically generated file for getdp simulation // file created with getdp.py on 20140415_120207 // Andre Antunes WM = 1; GM = 2; CSF = 3; Bone = 4; Scalp = 5; fn_dadt = "almi5_volume.msh"; Group { WM = Region[WM]; GM = Region[GM]; CSF = Region[CSF]; Bone = Region[Bone]; Scalp = Region[Scalp]; Omega = Region[{WM,GM,CSF,Bone,Scalp}]; } Function { dadt[] = VectorField[XYZ[]]{0}; sigma[WM] = 0.126; sigma[GM] = 0.275; sigma[CSF] = 1.654; sigma[Bone] = 0.01; sigma[Scalp] = 0.465; } Jacobian { { Name Volume; Case {{ Region Omega; Jacobian Vol; }} } } 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}]; 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; } } }} 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"]; }