Re: F9X twister & ADA (was: n-dim'l vectors)

["James Van Buskirk" <torsop@ix.netcom.com>]

bv wrote in message <38F28A85.53809F39@sdynamix.com>...

> vector function normalize(v)
> external vector
> vector u,v
>        real norm


>        vni = 1/norm(v)
> do i=1,len(v)
>    u(i) = v(i)*vni
> enddo
>        return u
>        end


>**      somewhere in your code  (sans interfaces & modules)
>        vector u,v, normalize
>                .
>                .
>                .
>        u = normalize(v)

module simple_1d_vector_stuff
   contains
      function normalize(v)
         real, intent(in) :: v(:)
         real normalize(size(v))
         real norm

         if(size(v) == 0) stop 'Can''t normalize zero length vector'
         norm = sqrt(dot_product(v,v)) ! OK for complex version, too
         if(norm == 0) stop 'Can''t normalize the zero vector'
         normalize = v/norm

      end function normalize
end module simple_1d_vector_stuff

**      somewhere in your code
        use simple_1d_vector_stuff
        integer, parameter :: vector_length = 3
        real u(vector_length),v(vector_length)
                .
                .
                .
        u = normalize(v)