Hi, I would like to use trilinos as the solver for laplacian examples provided in your tutorial. I have tried to start with step-33 example and replace the ConservationLaw by LaplaceProblem from step-3. I would like to keep the boundary conditions definition in some other class just like in EulerEquation.
I have a problem setting the dirichlet boundary conditions. It seems that there is no void MatrixTools::apply_boundary_values for Trilinos "matrices". How can I force dirichlet condition when using trilinos solver? Do you have examples regarding the heat equation problem with non-homogeneous neuman and robin boundary condtions? Thanks for your help _______________________________________________
