anchor.diff:         Wrk(i) = RR0(i,1)
      ENDDO
      Wrk33 = 0.0D0
  -   Wrk33 = BD_OuterProduct(Wrk,Wrk)
  +   Wrk33 = OuterProduct(Wrk,Wrk)
      C12 = C12 + cet*k1s*Wrk33
      C21 = C21 + cet*k1s*Wrk33
      C22 = C22 + cet*e1s*Wrk33
  
positive.code: 
  
     END SUBROUTINE BD_ComputeIniCoef
  
     FUNCTION BD_OuterProduct(vec1,vec2)
  
     REAL(ReKi),INTENT(IN):: vec1(:),vec2(:)
     REAL(ReKi)::BD_OuterProduct(SIZE(vec1),SIZE(vec2))
  
     INTEGER(IntKi)::i,j,n1,n2
  
     n1=SIZE(vec1)
     n2=SIZE(vec2)
  
     DO i=1,n1
         DO j=1,n2
             BD_OuterProduct(i,j) = vec1(i) * vec2(j)
         ENDDO
     ENDDO
  
     END FUNCTION BD_OuterProduct
  
     SUBROUTINE BD_Static(t,n,u,utimes,p,x,xd,z,OtherState,ErrStat,ErrMsg)
  
  
negative.code:    ENDDO
  
     e3 = 0.0D0
     e3 = Cross_Product(e1,e2)
     DO i=1,3
         Rr(i,3) = e3(i)
     ENDDO
  