anchor.diff:        MODULE PROCEDURE LAPACK_sgbsv
      END INTERFACE
   
  +   INTERFACE LAPACK_gemm   ! Computes scalar1*op( A )*op( B ) + scalar2*C where op(x) = x or op(x) = x**T for matrices A, B, and C
  +      MODULE PROCEDURE LAPACK_dgemm
  +      MODULE PROCEDURE LAPACK_sgemm
  +   END INTERFACE
  +   
  +   
      INTERFACE LAPACK_gesv ! Computes the solution to system of linear equations A * X = B for GE matrices
         MODULE PROCEDURE LAPACK_dgesv
         MODULE PROCEDURE LAPACK_sgesv
  
positive.code:    RETURN
     END SUBROUTINE LAPACK_SGBSV
  !=======================================================================
     SUBROUTINE LAPACK_DGESV ( N, A, IPIV, B, ErrStat, ErrMsg )
  
  
  
negative.code:    !      we need to check this somehow to make sure the right routines are called.
     ! (or define a directive that uses
  
     ! http://www.netlib.org/lapack/explore-html/ 
     
     
     IMPLICIT  NONE
  
     INTERFACE LAPACK_gbsv ! Computes the solution to system of linear equations A * X = B for GB matrices
  