Added plot

Author: brownbaerchen

pySDC/implementations/convergence_controller_classes/step_size_limiter.py

@@ -192,12 +192,7 @@ def _round_step_size(dt, fac, digits):
 
     def get_new_step_size(self, controller, S, **kwargs):
         """
-        Enforce an upper and lower limit to the step size here.
-        Be aware that this is only tested when a new step size has been determined. That means if you set an initial
-        value for the step size outside of the limits, and you don't do any further step size control, that value will
-        go through.
-        Also, the final step is adjusted such that we reach Tend as best as possible, which might give step sizes below
-        the lower limit set here.
+        Round step size here
 
         Args:
             controller (pySDC.Controller): The controller

pySDC/projects/GPU/analysis_scripts/plot_RBC_matrix.py

@@ -61,6 +61,36 @@ def plot_ultraspherical():
     plt.show()
 
 
+def plot_DCT():
+    fig, axs = plt.subplots(1, 3, figsize=figsize_by_journal('TUHH_thesis', 1, 0.28), sharey=True)
+
+    N = 8
+    color = 'black'
+
+    x = np.linspace(0, 3, N)
+    y = x**3 - 4 * x**2
+    axs[0].plot(y, marker='o', color=color)
+
+    y_m = np.append(y, y[::-1])
+    axs[1].scatter(np.arange(2 * N)[::2], y_m[::2], marker='<', color=color)
+    axs[1].scatter(np.arange(2 * N)[1::2], y_m[1::2], marker='>', color=color)
+    axs[1].plot(np.arange(2 * N), y_m, color=color)
+
+    v = y_m[::2]
+    axs[2].plot(np.arange(N), v, color=color, marker='x')
+
+    axs[0].set_title('original')
+    axs[1].set_title('mirrored')
+    axs[2].set_title('periodically reordered')
+
+    for ax in axs:
+        # ax.set_xlabel(r'$n$')
+        ax.set_yticks([])
+    fig.savefig('plots/DCT_via_FFT.pdf', bbox_inches='tight', dpi=300)
+
+
 if __name__ == '__main__':
     setup_mpl()
-    plot_ultraspherical()
+    plot_DCT()
+    # plot_ultraspherical()
+    plt.show()