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() (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
