anchor.diff:      real(real64) :: a
   
       a = exp(x)
  -    nrm_func = exp(lfwa(a))  * a
  +    nrm_func = exp(log_fwa(a))  * a
       
     end function nrm_func
   
  
positive.code:     real(real64) :: a
  
      a = exp(x)
      m2_func = exp(lfwa(a)) * hdir(a)**2 * a
      
    end function m2_func
  
  
negative.code:     real(real64) :: a
  
      a = exp(x)
      m_func = exp(log_fwa(a)) * hdir(a) * a
      
    end function m_func
  
  