Hi Denis, Thank you for your answer but I still need to further understand some things : * I think I need to pass the point (xo) during assembly since it will different for different cells ( dof coordinate or center of the cell), and I do not see how I can pass the point when declaring the function (ex. Enrichment <dim> (param1,param2, etc...) ) to create the FE_Enriched constructor, maybe I am missing something here.
* The thing is that during assembly the functions are already encapsulated in the “enrichments” vector and if I even create a member function like Enrichment<dim>::get_point(Point<dim> & ) how do I apply it to the functions that are already in the “enrichments” vector ? Maybe I can make the “enrichments” vector public?... Thank you for your help Edith > On Nov 1, 2017, at 5:29 PM, Denis Davydov <davyd...@gmail.com> wrote: > > Hi Edith, > > f(x-x0) is just your enrichment function derived from dealii::Function. So > define a class which takes that extra parameter (or any number of them) and > use it with FE_Enriched. > Assuming that your bilinear form needs gradients, then you would need to > implement gradients of your custom function so that FE_Enriched can work out > chain rule for you. > > Cheers, > Denis. > > On Wednesday, November 1, 2017 at 5:27:10 AM UTC+1, Edith Sotelo wrote: > Hello, > > I would like to implement enrichment functions of this type: f(x-xo) where xo > is a user defined point. It could be the center of the cell or the dof > coordinate. > I would like a suggestion about how to implement this, because I do not see > how to pass it in the constructor since it is not a static point.I think it > has to be during the assembly step, but FE_values call the enrich function > with a single argument (the quadrature point). > What sould be a way to do pass the addtional point? > > Thank you for you help. > > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <https://groups.google.com/d/forum/dealii?hl=en> > --- > You received this message because you are subscribed to a topic in the Google > Groups "deal.II User Group" group. > To unsubscribe from this topic, visit > https://groups.google.com/d/topic/dealii/fFbITWIMW5M/unsubscribe > <https://groups.google.com/d/topic/dealii/fFbITWIMW5M/unsubscribe>. > To unsubscribe from this group and all its topics, send an email to > dealii+unsubscr...@googlegroups.com > <mailto:dealii+unsubscr...@googlegroups.com>. > For more options, visit https://groups.google.com/d/optout > <https://groups.google.com/d/optout>. Edith -- 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.
