From: Vijay S. Mahadevan <vijay.m@gm...>  20090501 14:03:57

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. Vijay 