Thanks, Wolfgang. I'll dig into FE_Enriched! On Thursday, June 8, 2023 at 6:57:44 PM UTC-4 Wolfgang Bangerth wrote:
> On 6/8/23 09:40, Alex Quinlan wrote: > > > > I am looking to modify the FE_Q element to have non-polynomial shape > > functions. It looks like FE_Q_Base takes an argument of type > > ScalarPolynomialBase<dim> (fe_q_base.cc: line 418). > > > > I am looking at a 2D element and instead of 2D polynomial shape > functions like > > this: > > > > N = (1 - ξ)(1 - η) > > > > I would like to add some additional terms like this: > > > > N = (1 - ξ)(1 - η) + ζ > > > > where ζ changes at different quadrature points. > > > > It seems that I would need to access the values of ζ before computing > > fe_values, since the shape functions and their derivatives would be > dependent > > on ζ. > > > > Any thoughts on how feasible this would be to implement? > > Alex: > if your shape functions are not polynomial, then this is the wrong > approach -- > FE_Q_Base is derived from FE_Poly, which is specifically built to create > polynomial bases. > > I'm not sure we have a particularly good starting point for a completely > general choice of ζ(ξ,η), but you might want to take a look at what > FE_Enriched is doing. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/97405fd5-9732-4fb2-b1b7-ce6f2d1a155en%40googlegroups.com.
