anchor.diff:        
   c
         subroutine logit(theta,n,ltheta)
  -c Truncates theta to between 0 and 1.
  +c Maps (0,1) -> R.
   cf2py intent(hide) n
   cf2py intent(out) ltheta
         DOUBLE PRECISION theta(n), ltheta(n)
  
positive.code:       end do
        RETURN
        END      
  
negative.code:       end do
        RETURN
        END
  c 
        
  
        subroutine invlogit(ltheta,n,theta)
  c Maps R -> (0,1).
  cf2py intent(hide) n
  cf2py intent(out) theta
        DOUBLE PRECISION theta(n), ltheta(n)
  