Displaying 1 result from an estimated 1 matches for "testddot".
Did you mean:
testdot
2014 Jan 11
1
Fortran BLAS giving bad results
...he
problem predates Mavericks. The package used to build and check without
errors when I originally put it on r-forge. Would appreciate suggestions.
Thanks,
Emmanuel
Example: this program creates two double vectors, takes the dot product
with ddot, and prints the result:
test.f:
subroutine testddot(n,x,y)
integer i,n
double precision x,y,d1(n),d2(n)
do 100 i=1,n
d1(i)=x
d2(i)=y
100 continue
print *,ddot(n,d1,1,d2,1)
end
Compiling
R CMD SHLIB test.f
Running:
> dyn.load("test.so")
> test<-.Fortran("testddot",as.integ...