The output modules have been completely revised, as to ensure consistent output results across all simulation features:
Output files now have semantic names (output_ox, output_vmax), etc.
The time series output has been generalised, so that the output produced at the IC and IOT locations follow the same structure and are treated identically. Time series output from one location is stored in a file with the location index identifier (e.g. output_ot_1020 for location index 1020).
To reduce overhead on the time series output, the output from the location at which V = max(V) is stored in a separate file (ot_vmax).
Vmax and IOT output are now supported for MPI parallel simulations (#44)
Octave support was dropped, which could no longer be guaranteed (#47)
For 3D simulations, a separate snapshot output resolution was introduced for along-strike (NXOUT) and along-dip (NWOUT) directions (#48)
A user variable for the MPI path was added to the PyQDYN wrapper
For the time being, binary output has been disabled. The generalised output structure expedites the implementation of HDF5 binary output, which will be implemented in the future (#51).
3D tutorial notebook (#52)
To introduce new users into 3D simulations, a Jupyter notebook tutorial has been added which illustrates how to set-up a basic 3D simulation (similar to the single asperity tutorial).
Bug fixes since version 2.2.0
(#37) Fixed a small bug in pyqdyn.py which caused the grid nodes to be slightly offset (by half the grid spacing).
(#38) Previously slip was computed by forward Euler integration of slip rate, which has a (much) lower-order accuracy than the implemented ODE solvers. To prevent the accumulation of discrepancies in the slip output for long-term simulations (see #17), the slip computation was incorporated within the ODE solver routine.
(#38) The low precision output of the simulation time (in time series output) caused round-off issues for long-term simulations. The output precision has been increased for time quantities.
(#53) The Runge-Kutta solver is now compatible with MPI parallel simulations
Fixed various problems in the PyQDYN wrapper with running 3D simulations.
Fixed a NumPy (v1.18+) issue interpreting floating point numbers as integers, which would throw an error while running the test suite. As a result, the TravisCI automated tests would fail too.
Added a travis_wait command for the test suite to .travis.yml to prevent time-out during TravisCI testing.
What's new in version 2.3.0?
Better simulation output (#42, #45)
The output modules have been completely revised, as to ensure consistent output results across all simulation features:
output_ox,output_vmax), etc.ICandIOTlocations follow the same structure and are treated identically. Time series output from one location is stored in a file with the location index identifier (e.g.output_ot_1020for location index1020).ot_vmax).NXOUT) and along-dip (NWOUT) directions (#48)For the time being, binary output has been disabled. The generalised output structure expedites the implementation of HDF5 binary output, which will be implemented in the future (#51).
3D tutorial notebook (#52)
To introduce new users into 3D simulations, a Jupyter notebook tutorial has been added which illustrates how to set-up a basic 3D simulation (similar to the single asperity tutorial).
Bug fixes since version 2.2.0
pyqdyn.pywhich caused the grid nodes to be slightly offset (by half the grid spacing).travis_waitcommand for the test suite to.travis.ymlto prevent time-out during TravisCI testing.