Skip to content
Snippets Groups Projects
Commit 4f0e33d1 authored by Axel Modave's avatar Axel Modave
Browse files

AcademicWaves: update + more onelab (EM not checked)

parent dcfd65cc
No related branches found
No related tags found
No related merge requests found
Pipeline #3059 passed
DIR = "out/";
PBM_GAUSS_CUBOID = 1;
PBM_GAUSS_SQUARE = 2;
PBM_LOUDSPEAKER = 3;
PBM_SCATT_KITE = 4;
PBM_SCATT_PLATES = 5;
PBM_SCATT_SPHERE = 6;
RES_FREQ = 1;
RES_TIME = 2;
DefineConstant[
FLAG_PBM = {PBM_GAUSS_SQUARE,
Name "Input/0Type of problem", Highlight "Blue",
GmshOption "Reset", Autocheck 0,
Choices {PBM_GAUSS_SQUARE = "Gaussian pulse in square (2D, time)",
PBM_GAUSS_CUBOID = "Gaussian pulse in cuboid (3D, time)",
PBM_LOUDSPEAKER = "Loudspeaker (2D, freq)",
PBM_SCATT_KITE = "Scattering by Kite (2D, freq)",
PBM_SCATT_PLATES = "Scattering by Plates (2D, freq)",
PBM_SCATT_SPHERE = "Scattering by Sphere (3D, freq)"}}
];
If (FLAG_PBM==PBM_GAUSS_CUBOID)
LinkGeo = "gaussInCuboid.geo";
LinkPro = "gaussInCuboid.pro";
EndIf
If (FLAG_PBM==PBM_GAUSS_SQUARE)
LinkGeo = "gaussInSquare.geo";
LinkPro = "gaussInSquare.pro";
EndIf
If (FLAG_PBM==PBM_LOUDSPEAKER)
LinkGeo = "loudspeaker.geo";
LinkPro = "loudspeaker.pro";
EndIf
If (FLAG_PBM==PBM_SCATT_KITE)
LinkGeo = "scattKite.geo";
LinkPro = "scattKite.pro";
EndIf
If (FLAG_PBM==PBM_SCATT_PLATES)
LinkGeo = "scattPlates.geo";
LinkPro = "scattPlates.pro";
EndIf
If (FLAG_PBM==PBM_SCATT_SPHERE)
LinkGeo = "scattSphere.geo";
LinkPro = "scattSphere.pro";
EndIf
DefineConstant[
C_ = {"-solve -pos -bin -v2", Name "GetDP/9ComputeCommand", Visible 0 }
];
Include "main.dat" ;
ORDER = 1;
SetOrder ORDER;
Mesh.ElementOrder = ORDER;
Mesh.SecondOrderLinear = 0;
If (FLAG_PBM!=0)
Include Str[LinkGeo] ;
EndIf
Include "main.dat" ;
//===========================================
// JACOBIAN & INTEGRATION
//===========================================
Jacobian {
{ Name JVol; Case {{ Region All; Jacobian Vol; }}}
{ Name JSur; Case {{ Region All; Jacobian Sur; }}}
{ Name JLin; Case {{ Region All; Jacobian Lin; }}}
}
Integration {
{ Name I1;
Case {
{ Type Gauss;
Case {
{ GeoElement Point; NumberOfPoints 1; }
{ GeoElement Line; NumberOfPoints 4; }
{ GeoElement Line2; NumberOfPoints 6; }
{ GeoElement Triangle; NumberOfPoints 6; }
{ GeoElement Triangle2; NumberOfPoints 12; }
{ GeoElement Tetrahedron; NumberOfPoints 15; }
{ GeoElement Tetrahedron2; NumberOfPoints 15; }
}
}
}
}
}
//===========================================
// LOAD SPECIFIC .PRO
//===========================================
If (FLAG_PBM!=0)
Include Str[LinkPro] ;
EndIf
DefineConstant[
C_ = {"-solve -pos -bin -v2", Name "GetDP/9ComputeCommand", Visible 0 }
];
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment