Performs complex division in real arithmetic, avoiding unnecessary overflow.
Syntax
call sladiv( a, b, c, d, p, q )
call dladiv( a, b, c, d, p, q )
res = cladiv( x, y )
res = zladiv( x, y )
Include Files
- mkl.fi
Description
The routines sladiv/dladiv perform complex division in real arithmetic as
Complex functions cladiv/zladiv compute the result as
res = x/y,
where x and y are complex. The computation of x / y will not overflow on an intermediary step unless the results overflows.
The algorithm used is due to [Baudin12].
Input Parameters
- a, b, c, d
REAL for sladiv
DOUBLE PRECISION for dladiv
The scalars a, b, c, and d in the above expression (for real flavors only).
- x, y
COMPLEX for cladiv
DOUBLE COMPLEX for zladiv
The complex scalars x and y (for complex flavors only).
Output Parameters
- p, q
REAL for sladiv
DOUBLE PRECISION for dladiv
The scalars p and q in the above expression (for real flavors only).
- res
COMPLEX for cladiv
DOUBLE COMPLEX for zladiv
Contains the result of division x / y.