use INTENT(IN) to avoid violating Fortran aliasing rules

see https://github.com/Reference-LAPACK/lapack/issues/152 for details.

fixes issue #152

diff --git a/SRC/disnan.f b/SRC/disnan.f
index da89158..6d26254 100644 (file)
--- a/SRC/disnan.f
+++ b/SRC/disnan.f
@@ -65,7 +65,7 @@
 *     December 2016
 *
 *     .. Scalar Arguments ..
-      DOUBLE PRECISION   DIN
+      DOUBLE PRECISION, INTENT(IN) :: DIN
 *     ..
 *
 *  =====================================================================

diff --git a/SRC/dlaisnan.f b/SRC/dlaisnan.f
index 4b5ebb4..3886dc5 100644 (file)
--- a/SRC/dlaisnan.f
+++ b/SRC/dlaisnan.f
@@ -21,7 +21,7 @@
 *       LOGICAL FUNCTION DLAISNAN( DIN1, DIN2 )
 *
 *       .. Scalar Arguments ..
-*       DOUBLE PRECISION   DIN1, DIN2
+*       DOUBLE PRECISION, INTENT(IN) :: DIN1, DIN2
 *       ..
 *
 *