223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

trustworthy reference

This commit is contained in:
Martin Diehl 2022-05-08 19:57:25 +02:00
parent b376b10b7a
commit 53796fce7a
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/polynomials.f90
View File

@ -94,6 +94,7 @@ end function polynomial_from_dict
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief Evaluate a Polynomial. !> @brief Evaluate a Polynomial.
!> @details https://nvlpubs.nist.gov/nistpubs/jres/71b/jresv71bn1p11_a1b.pdf (eq. 1.2)
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
pure function eval(self,x) result(y) pure function eval(self,x) result(y)
@ -101,15 +102,15 @@ pure function eval(self,x) result(y)
real(pReal), intent(in) :: x real(pReal), intent(in) :: x
real(pReal) :: y real(pReal) :: y
integer :: i integer :: o
y = self%coef(ubound(self%coef,1)) y = self%coef(ubound(self%coef,1))
do i = ubound(self%coef,1)-1, 0, -1 do o = ubound(self%coef,1)-1, 0, -1
#ifndef __INTEL_LLVM_COMPILER #ifndef __INTEL_LLVM_COMPILER
y = y*(x-self%x_ref) +self%coef(i) y = y*(x-self%x_ref) +self%coef(o)
#else #else
y = IEEE_FMA(y,x-self%x_ref,self%coef(i)) y = IEEE_FMA(y,x-self%x_ref,self%coef(o))
#endif #endif
enddo enddo