Commit e5f110c8 authored by Martin Schneider's avatar Martin Schneider
Browse files

[wmpfa][nlmpfa] Update eps in flux calculation

parent 73d68e77
......@@ -250,7 +250,7 @@ public:
fluxIJ += (betaIJ - beta)*(sI - sJ);
fluxJI += (betaJI - beta)*(sJ - sI);
auto calcAbs = [](Scalar a){ return std::abs(a) + 1e-8; };
auto calcAbs = [](Scalar a){ return std::abs(a) + 1e-30; };
Scalar wIJ = (calcAbs(fluxJI))/(calcAbs(fluxIJ) + calcAbs(fluxJI));
Scalar wJI = (calcAbs(fluxIJ))/(calcAbs(fluxIJ) + calcAbs(fluxJI));
......
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