Hello, I'm trying to solve the 2D heat equation with the method of lines, instead of Rothe's method as explained in the tutorials, as I'd like to use ARKODE to march in time. Assuming that the *boundary data is 0* and the initial data is *x*y*(1-x)*(1-y)*, after writing the weak formulation, I have a large system of ODEs:
Mu' +Au=F u(0)=u_0 I started with modifying step-4 so that I assembled correctly the mass matrix. Before writing the lambdas that ARKODE needs, I need a way to enforce those Dirichlet boundary conditions for my system, since so far I am only using * MatrixTools::apply_boundary_values* applied to the stiffness matrix. I searched for a solution in several steps, but apparently I can't find a way. How can I impose them taking into account also the presence of the mass matrix? With kind regards, Bob -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/22786f10-2153-4cd7-8142-5fa542f0f662n%40googlegroups.com.
