From f0e610c178cf44ec4e647d68195ee7f7cc62d4c0 Mon Sep 17 00:00:00 2001
From: Paul Dechamps <paul.dechamps@uliege.be>
Date: Thu, 25 Apr 2024 12:17:11 +0200
Subject: [PATCH] (fix) Modified rae validation case for PARDISO linear solver
---
blast/validation/raeValidation.py | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/blast/validation/raeValidation.py b/blast/validation/raeValidation.py
index c99406c..855ccbf 100644
--- a/blast/validation/raeValidation.py
+++ b/blast/validation/raeValidation.py
@@ -104,7 +104,7 @@ def main():
else:
n = np.log10(1-(1-gr)*icfg['Pars']['xLgt']/icfg['Pars']['msLe'])/np.log10(gr)
icfg['Pars']['msF'] = icfg['Pars']['msLe']*gr**(n-1)
- print('msF =', icfg['Pars']['msF'])
+
tms['pre'].start()
coupler, isol, vsol = viscUtils.initBlast(icfg, vcfg)
tms['pre'].stop()
@@ -136,7 +136,10 @@ def main():
tests.add(CTest('Cd wake', vsol.Cdt, 0.0093, 1e-3, forceabs=True))
tests.add(CTest('Cd integral', isol.getCd() + vsol.Cdf, 0.0138, 1e-3, forceabs=True))
tests.add(CTest('Cdf', vsol.Cdf, 0.0069, 1e-3, forceabs=True))
- tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 42, 0, forceabs=True))
+ if icfg['LSolver'] == 'PARDISO':
+ tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 43, 0, forceabs=True))
+ elif icfg['LSolver'] == 'SparseLu':
+ tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 42, 0, forceabs=True))
tests.run()
expResults = np.loadtxt(os.path.dirname(os.path.dirname(os.path.abspath(__file__))) + '/models/references/rae2822_AR138_case6.dat')
--
GitLab