Timo, I'm making my way through unanswered email and yours seems to be one from a long time ago :-(
> Printing the polynomials I got -1 x^1 1 x^0 and 1 x^1 0 x^0. Outputting > the indices of the TensorProductPolynomials made me assume that the > following test functions are used: > 0 (-x+1)(-y+1) > 1 x(-y+1) > 2 (-x+1)y > 3 xy Yes, that's correct. > I then tried to reproduce the matrix of my actual problem and failed. So > I started integrating just one single test function and it was just > fine, same results using deal.II and doing it by hand. > But as soon as I compute a product of two functions I get differing > results. E.g. integrating the product 0*0 using deal.II on my mesh > results in 3.90625e-08, integrating by hand gives 6.9444444e-08. > Where the JxW value is the same and equals 6.25e-07. > > I'm pretty sure deal.II is working, so somehow I'm doing some things > wrong here. I would be very glad if you can help me out. Have you figured this one out? If not, what's the bilinear form and what's the quadrature formula you use? Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
