I'm having trouble understanding how to code 2-D Poisson's Equation with Dirichlet boundary conditions. What I have thus far is
- Constructed square mesh with triangular elements
- Assembled stiffness matrix. The function for this also gives the area of each element (all same, since square).
- Use Gaussian Quadrature to construct the right hand side.
edit: I may have been a little vague in my question. Hopefully this is a better rephrasing: Separate from FE, my code for gaussian quadrature works just fine. However, I don't understand how gaussian quadrature can even be relevant/involved with obtaining a solution to 2D Poisson; hence, I don't know how to continue coding.