Re: [deal.II] Re: imposing boundary conditions for a three-fields formulation

Dear Alberto, Now I understand your concern. The boolean flag that is passed as the final parameter to MatrixTools::apply_boundary_values determines whether column elimination is performe

Re: [deal.II] Re: imposing boundary conditions for a three-fields formulation

Dear Jean-Paul, thanks for your reply. I see your point. It was my understanding (wrong, I guess) that in imposing Dirichlet BC deal.II would modify the rows but also the columns, adjusting the rhs accordingly. In this way, the symmetry of the system matrix would have been preserved. I see now th

Re: [deal.II] Re: QProjector::project_to_face and Gauss Lobatto

GLL points include the endpoints of the domain...they are useful when you are trying to solve large-scale problems and have integration points same as the interpolation points...e.g. Lagrange interpolation with gll points etc, a very standard technique used in spectral element methods. I'm curious

Re: [deal.II] Building UMFPACK with VS 2015

On 04/19/2017 09:18 AM, Bill Greene wrote: I am continuing to try to use DEAL with Visual Studio. I have had some success but ran into a problem when I tried to include the bundled umfpack via the Cmake directive -D DEAL_II_WITH_UMFPACK=ON Deal II seemed to build without problems but when I tri

[deal.II] Building UMFPACK with VS 2015

I am continuing to try to use DEAL with Visual Studio. I have had some success but ran into a problem when I tried to include the bundled umfpack via the Cmake directive -D DEAL_II_WITH_UMFPACK=ON Deal II seemed to build without problems but when I tried to build example step-56, I got unsatisfie

[deal.II] Re: QProjector::project_to_face and Gauss Lobatto

vganpat, I am trying to solve a problem that has both the volume integral and a > surface integral. > You would normally just create a quadrature rule for the volume and the surface integral separately. > If I am using a higher order polynomial approximation and use Gauss > Lobatto quadratu