Dear Jesse, You can probably just provide an instance of the GrowingVectorMemory <https://www.dealii.org/8.5.0/doxygen/deal.II/classGrowingVectorMemory.html> class templated on BlockVector< double >. So this would look something like: GrowingVectorMemory< BlockVector<double> > memory; const SchurMatrix schur_matrix (A_inv, B, Dt, C, memory); ...
If its of any help to you, and if you have a C++11 compatible compiler, then you could alternatively use the (relatively) new LinearOperator <https://www.dealii.org/8.5.0/doxygen/deal.II/classLinearOperator.html> class and the schur_complement <https://www.dealii.org/8.5.0/doxygen/deal.II/group__LAOperators.html#ga76acca911f21089cd3bb385d20ccc995> operator. It provides a more flexible framework with which you could perform the same operation. Kind regards, Jean-Paul -- 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]. For more options, visit https://groups.google.com/d/optout.
