Feature/gifs

README.md

@@ -1,19 +1,14 @@
 # ODE Test Problems
 
 
+ <img src="images/logo.png" width="200" alt='ODE Test problems Logo' /> <img src="images/animations/nbody-OuterSolarSystem.gif" width="200" alt='KPR problem' loop=infinite  /> <img src="images/animations/lorenz63-Canonical.gif" width="200" alt='Pendulum problem' loop=infinite /> 
+
+ <img src="images/animations/bouncingball-RandomTerrain.gif" width="200" alt='Bouncing ball problem' loop=infinite  />  <img src="images/animations/pendulum-Canonical.gif" width="200" alt='Lorenz63 problem' loop=infinite  />  <img src="images/animations/kpr-Canonical.gif" width="200" alt='Lorenz63 problem' loop=infinite  />
+
+
+
+ ODE Test Problems (<code>OTP</code>) is an object-oriented Matlab/Octave package offering a broad range of initial value problems in the form of ordinary and differential-algebraic equations that can be used to test numerical methods such as time integration or data assimilation.  It includes problems that are linear and nonlinear, homogeneous and nonhomogeneous, autonomous and nonautonomous, scalar and high-dimensional, stiff and nonstiff, and chaotic and nonchaotic.  Many are real-world problems in fields such as chemistry, astrophysics, meteorology, and electrical engineering.  <code>OTP</code>  also supports partitioned ODEs for testing split, multirate, and other multimethods.  Functions for plotting solutions and creating movies are available for all problems, and exact solutions are included when available. <code>OTP</code> is designed for ease of use &#8212;meaning that working with and modifying problems is simple and intuitive. <code>OTP</code> is actively under development. We are currently writing full documentation in order to release version <code>1.0.0</code>. </p>
 
-<table border="0";>
-<tbody>
-    <tr>
-        <td> <img src="images/logo.png" width="5000" alt='ODE Test problems Logo' /></td>
-        <td>
-            <p> ODE Test Problems (<code>OTP</code>) is an object-oriented Matlab/Octave package offering a broad range of initial value problems in the form of ordinary and differential-algebraic equations that can be used to test numerical methods such as time integration or data assimilation.  It includes problems that are linear and nonlinear, homogeneous and nonhomogeneous, autonomous and nonautonomous, scalar and high-dimensional, stiff and nonstiff, and chaotic and nonchaotic.  Many are real-world problems in fields such as chemistry, astrophysics, meteorology, and electrical engineering.  <code>OTP</code>  also supports partitioned ODEs for testing split, multirate, and other multimethods.  Functions for plotting solutions and creating movies are available for all problems, and exact solutions are included when available. <code>OTP</code> is designed for ease of use &#8212;meaning that working with and modifying problems is simple and intuitive. 
-            <code>OTP</code> is actively under development. We
-            are currently writing full documentation in order to release version <code>1.0.0</code>. </p>
-        </td>
-    </tr>
-</tbody>
-</table>
 
 ![Version](https://img.shields.io/github/v/release/ComputationalScienceLaboratory/ODE-Test-Problems?label=Version)
 ![GitHub](https://img.shields.io/github/license/ComputationalScienceLaboratory/ODE-Test-Problems)
@@ -66,7 +61,7 @@ New feature requests, and bug reports can be made through
 We also accept pull requests that adhere to our
 [contributing guide](docs/contrib.rst). 
 
+
 An interactive tutorial on the main features of `OTP` is available in a Jupyter notebook [in the repository](/notebooks)
 
 Further discussion of the main contributions of this open-source package can be found in our [paper](/paper/paper.md). 
-

images/makeMovie.m

@@ -0,0 +1,68 @@
+close all
+
+problemName  = 'kpr';
+presetName   = 'Canonical';
+outputFormat = 'gif';
+targetTime   =  10;
+size         = [400,400];
+
+
+set(0,'DefaultAxesXMinorTick','on','DefaultAxesYMinorTick','on')
+set(0,'DefaultAxesLineWidth', 2,'DefaultAxesFontName','Times',...
+    'DefaultAxesFontSize',20,'DefaultAxesBox','on')
+
+ffmpegPath   = {'/opt/homebrew/bin/ffmpeg'};
+
+
+mkdir('animations')
+filename = strcat('animations/', problemName,'-', presetName);
+
+
+string   = strcat(['otp.', problemName,'.presets.', presetName] );
+problem  = eval(string);
+problem.TimeSpan = [0,30];
+
+% 
+
+sol      = problem.solve('RelTol', 1e-8);
+mov      = problem.movie(sol, 'Save', filename, ...
+                         'TargetDuration',targetTime, ...
+                         'Size', size, ...
+                         'FrameRate', 15, ...
+                         'view', [0,90]);
+
+
+% convert vido format if neccessary 
+switch outputFormat
+    case 'avi'
+
+    case 'gif'
+
+        args     = strcat(['-i ' ...
+        filename,'.avi ' ...
+        ' -vf "fps=10,scale=400:-1:flags=lanczos,split[s0][s1];[s0]palettegen[p];[s1][p]paletteuse" -loop 0' ...
+        ' '...
+        filename,'.gif']);
+
+
+    case 'webm'
+
+        args     = strcat(['-i ' ...
+        filename,'.avi ' ...
+        '-c:v libvpx -crf 10 -b:v 1M -c:a libvorbis '...
+        filename,'.webm']);
+
+end
+
+command = strjoin(strcat([ffmpegPath,' ' , args]));
+system(command);        
+
+
+
+
+
+
+
+
+
+

toolbox/+otp/+utils/+movie/Movie.m

@@ -66,7 +66,7 @@ function record(obj, t, y)
 
             fig = figure;
             if ~isempty(obj.Config.Size)
-                fig.Position = [0; 0: obj.Config.Size(:)];
+                fig.Position = [0; 0; obj.Config.Size(:)];
             end
 
             obj.init(fig, state);