UpsilonNumworks/Upsilon
1
0
Fork 0
mirror of https://github.com/UpsilonNumworks/Upsilon.git synced 2026-03-18 21:30:38 +01:00
Code Issues Packages Projects Releases Wiki Activity

[poincare] ApproximationHelper: fix precision threshold of

Browse Source 
NeglectRealOrImaginaryPartIfNeglectable
This commit is contained in:
Émilie Feral
2020-03-06 10:35:14 +01:00
committed by LeaNumworks
parent d431642e60
commit 879c1e3bf1
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
poincare/src/approximation_helper.cpp
View File

@@ -11,7 +11,7 @@ extern "C" {
namespace Poincare {
template <typename T> bool isNegligeable(T x, T precision, T norm1, T norm2) {
return x <= precision && x/norm1 <= precision && x/norm2 <= precision;
return x <= 10.0*precision && x/norm1 <= precision && x/norm2 <= precision;
}
template < typename T> T minimalNonNullMagnitudeOfParts(std::complex<T> c) {
@@ -47,7 +47,7 @@ template <typename T> std::complex<T> ApproximationHelper::NeglectRealOrImaginar
* small as 1E-308 (in double) and most results still be correct. */
T magnitude1 = minimalNonNullMagnitudeOfParts(input1);
T magnitude2 = minimalNonNullMagnitudeOfParts(input2);
T precision = 100.0*Expression::Epsilon<T>();
T precision = 10.0*Expression::Epsilon<T>();
if (isNegligeable(std::fabs(result.imag()), precision, magnitude1, magnitude2)) {
result.imag(0);
}