SCALAPACK 2.2.2
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches

◆ dcombssq()

subroutine dcombssq ( double precision, dimension( 2 )  v1,
double precision, dimension( 2 )  v2 
)

Definition at line 258 of file pdtreecomb.f.

259*
260* -- ScaLAPACK tools routine (version 1.7) --
261* University of Tennessee, Knoxville, Oak Ridge National Laboratory,
262* and University of California, Berkeley.
263* May 1, 1997
264*
265* .. Array Arguments ..
266 DOUBLE PRECISION V1( 2 ), V2( 2 )
267* ..
268*
269* Purpose
270* =======
271*
272* DCOMBSSQ does a scaled sum of squares on two scalars.
273*
274* Arguments
275* =========
276*
277* V1 (local input/local output) DOUBLE PRECISION array of
278* dimension 2. The first scaled sum. V1(1) = SCALE,
279* V1(2) = SUMSQ.
280*
281* V2 (local input) DOUBLE PRECISION array of dimension 2.
282* The second scaled sum. V2(1) = SCALE, V2(2) = SUMSQ.
283*
284* =====================================================================
285*
286* .. Parameters ..
287 DOUBLE PRECISION ZERO
288 parameter( zero = 0.0d+0 )
289* ..
290* .. Executable Statements ..
291*
292 IF( v1( 1 ).GE.v2( 1 ) ) THEN
293 IF( v1( 1 ).NE.zero )
294 $ v1( 2 ) = v1( 2 ) + ( v2( 1 ) / v1( 1 ) )**2 * v2( 2 )
295 ELSE
296 v1( 2 ) = v2( 2 ) + ( v1( 1 ) / v2( 1 ) )**2 * v1( 2 )
297 v1( 1 ) = v2( 1 )
298 END IF
299*
300 RETURN
301*
302* End of DCOMBSSQ
303*
Here is the caller graph for this function: