> 1. Does the package use legendre or lagrange polynomials or is there some > way to define which one to use?
The documentation to the FE_Q class describes these issues. > 2. I am trying to get the value at the end point of the neighbor element > for the flux. Right now I am trying to use Gauss Lobatto points and using > the last quadrature point to get the boundary value. Is there a simpler way > to access the values at the end points for each element? If you use a QTrapez quadrature formula, you get a quadrature formula that also has a quadrature point that sits at the end points of your elements. > 3. I am confused with the "reinit" function. Does it matter if I don't give > this command. I wanted to make sure because every program seems to use > this. There are several classes in deal.II that have a reinit() function. Which one do you mean? > 4. Below is the code I've written for assembling the system. When I try to > access get_function_values(), I get the folowing error: You try to evaluate the finite element field given by the nodal vector old_solution (or solution, I don't know which call produces the problem) on dof_handler. But dof_handler describes a mesh with 32 degrees of freedom while the vector you pass has length zero. I suspect you didn't initialize old_solution yet? W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
