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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
|
/*
-- 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/zpanel_to_q.cpp, normal z -> s, Wed Jan 22 14:39:05 2025
*/
#include "magma_internal.h"
/***************************************************************************//**
Put 0s in the upper triangular part of a panel and 1s on the diagonal.
Stores previous values in work array, to be restored later with magma_sq_to_panel().
@ingroup magma_panel2q
*******************************************************************************/
extern "C"
void magma_spanel_to_q(magma_uplo_t uplo, magma_int_t ib, float *A, magma_int_t lda, float *work)
{
magma_int_t i, j, k = 0;
float *col;
float c_zero = MAGMA_S_ZERO;
float c_one = MAGMA_S_ONE;
if (uplo == MagmaUpper) {
for (i = 0; i < ib; ++i) {
col = A + i*lda;
for (j = 0; j < i; ++j) {
work[k] = col[j];
col [j] = c_zero;
++k;
}
work[k] = col[i];
col [j] = c_one;
++k;
}
}
else {
for (i=0; i < ib; ++i) {
col = A + i*lda;
work[k] = col[i];
col [i] = c_one;
++k;
for (j=i+1; j < ib; ++j) {
work[k] = col[j];
col [j] = c_zero;
++k;
}
}
}
}
/***************************************************************************//**
Restores a panel, after call to magma_spanel_to_q().
@ingroup magma_panel2q
*******************************************************************************/
extern "C"
void magma_sq_to_panel(magma_uplo_t uplo, magma_int_t ib, float *A, magma_int_t lda, float *work)
{
magma_int_t i, j, k = 0;
float *col;
if (uplo == MagmaUpper) {
for (i = 0; i < ib; ++i) {
col = A + i*lda;
for (j = 0; j <= i; ++j) {
col[j] = work[k];
++k;
}
}
}
else {
for (i = 0; i < ib; ++i) {
col = A + i*lda;
for (j = i; j < ib; ++j) {
col[j] = work[k];
++k;
}
}
}
}
|
