On Wed, Sep 30, 2009 at 9:23 AM, Petry Stefan <[email protected]> wrote:
> unsigned int quadrature_points = 0; > double **Mx_on_qpoints, **My_on_qpoints; > double rVal = 0., phiVal = 0., MrTmp = 0., MphiTmp = 0.; > > MeshBase::const_element_iterator el = > mesh->elements_begin(); // prepare iterator for elements > const MeshBase::const_element_iterator end_el = > mesh->elements_end(); > > Mx_on_qpoints = new double *[mesh->n_elem()]; > My_on_qpoints = new double *[mesh->n_elem()]; > > FEType fe_type; // create and initialize finite element with > appropriate approximation order and quadrature rule > fe_type.order = (libMeshEnums::Order)VAL_ORDER_ES; > fe_type.family = (libMeshEnums::FEFamily)VAL_FAMILY; > AutoPtr<FEBase> fe (FEBase::build(dim, fe_type)); > QGauss qrule (dim, (libMeshEnums::Order)VAL_ORDER_GAUSS); > fe->attach_quadrature_rule(&qrule); > > /* > * we are computing all quadrature points first before we use > it, > * because then we have the right number for e.g. using > triangles and quads > */ > Compute_Quadrature_Points(quadrature_points); Is the value of "quadrature_points" still zero when you call this function? Is it passed by reference? What does this Compute_Quadrature_Points function do? -- John ------------------------------------------------------------------------------ Come build with us! The BlackBerry® Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9-12, 2009. Register now! http://p.sf.net/sfu/devconf _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
