Merge pull request #3 from jchristopherson/v1.2.0

Bug fix in ODE_AUTO
jchristopherson · Jul 7, 2022 · 0ac9900 · 0ac9900
2 parents 3c90c31 + fd06ec2
commit 0ac9900
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 1 deletion.
diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -4,7 +4,7 @@ project(integral C CXX Fortran)
 
 # Define version information
 set(INTEGRAL_MAJOR_VERSION 1)
-set(INTEGRAL_MINOR_VERSION 1)
+set(INTEGRAL_MINOR_VERSION 2)
 set(INTEGRAL_PATCH_VERSION 0)
 set(INTEGRAL_VERSION ${INTEGRAL_MAJOR_VERSION}.${INTEGRAL_MINOR_VERSION}.${INTEGRAL_PATCH_VERSION})
 

diff --git a/src/integral_ode_auto.f90 b/src/integral_ode_auto.f90
@@ -188,19 +188,27 @@ module subroutine oa_init_workspace(this, liw, lrw, ncnsts, err)
             if (size(this%m_iwork) < liw) then
                 deallocate(this%m_iwork)
                 allocate(this%m_iwork(liw), stat = flag)
+                this%m_iwork = 0
+                this%m_iwork(8) = 12
+                this%m_iwork(9) = 5
             end if
         else
             allocate(this%m_iwork(liw), stat = flag)
+            this%m_iwork = 0
+            this%m_iwork(8) = 12
+            this%m_iwork(9) = 5
         end if
 
         if (flag == 0) then
             if (allocated(this%m_rwork)) then
                 if (size(this%m_rwork) < lrw) then
                     deallocate(this%m_rwork)
                     allocate(this%m_rwork(lrw), stat = flag)
+                    this%m_rwork = 0.0d0
                 end if
             else
                 allocate(this%m_rwork(lrw), stat = flag)
+                this%m_rwork = 0.0d0
             end if
         end if