From mboxrd@z Thu Jan  1 00:00:00 1970
X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on polar.synack.me
X-Spam-Level: 
X-Spam-Status: No, score=-1.3 required=5.0 tests=BAYES_00,INVALID_MSGID
	autolearn=no autolearn_force=no version=3.4.4
X-Google-Language: ENGLISH,ASCII-7-bit
X-Google-Thread: 103376,fc050a66c3b5d87d
X-Google-Attributes: gid103376,public
X-Google-Thread: 1094ba,86e8c626be2471ae
X-Google-Attributes: gid1094ba,public
From: "James Van Buskirk" <torsop@ix.netcom.com>
Subject: Re: F9X twister & ADA (was: n-dim'l vectors)
Date: 2000/04/10
Message-ID: <8cu8fj$eet$1@nntp9.atl.mindspring.net>#1/1
X-Deja-AN: 609403826
References: <8cctts$ujr$1@nnrp1.deja.com> <38EA0440.1ECBC158@ncep.noaa.gov>
 <38ED4ECA.ADB698C9@sdynamix.com> <m2n1n6blie.fsf@vega.qnet.com>
 <38F28A85.53809F39@sdynamix.com>
X-Server-Date: 11 Apr 2000 04:08:19 GMT
X-MimeOLE: Produced By Microsoft MimeOLE V4.72.3110.3
Organization: MindSpring Enterprises
Newsgroups: comp.lang.fortran,comp.lang.ada
Date: 2000-04-11T04:08:19+00:00
List-Id: <comp.lang.ada>


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)