dolfin team mailing list archive

[Harish Narayanan <harish.mlists@xxxxxxxxx>]

On 5/18/11 3:23 PM, Anders Logg wrote:
> On Wed, May 18, 2011 at 03:03:18PM +0200, Marie E. Rognes wrote:
>> On 05/18/2011 02:43 PM, Harish Narayanan wrote:
>>> Hello all,
>>>
>>> I used to be able pass marked exterior facet domains like so in Python.
>>>
>>> VariationalProblem(a, L, exterior_facet_domains=boundary)
>>>
>>> This was useful in setting boundary conditions and doesn't seem to work
>>> any more. Is there a new syntax for doing this?
>>>
>>
>> Yes, attach the mesh function to the relevant form(s):
>>
>> 	L.exterior_facet_domains = boundary

Cool, I will try this.

> The old syntax is still supported.

I get the following when I try to run a cbc.twist example:

Traceback (most recent call last):
  File "twist.py", line 57, in <module>
    u = twist.solve()
  File "cbc.solve/cbc/twist/problem_definitions.py", line 26, in solve
    self.solver = StaticMomentumBalanceSolver(self,
self.parameters["solver_parameters"])
  File "cbc.solve/cbc/twist/solution_algorithms.py", line 85, in __init__
    equation = VariationalProblem(L, a, bcu,
exterior_facet_domains=boundary)
TypeError: __init__() got an unexpected keyword argument
'exterior_facet_domains

Harish