[Getdp] questions about higher order BF
gdemesy at physics.utoronto.ca
gdemesy at physics.utoronto.ca
Fri Mar 19 16:01:28 CET 2010
Dear Getdp Developpers,
I've successfully solved a vector Helmholtz equation with
quasi-periodic (or Bloch) conditions in two directions with Getdp
(e.g. a crossed-grating 3D vector harmonic problem).
Now I would like to add some BF to my FunctionSpace for a faster
convergence as a function of the mesh refinement. I have some
questions about these additionnal BF (I would like to have a full
second order). Maybe I simply don't get what the X in BF_Edge_XE or
BF_Edge_XF stands for.
As far as I understand, for the curl-conforming tetrahedral element,
BF_edge, BF_edge_2E and BF_edge_2F should somehow be related to
something like :
BF_Edge : li grad(lj) ? lj grad(li)
BF_Edge_2E : grad(li lj)
BF_Edge_2F : lk lj grad(li) and lk li grad(lj)
where li is the barycentric coordinate of node i.
When I try to add BF_Edge_2F or BF_Edge_3E, I get the message "You
should never end up here!". Why?
BF_Edge_3F seems to be implemented in BF_Edge_3.cpp but I get the
message "Unknown Function for BasisFunction: BF_Edge_3F" when the
"Problem Definition is being loaded".
Finally, BF_Edge_4E is implemented, but I should first add some Basis
Functions related to faces to my set of BF, shouldn't I ?
(I am using getdp_2.0.0)
Could you please give me some explanations or references ?
Thanks a lot for you time.
Best Regards,
Guillaume Demésy
________________________________________
FunctionSpace {
{ Name Hcurl; Type Form1;
BasisFunction {
{ Name sn; NameOfCoef un; Function BF_Edge;
Support Region[Omega]; Entity EdgesOf[Omega]; }
{ Name sn2; NameOfCoef un2; Function BF_Edge_2E;
Support Region[Omega]; Entity EdgesOf[Omega]; }
// { Name sn3; NameOfCoef un3; Function BF_Edge_2F;
// Support Region[Omega]; Entity EdgesOf[Omega]; }
}
Constraint {
{ NameOfCoef un; EntityType EdgesOf ; NameOfConstraint BlochX; }
{ NameOfCoef un; EntityType EdgesOf ; NameOfConstraint BlochY; }
{ NameOfCoef un; EntityType EdgesOf ; NameOfConstraint Dirichlet; }
{ NameOfCoef un2; EntityType EdgesOf ; NameOfConstraint BlochX; }
{ NameOfCoef un2; EntityType EdgesOf ; NameOfConstraint BlochY; }
{ NameOfCoef un2; EntityType EdgesOf ; NameOfConstraint Dirichlet; }
// { NameOfCoef un3; EntityType EdgesOf ; NameOfConstraint BlochX; }
// { NameOfCoef un3; EntityType EdgesOf ; NameOfConstraint BlochY; }
// { NameOfCoef un3; EntityType EdgesOf ; NameOfConstraint Dirichlet; }
}
}
}
________________________________________
----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.