On Fri, May 1, 2009 at 9:03 AM, Vijay S. Mahadevan <vija...@gmail.com> wrote: > Hi, > > This might sound ridiculous but I was wondering if there is a scenario > when the quadrature rule you created for say, 2 dimensions and SIXTH > order returns zero points ? Is there ever a need for initing > quadrature rules explicitly ? > > To give you a background, I have an ImplicitSystem and an > ExplicitSystem in my EquationSystems. I want to update the solution in > ExplicitSystem (solution_exp) based on the solution in the > ImplicitSystem (solution_imp). Hence, I loop over the local elements > in the mesh, perform user code and compute solution_exp. But at the > beginning of this loop, I attach the quadrature rule to the FEType of > variable in ImplicitSystem which returns zero points (N_Q_Points=0). > The order and dimension variables in the rule are correct but just > that the points have not been computed. > > I have tried to debug this code for a while but I cannot get past the > zero quadrature points. Any ideas or suggestions to get further on > this will be extremely helpful. > > I am terribly sorry if I was vague or confusing or if this is > downright absurd. I spent an entire day trying to get to the root of > this problem but got nowhere. Please feel free to let me know if you > need more details. Thanks a bunch in advance.
Which element type, triangles? We may need you to post your code or at least email it to one of us to take a closer look... -- John ------------------------------------------------------------------------------ Register Now & Save for Velocity, the Web Performance & Operations Conference from O'Reilly Media. Velocity features a full day of expert-led, hands-on workshops and two days of sessions from industry leaders in dedicated Performance & Operations tracks. Use code vel09scf and Save an extra 15% before 5/3. http://p.sf.net/sfu/velocityconf _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
