Added and used functions deq() and seq() to compare floating-point variables

* Added and used functions deq() and seq() to compare floating-point variables
  * deq(x,y) = true if abs(x-y).le.2*tiny(x)
  * seq(): same, but using single rather than double precision