anchor.diff:    use, intrinsic :: iso_c_binding, only: c_bool, c_double
     use surface, only: &
          de_casteljau_one_round, evaluate_barycentric, &
  -       evaluate_barycentric_multi, evaluate_cartesian_multi
  +       evaluate_barycentric_multi, evaluate_cartesian_multi, jacobian_both, &
  +       jacobian_det
     use types, only: dp
     use unit_test_helpers, only: print_status, get_random_nodes
     implicit none
     private &
          test_de_casteljau_one_round, test_evaluate_barycentric, &
  -       test_evaluate_barycentric_multi, test_evaluate_cartesian_multi
  +       test_evaluate_barycentric_multi, test_evaluate_cartesian_multi, &
  +       test_jacobian_both, test_jacobian_det
     public surface_all_tests
   
   contains
  
positive.code: 
    end subroutine test_evaluate_cartesian_multi
  
  end module test_surface
  
negative.code:     call test_evaluate_barycentric(success)
      call test_evaluate_barycentric_multi(success)
      call test_evaluate_cartesian_multi(success)
      call test_jacobian_both(success)
      call test_jacobian_det(success)
  
    end subroutine surface_all_tests
  
  