From f92078c518eac9f357a735e555a659180cdbcb4a Mon Sep 17 00:00:00 2001
From: Markus Wolff <markus.wolff@twt-gmbh.de>
Date: Mon, 18 Mar 2013 08:18:13 +0000
Subject: [PATCH] changed criterion for choice of transmissibility to newer
version
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@10429 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
.../diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2p.hh | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2p.hh b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2p.hh
index d2fd8ff2e6..d2c6fe00f3 100644
--- a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2p.hh
+++ b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2p.hh
@@ -2242,9 +2242,11 @@ bool FVMPFAL2PFABoundPressure2P<TypeTag>::calculateTransmissibility(
D += B.leftmultiply(C.rightmultiply(A));
Dune::FieldMatrix<Scalar,dim,2*dim-dim+1> TL1(D);
+ Scalar sR = std::abs(TR1[1][2] - TR1[1][0]);
+ Scalar sL = std::abs(TL1[1][0] - TL1[1][2]);
// 3.decide which triangle (which transmissibility coefficients) to use
- if (std::abs(TR1[1][0]) < std::abs(TL1[1][0]))
+ if (sR <= sL)
{
transmissibility = TR1;
// if (isnan(transmissibility.frobenius_norm()))
@@ -2273,7 +2275,7 @@ bool FVMPFAL2PFABoundPressure2P<TypeTag>::calculateTransmissibility(
// std::cout<<"\n";
// }
// std::cout<<"transmissibility = "<<transmissibility<<"\n";
- return true;
+// return true;
}
else
{
--
GitLab