[Getdp] Pyramids not working ?
Kubicek Bernhard
Bernhard.Kubicek at arsenal.ac.at
Fri Jul 28 11:46:58 CEST 2006
Sorry, I just found out that I only replied to single persons, instead of the whole mailling-list:
---mail1---
for scalar, node attached Dofs pyramid cells seem to work, with my #define line .. I prepared some nice "proof" images:
http://wiki.kariert.org/index.php/GetDP
nice greetings,
Bernhard
---mail2---
I am using Pyramidal cells, however i had to enable them in the source code: file include/GetDP.h add a line containing #define NEW_PYRAMIDS
However I did not verify 100% wether the calculations then are really correct. Still, i found the same line as you: { GeoElement Pyramid ; NumberOfPoints 8 ; } But I am still wondering how a pyramid can have 8 points...
nice greetings, hope it works.
Bernhard Kubicek
-----Ursprüngliche Nachricht-----
Von: getdp-bounces at geuz.org [mailto:getdp-bounces at geuz.org] Im Auftrag von Ruth V. Sabariego
Gesendet: Freitag, 28. Juli 2006 10:15
An: Thomas Jung
Cc: getdp at geuz.org
Betreff: Re: [Getdp] Pyramids not working ?
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/
_______________________________________________
getdp mailing list
getdp at geuz.org
http://www.geuz.org/mailman/listinfo/getdp