Hi,
I want to obtain information of dphi/dxi on the elements at the boundary,
using
fe_face_ptr->get_dphidxi(); but this returns a null vector.
However, fe_face_ptr->get_phi() can give me values
fe_face_ptr = fe_face.get(); fe_face is one element on the boundary, as
shown below:
// define fe_face
AutoPtr<FEBase> fe_face (FEBase::build(dim, fe_type));
QGauss qface(dim-1, fe_type.default_quadrature_order());
fe_face->attach_quadrature_rule (&qface);
// at the boundary
if (elem->neighbor(side) == NULL)
{ fe_face->reinit(elem, side);
fe_face_ptr = fe_face.get(); // some functions define the argument as a
normal pointer
fe_face_ptr->get_phi(); // OK gives the number
const vector <vector <Real > >& dphi_i_q = fe_face_ptr->get_dphidxi();
/// not OK, dphi_i_q.size() =0 ??
Thanks,
Walter
------------------------------------------------------------------------------
"Accelerate Dev Cycles with Automated Cross-Browser Testing - For FREE
Instantly run your Selenium tests across 300+ browser/OS combos. Get
unparalleled scalability from the best Selenium testing platform available.
Simple to use. Nothing to install. Get started now for free."
http://p.sf.net/sfu/SauceLabs
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
