dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Wed, May 18, 2011 at 03:43:19PM +0200, Harish Narayanan wrote:
> 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

Ah, yes. It's supported by the low-level assemble function, but not by
VariationalProblem.

--
Anders