I'm trying to compute the xyz, JxW, and other values at arbitrary points for a FEBase object during the evaluation of a MOOSE UserObject. I implemented the following in code:
``` //Create a vector of points std::vector< Point > cell_points; ...things to fill cell_points... //Create a vector of ones (the weights) ones = std::vector< double > (cell_points.size(), 1); std::unique_ptr< libMesh::FEBase > fe( libMesh::FEBase::build(_mesh.dimension(), libMesh::FEType(_current_elem->default_order()))); fe->reinit(_current_elem, &cell_points, &ones); std::vector< Point > xyz = fe->get_xyz(); ``` following basically what I found in the third example. When I run this code I hit the assertion: "Assertion `!calculations_started || calculate_xyz' failed." I have also tried substituting FEBase for FEGenericBase<Real> to identical results. Is there something that I'm missing in my code? _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
