Commit 4364d986 authored by wenqing's avatar wenqing
Browse files

[ProcessLib/TM] Removed the global vector of _previous_T and its associated function

parent 01519054
......@@ -353,8 +353,6 @@ void ThermoMechanicsProcess<DisplacementDim>::
dof_tables.emplace_back(
*_local_to_global_index_map_single_component);
}
setCoupledSolutionsOfPreviousTimeStep();
}
ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
......@@ -410,21 +408,6 @@ void ThermoMechanicsProcess<DisplacementDim>::preTimestepConcreteProcess(
*x[process_id], t, dt);
return;
}
// For the staggered scheme.
if (!_previous_T)
{
_previous_T = MathLib::MatrixVectorTraits<GlobalVector>::newInstance(
*x[process_id]);
}
else
{
auto& x0 = *_previous_T;
MathLib::LinAlg::copy(*x[process_id], x0);
}
auto& x0 = *_previous_T;
MathLib::LinAlg::setLocalAccessibleVector(x0);
}
template <int DisplacementDim>
......@@ -447,14 +430,6 @@ void ThermoMechanicsProcess<DisplacementDim>::postTimestepConcreteProcess(
t, dt);
}
template <int DisplacementDim>
void ThermoMechanicsProcess<
DisplacementDim>::setCoupledSolutionsOfPreviousTimeStep()
{
_coupled_solutions->coupled_xs_t0.resize(1);
_coupled_solutions->coupled_xs_t0[0] = _previous_T.get();
}
template <int DisplacementDim>
NumLib::LocalToGlobalIndexMap const&
ThermoMechanicsProcess<DisplacementDim>::getDOFTable(const int process_id) const
......
......@@ -104,13 +104,6 @@ private:
MeshLib::PropertyVector<double>* _nodal_forces = nullptr;
MeshLib::PropertyVector<double>* _heat_flux = nullptr;
/// Temperature of the previous time step for staggered scheme.
std::unique_ptr<GlobalVector> _previous_T;
/// Set the increment solutions of the present time step to the coupled
/// term.
void setCoupledSolutionsOfPreviousTimeStep();
};
extern template class ThermoMechanicsProcess<2>;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment