Hi All,
I have sent the following fix to John to test:
> vlength = 1.0d0/sqrt(vlength) !DIVIDE BY ZERO
> if (vlength<1000.0d0) then
> vector(1:dim) = vector(1:dim)*vlength
> else
> vector(1:dim) = 0.d0
> end if
becomes:
if (vlength>1.0d-6) then
vlength = 1.0d0/sqrt(vlength)
vector(1:dim) = vector(1:dim)*vlength
else
vector(1:dim) = 0.d0
end if
All the best
Ross
/\
\/
|\oss Walker
| HPC Consultant and Staff Scientist |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- ross.rosswalker.co.uk |
|
http://www.rosswalker.co.uk | PGP Key available on request |
Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.
Received on Wed Apr 05 2006 - 23:49:35 PDT
