← Back to team overview

ffc team mailing list archive

Re: [HG FFC] Several updates for handling subdomain integrals:

 

On Thu, Jun 12, 2008 at 12:04:35PM -0600, Jake Ostien wrote:
> 
> 
> Anders Logg wrote:
> > On Thu, May 29, 2008 at 07:18:13PM -0400, Jake Ostien wrote:
> >   
> >> Is the update mentioned below also supposed to work for interior 
> >> facets?  If I just specify dS(2), I get
> >>
> >>     *** list index out of range
> >>     
> >
> > I checked this and it works fine for me. I can do dS(2), even
> > dS(999999999999999).
> >
> > What exact code did you try?
> >
> >   
> 
> Hi,
> 
> Weird.  I use the following code:
> 
> element = FiniteElement("CG", "triangle", 1)
> 
> v = TestFunction(element)
> u = TrialFunction(element)
> 
> a = dot(avg(v),jump(u))*dS(3)
> 
> 
> Then I simple call ffc as:
> 
> % ffc -l dolfin test.form
> 
> and I get the following error:
> 
> ...
> Compiler phase 4: Generating code
> ---------------------------------
> 
>   Generating code for finite elements... done
>   Generating code for finite dof maps... done
>   Generating code for cell integrals... done
>   Generating code for exterior facet integrals... done
>   Generating code for interior facet integrals... *** list index out of 
> range
> 
> Something might be fishy over here, but I pulled code this morning just 
> in case. 
> 
> Can you think of anything else that I might need to update?

I never tried to compile a form with dS(n). I just checked that
creating it the integral dS(n) worked. Anyway, it should be fixed now
so try again.

-- 
Anders

Attachment: signature.asc
Description: Digital signature


References