Update to version 2.1

munstermonster · Aug 11, 2019 · 79c2afa · 79c2afa
1 parent 1a4b7be
commit 79c2afa
Show file tree

Hide file tree

Showing 56 changed files with 3,084 additions and 596 deletions.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -6,6 +6,7 @@ CHANGELOG
 
 | Contents               |
 | :--------------------- |
+| [cuSten 2.1](#cuSten-21) |
 | [cuSten 2.0](#cuSten-20) |
 | [cuSten 1.5](#cuSten-15) |
 | [cuSten 1.4](#cuSten-14) |
@@ -18,6 +19,15 @@ CHANGELOG
 
 </details>
 
+cuSten 2.1
+---------
+
+* Updated naming convention to consistent cuSten for all functions and types
+* Replaced the 0 and 1 for loading/unloading data with MACROS HOST and DEVICE to specify where the programmer desires the data to be stored
+* Restructured cuPentCahnADI to use a better guess for the initial n - 1 time step, scheme now more stable.
+* Added serial example of the same for GPU versus CPU timing benchmarks.
+* doxygen updated to reflect above changes
+
 cuSten 2.0
 ---------
 

diff --git a/README.md b/README.md
@@ -29,7 +29,7 @@ Run the command make doc in the cuSten directory.
 # Where is cuSten used
 cuSten is currently used in the following academic papers:
 
-- cuPentBatch -- A batched pentadiagonal solver for NVIDIA GPUs (preprint: https://arxiv.org/abs/1807.07382)
+- cuPentBatch -- A batched pentadiagonal solver for NVIDIA GPUs (preprint: https://arxiv.org/abs/1807.07382, published article: https://www.sciencedirect.com/science/article/pii/S0010465519300979)
 - cuSten -- CUDA Finite Difference and Stencil Library (preprint: https://arxiv.org/abs/1902.09931)
 
 

diff --git a/cuPentCahnADI/movie.py b/cuPentCahnADI/movie.py
@@ -25,7 +25,7 @@
 
 output = 'output/'
 
-nx = 1024
+nx = 512
 
 # Get filenames and times
 hdfiles = os.listdir(output)

diff --git a/cuPentCahnADI/plotting.py b/cuPentCahnADI/plotting.py
@@ -24,7 +24,7 @@
 
 output = 'output/'
 
-nx = 1024
+nx = 512
 
 s_t = []
 k1 = []
@@ -60,8 +60,8 @@
 	a_group_key = f.keys()[0]
 	data = np.array(f[a_group_key])
 
-	x = np.linspace(0, 2.0 * math.pi, 2 ** 10)
-	y = np.linspace(0, 2.0 * math.pi, 2 ** 10)
+	x = np.linspace(0, 2.0 * math.pi, nx)
+	y = np.linspace(0, 2.0 * math.pi, nx)
 	avg = (1.0 / ((2.0 * math.pi) ** 2)) * simps(simps(np.square(data), y), x)
 	s_t.append(1.0 / (1.0 - avg))
 
@@ -78,7 +78,7 @@
 plt.loglog(t, k1, label = '1 / k_1')
 plt.legend(loc='upper left')
 plt.xlabel('t')
-plt.savefig('analysis.png')
+plt.savefig('analysis.png', dpi = 300)
 
 
 

diff --git a/cuPentCahnADI/src/CHANGELOG.md b/cuPentCahnADI/src/CHANGELOG.md