hi all,
I am getting mad with the results. I executed same program on compaq and intel compilers. I am getting differnt results. I am calculating function values for different variables.
Da=y(5)
y1=y(1)
y2=y(2)
x1=y(3)
x2=y(4)
call dataf(Da,gamma,beta,daph,dapm,km,kh)
xf=Da*(x1+x2)/(2.d0*(Da+Dapm))
yf=Da*(y1+y2)/(2.d0*(Da+Daph))
f(1)=km*(x2-x1)+(xf-x1)/Dapm+(1.d0-x1)*dexp(gamma*y1/(1.d0+y1))
f(2)=km*(x1-x2)+(xf-x2)/Dapm+(1.d0-x2)*dexp(gamma*y2/(1.d0+y2))
f(3)=kh*(y2-y1)+(yf-y1)/Daph+beta*(1.d0-x1)*dexp(gamma*y1/(1.d0+y1))
f(4)=kh*(y1-y2)+(yf-y2)/Daph+beta*(1.d0-x2)*dexp(gamma*y2/(1.d0+y2))
for
y(:,1)=(/ 0.01410494974141788d0,0.01410494974141788d0 ,0.04402030190585460d0,0.04402030190585460d0,0.001d0/)
I got the same f(1) to f(5) on both
but on intel when I incremented y with samll value, I get different results. I get f(2) and f(4) same. but f(1) and f(3) are correct upto 11 digits only(compaq compiler results exactly match upto 15 digits with literature).
can anybody tell me what is going wrong. I
