34
34
namespace mkl_blas = oneapi::mkl::blas::row_major;
35
35
namespace mkl_lapack = oneapi::mkl::lapack;
36
36
37
-
38
37
template <typename _DataType>
39
38
void dpnp_cholesky_c (void * array1_in, void * result1, const size_t size, const size_t data_size)
40
39
{
@@ -53,7 +52,6 @@ void dpnp_cholesky_c(void* array1_in, void* result1, const size_t size, const si
53
52
for (size_t t = 0 ; t < data_size * data_size; ++t)
54
53
{
55
54
matrix[t] = in_array[k * (data_size * data_size) + t];
56
-
57
55
}
58
56
59
57
for (size_t it = 0 ; it < data_size * data_size; ++it)
@@ -65,18 +63,12 @@ void dpnp_cholesky_c(void* array1_in, void* result1, const size_t size, const si
65
63
66
64
const std::int64_t lda = std::max<size_t >(1UL , n);
67
65
68
- const std::int64_t scratchpad_size = mkl_lapack::potrf_scratchpad_size<_DataType>(
69
- DPNP_QUEUE, oneapi::mkl::uplo::upper, n, lda);
66
+ const std::int64_t scratchpad_size =
67
+ mkl_lapack::potrf_scratchpad_size<_DataType>( DPNP_QUEUE, oneapi::mkl::uplo::upper, n, lda);
70
68
71
69
_DataType* scratchpad = reinterpret_cast <_DataType*>(dpnp_memory_alloc_c (scratchpad_size * sizeof (_DataType)));
72
70
73
- event = mkl_lapack::potrf (DPNP_QUEUE,
74
- oneapi::mkl::uplo::upper,
75
- n,
76
- result_,
77
- lda,
78
- scratchpad,
79
- scratchpad_size);
71
+ event = mkl_lapack::potrf (DPNP_QUEUE, oneapi::mkl::uplo::upper, n, result_, lda, scratchpad, scratchpad_size);
80
72
81
73
event.wait ();
82
74
@@ -101,11 +93,8 @@ void dpnp_cholesky_c(void* array1_in, void* result1, const size_t size, const si
101
93
for (size_t t = 0 ; t < data_size * data_size; ++t)
102
94
{
103
95
result[k * (data_size * data_size) + t] = result_[t];
104
-
105
96
}
106
-
107
97
}
108
-
109
98
}
110
99
111
100
template <typename _DataType>
0 commit comments