fast evaluation for (most common) case of constant

This commit is contained in:
Martin Diehl 2022-05-08 16:49:30 +02:00
parent b2052cb3c7
commit d713026f7e
1 changed files with 2 additions and 2 deletions

View File

@ -105,8 +105,8 @@ pure function eval(self,x) result(y)
integer :: i
y = 0.0_pReal
do i = ubound(self%coef,1), 0, -1
y = self%coef(ubound(self%coef,1))
do i = ubound(self%coef,1)-1, 0, -1
#ifndef __INTEL_COMPILER
y = y*(x-self%x_ref) +self%coef(i)
#else