anchor.diff:        do h1=1,h1d
         do h2=1,h2d
         do h3=1,h3d
  -        triplesx(h3,h2,h1,p6,p5,p4) = o(h3,h2,h1)
  +        triplesx(h3,h2,h1,p6,p5,p4) = o(h1,h3,h2)
         enddo
         enddo
         enddo
  
positive.code:       enddo
        return
        end
  c
        subroutine blas_sd_t_d1_2(h3d,h2d,h1d,p6d,p5d,p4d,h7d,
       2               triplesx,t2sub,v2sub)
  
negative.code:       return
        end
  c
  #if BROKEN
        subroutine blas_sd_t_d1_1(h3d,h2d,h1d,p6d,p5d,p4d,h7d,
       2               triplesx,t2sub,v2sub)
        IMPLICIT NONE
  