Oh, this is because the interface counts layers as being the number of top/bottom surfaces, not the number of elements. It really needs numbering as we kept getting bugs because of it!

--cjc
________________________________
From: noreply@github.com [noreply@github.com] on behalf of Lawrence Mitchell [notifications@github.com]
Sent: 10 December 2013 09:55
To: firedrakeproject/firedrake
Subject: Re: [firedrake] Extrusion side bcs new (#113)


In python/firedrake/core_types.pyx:

> @@ -285,7 +296,7 @@ cdef ft.element_t as_element(object fiat_element):
>
>  class _Facets(object):
>      """Wrapper class for facet interation information on a Mesh"""
> -    def __init__(self, mesh, count, kind, facet_cell, local_facet_number, markers=None):
> +    def __init__(self, mesh, count, kind, facet_cell, local_facet_number, markers=None, layers=1):


So I still don't like layers=1 meaning not extruded, because it doesn't right? Can we have layers=None for the non-extruded case?

—
Reply to this email directly or view it on GitHub<https://github.com/firedrakeproject/firedrake/pull/113/files#r8226923>.