← Back to team overview

fenics team mailing list archive

Re: FEM for Electromagnetics

 

On Wed, Mar 26, 2008 at 11:46:44AM -0500, Andy Ray Terrel wrote:
> Anders Logg wrote:
> > On Wed, Mar 26, 2008 at 08:08:10AM -0500, Andy Ray Terrel wrote:
> >   
> >> Hello Evan,
> >>
> >> I would look at start at the demos (dolfin/src/demo/pde). The 
> >> mixedpoisson and stokes demos have vector elements, unfortunately in my 
> >> quick test it looks like Nedelec hasn't been fully integrated in dolfin 
> >> just yet. The convectiondiffusion demo gives some time stepping routines 
> >> that might be helpful.
> >>
> >> -- Andy
> >>     
> >
> > What's the problem with Nedelec? I thought it worked.
> >
> >   
>  
> First time I tried using it I got:
> 
> terminate called after throwing an instance of 'std::runtime_error'
>   what():  evaluate_dof not implemented for this type of element
> Abort trap
> 
> Which comments in ffc claim that nedelec is not supported for the 
> evaluate_basis functions.
> 
> If I just stick to the mixed poisson demo it works fine, but it is 
> producing a different picture than the BDM element.  Of course this 
> could just be that the plotting is not really doing the right thing.
> 
> I guess it really depends on what you mean by working, although the 
> issues I bring up could not be important for a specific user.
> 
> -- Andy

Yes, evaluate_dof is needed for setting BCs and this function was long
missing, but it should be there now. (Marie and Rob added the missing
pieces in FIAT with the class Functionaltype.)

-- 
Anders



> ====================
> Andy Terrel
> Computer Science Dept
> University of Chicago
> aterrel@xxxxxxxxxxxx
> ---------------------
> 
> Failure has a thousand explanations,
> Success doesn't need one.
>                    -Sir Alec Guiness
> 
> _______________________________________________
> FEniCS-dev mailing list
> FEniCS-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/fenics-dev


References