// constructors ParallelEvolver(const operator_type& L, const bc_set& bcs) { evolvers_.reserve(L.size()); for (Size i=0; i < L.size(); i++) { evolvers_.push_back(boost::shared_ptr<Evolver>(new Evolver(L[i], bcs[i]))); } }
// constructors TRBDF2(const operator_type& L, const bc_set& bcs) : L_(L), I_(operator_type::identity(L.size())), dt_(0.0), bcs_(bcs), alpha_(2.0-sqrt(2.0)) {}
// constructors MixedScheme(const operator_type& L, Real theta, const bc_set& bcs) : L_(L), I_(operator_type::identity(L.size())), dt_(0.0), theta_(theta) , bcs_(bcs) {}