Berwin,
Many thanks, adding a statement:
double precision ddot
at the top fixed the problem! Once I add that the printing to screen works
well too.
Based on what you're telling me about the expected types when not declared,
I'm surprised the code ever worked at all!
Regards,
Emmanuel
On Sat
On 12-01-2014, at 08:26, Berwin A Turlach wrote:
> G'day Emmanuel,
>
> On Sat, 11 Jan 2014 14:08:41 -0800
> Emmanuel Sharef wrote:
>
>> Example: this program creates two double vectors, takes the dot
>> product with ddot, and prints the result:
>>
>> test.f:
>> subroutine testddot(n,x,y
