> That looks good, but `P` is already used as a name for the P
> elements in the latest draft so this does not seem to work out.
Would it work if the syntax were
("P-", cell, r, k=0)
("P", cell, r, k=0)
("Q-", cell, r, k=0)
("S", cell, r, k=0)
i.e., if the fourth argument is left out it defaults to 0.
That way ("P", cell, r) would continue to have the same meaning.
> Are there other options? Something like
>
> FEEC("P", cell, r, k)
Of course I am happy to have FEEC built into the syntax
if you and others think that is appropriate.
On 03/03/2014 12:51 PM, Anders Logg wrote:
On Mon, Mar 03, 2014 at 12:21:11PM -0600, Douglas N Arnold wrote:
Are you happy with these names? Is `P Lambda` the best choice?
It seems a bit odd with a space in the name but I can't think of
anything better.
Actually my choice would be
("P-", cell, r, k)
("P", cell, r, k)
("Q-", cell, r, k)
("S", cell, r, k)
The Lambda is always present so conveys no information.
That looks good, but `P` is already used as a name for the P
elements in the latest draft so this does not seem to work out.
Are there other options? Something like
FEEC("P", cell, r, k)
?
--
Anders
Not that I'm aware of, but we could easily add a function `simplex`
which would allow this:
FiniteElement("P Lambda", simplex(n), r, k)
Seems like a good idea to me.
On 03/03/2014 11:54 AM, Anders Logg wrote:
On Mon, Mar 03, 2014 at 09:49:39AM -0600, Douglas N Arnold wrote:
Let me recall that there is one way to specify
elements that is truly consistent and leads to
dimension free code:
(name, cell, r, k)
where name is 'P- Lambda' or 'P Lambda' for cell
= 'interval', 'triangle', or 'tetrahedron',
and name is 'Q- Lambda' or 'S Lambda' for
cell = 'interval', 'quadrilateral', or
'hexahedron'. From one point of view
(the one conveyed by the periodic table), everything
else is a convenience alias.
Are you happy with these names? Is `P Lambda` the best choice?
It seems a bit odd with a space in the name but I can't think of
anything better.
BTW, is there a dimension-independent way
to refer to a simplex or cube, (e.g., (simplex, 2)
instead of triangle) in UFL? If you want code
that changes dimension by just changing the one
line "n = 2", I guess you need this.
Not that I'm aware of, but we could easily add a function `simplex`
which would allow this:
FiniteElement("P Lambda", simplex(n), r, k)
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
