diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2p.hh b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2p.hh index d2fd8ff2e6c4b3328b1169a4bf0331f85bc3a23a..d2c6fe00f3d820d418157a114d5d13ff1cdebe56 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 {