From 3b057f1ca1564dfd90395aa67b711b4a7a142a2d Mon Sep 17 00:00:00 2001
From: Norihiro Watanabe <norihiro.watanabe@aist.go.jp>
Date: Mon, 10 Jul 2017 09:18:46 +0900
Subject: [PATCH] [PCS/LIE/HM] use LinAlg::setLocalAccessibleVector()
---
.../HydroMechanics/HydroMechanicsProcess.cpp | 17 ++---------------
1 file changed, 2 insertions(+), 15 deletions(-)
diff --git a/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcess.cpp b/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcess.cpp
index e20e945249f..814aa30eab8 100644
--- a/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcess.cpp
+++ b/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcess.cpp
@@ -373,15 +373,7 @@ void HydroMechanicsProcess<GlobalDim>::computeSecondaryVariableConcrete(
}
}
-#ifdef USE_PETSC
- // TODO It is also possible directly to copy the data for single process
- // variable to a mesh property. It needs a vector of global indices and
- // some PETSc magic to do so.
- std::vector<double> x_copy(x.getLocalSize() + x.getGhostSize());
-#else
- std::vector<double> x_copy(x.size());
-#endif
- x.copyValues(x_copy);
+ MathLib::LinAlg::setLocalAccessibleVector(x);
ProcessVariable& pv_g = this->getProcessVariables()[g_variable_id];
auto& mesh_prop_g = pv_g.getOrCreateMeshProperty();
@@ -400,13 +392,8 @@ void HydroMechanicsProcess<GlobalDim>::computeSecondaryVariableConcrete(
mesh_id, MeshLib::MeshItemType::Node, node->getID());
auto const global_component_id = g_global_component_offset + component_id;
- auto const index =
- _local_to_global_index_map->getLocalIndex(
- l, global_component_id, x.getRangeBegin(),
- x.getRangeEnd());
-
mesh_prop_g[node->getID() * num_comp + component_id] =
- x_copy[index];
+ x[global_component_id];
}
}
}
--
GitLab