anchor.diff:      end function
   
       ! Non-parallel version.
  -    function all_rots(p) result(ps)
  +    pure function all_rots(p) result(ps)
           integer, intent(in) :: p(:, :)
           integer             :: p1(size(p, 1), size(p, 2))
           integer             :: d(3, 3, 3)
  
positive.code:         end forall
      end function
  
      function contains_p (ps, p) result(b)
          type(list), intent(in) :: ps
          integer, intent(in)    :: p(:, :)
          logical                :: b
  
negative.code:         end forall
      end subroutine
  
      pure function rotate(p1, axis, theta) result(p2)
          integer, intent(in) :: p1(:, :), axis, theta
          integer             :: p2(size(p1, 1), size(p1, 2))
  
  