Andreas,
> > How does GetDP do that? Is it a global flag that is set once? > > In GetDP the user defines in the input script an object called "Jacobian", > which can be different for each domain region. When setting up an integral > term in the weak formulation then, you have to specify which jacobian method > will be used to evaluate the integrals. In analogy to Deal.II terminology, > you > define a certain mapping or FEValues object (I am not exactly sure which > class > calculates the jacobian in DealII), which applies to different parts of the > triangulation (e.g. separated by material_id). Then whenever evaluating > integrals in a certain part of the triangulation, the corresponding > mapping/FEValue object would invoke the corresponding appropriate JxW that > corresponds to that triangulartion object. Interesting. That would be some kind of metadata that the triangulation stores -- "I am a cross section of a cylindrical domain". The FEValues class would then query this flag when computing the JxW values and multiply them by 2*pi*r. It's an interesting concept. I'm not sure it's ever going to be particularly high on my TODO list, but it's something I wouldn't be opposed to either if anyone would be interested in implementing it! Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.colostate.edu/~bangerth/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/6b72fc92-7c63-319f-2e21-d78876e9c8a5%40colostate.edu.
