anchor.diff:  contains
   
       real(rk) function dot_product_2(v, w) result (dot)
  -        real(rk), intent(in), dimension(d + 1, num_points) :: v, w
  +        real(rk), intent(in), contiguous, dimension(d + 1, num_points) :: v, w
   
   #if defined(PCREO_USE_MKL) .and. defined(PCREO_SINGLE_PREC)
           dot = sdot(num_vars, v, 1, w, 1)
  
positive.code: 
  
      subroutine matrix_multiply_42(c, A, b)
          real(rk), intent(out) :: c(d + 1, num_points)
          real(rk), intent(in) :: A(d + 1, num_points, d + 1, num_points)
          real(rk), intent(in) :: b(d + 1, num_points)
  
  #if defined(PCREO_USE_MKL) .and. defined(PCREO_SINGLE_PREC)
          call ssymv('U', num_vars, 1.0_rk, A, num_vars, b, 1, 0.0_rk, c, 1)
  
negative.code:     end subroutine matrix_multiply_42
  
  
      subroutine symmetric_update_4(A, c, x, y)
          real(rk), intent(inout), contiguous :: &
                  & A(d + 1, num_points, d + 1, num_points)
          real(rk), intent(in) :: c
          real(rk), intent(in), contiguous :: &
                  & x(d + 1, num_points), y(d + 1, num_points)
  
  #if defined(PCREO_USE_MKL) .and. defined(PCREO_SINGLE_PREC)
          call ssyr2('U', num_vars, c, x, 1, y, 1, A, num_vars)
  