anchor.diff:        return
         end
   c
  +#if BROKEN
         subroutine blas_sd_t_d1_1(h3d,h2d,h1d,p6d,p5d,p4d,h7d,
        2               triplesx,t2sub,v2sub)
         IMPLICIT NONE
  
positive.code:       do h1=1,h1d
        do h2=1,h2d
        do h3=1,h3d
          triplesx(h3,h2,h1,p6,p5,p4) = o(h3,h2,h1)
        enddo
        enddo
        enddo
  
negative.code:       enddo
        return
        end
  c
  #else
        subroutine blas_sd_t_d1_1(h3d,h2d,h1d,p6d,p5d,p4d,h7d,
       2               triplesx,t2sub,v2sub)
        IMPLICIT NONE
        integer h3d,h2d,h1d,p6d,p5d,p4d,h7d
        integer h3,h2,h1,p6,p5,p4,h7
        double precision triplesx(h3d,h2d,h1d,p6d,p5d,p4d)
        double precision t2sub(h7d,p4d,p5d,h1d)
        double precision v2sub(h3d,h2d,p6d,h7d)
  !     TEMP
        double precision o(h2d,h1d)
        double precision t(h7d,h1d)
        double precision v(h2d,h7d)
  !     GOAL
  !     for all 654:
  !       2,1 += 7,1 * 2,7 (TT-DGEMM)
        do p4=1,p4d
        do p5=1,p5d
        do p6=1,p6d
        do h3=1,h3d
  !     CREATE TEMP (should be eliminated)
        do h1=1,h1d
        do h2=1,h2d
          o(h2,h1) = triplesx(h3,h2,h1,p6,p5,p4)
        enddo
        enddo
  !     TRANSPOSE T
        do h1=1,h1d
        do h7=1,h7d
          t(h1,h7) = t2sub(h7,p4,p5,h1)
        enddo
        enddo
  !     TRANSPOSE V
        do h7=1,h7d
        do h2=1,h2d
          v(h7,h2) = v2sub(h3,h2,p6,h7)
        enddo
        enddo
  !     DGEMM
  !        C : h2d rows x h1d cols
          call dgemm('N','N',h2d,h1d,h7d,     ! m,n,k
       &             -1.0d0,t,h2d,            ! alpha,a,lda=m
       &             v,h7d,1.0d0,             ! b,ldb=k,beta
       &             o,h2d)                   ! c,ldc=m
  !     END DGEMM
  !     OUTPUT
        do h1=1,h1d
        do h2=1,h2d
          triplesx(h3,h2,h1,p6,p5,p4) = o(h1,h2)
        enddo
        enddo
  !     END OUTPUT
        enddo
        enddo
        enddo
        enddo
        return
        end
  #endif
  c
        subroutine blas_sd_t_d1_2(h3d,h2d,h1d,p6d,p5d,p4d,h7d,
       2               triplesx,t2sub,v2sub)
  