02a78c09404936e509bcbb3b10d9f9c6129c2263,geomstats/symmetric_positive_definite_matrices.py,,vector_to_matrix,#Any#,98

Before Change

    for j in range(dim_mat):
        for i in range(j + 1, dim_mat):
            matrix[i, j] = matrix[j, i]

    return matrix

After Change

    dim_mat = int((np.sqrt(8 * dim_vec + 1) - 1) / 2)
    matrix = np.zeros((dim_mat, dim_mat))

    lower_triangle_indices = np.tril_indices(dim_mat)
    diag_indices = np.diag_indices(dim_mat)

    matrix[lower_triangle_indices] = 2 * vector / np.sqrt(2)
    matrix[diag_indices] = vector

Instances