[Getdp] Pyramids not working ?

Ruth V. Sabariego r.sabariego at ulg.ac.be
Fri Jul 28 10:15:19 CEST 2006


Indeed, it is not yet implemented. You are using BF_Edge for {a}, what 
makes BF_CurlEdge for {d a}.
For the time being, your mesh should contained either only hexahedrons 
or only tetrahedrons.

Ruth


Thomas Jung wrote:

>Hi, 
>
>I am using a 3D a-v-formulation, up to now I used only tetrahedrons, and 
>everything was fine.
>I have a different grid now, containing also hexahedrons and pyramids, and get 
>the message:
>
>===========================================================
>P r e - P r o c e s s i n g . . .
>Operation : Treatment Formulation 'MagnetoDynamics_a_3D'
>Info      :   Generate ExtendedGroup '_CO_Entity_14' (EdgesOf)
>Info      :   Generate ExtendedGroup '_CO_Entity_15' (EdgesOfTreeIn)
>Info      :   Generate ExtendedGroup '_CO_Entity_11' (NodesOf)
>Info      :   Generate ExtendedGroup '_CO_Entity_12' (NodesOf)
>Resources : cpu 2.02813 s
>E n d   P r e - P r o c e s s i n g
>P r o c e s s i n g . . .
>Operation : Generate[A2]
>Solver    : Loading parameter file 'solver.par'
>Info      : Setting System {A,b} to zero
>Error     : Unkown type of Element in BF_CurlEdge
>
>============================================================
>
>Here are the Integration, Functionspace, and formulation:
>
>Integration {
>    { Name CurlCurl ;
>    Case { {Type Gauss ;
>    Case { { GeoElement Triangle    ; NumberOfPoints  4 ; }
>    { GeoElement Quadrangle  ; NumberOfPoints  4 ; }
>    { GeoElement Tetrahedron ; NumberOfPoints  4 ; }
>    { GeoElement Hexahedron  ; NumberOfPoints  6 ; }
>    { GeoElement Prism       ; NumberOfPoints  9 ; } 
>    { GeoElement Pyramid     ; NumberOfPoints  8 ; } }
>    }
>    }
>}
>}
>
>FunctionSpace {
>  { Name v_Ele; Type Form0;
>    BasisFunction {
>      { Name sn; NameOfCoef vn; Function BF_Node;
>        Support DomainC_Mag; Entity NodesOf[ All ]; }
>    }
>    Constraint {
>      { NameOfCoef vn; EntityType NodesOf; 
>        NameOfConstraint ElectricScalarPotential; }
>    }
>  }
>
>  // Magnetic vector potential a (b = curl a)
>  { Name Hcurl_a_Mag_3D; Type Form1;
>    BasisFunction {
>	{ Name se; NameOfCoef ae; Function BF_Edge;
>        Support Domain_Mag; Entity EdgesOf[ All ]; }
>    }
>    Constraint {
>	{ NameOfCoef ae; EntityType EdgesOf;NameOfConstraint a; }
>	{ NameOfCoef ae; EntityType EdgesOfTreeIn; EntitySubType StartingOn; 
>NameOfConstraint Gauge; }
>    }
>  }
>
>}
>
>Formulation {
>
>  { Name MagnetoDynamics_a_3D; Type FemEquation;
>    Quantity {
>      { Name a ; Type Local; NameOfSpace Hcurl_a_Mag_3D; }
>      { Name v; Type Local; NameOfSpace v_Ele; }
>    }
>    Equation {
>	Galerkin { [ nu[] * Dof{d a} , {d a} ]; In Domain_Mag; Jacobian Vol; 
>Integration CurlCurl; }
>
>	Galerkin { DtDof [ sigma[] * Dof{a} , {a} ]; In DomainC_Mag; Jacobian Vol; 
>Integration CurlCurl; }
>
>	Galerkin { [ sigma[] * Dof{d v} , {a} ]; In DomainC_Mag; Jacobian Vol; 
>Integration CurlCurl; }
>
>	Galerkin { DtDof [ sigma[] * Dof{a} , {d v} ]; In DomainC_Mag;	Jacobian Vol; 
>Integration CurlCurl; }
>
>	Galerkin { [ sigma[] * Dof{d v} , {d v} ]; In DomainC_Mag; Jacobian Vol; 
>Integration CurlCurl; }
>
>
>    }
>  }
>
>}
>
>
>
>Thank you very much for any hint !
>
>
>
> Thomas Jung
>_______________________________________________
>getdp mailing list
>getdp at geuz.org
>http://www.geuz.org/mailman/listinfo/getdp
>
>
>  
>


-- 
Dr. Ir. Ruth V. Sabariego
University of Liege
Dept. of Electrical Engineering (ELAP), Institut Montefiore
Sart Tilman Campus, Grande Traverse, 10 (B28), B-4000 LIEGE, Belgium
phone: +32-4-3663737 -- fax: +32-4-3662910 -- http://elap.montefiore.ulg.ac.be/