From 557e7ebda848cfbeaa3c9bd6db45748187dfde9e Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Tue, 8 Apr 2014 11:59:41 +0100
Subject: [PATCH] ENH: forces FO - porosity calc updated for mesh changes
---
src/postProcessing/functionObjects/forces/forces/forces.C | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/postProcessing/functionObjects/forces/forces/forces.C b/src/postProcessing/functionObjects/forces/forces/forces.C
index 488a71a8818..539e14b3093 100644
--- a/src/postProcessing/functionObjects/forces/forces/forces.C
+++ b/src/postProcessing/functionObjects/forces/forces/forces.C
@@ -922,7 +922,8 @@ void Foam::forces::calcForcesMoment()
forAllConstIter(HashTable<const porosityModel*>, models, iter)
{
- const porosityModel& pm = *iter();
+ // non-const access required if mesh is changing
+ porosityModel& pm = const_cast<porosityModel&>(*iter());
vectorField fPTot(pm.force(U, rho, mu));
--
GitLab