Wolfgang Bangerth
Wolfgang Bangerth
Can you get a backtrace for where this error comes from?
A direct solver solving a 2d problem with more than 1M unknowns will need a *very* large amount of memory (I'd estimate several GB to several 10 GB). How much...
What do you get with ``` std::make_unique(FE_BDM(degree), 1, FE_DGP(degree), // same 'degree' 1); ``` ? Some elements have funny definitions of what 'degree' means. For example, for RT, you need...
How do we move forward with this one?
I'm not a huge fan of merging patches we know yield the wrong result. I think you probably want to come up with (a separate set of) tests that check...
Do you think you can come up with a small test case that shows this? The way we often deal with this is that we keep pointers to objects such...
I think that's the key question: all of the `TrilinosWrappers` classes build on Epetra, but we don't have a complete replacement for these classes that would build on Tpetra. We...
@mayrmt Thanks for working on this! As for template parameters: If you just want to convert the existing Epetra wrappers, then think `Scalar=double`, `LocalOrdinal=unsigned int` and `GlobalOrdinal=types::global_dof_index` are appropriate, as...
Yes, the idea is that there is only one or two (deal.II) classes that implement the `ReadWriteVector` interface, and that all other vector types are `VectorSpaceVector`-type classes that import/export their...
Martin -- nice plan, I'm looking forward to the implementation! > The main characteristics and changes would be as follows: > - Indices are stored in MPI-local index space. DoFHandler...