tadded information to maxNormResNS nan test - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
(HTM) git clone git://src.adamsgaard.dk/sphere
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) LICENSE
---
(DIR) commit a1210156b64db147b3827d400a4f0ce3dd2059cb
(DIR) parent d843fd95007da8e532a444c0e3aca8df9010fd6f
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Fri, 23 May 2014 11:18:43 +0200
added information to maxNormResNS nan test
Diffstat:
M src/debug.h | 4 ++--
M src/navierstokes.cpp | 2 ++
2 files changed, 4 insertions(+), 2 deletions(-)
---
(DIR) diff --git a/src/debug.h b/src/debug.h
t@@ -51,7 +51,7 @@ const int conv_log_interval = 10;
// simulation of particulate systems: Theoretical developments".
// SET_2 corresponds approximately to Model A in Zhu et al. 2007.
// Choose exactly one.
-//#define SET_1
-#define SET_2
+#define SET_1
+//#define SET_2
#endif
(DIR) diff --git a/src/navierstokes.cpp b/src/navierstokes.cpp
t@@ -277,6 +277,8 @@ double DEM::maxNormResNS()
<< x << "," << y << "," << z << std::endl;
std::cerr << "\tt = " << time.current << ", iter = "
<< int(time.current/time.dt) << std::endl;
+ std::cerr << "This often happens if the system has become "
+ "unstable." << std::endl;
exit(1);
}
if (norm_res > max_norm_res)