[deal.II] On using TrilinosWrappers::SolverDirect

Dear all I am using the new TrilinosWrappers::SolverDirect which allows two step solution: initialize and solve. Since I want to calculate LU decomposition only once, I do something like this static int first_time = 1; static TrilinosWrappers::SolverDirect::AdditionalData data (fals

[deal.II] TrilinosWrappers::SolverDirect

Dear all I am using the new TrilinosWrappers::SolverDirect which allows two step solution: initialize and solve. Since I want to calculate LU decomposition only once, I do something like this static int first_time = 1; static TrilinosWrappers::SolverDirect::AdditionalData data (fals

Re: [deal.II] Adding dirichlet values to ConstraintMatrix

On 06/23/2016 02:46 AM, Praveen C wrote: void add_dirichlet_constraints (const std::map &values, ConstraintMatrix &constraints) { for (const auto &pair : values) { constraints.add_line (pair.first); constraints.se

Re: [deal.II] Fast FE evaluation over several points and several functions

Thanks for the fast reply, this is exactly what I needed! El miércoles, 22 de junio de 2016, 19:38:35 (UTC+2), bangerth escribió: > > > > I am solving an eigenvalue problem similar than step-36. After solving > for the > > eigenpairs, I evaluate the eigenfunctions in the standard way: > > > >

Re: [deal.II] Adding dirichlet values to ConstraintMatrix

Thanks Wolfgang. I will test this. Following function should work. best praveen void add_dirichlet_constraints (const std::map &values, ConstraintMatrix &constraints) { for (const auto &pair : values) { constraints.add_line (pair