diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2011-02-18 15:11:31 +0100 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2011-02-18 15:11:31 +0100 |
| commit | f7cd63b964b272fecf8d3411393f60803432fc70 (patch) | |
| tree | 44eb7b949b29c174aa168b1f55b1f355076c79e2 /blas/complexdots.f | |
| parent | 69cecc45e59acbc4d050b9bfc84179f57e6a3ee2 (diff) | |
fix bug #189 (issue with fortran concentions to return COMPLEX values)
Diffstat (limited to 'blas/complexdots.f')
| -rw-r--r-- | blas/complexdots.f | 43 |
1 files changed, 43 insertions, 0 deletions
diff --git a/blas/complexdots.f b/blas/complexdots.f new file mode 100644 index 000000000..a7da51d16 --- /dev/null +++ b/blas/complexdots.f @@ -0,0 +1,43 @@ + COMPLEX FUNCTION CDOTC(N,CX,INCX,CY,INCY) + INTEGER INCX,INCY,N + COMPLEX CX(*),CY(*) + COMPLEX RES + EXTERNAL CDOTCW + + CALL CDOTCW(N,CX,INCX,CY,INCY,RES) + CDOTC = RES + RETURN + END + + COMPLEX FUNCTION CDOTU(N,CX,INCX,CY,INCY) + INTEGER INCX,INCY,N + COMPLEX CX(*),CY(*) + COMPLEX RES + EXTERNAL CDOTUW + + CALL CDOTUW(N,CX,INCX,CY,INCY,RES) + CDOTU = RES + RETURN + END + + DOUBLE COMPLEX FUNCTION ZDOTC(N,CX,INCX,CY,INCY) + INTEGER INCX,INCY,N + DOUBLE COMPLEX CX(*),CY(*) + DOUBLE COMPLEX RES + EXTERNAL ZDOTCW + + CALL ZDOTCW(N,CX,INCX,CY,INCY,RES) + ZDOTC = RES + RETURN + END + + DOUBLE COMPLEX FUNCTION ZDOTU(N,CX,INCX,CY,INCY) + INTEGER INCX,INCY,N + DOUBLE COMPLEX CX(*),CY(*) + DOUBLE COMPLEX RES + EXTERNAL ZDOTUW + + CALL ZDOTUW(N,CX,INCX,CY,INCY,RES) + ZDOTU = RES + RETURN + END |