File: xdot.F

package info (click to toggle)
aces3 3.0.6-7
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 82,460 kB
  • sloc: fortran: 225,647; ansic: 20,413; cpp: 4,349; makefile: 953; sh: 137
file content (27 lines) | stat: -rw-r--r-- 507 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
 

      double precision function xdot(n,dx,incx,dy,incy)

      double precision dx(*),dy(*),ddot
      integer incx, incy, n

#ifdef _32BIT_BLAS
      integer*4 int_n, int_incx, int_incy
#else
      integer   int_n, int_incx, int_incy
#endif /* _32BIT_BLAS */

      if (n.le.0) return

      int_incx = incx
      int_incy = incy
      int_n    = n

#ifdef USE_SP_BLAS
      xdot = sdot(int_n,dx,int_incx,dy,int_incy)
#else
      xdot = ddot(int_n,dx,int_incx,dy,int_incy)
#endif

      return
      end