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 _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
