Fine by me. That matches well how I envision the domains attached to ufl forms as well.
Martin Den 7. feb. 2013 09:17 skrev "Garth N. Wells" <[email protected]> følgende: > I'm about to start cleaning up the assemblers, and fixing some bugs > will necessitate some interface changes. For example we have: > > /// Assemble tensor > void assemble(GenericTensor& A, > const Form& a, > bool reset_sparsity=true, > bool add_values=false, > bool finalize_tensor=true, > bool keep_diagonal=false); > > /// Assemble tensor on sub domains > void assemble(GenericTensor& A, > const Form& a, > const MeshFunction<std::size_t>* cell_domains, > const MeshFunction<std::size_t>* exterior_facet_domains, > const MeshFunction<std::size_t>* interior_facet_domains, > bool reset_sparsity=true, > bool add_values=false, > bool finalize_tensor=true, > bool keep_diagonal=false); > > but calling > > assemble(A, a, 0, 0, 0); > > is to me ambiguous in terms which function will be called. I think the > best solution, which we've discussed, it to always attach the domain > data to forms. > > Garth > > _______________________________________________ > Mailing list: https://launchpad.net/~dolfin > Post to : [email protected] > Unsubscribe : https://launchpad.net/~dolfin > More help : https://help.launchpad.net/ListHelp >
_______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp

