18 template <
typename scalar>
24 graph = Teuchos::rcp(
new Graph(_map, domain->elementsList));
39 template <
typename scalar>
40 Teuchos::RCP<Tpetra::Vector<scalar, Map::local_ordinal_type, Map::global_ordinal_type>>
46 std::ifstream inputfile;
50 RCP<Tpetra::Vector<scalar, Map::local_ordinal_type, Map::global_ordinal_type>> x_0 =
51 rcp(
new Tpetra::Vector<scalar, Map::local_ordinal_type, Map::global_ordinal_type>(map->mapOutput,
true));
52 RCP<Tpetra::Vector<scalar, Map::local_ordinal_type, Map::global_ordinal_type>> x =
53 rcp(
new Tpetra::Vector<scalar, Map::local_ordinal_type, Map::global_ordinal_type>(map->mapDofs,
true));
57 auto line_size = map->mapOutput->getGlobalNumElements();
60 const int ensemble_size = ET::size;
65 for (
auto j = 0; j < ensemble_size; ++j)
66 inputfile >> ET::coeff(tmp, j);
68 x_0->replaceGlobalValue(global_index, tmp);
73 Tpetra::Export<> exportDof(map->mapOutput, map->mapDofs);
75 x->doExport(*x_0, exportDof, Tpetra::ADD);