On Tue, Jun 6, 2017 at 8:52 AM, Barna Becsek <barnabec...@gmail.com> wrote:
> Hello, > > I am in a situation where I need to setup test functions for an > L2projection. Essentially, all I have to access in the end is get_phi() and > get_JxW(). > Now, I need to call fe->reinit() for every element that I add. Since I do > not need anything else than get_phi() and get_JxW(), I would like to lean > the computation as much as possible. I saw that there are a lot of boolean > variables that control what is executed in FE and FE_map classes, e.g. > calculate_phi, calculate_dphi, calculate_xyz, calculate_dxyz,… I saw that > if the user does not set anything then determine_calculations() in reinit() > will just set everything to true. > What is the correct way to do lean the computation? Do I call get_phi() > and get_JxW() before calling reinit() Yes, this is correct. That's exactly why we do this in the examples, e.g. https://github.com/libMesh/libmesh/blob/master/examples/introduction/introduction_ex3/introduction_ex3.C#L236 Please let us know if you have any problems. > (the boolean flags are protected)? How can I make sure that all the second > derivatives are not calculated in FE’s _fe_map? > > Sincerely, > > Barna > ------------------------------------------------------------ > ------------------ > Check out the vibrant tech community on one of the world's most > engaging tech sites, Slashdot.org! http://sdm.link/slashdot > _______________________________________________ > Libmesh-users mailing list > Libmesh-users@lists.sourceforge.net > https://lists.sourceforge.net/lists/listinfo/libmesh-users > ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
