anchor.diff:  
   
      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)
   
  
positive.code:    ENDDO
  
     e3 = 0.0D0
     e3 = BD_CrossProduct(e1,e2)
     DO i=1,3
         Rr(i,3) = e3(i)
     ENDDO
  
negative.code:    PUBLIC :: BD_CrvCompose                ! Tight coupling routine for updating discrete states
     PUBLIC :: BD_CrvExtractCrv                ! Tight coupling routine for updating discrete states
     PUBLIC :: BD_Tilde
     PUBLIC :: ludcmp
     PUBLIC :: lubksb
     PUBLIC :: BD_MotionTensor
  