On Thu, Mar 21, 2019 at 8:17 AM Povolotskyi, Mykhailo <mpovo...@purdue.edu> wrote:
> Thank you, John, > > at least now we are on the same page. > > Let me tell you what I'm doing, and hopefully you can suggest an effective > way without using protected methods. > > I'm solving a Newmark system with vector finite elements (Nedelec > elements). > > I need to compute solution, curl of the solution and derivative over time > of the solution at a specific set of points for every time step. > > To do so I pre-compute local coordinates of the points, then I do > something like this: > > FEGenericBase<RealGradient>* fe_new = > fecontext.cached_fe<RealGradient>(elem->dim(), fe_type); > fe_new->reinit (elem, &(local_points)); > > After this I have access to fe_new->get_phi() and can compute all what I > need. > > What would you suggest? > Would one of the various public get_element_fe() APIs in FEMContext not work for your case? -- John _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
