diff --git a/auto_diff/private/auto_diff_real_15var_order1_module.f90 b/auto_diff/private/auto_diff_real_15var_order1_module.f90
index abd3a253a..06da648a3 100644
--- a/auto_diff/private/auto_diff_real_15var_order1_module.f90
+++ b/auto_diff/private/auto_diff_real_15var_order1_module.f90
@@ -306,14 +306,14 @@ subroutine assign_from_real_dp(this, other)
       type(auto_diff_real_15var_order1), intent(out) :: this
       real(dp), intent(in) :: other
       this%val = other
-      this%d1Array = 0_dp
+      this%d1Array = 0.0_dp
    end subroutine assign_from_real_dp
    
    subroutine assign_from_int(this, other)
       type(auto_diff_real_15var_order1), intent(out) :: this
       integer, intent(in) :: other
       this%val = other
-      this%d1Array = 0_dp
+      this%d1Array = 0.0_dp
    end subroutine assign_from_int
    
    function equal_self(this, other) result(z)
diff --git a/auto_diff/private/auto_diff_real_1var_order1_module.f90 b/auto_diff/private/auto_diff_real_1var_order1_module.f90
index 57b87a844..9b86f2efe 100644
--- a/auto_diff/private/auto_diff_real_1var_order1_module.f90
+++ b/auto_diff/private/auto_diff_real_1var_order1_module.f90
@@ -372,14 +372,14 @@ subroutine assign_from_real_dp(this, other)
       type(auto_diff_real_1var_order1), intent(out) :: this
       real(dp), intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
+      this%d1val1 = 0.0_dp
    end subroutine assign_from_real_dp
    
    subroutine assign_from_int(this, other)
       type(auto_diff_real_1var_order1), intent(out) :: this
       integer, intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
+      this%d1val1 = 0.0_dp
    end subroutine assign_from_int
    
    function equal_self(this, other) result(z)
@@ -1294,7 +1294,7 @@ function differentiate_auto_diff_real_1var_order1_1(this) result(derivative)
       type(auto_diff_real_1var_order1), intent(in) :: this
       type(auto_diff_real_1var_order1) :: derivative
       derivative%val = this%d1val1
-      derivative%d1val1 = 0_dp
+      derivative%d1val1 = 0.0_dp
    end function differentiate_auto_diff_real_1var_order1_1
    
 end module auto_diff_real_1var_order1_module
\ No newline at end of file
diff --git a/auto_diff/private/auto_diff_real_2var_order1_module.f90 b/auto_diff/private/auto_diff_real_2var_order1_module.f90
index f83f2a920..257939d9f 100644
--- a/auto_diff/private/auto_diff_real_2var_order1_module.f90
+++ b/auto_diff/private/auto_diff_real_2var_order1_module.f90
@@ -379,16 +379,16 @@ subroutine assign_from_real_dp(this, other)
       type(auto_diff_real_2var_order1), intent(out) :: this
       real(dp), intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1val2 = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1val2 = 0.0_dp
    end subroutine assign_from_real_dp
    
    subroutine assign_from_int(this, other)
       type(auto_diff_real_2var_order1), intent(out) :: this
       integer, intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1val2 = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1val2 = 0.0_dp
    end subroutine assign_from_int
    
    function equal_self(this, other) result(z)
@@ -1511,16 +1511,16 @@ function differentiate_auto_diff_real_2var_order1_1(this) result(derivative)
       type(auto_diff_real_2var_order1), intent(in) :: this
       type(auto_diff_real_2var_order1) :: derivative
       derivative%val = this%d1val1
-      derivative%d1val1 = 0_dp
-      derivative%d1val2 = 0_dp
+      derivative%d1val1 = 0.0_dp
+      derivative%d1val2 = 0.0_dp
    end function differentiate_auto_diff_real_2var_order1_1
    
    function differentiate_auto_diff_real_2var_order1_2(this) result(derivative)
       type(auto_diff_real_2var_order1), intent(in) :: this
       type(auto_diff_real_2var_order1) :: derivative
       derivative%val = this%d1val2
-      derivative%d1val1 = 0_dp
-      derivative%d1val2 = 0_dp
+      derivative%d1val1 = 0.0_dp
+      derivative%d1val2 = 0.0_dp
    end function differentiate_auto_diff_real_2var_order1_2
    
 end module auto_diff_real_2var_order1_module
\ No newline at end of file
diff --git a/auto_diff/private/auto_diff_real_2var_order2_module.f90 b/auto_diff/private/auto_diff_real_2var_order2_module.f90
index 58534f842..ebd151bf3 100644
--- a/auto_diff/private/auto_diff_real_2var_order2_module.f90
+++ b/auto_diff/private/auto_diff_real_2var_order2_module.f90
@@ -385,22 +385,22 @@ subroutine assign_from_real_dp(this, other)
       type(auto_diff_real_2var_order2), intent(out) :: this
       real(dp), intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1val2 = 0_dp
-      this%d2val1 = 0_dp
-      this%d1val1_d1val2 = 0_dp
-      this%d2val2 = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1val2 = 0.0_dp
+      this%d2val1 = 0.0_dp
+      this%d1val1_d1val2 = 0.0_dp
+      this%d2val2 = 0.0_dp
    end subroutine assign_from_real_dp
    
    subroutine assign_from_int(this, other)
       type(auto_diff_real_2var_order2), intent(out) :: this
       integer, intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1val2 = 0_dp
-      this%d2val1 = 0_dp
-      this%d1val1_d1val2 = 0_dp
-      this%d2val2 = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1val2 = 0.0_dp
+      this%d2val1 = 0.0_dp
+      this%d1val1_d1val2 = 0.0_dp
+      this%d2val2 = 0.0_dp
    end subroutine assign_from_int
    
    function equal_self(this, other) result(z)
@@ -1969,9 +1969,9 @@ function differentiate_auto_diff_real_2var_order2_1(this) result(derivative)
       derivative%val = this%d1val1
       derivative%d1val1 = this%d2val1
       derivative%d1val2 = this%d1val1_d1val2
-      derivative%d2val1 = 0_dp
-      derivative%d1val1_d1val2 = 0_dp
-      derivative%d2val2 = 0_dp
+      derivative%d2val1 = 0.0_dp
+      derivative%d1val1_d1val2 = 0.0_dp
+      derivative%d2val2 = 0.0_dp
    end function differentiate_auto_diff_real_2var_order2_1
    
    function differentiate_auto_diff_real_2var_order2_2(this) result(derivative)
@@ -1980,9 +1980,9 @@ function differentiate_auto_diff_real_2var_order2_2(this) result(derivative)
       derivative%val = this%d1val2
       derivative%d1val1 = this%d1val1_d1val2
       derivative%d1val2 = this%d2val2
-      derivative%d2val1 = 0_dp
-      derivative%d1val1_d1val2 = 0_dp
-      derivative%d2val2 = 0_dp
+      derivative%d2val1 = 0.0_dp
+      derivative%d1val1_d1val2 = 0.0_dp
+      derivative%d2val2 = 0.0_dp
    end function differentiate_auto_diff_real_2var_order2_2
    
 end module auto_diff_real_2var_order2_module
\ No newline at end of file
diff --git a/auto_diff/private/auto_diff_real_2var_order3_module.f90 b/auto_diff/private/auto_diff_real_2var_order3_module.f90
index 92a7d7ef0..1cc5ff8e2 100644
--- a/auto_diff/private/auto_diff_real_2var_order3_module.f90
+++ b/auto_diff/private/auto_diff_real_2var_order3_module.f90
@@ -393,30 +393,30 @@ subroutine assign_from_real_dp(this, other)
       type(auto_diff_real_2var_order3), intent(out) :: this
       real(dp), intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1val2 = 0_dp
-      this%d2val1 = 0_dp
-      this%d1val1_d1val2 = 0_dp
-      this%d2val2 = 0_dp
-      this%d3val1 = 0_dp
-      this%d2val1_d1val2 = 0_dp
-      this%d1val1_d2val2 = 0_dp
-      this%d3val2 = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1val2 = 0.0_dp
+      this%d2val1 = 0.0_dp
+      this%d1val1_d1val2 = 0.0_dp
+      this%d2val2 = 0.0_dp
+      this%d3val1 = 0.0_dp
+      this%d2val1_d1val2 = 0.0_dp
+      this%d1val1_d2val2 = 0.0_dp
+      this%d3val2 = 0.0_dp
    end subroutine assign_from_real_dp
    
    subroutine assign_from_int(this, other)
       type(auto_diff_real_2var_order3), intent(out) :: this
       integer, intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1val2 = 0_dp
-      this%d2val1 = 0_dp
-      this%d1val1_d1val2 = 0_dp
-      this%d2val2 = 0_dp
-      this%d3val1 = 0_dp
-      this%d2val1_d1val2 = 0_dp
-      this%d1val1_d2val2 = 0_dp
-      this%d3val2 = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1val2 = 0.0_dp
+      this%d2val1 = 0.0_dp
+      this%d1val1_d1val2 = 0.0_dp
+      this%d2val2 = 0.0_dp
+      this%d3val1 = 0.0_dp
+      this%d2val1_d1val2 = 0.0_dp
+      this%d1val1_d2val2 = 0.0_dp
+      this%d3val2 = 0.0_dp
    end subroutine assign_from_int
    
    function equal_self(this, other) result(z)
@@ -3317,10 +3317,10 @@ function differentiate_auto_diff_real_2var_order3_1(this) result(derivative)
       derivative%d2val1 = this%d3val1
       derivative%d1val1_d1val2 = this%d2val1_d1val2
       derivative%d2val2 = this%d1val1_d2val2
-      derivative%d3val1 = 0_dp
-      derivative%d2val1_d1val2 = 0_dp
-      derivative%d1val1_d2val2 = 0_dp
-      derivative%d3val2 = 0_dp
+      derivative%d3val1 = 0.0_dp
+      derivative%d2val1_d1val2 = 0.0_dp
+      derivative%d1val1_d2val2 = 0.0_dp
+      derivative%d3val2 = 0.0_dp
    end function differentiate_auto_diff_real_2var_order3_1
    
    function differentiate_auto_diff_real_2var_order3_2(this) result(derivative)
@@ -3332,10 +3332,10 @@ function differentiate_auto_diff_real_2var_order3_2(this) result(derivative)
       derivative%d2val1 = this%d2val1_d1val2
       derivative%d1val1_d1val2 = this%d1val1_d2val2
       derivative%d2val2 = this%d3val2
-      derivative%d3val1 = 0_dp
-      derivative%d2val1_d1val2 = 0_dp
-      derivative%d1val1_d2val2 = 0_dp
-      derivative%d3val2 = 0_dp
+      derivative%d3val1 = 0.0_dp
+      derivative%d2val1_d1val2 = 0.0_dp
+      derivative%d1val1_d2val2 = 0.0_dp
+      derivative%d3val2 = 0.0_dp
    end function differentiate_auto_diff_real_2var_order3_2
    
 end module auto_diff_real_2var_order3_module
\ No newline at end of file
diff --git a/auto_diff/private/auto_diff_real_4var_order1_module.f90 b/auto_diff/private/auto_diff_real_4var_order1_module.f90
index 1c8c97be2..78001364f 100644
--- a/auto_diff/private/auto_diff_real_4var_order1_module.f90
+++ b/auto_diff/private/auto_diff_real_4var_order1_module.f90
@@ -393,20 +393,20 @@ subroutine assign_from_real_dp(this, other)
       type(auto_diff_real_4var_order1), intent(out) :: this
       real(dp), intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1val2 = 0_dp
-      this%d1val3 = 0_dp
-      this%d1val4 = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1val2 = 0.0_dp
+      this%d1val3 = 0.0_dp
+      this%d1val4 = 0.0_dp
    end subroutine assign_from_real_dp
    
    subroutine assign_from_int(this, other)
       type(auto_diff_real_4var_order1), intent(out) :: this
       integer, intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1val2 = 0_dp
-      this%d1val3 = 0_dp
-      this%d1val4 = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1val2 = 0.0_dp
+      this%d1val3 = 0.0_dp
+      this%d1val4 = 0.0_dp
    end subroutine assign_from_int
    
    function equal_self(this, other) result(z)
@@ -1693,40 +1693,40 @@ function differentiate_auto_diff_real_4var_order1_1(this) result(derivative)
       type(auto_diff_real_4var_order1), intent(in) :: this
       type(auto_diff_real_4var_order1) :: derivative
       derivative%val = this%d1val1
-      derivative%d1val1 = 0_dp
-      derivative%d1val2 = 0_dp
-      derivative%d1val3 = 0_dp
-      derivative%d1val4 = 0_dp
+      derivative%d1val1 = 0.0_dp
+      derivative%d1val2 = 0.0_dp
+      derivative%d1val3 = 0.0_dp
+      derivative%d1val4 = 0.0_dp
    end function differentiate_auto_diff_real_4var_order1_1
    
    function differentiate_auto_diff_real_4var_order1_2(this) result(derivative)
       type(auto_diff_real_4var_order1), intent(in) :: this
       type(auto_diff_real_4var_order1) :: derivative
       derivative%val = this%d1val2
-      derivative%d1val1 = 0_dp
-      derivative%d1val2 = 0_dp
-      derivative%d1val3 = 0_dp
-      derivative%d1val4 = 0_dp
+      derivative%d1val1 = 0.0_dp
+      derivative%d1val2 = 0.0_dp
+      derivative%d1val3 = 0.0_dp
+      derivative%d1val4 = 0.0_dp
    end function differentiate_auto_diff_real_4var_order1_2
    
    function differentiate_auto_diff_real_4var_order1_3(this) result(derivative)
       type(auto_diff_real_4var_order1), intent(in) :: this
       type(auto_diff_real_4var_order1) :: derivative
       derivative%val = this%d1val3
-      derivative%d1val1 = 0_dp
-      derivative%d1val2 = 0_dp
-      derivative%d1val3 = 0_dp
-      derivative%d1val4 = 0_dp
+      derivative%d1val1 = 0.0_dp
+      derivative%d1val2 = 0.0_dp
+      derivative%d1val3 = 0.0_dp
+      derivative%d1val4 = 0.0_dp
    end function differentiate_auto_diff_real_4var_order1_3
    
    function differentiate_auto_diff_real_4var_order1_4(this) result(derivative)
       type(auto_diff_real_4var_order1), intent(in) :: this
       type(auto_diff_real_4var_order1) :: derivative
       derivative%val = this%d1val4
-      derivative%d1val1 = 0_dp
-      derivative%d1val2 = 0_dp
-      derivative%d1val3 = 0_dp
-      derivative%d1val4 = 0_dp
+      derivative%d1val1 = 0.0_dp
+      derivative%d1val2 = 0.0_dp
+      derivative%d1val3 = 0.0_dp
+      derivative%d1val4 = 0.0_dp
    end function differentiate_auto_diff_real_4var_order1_4
    
 end module auto_diff_real_4var_order1_module
\ No newline at end of file
diff --git a/auto_diff/private/auto_diff_real_star_order1_module.f90 b/auto_diff/private/auto_diff_real_star_order1_module.f90
index dad0e88c0..cc9df242d 100644
--- a/auto_diff/private/auto_diff_real_star_order1_module.f90
+++ b/auto_diff/private/auto_diff_real_star_order1_module.f90
@@ -367,14 +367,14 @@ subroutine assign_from_real_dp(this, other)
       type(auto_diff_real_star_order1), intent(out) :: this
       real(dp), intent(in) :: other
       this%val = other
-      this%d1Array = 0_dp
+      this%d1Array = 0.0_dp
    end subroutine assign_from_real_dp
    
    subroutine assign_from_int(this, other)
       type(auto_diff_real_star_order1), intent(out) :: this
       integer, intent(in) :: other
       this%val = other
-      this%d1Array = 0_dp
+      this%d1Array = 0.0_dp
    end subroutine assign_from_int
    
    function equal_self(this, other) result(z)
diff --git a/auto_diff/private/auto_diff_real_tdc_module.f90 b/auto_diff/private/auto_diff_real_tdc_module.f90
index 6995f681f..39a11feae 100644
--- a/auto_diff/private/auto_diff_real_tdc_module.f90
+++ b/auto_diff/private/auto_diff_real_tdc_module.f90
@@ -376,18 +376,18 @@ subroutine assign_from_real_dp(this, other)
       type(auto_diff_real_tdc), intent(out) :: this
       real(dp), intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1Array = 0_dp
-      this%d1val1_d1Array = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1Array = 0.0_dp
+      this%d1val1_d1Array = 0.0_dp
    end subroutine assign_from_real_dp
    
    subroutine assign_from_int(this, other)
       type(auto_diff_real_tdc), intent(out) :: this
       integer, intent(in) :: other
       this%val = other
-      this%d1val1 = 0_dp
-      this%d1Array = 0_dp
-      this%d1val1_d1Array = 0_dp
+      this%d1val1 = 0.0_dp
+      this%d1Array = 0.0_dp
+      this%d1val1_d1Array = 0.0_dp
    end subroutine assign_from_int
    
    function equal_self(this, other) result(z)
@@ -1678,9 +1678,9 @@ function differentiate_auto_diff_real_tdc_1(this) result(derivative)
       type(auto_diff_real_tdc), intent(in) :: this
       type(auto_diff_real_tdc) :: derivative
       derivative%val = this%d1val1
-      derivative%d1val1 = 0_dp
+      derivative%d1val1 = 0.0_dp
       derivative%d1Array = this%d1val1_d1Array
-      derivative%d1val1_d1Array = 0_dp
+      derivative%d1val1_d1Array = 0.0_dp
    end function differentiate_auto_diff_real_tdc_1
    
 end module auto_diff_real_tdc_module