dolfin team mailing list archive

[Dag Lindbo <dag@xxxxxxxxxx>]

Anders Logg wrote:
> On Tue, Jun 17, 2008 at 10:52:06PM +0200, Dag Lindbo wrote:
>> Hello again!
>>
>> Is it correct that Periodic BCs won't work for vector-valued systems?
>> There is a comment to suggest this in PeriodicBC.cpp. Is there a
>> reasonable workaround for the needy?
>>
>> /Dag
> 
> I don't know of an easy workaround.
> 
> But maybe you can take a look at the code in PeriodicBC.cpp to see if
> it's possible to understand what it does and then extend it? I guess
> the extension will be nontrivial, but the code is fairly well
> documented so it should not be impossible.

Paper submission deadline approaching fast so I don't really have a
choice but to fix it!

It currently works for "simple" forms (i.e. scalar non-mixed), so it
seems natural to have
PeriodicBC::apply(GenericMatrix& A, GenericVector& b, const Form& form)
pick out the "simple" finite elements and dof-maps from the form and
call apply(A, b, dof_map, form) for each. Any thoughts?

/Dag

> 
> Otherwise, file a bug report.
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev