norm2 is (yet) not available in PGI

src/plastic_nonlocal.f90

@@ -1624,6 +1624,10 @@ use debug,    only: debug_level, &
                     debug_e
 #endif
 use math,     only: math_mul3x3, &
+#ifdef __PGI
+   norm2, &
+#endif
+
                     math_mul33x3, &
                     math_mul33xx33, &
                     math_mul33x33, &

src/quaternions.f90

@@ -354,6 +354,10 @@ end function pow_quat__
 !> ToDo: Lacks any check for invalid operations
 !---------------------------------------------------------------------------------------------------
 type(quaternion) elemental function exp__(self)
+#ifdef __PGI
+ use math, only: &
+   norm2
+#endif
 
  implicit none
  class(quaternion), intent(in) :: self
@@ -374,6 +378,10 @@ end function exp__
 !> ToDo: Lacks any check for invalid operations
 !---------------------------------------------------------------------------------------------------
 type(quaternion) elemental function log__(self)
+#ifdef __PGI
+ use math, only: &
+   norm2
+#endif
 
  implicit none
  class(quaternion), intent(in) :: self
@@ -393,6 +401,10 @@ end function log__
 !> norm of a quaternion
 !---------------------------------------------------------------------------------------------------
 real(pReal) elemental function abs__(a)
+#ifdef __PGI
+ use math, only: &
+   norm2
+#endif
 
  implicit none
  class(quaternion), intent(in) :: a

src/rotations.f90

@@ -157,6 +157,10 @@ end subroutine
 function rotVector(self,v,active)
  use prec, only: &
    dEq
+#ifdef __PGI
+ use math, only: &
+   norm2
+#endif
    
  implicit none
  real(pReal),                 dimension(3) :: rotVector
@@ -573,6 +577,9 @@ pure function ro2ax(ro) result(ax)
  use prec, only: &
    dEq0
  use math, only: &
+#ifdef __PGI
+   norm2, &
+#endif
    PI
 
  implicit none 
@@ -662,6 +669,9 @@ pure function ro2ho(ro) result(ho)
  use prec, only: &
    dEq0
  use math, only: &
+#ifdef __PGI
+   norm2, &
+#endif
    PI
 
  implicit none
@@ -718,6 +728,10 @@ end function qu2om
 function om2qu(om) result(qu)
  use prec, only: &
    dEq
+#ifdef __PGI
+ use math, only: &
+   norm2
+#endif
 
  implicit none
  real(pReal),      intent(in), dimension(3,3) :: om
@@ -791,6 +805,9 @@ pure function qu2ro(qu) result(ro)
  use prec, only: &
    dEq0
  use math, only: &
+#ifdef __PGI
+   norm2, &
+#endif
    math_clip
  
  type(quaternion), intent(in)              :: qu
@@ -819,6 +836,9 @@ pure function qu2ho(qu) result(ho)
  use prec, only: &
    dEq0
  use math, only: &
+#ifdef __PGI
+   norm2, &
+#endif
    math_clip
 
  implicit none