> I'm coming unstuck at correctly setting up the Trilinos matrices and > vectors. I've tried numerous permutations of the Trilinos::SparseMatrix > and Trilinos::MPI::Vector setup and initialisations using deal.II > matrices and either end up with some sort of segmentation fault, a > Trilinos error "-1" or, for my best attempt, a hung solver. I can get > the local dofs using the > "DoFTools::dof_indices_with_subdomain_association" function (it appears > that ghost cells might be included in the filled IndexSet),
Yes. I figure what you probably need is a *disjoint* partitioning of the degrees of freedom. DoFTools should have functions for this as well, and I think we use them in step-17/18. W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
