Thanks for this suggestion - I will have a go at this as well. On Tuesday, July 31, 2018 at 8:05:11 PM UTC+1, Jean-Paul Pelteret wrote: > > Hi Jane, > > I’m not quite sure that I understand exactly what you’re trying to do, but > as an additional consideration to Bruno’s suggestions is the > FEFieldFunction > <https://www.dealii.org/9.0.0/doxygen/deal.II/classFunctions_1_1FEFieldFunction.html> > class > possibly of any use to you in this situation? > > Best, > Jean-Paul > > On 31 Jul 2018, at 18:08, Bruno Turcksin <bruno.t...@gmail.com > <javascript:>> wrote: > > Le mar. 31 juil. 2018 à 11:49, Jane Lee <jane...@maths.ox.ac.uk > <javascript:>> a écrit : > >> What i then have got is an extra bit of domain added onto my domain 1. To >> solve for another variable in this new domain, I need to know f in my >> domain 1 (which I have solved for), and for the rest of the added domain >> (let's call it 2, where the new domain is the union of 1 and 2) I know what >> the variable f is, but basically I don't know how to tell the next solver >> what f is in domain 2. For the bits in domain 1, i could just do something >> like: fe_values.get_function_values if my domain remained as domain 1. But >> since my domain has changed, I'm not sure what to do. >> > So if I understand correctly you have f on domain 1 which you know > numerically, g on domain 2 that you know analytically, and h which is equal > to f on domain 1 and g on domain 2 and your problem is that you can't use > VectorTools::interpolate because you don't have an analytical formula for > f. Can you say that alpha is a function that is equal to f on domain 1 and > zero on domain 2 and use fe_values.get_function_values. Similarly say that > beta is equal to g on domain 2 and zero on domain 1 and use > VectorTools::interpolate. Now you have two vectors one associate to alpha > and to beta and you can simply add them to get the evaluation of h. Does > that do what you want? > > Best, > > Bruno > > -- > 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 dealii+un...@googlegroups.com <javascript:>. > For more options, visit https://groups.google.com/d/optout. > > >
-- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
