program array_with_array_eval

use iso_fortran_env

use FEQParse

implicit none

integer,parameter :: N = 10000

type(EquationParser) :: f

character(LEN=1),dimension(1) :: independentVars

character(LEN=30) :: eqChar

real :: x(1:N,1)

real :: feval(1:N)

integer :: i

real :: t1,t2

! Specify the independent variables

independentVars = (/'x'/)

! Specify an equation string that we want to evaluate

eqChar = 'f = exp( -(x^2) )'

! Create the EquationParser object

f = EquationParser(eqChar,independentVars)

! Evaluate the equation

call cpu_time(t1)

do i = 1,N

x(i,1) = -1.0_real32+(2.0_real32)/real(N,real32)*real(i-1,real32)

enddo

feval = f%evaluate(x)

call cpu_time(t2)

print*,"runtime :",(t2-t1)," s"

endprogram array_with_array_eval