Hi David,

I have updated the whole toolchain. And I have constructed a test example with geometric BC on a DG space with an extruded mesh (see attached code). This seems to work.

Unfortunately my code now breaks with an error (see 2) below). During PyOP2 installation I also got a warning (see 1) below)

Maybe something (unrelated to the geometric BC?) broke while updating PyOP2?

Henrik

1) Warning during installation of pyop2:

/rwthfs/rz/SW/UTIL.common/Python/2.7.9/x86_64/lib/python2.7/distutils/dist.py:267: UserWarning: Unknown distribution option: 'test_requires'

  warnings.warn(msg)

/rwthfs/rz/SW/UTIL.common/Python/2.7.9/x86_64/lib/python2.7/site-packages/setuptools/dist.py:291: UserWarning: The version specified ('0.11.0-415-g450bf5b') is an invalid version, this may not work as expected with newer versions of setuptools, pip, and PyPI. Please see PEP 440 for more details.

  "details." % self.metadata.version

running build_ext

2) Error during execution of test_geometric_BC.py

Traceback (most recent call last):                                       

  File "test_geometric_BC.py", line 37, in <module>                      

    outfile << project(pw,V_out, name="Pressure")                         

  File "/home/hb111949/.local/lib/python2.7/site-packages/firedrake/projection.py", line 97, in project                                                                                        

    form_compiler_parameters=form_compiler_parameters)                                          

  File "/home/hb111949/.local/lib/python2.7/site-packages/firedrake/solving.py", line 120, in solve                                                                                            

    _solve_varproblem(*args, **kwargs)                                                         

  File "/home/hb111949/.local/lib/python2.7/site-packages/firedrake/solving.py", line 147, in _solve_varproblem                                                                                 

    solver.solve()                                                                             

  File "<string>", line 2, in solve                                                            

  File "/home/hb111949/Code/PyOP2/pyop2/profiling.py", line 203, in wrapper                    

    return f(*args, **kwargs)                                                                  

  File "/home/hb111949/.local/lib/python2.7/site-packages/firedrake/variational_solver.py", line 173, in solve                                                                                  

    self.snes.solve(None, v)                                                                   

  File "PETSc/SNES.pyx", line 520, in petsc4py.PETSc.SNES.solve (src/petsc4py.PETSc.c:167382)  

  File "PETSc/petscsnes.pxi", line 253, in petsc4py.PETSc.SNES_Function (src/petsc4py.PETSc.c:32821)                                                                                           

  File "/home/hb111949/.local/lib/python2.7/site-packages/firedrake/solving_utils.py", line 193, in form_function                                                                              

    with ctx._Fs[lvl].dat.vec_ro as v:

  File "/rwthfs/rz/SW/UTIL.common/Python/2.7.9/x86_64/lib/python2.7/contextlib.py", line 17, in __enter__

    return self.gen.next()

  File "/home/hb111949/Code/PyOP2/pyop2/petsc_base.py", line 230, in vec_context

    self._force_evaluation()

  File "/home/hb111949/Code/PyOP2/pyop2/base.py", line 1566, in _force_evaluation

    _trace.evaluate(reads, writes)

  File "/home/hb111949/Code/PyOP2/pyop2/base.py", line 170, in evaluate

    comp._run()

  File "/home/hb111949/Code/PyOP2/pyop2/base.py", line 4042, in _run

    return self.compute()

  File "/home/hb111949/Code/PyOP2/pyop2/base.py", line 4083, in compute

    fun = self._jitmodule

  File "/home/hb111949/Code/PyOP2/pyop2/utils.py", line 64, in __get__

    obj.__dict__[self.__name__] = result = self.fget(obj)

  File "/home/hb111949/Code/PyOP2/pyop2/sequential.py", line 151, in _jitmodule

    direct=self.is_direct, iterate=self.iteration_region)

  File "/home/hb111949/Code/PyOP2/pyop2/caching.py", line 203, in __new__

    obj = make_obj()

  File "/home/hb111949/Code/PyOP2/pyop2/caching.py", line 193, in make_obj

    obj.__init__(*args, **kwargs)

  File "/home/hb111949/Code/PyOP2/pyop2/host.py", line 700, in __init__

    self.compile()

  File "/home/hb111949/Code/PyOP2/pyop2/host.py", line 725, in compile

    if self._kernel._cpp:

AttributeError: 'Kernel' object has no attribute '_cpp'