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
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
/*
-- MAGMA (version 2.9.0) --
Univ. of Tennessee, Knoxville
Univ. of California, Berkeley
Univ. of Colorado, Denver
@date January 2025
@author Mark Gates
@generated from control/magma_zauxiliary.cpp, normal z -> c, Wed Jan 22 14:39:04 2025
*/
#include "magma_internal.h"
#define PRECISION_c
/***************************************************************************//**
Purpose
-------
This deals with a subtle bug with returning lwork as a Float.
If lwork > 2**24, then it will get rounded as a Float;
we need to ensure it is rounded up instead of down,
by multiplying by 1.+eps in Double precision:
float( 16777217 ) == 16777216
float( 16777217 * (1.+eps) ) == 16777218
where eps is Single precision machine epsilon.
(Could use 1+2*eps in Single precision, but that can add more than necesary.)
If lwork > 2**53, rounding would happen in Double, too, but that's 94M x 94M!
@param[in]
lwork Workspace size.
@return lwork, converted to magmaFloatComplex and rounded up slightly
if necesary so that returned lwork >= input lwork.
@ingroup magma_make_lwork
*******************************************************************************/
magmaFloatComplex magma_cmake_lwork( magma_int_t lwork )
{
#if defined(PRECISION_s) || defined(PRECISION_c)
real_Double_t one_eps = 1. + lapackf77_slamch("Epsilon");
return MAGMA_C_MAKE( float(lwork*one_eps), 0 );
#else
return MAGMA_C_MAKE( float(lwork), 0 );
#endif
}
|
