...solved all of them but the double precision product. The following piece of code is not working as it should: program euclidean_norm double precision x(2) As highegg already said, you should add here: double precision,external:: norm The "external" attribute ... and 1.52316362648 are REAL (not DOUBLE PRECISION) constants. Use 2.70753335146d0 and ...

...> v2.21). I solved all of them but the double precision product. The following piece of code...working as it should: program euclidean_norm double precision x(2) x(1) ...> write(*,*)x write(*,*)norm(x) end double precision function norm(x) ...different way to get it using double precision? Many thanks in advance. ...

... says that the relative ordering of a couple priorities is not the same as for an optimizing production code compiler. It is not a crticism. In fact it is a complement when there are several .... A student oriented compiler had better not miss any features as might be tolerated on a production compiler if it promised great optimization. I can no longer check the Polyhedron benchmarks for ...

...+x(2)**2) WARNING - Variable NORMA has been given a value but never used 0009) double precision function norm(x) WARNING - The result of FUNCTION NORM has never been set That seems pretty...+x(2)**2) WARNING - Variable NORMA has been given a value but never used 0009) double precision function norm(x) WARNING - The result of FUNCTION NORM has never been set 0012)...

... Snouck <dsno...@yahoo.nospam> wrote: Jorge Lopez wrote: double precision function norm(x) double precision x(2) norma = sqrt(x(1)**2+x(2)**2) ... been given a value but never used 0009) double precision function norm(x) WARNING - The result ... been given a value but never used 0009) double precision function norm(x) WARNING - The ...

...! The following piece of code is not working as it should: program euclidean_norm double precision x(2) x(1) = 2.70753335146 x(2) = 1.52316362648 write(*,*)x write(*,*)norm(x) end double precision function norm(x) double precision x(2) norma = sqrt(x(1)**2+x(2)**2) end

On 25 mei, 14:08, Daniel Snouck <dsno...@yahoo.nospam> wrote: Jorge Lopez wrote: double precision function norm(x) double precision x(2) norma = sqrt(x(1)**2+x(2)**2) end Others already mentioned a few problems. I would like to add that your...

...$si3$1@nsnmpen2-gest.nuria.telefonica-data.net>, "Jorge Lopez" <not@telling.com> wrote: Is there any intrinsinc funcion to compute the norm, or a different way to get it using double precision? Normally for vectors of arbitrary length you might write something like sqrt(dot_product(x,x)) but for vectors of length 2 it isn't clear if this is simpler or more complicated. ...

Michel Olagnon wrote: Jorge Lopez wrote: Is there any intrinsinc funcion to compute the norm, or a different way to get it using double precision? For 2-D, a simple trick is Abs (Cmplx ( x(1), x(2), Kind(0.0d0))) That should avoid the overflow problem for larger x(1) and x(2), also. It should implement something similar to: a=max(x(1),x(2)) b=min(x(2),x(2)) norm=...

Jorge Lopez wrote: double precision function norm(x) double precision x(2) norma = sqrt(x(1)**2+x(2)**2) end Others already mentioned a few problems. I would like to add that your function is named "norm" but then you assign a value to "norma". Always using "implicit none" seems useful to me.