elementwise_product#

Calculates the elementwise product / Hadamard product. Multiplies every element of \(\boldsymbol{A}\) with the respective element of \(\boldsymbol{B}\) if dimensions of \(\boldsymbol{A}\), \(\boldsymbol{B}\) and \(\boldsymbol{C}\) are equal.

\[\boldsymbol{C}=\boldsymbol{A} \odot \boldsymbol{B}\]
\[\begin{split}\begin{bmatrix} a_{11} & a_{12} & a_{1n} \\ a_{21} & a_{22} & a_{2n} \\ a_{m1} & a_{m2} & a_{mn} \end{bmatrix} \begin{bmatrix} b_{11} & b_{12} & b_{1n} \\ b_{21} & b_{22} & b_{2n} \\ b_{m1} & b_{m2} & b_{mn} \end{bmatrix} = \begin{bmatrix} a_{11} \cdot b_{11} & a_{12} \cdot b_{12} & a_{1n} \cdot b_{1n} \\ a_{21} \cdot b_{21} & a_{22} \cdot b_{22} & a_{2n} \cdot b_{2n} \\ a_{m1} \cdot b_{m1} & a_{m2} \cdot b_{m2} & a_{mn} \cdot b_{mn} \end{bmatrix}\end{split}\]
void uz_matrix_elementwise_product(uz_matrix_t const *const A, uz_matrix_t const *const B, uz_matrix_t *const C_out)#

Calculates the elementwise product C_out= A .* B of all elements of the matrix A and B (also called Hadamard-Product)

Parameters:

  • A – Pointer to a uz_matrix_t instance

  • B – Pointer to a uz_matrix_t instance

  • C_out – Result of the elementwise product is written to C_out

Column vector Example#

\[\begin{split}\begin{bmatrix} 1 \\ 2 \\ 3 \end{bmatrix} \odot \begin{bmatrix} 1 \\ 2 \\ 3 \end{bmatrix} = \begin{bmatrix} 1 \\ 4 \\ 9 \end{bmatrix}\end{split}\]
 1
 2void test_uz_matrix_vector_elementwise_multiplication_3_times_1(void){
 3    uint32_t rows=3;
 4    uint32_t columns=1;
 5    float A_data[3]={1,2,3};
 6    float B_data[3]={1,2,3};
 7    float C_data[3]={134513};
 8    uz_matrix_t* A=init_array_test_helper(A_data,UZ_MATRIX_SIZE(A_data),rows, columns );
 9    uz_matrix_t* B=init_array_test_helper(B_data,UZ_MATRIX_SIZE(B_data),rows, columns );
10    uz_matrix_t* C=init_array_test_helper(C_data,UZ_MATRIX_SIZE(C_data),rows, columns );
11    // C=A .* B
12    uz_matrix_elementwise_product(A,B, C);
13    TEST_ASSERT_EQUAL_FLOAT(1,uz_matrix_get_element_zero_based(C,0,0));
14    TEST_ASSERT_EQUAL_FLOAT(4,uz_matrix_get_element_zero_based(C,1,0));
15    TEST_ASSERT_EQUAL_FLOAT(9,uz_matrix_get_element_zero_based(C,2,0));

Row vector Example#

\[\begin{bmatrix} 1 & 2 & 3 & 4 \end{bmatrix} \odot \begin{bmatrix} 1 & 2 & 3 & 4 \end{bmatrix} = \begin{bmatrix} 1 & 4 & 9 & 16 \end{bmatrix}\]
 1
 2void test_uz_matrix_matrix_elementwise_multiplication_1_times_4(void){
 3    uint32_t rows=1;
 4    uint32_t columns=4;
 5    float A_data[4]={1,2,3,4};
 6    float B_data[4]={1,2,3,4};
 7    float C_data[4]={21345};
 8    uz_matrix_t* A=init_array_test_helper(A_data,UZ_MATRIX_SIZE(A_data),rows, columns );
 9    uz_matrix_t* B=init_array_test_helper(B_data,UZ_MATRIX_SIZE(B_data),rows, columns );
10    uz_matrix_t* C=init_array_test_helper(C_data,UZ_MATRIX_SIZE(C_data),rows, columns );
11    // C=A .* B
12    uz_matrix_elementwise_product(A,B, C);
13    TEST_ASSERT_EQUAL_FLOAT(1, uz_matrix_get_element_zero_based(C,0,0) );
14    TEST_ASSERT_EQUAL_FLOAT(4, uz_matrix_get_element_zero_based(C,0,1) );
15    TEST_ASSERT_EQUAL_FLOAT(9, uz_matrix_get_element_zero_based(C,0,2) );
16    TEST_ASSERT_EQUAL_FLOAT(16,uz_matrix_get_element_zero_based(C,0,3) );

Matrix Example#

\[\begin{split}\begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} \odot \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} = \begin{bmatrix} 1 & 4 \\ 9 & 16 \end{bmatrix}\end{split}\]
 1
 2void test_uz_matrix_matrix_elementwise_multiplication_2_times_2(void){
 3    uint32_t rows=2;
 4    uint32_t columns=2;
 5    // | 1 2 | .* | 1 2 | = | 1 4  |
 6    // | 3 4 |    | 3 4 |   | 9 16 |
 7    float A_data[4]={1,2,3,4};
 8    float B_data[4]={1,2,3,4};
 9    float C_data[4]={63};
10    uz_matrix_t* A=init_array_test_helper(A_data,UZ_MATRIX_SIZE(A_data),rows, columns );
11    uz_matrix_t* B=init_array_test_helper(B_data,UZ_MATRIX_SIZE(B_data),rows, columns );
12    uz_matrix_t* C=init_array_test_helper(C_data,UZ_MATRIX_SIZE(C_data),rows, columns );
13    // C=A .* B
14    uz_matrix_elementwise_product(A,B, C);
15    TEST_ASSERT_EQUAL_FLOAT(1, uz_matrix_get_element_zero_based(C,0,0));
16    TEST_ASSERT_EQUAL_FLOAT(4, uz_matrix_get_element_zero_based(C,0,1));
17    TEST_ASSERT_EQUAL_FLOAT(9, uz_matrix_get_element_zero_based(C,1,0));
18    TEST_ASSERT_EQUAL_FLOAT(16,uz_matrix_get_element_zero_based(C,1,1));