1. Fix banded matric derivative checking

2. Standardize the output generation.

Author: sriharikrishna

BLAS/test/test_caxpy.f90

@@ -47,14 +47,14 @@ subroutine run_test_for_size(n, passed)
     integer :: incy
 
     ! Derivative variables
+    complex(4), dimension(n) :: cx_d
     complex(4), dimension(n) :: cy_d
     complex(4) :: ca_d
-    complex(4), dimension(n) :: cx_d
 
     ! Array restoration and derivative storage
+    complex(4), dimension(n) :: cx_orig, cx_d_orig
     complex(4), dimension(n) :: cy_orig, cy_d_orig
     complex(4) :: ca_orig, ca_d_orig
-    complex(4), dimension(n) :: cx_orig, cx_d_orig
     real(4) :: temp_re, temp_im  ! For complex random init
     integer :: i, j
 
@@ -80,24 +80,24 @@ subroutine run_test_for_size(n, passed)
     do i = 1, n
       call random_number(temp_re)
       call random_number(temp_im)
-      cy_d(i) = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
+      cx_d(i) = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
     end do
-    call random_number(temp_re)
-    call random_number(temp_im)
-    ca_d = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
     do i = 1, n
       call random_number(temp_re)
       call random_number(temp_im)
-      cx_d(i) = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
+      cy_d(i) = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
     end do
+    call random_number(temp_re)
+    call random_number(temp_im)
+    ca_d = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
 
     ! Store _orig and _d_orig
+    cx_d_orig = cx_d
     cy_d_orig = cy_d
     ca_d_orig = ca_d
-    cx_d_orig = cx_d
+    cx_orig = cx
     cy_orig = cy
     ca_orig = ca
-    cx_orig = cx
 
     write(*,*) 'Testing CAXPY (n =', n, ')'
     cy_orig = cy
@@ -108,17 +108,17 @@ subroutine run_test_for_size(n, passed)
     write(*,*) 'Function calls completed successfully'
 
     ! Numerical differentiation check
-    call check_derivatives_numerically(n, nsize, cy_orig, ca_orig, cx_orig, cy_d_orig, ca_d_orig, cx_d_orig, cy_d, passed)
+    call check_derivatives_numerically(n, nsize, cx_orig, cy_orig, ca_orig, cx_d_orig, cy_d_orig, ca_d_orig, cy_d, passed)
 
   end subroutine run_test_for_size
 
-  subroutine check_derivatives_numerically(n, nsize, cy_orig, ca_orig, cx_orig, cy_d_orig, ca_d_orig, cx_d_orig, cy_d, passed)
+  subroutine check_derivatives_numerically(n, nsize, cx_orig, cy_orig, ca_orig, cx_d_orig, cy_d_orig, ca_d_orig, cy_d, passed)
     implicit none
     integer, intent(in) :: n
     integer, intent(in) :: nsize
+    complex(4), intent(in) :: cx_orig(n), cx_d_orig(n)
     complex(4), intent(in) :: cy_orig(n), cy_d_orig(n)
     complex(4), intent(in) :: ca_orig, ca_d_orig
-    complex(4), intent(in) :: cx_orig(n), cx_d_orig(n)
     complex(4), intent(in) :: cy_d(n)
     logical, intent(out) :: passed
 
@@ -129,9 +129,9 @@ subroutine check_derivatives_numerically(n, nsize, cy_orig, ca_orig, cx_orig, cy
     logical :: has_large_errors
     complex(4), dimension(n) :: cy_forward, cy_backward
     integer :: i, j
+    complex(4), dimension(n) :: cx
     complex(4), dimension(n) :: cy
     complex(4) :: ca
-    complex(4), dimension(n) :: cx
 
     max_error = 0.0e0
     has_large_errors = .false.
@@ -140,16 +140,16 @@ subroutine check_derivatives_numerically(n, nsize, cy_orig, ca_orig, cx_orig, cy
     write(*,*) 'Step size h =', h
 
     ! Forward perturbation: f(x + h)
+    cx = cx_orig + h * cx_d_orig
     cy = cy_orig + h * cy_d_orig
     ca = ca_orig + h * ca_d_orig
-    cx = cx_orig + h * cx_d_orig
     call caxpy(nsize, ca, cx, 1, cy, 1)
     cy_forward = cy
 
     ! Backward perturbation: f(x - h)
+    cx = cx_orig - h * cx_d_orig
     cy = cy_orig - h * cy_d_orig
     ca = ca_orig - h * ca_d_orig
-    cx = cx_orig - h * cx_d_orig
     call caxpy(nsize, ca, cx, 1, cy, 1)
     cy_backward = cy
 
@@ -178,7 +178,7 @@ subroutine check_derivatives_numerically(n, nsize, cy_orig, ca_orig, cx_orig, cy
     write(*,*) 'Tolerance thresholds: rtol=1.0e-3, atol=1.0e-3'
     passed = .not. has_large_errors
     if (has_large_errors) then
-      write(*,*) 'FAIL: Large errors detected in derivatives (outside tolerance)'
+      write(*,*) 'FAIL: Derivatives are outside tolerance'
     else
       write(*,*) 'PASS: Derivatives are within tolerance (rtol + atol)'
     end if

BLAS/test/test_caxpy_reverse.f90

@@ -205,13 +205,11 @@ subroutine check_vjp_numerically(n, nsize, incx_val, incy_val, ca_orig, cx_orig,
       relative_error = abs_error
     end if
     max_error = relative_error
-
-    write(*,*) ''
     write(*,*) 'Maximum relative error:', max_error
     write(*,*) 'Tolerance thresholds: rtol=1.0e-3, atol=1.0e-3'
     passed = .not. has_large_errors
     if (has_large_errors) then
-      write(*,*) 'FAIL: Large errors detected in derivatives (outside tolerance)'
+      write(*,*) 'FAIL: Derivatives are outside tolerance'
     else
       write(*,*) 'PASS: Derivatives are within tolerance (rtol + atol)'
     end if

BLAS/test/test_caxpy_vector_forward.f90

@@ -29,9 +29,9 @@ program test_caxpy_vector_forward
     all_passed = all_passed .and. passed
   end do
   if (all_passed) then
-    write(*,*) 'PASS: Vector forward mode - all sizes completed successfully'
+    write(*,*) 'PASS: All sizes completed successfully'
   else
-    write(*,*) 'FAIL: Vector forward mode - one or more sizes had derivative errors'
+    write(*,*) 'FAIL: One or more sizes had derivative errors'
   end if
 
 contains
@@ -126,7 +126,7 @@ subroutine check_derivatives_numerically(n, nbdirs, nsize, incx_val, incy_val, a
     max_error = 0.0e0
     has_large_errors = .false.
 
-    write(*,*) 'Checking vector derivatives against numerical differentiation:'
+    write(*,*) 'Checking derivatives against numerical differentiation:'
     write(*,*) 'Step size h =', h
 
     do idir = 1, nbdirs
@@ -152,13 +152,13 @@ subroutine check_derivatives_numerically(n, nbdirs, nsize, incx_val, incy_val, a
       end do
     end do
 
-    write(*,*) 'Maximum relative error across all directions:', max_error
+    write(*,*) 'Maximum relative error:', max_error
     write(*,*) 'Tolerance thresholds: rtol=1.0e-3, atol=1.0e-3'
     passed = .not. has_large_errors
     if (has_large_errors) then
-      write(*,*) 'FAIL: Large errors detected in vector derivatives (outside tolerance)'
+      write(*,*) 'FAIL: Derivatives are outside tolerance'
     else
-      write(*,*) 'PASS: Vector derivatives are within tolerance (rtol + atol)'
+      write(*,*) 'PASS: Derivatives are within tolerance (rtol + atol)'
     end if
 
   end subroutine check_derivatives_numerically

BLAS/test/test_caxpy_vector_reverse.f90

@@ -29,9 +29,9 @@ program test_caxpy_vector_reverse
     all_passed = all_passed .and. passed
   end do
   if (all_passed) then
-    write(*,*) 'PASS: Vector reverse mode - all sizes completed successfully'
+    write(*,*) 'PASS: All sizes completed successfully'
   else
-    write(*,*) 'FAIL: Vector reverse mode - one or more sizes had derivative errors'
+    write(*,*) 'FAIL: One or more sizes had derivative errors'
   end if
 
 contains
@@ -173,13 +173,11 @@ subroutine check_vjp_numerically(n, nbdirs, nsize, incx_val, incy_val, alpha_ori
       end if
       if (relative_error > max_error) max_error = relative_error
     end do
-
-    write(*,*) ''
     write(*,*) 'Maximum relative error:', max_error
     write(*,*) 'Tolerance thresholds: rtol=1.0e-3, atol=1.0e-3'
     passed = .not. has_large_errors
     if (has_large_errors) then
-      write(*,*) 'FAIL: Large errors detected in derivatives (outside tolerance)'
+      write(*,*) 'FAIL: Derivatives are outside tolerance'
     else
       write(*,*) 'PASS: Derivatives are within tolerance (rtol + atol)'
     end if

BLAS/test/test_ccopy.f90

@@ -46,12 +46,12 @@ subroutine run_test_for_size(n, passed)
     integer :: incy
 
     ! Derivative variables
-    complex(4), dimension(n) :: cy_d
     complex(4), dimension(n) :: cx_d
+    complex(4), dimension(n) :: cy_d
 
     ! Array restoration and derivative storage
-    complex(4), dimension(n) :: cy_orig, cy_d_orig
     complex(4), dimension(n) :: cx_orig, cx_d_orig
+    complex(4), dimension(n) :: cy_orig, cy_d_orig
     real(4) :: temp_re, temp_im  ! For complex random init
     integer :: i, j
 
@@ -74,19 +74,19 @@ subroutine run_test_for_size(n, passed)
     do i = 1, n
       call random_number(temp_re)
       call random_number(temp_im)
-      cy_d(i) = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
+      cx_d(i) = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
     end do
     do i = 1, n
       call random_number(temp_re)
       call random_number(temp_im)
-      cx_d(i) = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
+      cy_d(i) = cmplx(temp_re * 2.0 - 1.0, temp_im * 2.0 - 1.0, kind=4)
     end do
 
     ! Store _orig and _d_orig
-    cy_d_orig = cy_d
     cx_d_orig = cx_d
-    cy_orig = cy
+    cy_d_orig = cy_d
     cx_orig = cx
+    cy_orig = cy
 
     write(*,*) 'Testing CCOPY (n =', n, ')'
 
@@ -169,7 +169,7 @@ subroutine check_derivatives_numerically(n, nsize, cx_orig, cy_orig, cx_d_orig,
     write(*,*) 'Tolerance thresholds: rtol=1.0e-3, atol=1.0e-3'
     passed = .not. has_large_errors
     if (has_large_errors) then
-      write(*,*) 'FAIL: Large errors detected in derivatives (outside tolerance)'
+      write(*,*) 'FAIL: Derivatives are outside tolerance'
     else
       write(*,*) 'PASS: Derivatives are within tolerance (rtol + atol)'
     end if

BLAS/test/test_ccopy_reverse.f90

@@ -187,13 +187,11 @@ subroutine check_vjp_numerically(n, nsize, incx_val, incy_val, cx_orig, cy_orig,
       relative_error = abs_error
     end if
     max_error = relative_error
-
-    write(*,*) ''
     write(*,*) 'Maximum relative error:', max_error
     write(*,*) 'Tolerance thresholds: rtol=1.0e-3, atol=1.0e-3'
     passed = .not. has_large_errors
     if (has_large_errors) then
-      write(*,*) 'FAIL: Large errors detected in derivatives (outside tolerance)'
+      write(*,*) 'FAIL: Derivatives are outside tolerance'
     else
       write(*,*) 'PASS: Derivatives are within tolerance (rtol + atol)'
     end if

BLAS/test/test_ccopy_vector_forward.f90

@@ -29,9 +29,9 @@ program test_ccopy_vector_forward
     all_passed = all_passed .and. passed
   end do
   if (all_passed) then
-    write(*,*) 'PASS: Vector forward mode - all sizes completed successfully'
+    write(*,*) 'PASS: All sizes completed successfully'
   else
-    write(*,*) 'FAIL: Vector forward mode - one or more sizes had derivative errors'
+    write(*,*) 'FAIL: One or more sizes had derivative errors'
   end if
 
 contains
@@ -113,7 +113,7 @@ subroutine check_derivatives_numerically(n, nbdirs, nsize, incx_val, incy_val, x
     max_error = 0.0e0
     has_large_errors = .false.
 
-    write(*,*) 'Checking vector derivatives against numerical differentiation:'
+    write(*,*) 'Checking derivatives against numerical differentiation:'
     write(*,*) 'Step size h =', h
 
     do idir = 1, nbdirs
@@ -137,13 +137,13 @@ subroutine check_derivatives_numerically(n, nbdirs, nsize, incx_val, incy_val, x
       end do
     end do
 
-    write(*,*) 'Maximum relative error across all directions:', max_error
+    write(*,*) 'Maximum relative error:', max_error
     write(*,*) 'Tolerance thresholds: rtol=1.0e-3, atol=1.0e-3'
     passed = .not. has_large_errors
     if (has_large_errors) then
-      write(*,*) 'FAIL: Large errors detected in vector derivatives (outside tolerance)'
+      write(*,*) 'FAIL: Derivatives are outside tolerance'
     else
-      write(*,*) 'PASS: Vector derivatives are within tolerance (rtol + atol)'
+      write(*,*) 'PASS: Derivatives are within tolerance (rtol + atol)'
     end if
 
   end subroutine check_derivatives_numerically

BLAS/test/test_ccopy_vector_reverse.f90

@@ -29,9 +29,9 @@ program test_ccopy_vector_reverse
     all_passed = all_passed .and. passed
   end do
   if (all_passed) then
-    write(*,*) 'PASS: Vector reverse mode - all sizes completed successfully'
+    write(*,*) 'PASS: All sizes completed successfully'
   else
-    write(*,*) 'FAIL: Vector reverse mode - one or more sizes had derivative errors'
+    write(*,*) 'FAIL: One or more sizes had derivative errors'
   end if
 
 contains
@@ -154,13 +154,11 @@ subroutine check_vjp_numerically(n, nbdirs, nsize, incx_val, incy_val, x_orig, y
       end if
       if (relative_error > max_error) max_error = relative_error
     end do
-
-    write(*,*) ''
     write(*,*) 'Maximum relative error:', max_error
     write(*,*) 'Tolerance thresholds: rtol=1.0e-3, atol=1.0e-3'
     passed = .not. has_large_errors
     if (has_large_errors) then
-      write(*,*) 'FAIL: Large errors detected in derivatives (outside tolerance)'
+      write(*,*) 'FAIL: Derivatives are outside tolerance'
     else
       write(*,*) 'PASS: Derivatives are within tolerance (rtol + atol)'
     end if