From b7c8b442bedb814bd9e720c5992a57d318522534 Mon Sep 17 00:00:00 2001 From: Christoph Lehmann <christoph.lehmann@ufz.de> Date: Tue, 25 Jul 2017 20:37:27 +0200 Subject: [PATCH] [PL] adapted LIE to new getIntPt...() signature --- .../SmallDeformationLocalAssemblerFracture.h | 18 +++++++++++ .../SmallDeformationLocalAssemblerInterface.h | 18 +++++++++++ .../SmallDeformationLocalAssemblerMatrix.h | 18 +++++++++++ ...ormationLocalAssemblerMatrixNearFracture.h | 18 +++++++++++ .../SmallDeformationProcess.cpp | 32 +++++++++---------- 5 files changed, 88 insertions(+), 16 deletions(-) diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture.h index cf2058304a4..6a4a398b18b 100644 --- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture.h +++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture.h @@ -158,30 +158,45 @@ public: } std::vector<double> const& getIntPtEpsilonXX( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 0); } std::vector<double> const& getIntPtEpsilonYY( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 1); } std::vector<double> const& getIntPtEpsilonZZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 2); } std::vector<double> const& getIntPtEpsilonXY( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 3); } std::vector<double> const& getIntPtEpsilonXZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { assert(DisplacementDim == 3); @@ -189,6 +204,9 @@ public: } std::vector<double> const& getIntPtEpsilonYZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { assert(DisplacementDim == 3); diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerInterface.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerInterface.h index 7d5e91fcbb0..907b70fa56c 100644 --- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerInterface.h +++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerInterface.h @@ -119,21 +119,39 @@ public: std::vector<double>& cache) const = 0; virtual std::vector<double> const& getIntPtEpsilonXX( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const = 0; virtual std::vector<double> const& getIntPtEpsilonYY( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const = 0; virtual std::vector<double> const& getIntPtEpsilonZZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const = 0; virtual std::vector<double> const& getIntPtEpsilonXY( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const = 0; virtual std::vector<double> const& getIntPtEpsilonXZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const = 0; virtual std::vector<double> const& getIntPtEpsilonYZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const = 0; private: diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix.h index 132982937db..b8eb7ff4cb9 100644 --- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix.h +++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix.h @@ -158,30 +158,45 @@ public: } std::vector<double> const& getIntPtEpsilonXX( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 0); } std::vector<double> const& getIntPtEpsilonYY( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 1); } std::vector<double> const& getIntPtEpsilonZZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 2); } std::vector<double> const& getIntPtEpsilonXY( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 3); } std::vector<double> const& getIntPtEpsilonXZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { assert(DisplacementDim == 3); @@ -189,6 +204,9 @@ public: } std::vector<double> const& getIntPtEpsilonYZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { assert(DisplacementDim == 3); diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrixNearFracture.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrixNearFracture.h index c78937f96c1..a91d8721565 100644 --- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrixNearFracture.h +++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrixNearFracture.h @@ -160,30 +160,45 @@ public: } std::vector<double> const& getIntPtEpsilonXX( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 0); } std::vector<double> const& getIntPtEpsilonYY( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 1); } std::vector<double> const& getIntPtEpsilonZZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 2); } std::vector<double> const& getIntPtEpsilonXY( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { return getIntPtEpsilon(cache, 3); } std::vector<double> const& getIntPtEpsilonXZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { assert(DisplacementDim == 3); @@ -191,6 +206,9 @@ public: } std::vector<double> const& getIntPtEpsilonYZ( + const double /*t*/, + GlobalVector const& /*current_solution*/, + NumLib::LocalToGlobalIndexMap const& /*dof_table*/, std::vector<double>& cache) const override { assert(DisplacementDim == 3); diff --git a/ProcessLib/LIE/SmallDeformation/SmallDeformationProcess.cpp b/ProcessLib/LIE/SmallDeformation/SmallDeformationProcess.cpp index 6eb8c4598b5..1035c8e2214 100644 --- a/ProcessLib/LIE/SmallDeformation/SmallDeformationProcess.cpp +++ b/ProcessLib/LIE/SmallDeformation/SmallDeformationProcess.cpp @@ -212,54 +212,54 @@ void SmallDeformationProcess<DisplacementDim>::initializeConcreteProcess( if (DisplacementDim == 3) { Base::_secondary_variables.addSecondaryVariable( - "sigma_xz", 1, + "sigma_xz", makeExtrapolator( - getExtrapolator(), _local_assemblers, + 1, getExtrapolator(), _local_assemblers, &SmallDeformationLocalAssemblerInterface::getIntPtSigmaXZ)); Base::_secondary_variables.addSecondaryVariable( - "sigma_yz", 1, + "sigma_yz", makeExtrapolator( - getExtrapolator(), _local_assemblers, + 1, getExtrapolator(), _local_assemblers, &SmallDeformationLocalAssemblerInterface::getIntPtSigmaYZ)); } Base::_secondary_variables.addSecondaryVariable( - "epsilon_xx", 1, + "epsilon_xx", makeExtrapolator( - getExtrapolator(), _local_assemblers, + 1, getExtrapolator(), _local_assemblers, &SmallDeformationLocalAssemblerInterface::getIntPtEpsilonXX)); Base::_secondary_variables.addSecondaryVariable( - "epsilon_yy", 1, + "epsilon_yy", makeExtrapolator( - getExtrapolator(), _local_assemblers, + 1, getExtrapolator(), _local_assemblers, &SmallDeformationLocalAssemblerInterface::getIntPtEpsilonYY)); Base::_secondary_variables.addSecondaryVariable( - "epsilon_zz", 1, + "epsilon_zz", makeExtrapolator( - getExtrapolator(), _local_assemblers, + 1, getExtrapolator(), _local_assemblers, &SmallDeformationLocalAssemblerInterface::getIntPtEpsilonZZ)); Base::_secondary_variables.addSecondaryVariable( - "epsilon_xy", 1, + "epsilon_xy", makeExtrapolator( - getExtrapolator(), _local_assemblers, + 1, getExtrapolator(), _local_assemblers, &SmallDeformationLocalAssemblerInterface::getIntPtEpsilonXY)); if (DisplacementDim == 3) { Base::_secondary_variables.addSecondaryVariable( - "epsilon_xz", 1, + "epsilon_xz", makeExtrapolator( - getExtrapolator(), _local_assemblers, + 1, getExtrapolator(), _local_assemblers, &SmallDeformationLocalAssemblerInterface::getIntPtEpsilonXZ)); Base::_secondary_variables.addSecondaryVariable( - "epsilon_yz", 1, + "epsilon_yz", makeExtrapolator( - getExtrapolator(), _local_assemblers, + 1, getExtrapolator(), _local_assemblers, &SmallDeformationLocalAssemblerInterface::getIntPtEpsilonYZ)); } -- GitLab